Catalog footprint

What is connected

90works
49topics
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

90 published item(s)

preprint2026arXiv

A simple, flexible method for timing cross-calibration of space missions

The timing (cross-)calibration of astronomical instruments is often done by comparing pulsar times-of-arrival (TOAs) to a reference timing model. In high-energy astronomy, the choice of solar system ephemerides and source positions used to barycenter the photon arrival times has a significant impact on the procedure, requiring a full reprocessing of the data each time a new convention is used. Our method, developed as part of the activities of the International Astronomical Consortium for High Energy Calibration (IACHEC), adapts an existing pulsar solution to arbitrary JPL ephemerides and source positions by simulating geocentric TOAs and refitting timing models (implemented with PINT). We validate the procedure and apply it to thousands of observations of the Crab pulsar from 15 missions spanning 1996--2025, demonstrating inter-ephemeris TOA consistency at the $\lesssim5 μ$s level, using the DE200/FK5-based Jodrell Bank Monthly Ephemeris as a common reference. We release the TOAExtractor open-source tool and a TOA database to support future calibration and scientific studies. Instrument timing performance is broadly consistent with mission specifications; the X-ray-to-radio phase offset varies with energy and time at a level that is marginally consistent with the uncertainties of the radio ephemeris, motivating coordinated multiwavelength follow-up.

preprint2026arXiv

Digital Twin AI: Opportunities and Challenges from Large Language Models to World Models

Digital twins, as precise digital representations of physical systems, have evolved from passive simulation tools into intelligent and autonomous entities through the integration of artificial intelligence technologies. This paper presents a unified four-stage framework that systematically characterizes AI integration across the digital twin lifecycle, spanning modeling, mirroring, intervention, and autonomous management. By synthesizing existing technologies and practices, we distill a unified four-stage framework that systematically characterizes how AI methodologies are embedded across the digital twin lifecycle: (1) modeling the physical twin through physics-based and physics-informed AI approaches, (2) mirroring the physical system into a digital twin with real-time synchronization, (3) intervening in the physical twin through predictive modeling, anomaly detection, and optimization strategies, and (4) achieving autonomous management through large language models, foundation models, and intelligent agents. We analyze the synergy between physics-based modeling and data-driven learning, highlighting the shift from traditional numerical solvers to physics-informed and foundation models for physical systems. Furthermore, we examine how generative AI technologies, including large language models and generative world models, transform digital twins into proactive and self-improving cognitive systems capable of reasoning, communication, and creative scenario generation. Through a cross-domain review spanning eleven application domains, including healthcare, aerospace, smart manufacturing, robotics, and smart cities, we identify common challenges related to scalability, explainability, and trustworthiness, and outline directions for responsible AI-driven digital twin systems.

preprint2026arXiv

Disentangling Shared and Task-Specific Representations from Multi-Modal Clinical Data

Real-world clinical data is inherently multimodal, providing complementary evidence that mirrors the practical necessity of jointly assessing multiple related outcomes. Although multi-task learning can improve efficiency by sharing information across outcomes, existing approaches often fail to balance shared representation learning with outcome-specific modeling. Hard parameter sharing can trigger negative transfer when task gradients conflict, while flexible sharing may still entangle shared and task-specific signals. To address this, we propose a multi-task framework built on a unified Transformer for multimodal fusion, augmented with Orthogonal Task Decomposition (OrthTD) to split patient representations into shared and task-specific subspaces and impose a geometric orthogonality constraint to reduce redundancy and isolate task-specific signals. We evaluated OrthTD on a real-world cohort of 12,430 surgical patients for predicting four outcomes. OrthTD achieved average AUC (area under the receiver operating characteristic curve) of 87.5% and average AUPRC (area under the precision-recall curve) of 37.2%, consistently outperformed advanced tabular and multi-task methods. Notably, OrthTD achieves substantial gains in AUPRC, indicating superior performance in identifying rare events within imbalanced clinical data. These results suggest that enforcing non-redundant shared and task-specific representations can improve multi-outcome prediction from multimodal clinical data.

preprint2026arXiv

DRL-STAF: A Deep Reinforcement Learning Framework for State-Aware Forecasting of Complex Multivariate Hidden Markov Processes

Forecasting multivariate hidden Markov processes is challenging due to nonlinear and nonstationary observations, latent state transitions, and cross-sequence dependencies. While deep learning methods achieve strong predictive accuracy, they typically lack explicit state modeling, whereas Hidden Markov Models (HMMs) provide interpretable latent states but struggle with complex nonlinear emissions and scalability. To address these limitations, we propose DRL-STAF, a Deep Reinforcement Learning based STate-Aware Forecasting framework that jointly predicts next-step observations and estimates the corresponding hidden states for complex multivariate hidden Markov processes. Specifically, DRL-STAF models complex nonlinear emissions using deep neural networks and estimates discrete hidden states using reinforcement learning, reducing the reliance on predefined transition structures and enabling flexible adaptation to diverse temporal dynamics. In particular, DRL-STAF mitigates the state-space explosion encountered by typical multivariate HMM-based methods. Extensive experiments demonstrate that DRL-STAF outperforms HMM variants, standalone deep learning models, and existing DL-HMM hybrids in most cases, while also providing reliable hidden-state estimates.

preprint2026arXiv

DynaDebate: Breaking Homogeneity in Multi-Agent Debate with Dynamic Path Generation

Recent years have witnessed the rapid development of Large Language Model-based Multi-Agent Systems (MAS), which excel at collaborative decision-making and complex problem-solving. Recently, researchers have further investigated Multi-Agent Debate (MAD) frameworks, which enhance the reasoning and collaboration capabilities of MAS through information exchange and debate among multiple agents. However, existing approaches often rely on unguided initialization, causing agents to adopt identical reasoning paths that lead to the same errors. As a result, effective debate among agents is hindered, and the final outcome frequently degenerates into simple majority voting. To solve the above problem, in this paper, we introduce Dynamic Multi-Agent Debate (DynaDebate), which enhances the effectiveness of multi-agent debate through three key mechanisms: (1) Dynamic Path Generation and Allocation, which employs a dedicated Path Generation Agent to generate diverse and logical solution paths with adaptive redundancy; (2) Process-Centric Debate, which shifts the focus from surface-level outcome voting to rigorous step-by-step logic critique to ensure process correctness; (3) A Trigger-Based Verification Agent, which is activated upon disagreement and uses external tools to objectively resolve deadlocks. Extensive experiments demonstrate that DynaDebate achieves superior performance across various benchmarks, surpassing existing state-of-the-art MAD methods.

preprint2026arXiv

Foundation Models to Unlock Real-World Evidence from Nationwide Medical Claims

Evidence derived from large-scale real-world data (RWD) is increasingly informing regulatory evaluation and healthcare decision-making. Administrative claims provide population-scale, longitudinal records of healthcare utilization, expenditure, and detailed coding of diagnoses, procedures, and medications, yet their potential as a substrate for healthcare foundation models remains largely unexplored. Here we present ReClaim, a generative transformer trained from scratch on 43.8 billion medical events from more than 200 million enrollees in the MarketScan claims data spanning 2008-2022. ReClaim models longitudinal trajectories across diagnoses, procedures, medications, and expenditure, and was scaled to 140 million, 700 million, and 1.7 billion parameters. Across over 1,000 disease-onset prediction tasks, ReClaim achieved a mean AUC of 75.6%, substantially outperforming disease-specific LightGBM (66.3%) and the transformer-based Delphi model (69.4%), with the largest gains for rare diseases. These advantages held across retrospective and prospective evaluations and in external validation on two independent datasets. Performance improved monotonically with scale, and post-training added 13.8 percentage points over pre-training alone. Beyond disease prediction, ReClaim captured financial outcomes and improved real-world evidence (RWE) analyses: for healthcare expenditure forecasting it increased explained variance from 0.28 to 0.37 relative to LightGBM, and in a target trial emulation it reduced systematic bias by 72% on average relative to Delphi. Together, these results establish administrative claims as a scalable substrate for healthcare foundation models and show that learned representations generalize across time periods and data sources, supporting disease surveillance, expenditure forecasting, and RWE generation.

preprint2026arXiv

VIGIL: Defending LLM Agents Against Tool Stream Injection via Verify-Before-Commit

LLM agents operating in open environments face escalating risks from indirect prompt injection, particularly within the tool stream where manipulated metadata and runtime feedback hijack execution flow. Existing defenses encounter a critical dilemma as advanced models prioritize injected rules due to strict alignment while static protection mechanisms sever the feedback loop required for adaptive reasoning. To reconcile this conflict, we propose \textbf{VIGIL}, a framework that shifts the paradigm from restrictive isolation to a verify-before-commit protocol. By facilitating speculative hypothesis generation and enforcing safety through intent-grounded verification, \textbf{VIGIL} preserves reasoning flexibility while ensuring robust control. We further introduce \textbf{SIREN}, a benchmark comprising 959 tool stream injection cases designed to simulate pervasive threats characterized by dynamic dependencies. Extensive experiments demonstrate that \textbf{VIGIL} outperforms state-of-the-art dynamic defenses by reducing the attack success rate by over 22\% while more than doubling the utility under attack compared to static baselines, thereby achieving an optimal balance between security and utility.

preprint2025arXiv

PDA in Action: Ten Principles for High-Quality Multi-Site Clinical Evidence Generation

Background: Distributed Research Networks (DRNs) offer significant opportunities for collaborative multi-site research and have significantly advanced healthcare research based on clinical observational data. However, generating high-quality real-world evidence using fit-for-use data from multi-site studies faces important challenges, including biases associated with various types of heterogeneity within and across sites and data sharing difficulties. Over the last ten years, Privacy-Preserving Distributed Algorithms (PDA) have been developed and utilized in numerous national and international real-world studies spanning diverse domains, from comparative effectiveness research, target trial emulation, to healthcare delivery, policy evaluation, and system performance assessment. Despite these advances, there remains a lack of comprehensive and clear guiding principles for generating high-quality real-world evidence through collaborative studies leveraging the methods under PDA. Objective: The paper aims to establish ten principles of best practice for conducting high-quality multi-site studies using PDA. These principles cover all phases of research, including study preparation, protocol development, analysis, and final reporting. Discussion: The ten principles for conducting a PDA study outline a principled, efficient, and transparent framework for employing distributed learning algorithms within DRNs to generate reliable and reproducible real-world evidence.

preprint2023arXiv

Unconventional ferroelectricity in half-filling states of antiparallel stacking of twisted WSe2

Abstract: We report on emergence of an abnormal electronic polarization in twisted double bilayer WSe2 in antiparallel interface stacking geometry, where local centrosymmetry of atomic registries at the twist interface does not favor the spontaneous electronic polarizations as recently observed in the parallel interface stacking geometry. The unconventional ferroelectric behaviors probed by electronic transport measurement occur at half filling insulating states at 1.5 K and gradually disappear at about 40 K. Single band Hubbard model based on the triangular moiré lattice and the interlayer charge transfer controlled by insulating phase transition are proposed to interpret the formation of electronic polarization states near half filling in twisted WSe2 devices. Our work highlights the prominent role of many-body electronic interaction in fostering novel quantum states in moiré-structured systems.

preprint2023arXiv

Workload Failure Prediction for Data Centers

Failed workloads that consumed significant computational resources in time and space affect the efficiency of data centers significantly and thus limit the amount of scientific work that can be achieved. While the computational power has increased significantly over the years, detection and prediction of workload failures have lagged far behind and will become increasingly critical as the system scale and complexity further increase. In this study, we analyze workload traces collected from a production cluster and train machine learning models on a large amount of data sets to predict workload failures. Our prediction models consist of a queue-time model that estimates the probability of workload failures before execution and a runtime model that predicts failures at runtime. Evaluation results show that the queue-time model and runtime model can predict workload failures with a maximum precision score of 90.61% and 97.75%, respectively. By integrating the runtime model with the job scheduler, it helps reduce CPU time, and memory usage by up to 16.7% and 14.53%, respectively.

preprint2022arXiv

Data-driven vector localized waves and parameters discovery for Manakov system using deep learning approach

An improved physics-informed neural network (IPINN) algorithm with four output functions and four physics constraints, which possesses neuron-wise locally adaptive activation function and slope recovery term, is appropriately proposed to obtain the data-driven vector localized waves, including vector solitons, breathers and rogue waves (RWs) for the Manakov system with initial and boundary conditions, as well as data-driven parameters discovery for Manakov system with unknown parameters. The data-driven vector RWs which also contain interaction waves of RWs and bright-dark solitons, interaction waves of RWs and breathers, as well as RWs evolved from bright-dark solitons are learned to verify the capability of the IPINN algorithm in training complex localized wave. In the process of parameter discovery, routine IPINN can not accurately train unknown parameters whether using clean data or noisy data. Thus we introduce parameter regularization strategy with adjustable weight coefficients into IPINN to effectively and accurately train prediction parameters, then find that once setting the appropriate weight coefficients, the training effect is better as using noisy data. Numerical results show that IPINN with parameter regularization shows superior noise immunity in parameters discovery problem.

preprint2022arXiv

Double and Triple-Pole Solutions for the Third-Order Flow Equation of the Kaup-Newell System with Zero/Nonzero Boundary Conditions

In this work, the double and triple-pole solutions for the third-order flow equation of Kaup-Newell system (TOFKN) with zero boundary conditions (ZBCs) and non-zero boundary conditions (NZBCs) are investigated by means of the Riemann-Hilbert (RH) approach stemming from the inverse scattering transformation. Starting from spectral problem of the TOFKN, the analyticity, symmetries, asymptotic behavior of the Jost function and scattering matrix, the matrix RH problem with ZBCs and NZBCs are constructed. Then the obtained RH problem with ZBCs and NZBCs can be solved in the case of scattering coefficients with double or triple zeros, and the reconstruction formula of potential, trace formula as well as theta condition are also derived correspondingly. Specifically, the general formulas of $N$-double and $N$-triple poles solutions with ZBCs and NZBCs are derived systematically by means of determinants. The vivid plots and dynamics analysis for double and triple-pole soliton solutions with the ZBCs as well as double and triple-pole interaction solutions with the NZBCs are exhibited in details. Compared with the most classical second-order flow Kaup-Newell system, we find the third-order dispersion and quintic nonlinear term of the Kaup-Newell system change the trajectory and velocity of solutions. Furthermore, the asymptotic states of the 1-double poles soliton solution and the 1-triple poles soliton solution are analyzed when $t$ tends to infinity.

preprint2022arXiv

DV-QKD Coexistence With 1.6 Tbps Classical Channels Over Hollow Core Fibre

The feasibility of coexisting a quantum channel with carrier-grade classical optical channels over Hollow Core Nested Antiresonant Nodeless Fibre (HC-NANF) is experimentally explored for the first time in terms of achievable quantum bit error rate (QBER), secret key rate (SKR) as well as classical signal bit error rates (BER). A coexistence transmission of 1.6 Tbps is achieved for the classical channels simultaneously with a quantum channel over a 2 km-long HC-NANF with a total coexistence power of 0 dBm. To find the best and worst wavelength position for the classical channels, we simulated different classical channels bands with different spacing between the quantum and classical channels considering the crosstalk generated from both Raman scattering and four-wave-mixing (FWM) on the quantum channel. Following our simulation, we numerically estimate the best (Raman spectrum dip) and worst locations (Raman spectrum peak) of the classical channel with respect to its impact on the performance on the quantum channel in terms of SKR and QBER. We further implemented a testbed to experimentally test both single mode fibre (SMF) and HC-NANF in the best and worst-case scenarios. In the best-case scenario, the spacing between quantum and classical is 200 GHz (1.6 nm) with 50 GHz (0.4 nm) spacing between each classical channel. The SKR was preserved without any noticeable changes when coexisting the quantum channel with eight classical channels at 0 dBm total coexistence power in HC-NANF compared to a significant drop of 73% when using SMF at -24 dBm total coexistence power which is 250 times lower than the power used in HC-NANF. In the worst-case scenario using the same powers, and with 1 THz (8 nm) spacing between quantum and classical channels, the SKR dropped 10% using the HC-NANF, whereas in the SMF the SKR plummeted to zero.

preprint2022arXiv

Federated Learning Algorithms for Generalized Mixed-effects Model (GLMM) on Horizontally Partitioned Data from Distributed Sources

Objectives: This paper develops two algorithms to achieve federated generalized linear mixed effect models (GLMM), and compares the developed model's outcomes with each other, as well as that from the standard R package (`lme4'). Methods: The log-likelihood function of GLMM is approximated by two numerical methods (Laplace approximation and Gaussian Hermite approximation), which supports federated decomposition of GLMM to bring computation to data. Results: Our developed method can handle GLMM to accommodate hierarchical data with multiple non-independent levels of observations in a federated setting. The experiment results demonstrate comparable (Laplace) and superior (Gaussian-Hermite) performances with simulated and real-world data. Conclusion: We developed and compared federated GLMMs with different approximations, which can support researchers in analyzing biomedical data to accommodate mixed effects and address non-independence due to hierarchical structures (i.e., institutes, region, country, etc.).

preprint2022arXiv

Kernel representation formula from complex to real Wiener-Ito integrals and vice versa

We clearly characterize the relation between real and complex Wiener-Ito integrals. Given a complex multiple Wiener-Ito integral, we get explicit expressions for two kernels of its real and imaginary parts. Conversely, consider a two-dimensional real Wiener-Ito integral, we obtain the representation formula by a finite sum of complex Wiener-Ito integrals. The main tools are a recursion technique and Malliavin derivative operators. We build a bridge between real and complex Wiener-Ito integrals.

preprint2022arXiv

Long time asymptotic analysis for a nonlocal Hirota equation via the Dbar steepest descent method

In this paper, we mainly focus on the Cauchy problem of an integrable nonlocal Hirota equation with initial value in weighted Sobolev space. Through the spectral analysis of Lax pairs, we successfully transform the Cauchy problem of the nonlocal Hirota equation into a solvable Riemann-Hilbert problem. Furthermore, in the absence of discrete spectrum, the long-time asymptotic behavior of the solution for the nonlocal Hirota equation is obtained through the Dbar steepest descent method. Different from the local Hirota equation, the leading order term on the continuous spectrum and residual error term of $q(x,t)$ are affected by the function $Imν(z_j)$.

preprint2022arXiv

Multiple-high-order pole solutions for the NLS equation with quartic terms

The aim of this article is to investigate the multiple-high-order pole solutions to the focusing NLS equation with quartic terms(QNLS) under the non-vanishing boundary conditions(NVBC) via the Riemann-Hilbert(RH) method. The determinant formula of multiple-high-order pole soliton solutions for NVBC is given. Further the double 1nd-order, mixed 2nd- and 1nd-order pole solutions are obtained.

preprint2022arXiv

Physics-informed neural network methods based on Miura transformations and discovery of new localized wave solutions

We put forth two physics-informed neural network (PINN) schemes based on Miura transformations and the novelty of this research is the incorporation of Miura transformation constraints into neural networks to solve nonlinear PDEs. The most noteworthy advantage of our method is that we can simply exploit the initial-boundary data of a solution of a certain nonlinear equation to obtain the data-driven solution of another evolution equation with the aid of PINNs and during the process, the Miura transformation plays an indispensable role of a bridge between solutions of two separate equations. It is tailored to the inverse process of the Miura transformation and can overcome the difficulties in solving solutions based on the implicit expression. Moreover, two schemes are applied to perform abundant computational experiments to effectively reproduce dynamic behaviors of solutions for the well-known KdV equation and mKdV equation. Significantly, new data-driven solutions are successfully simulated and one of the most important results is the discovery of a new localized wave solution: kink-bell type solution of the defocusing mKdV equation and it has not been previously observed and reported to our knowledge. It provides a possibility for new types of numerical solutions by fully leveraging the many-to-one relationship between solutions before and after Miura transformations. Performance comparisons in different cases as well as advantages and disadvantages analysis of two schemes are also discussed. On the basis of the performance of two schemes and no free lunch theorem, they both have their own merits and thus more appropriate one should be chosen according to specific cases.

preprint2022arXiv

Quasi-periodic oscillations of the X-ray burst from the magnetar SGR J1935+2154 and associated with the fast radio burst FRB 200428

The origin(s) and mechanism(s) of fast radio bursts (FRBs), which are short radio pulses from cosmological distances, have remained a major puzzle since their discovery. We report a strong Quasi-Periodic Oscillation(QPO) of 40 Hz in the X-ray burst from the magnetar SGR J1935+2154 and associated with FRB 200428, significantly detected with the Hard X-ray Modulation Telescope (Insight-HXMT) and also hinted by the Konus-Wind data. QPOs from magnetar bursts have only been rarely detected; our 3.4 sigma (p-value is 2.9e-4) detection of the QPO reported here reveals the strongest QPO signal observed from magnetars (except in some very rare giant flares), making this X-ray burst unique among magnetar bursts. The two X-ray spikes coinciding with the two FRB pulses are also among the peaks of the QPO. Our results suggest that at least some FRBs are related to strong oscillation processes of neutron stars. We also show that we may overestimate the significance of the QPO signal and underestimate the errors of QPO parameters if QPO exists only in a fraction of the time series of a X-ray burst which we use to calculate the Leahy-normalized periodogram.

preprint2022arXiv

The Berry-Esséen Upper Bounds of Vasicek Model Estimators

The Berry-Esséen upper bounds of moment estimators and least squares estimators of the mean and drift coefficients in Vasicek models driven by general Gaussian processes are studied. When studying the parameter estimation problem of Ornstein-Uhlenbeck (OU) process driven by fractional Brownian motion, the commonly used methods are mainly given by Kim and Park, they show the upper bound of Kolmogorov distance between the distribution of the ratio of two double Wiener-Itô stochastic integrals and the Normal distribution. The main innovation in this paper is extending the above ratio process, that is to say, the numerator and denominator respectively contain triple Wiener-Itô stochastic integrals at most. As far as we know, the upper bounds between the distribution of above estimators and the Normal distribution are novel.

preprint2022arXiv

Two-stage Hypothesis Tests for Variable Interactions with FDR Control

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

preprint2021arXiv

$N$-double poles solutions for nonlocal Hirota equation with nonzero boundary conditions using Riemann-Hilbert method and PINN algorithm

We systematically investigate the nonlocal Hirota equation with nonzero boundary conditions via Riemann-Hilbert method and multi-layer physics-informed neural networks algorithm. Starting from the Lax pair of nonzero nonlocal Hirota equation, we first give out the Jost function, scattering matrix, their symmetry and asymptotic behavior. Then, the Riemann-Hilbert problem with nonzero boundary conditions are constructed and the precise formulaes of $N$-double poles solutions and $N$-simple poles solutions are written by determinants. Different from the local Hirota equation, the symmetry of scattering data for nonlocal Hirota equation is completely different, which results in disparate discrete spectral distribution. In particular, it could be more complicated and difficult to obtain the symmetry of scattering data under the circumstance of double poles. Besides, we also analyse the asymptotic state of one-double poles solution as $t\rightarrow \infty$. Whereafter, the PINN algorithm is applied to research the data-driven soliton solutions of the nonzero nonlocal Hirota equation by using the training data obtained from the Riemann-Hilbert method. Most strikingly, the integrable nonlocal equation is firstly solved via PINN algorithm. As we all know, the nonlocal equations contain the $\mathcal{PT}$ symmetry $\mathcal{P}:x\rightarrow -x,$ or $\mathcal{T}:t\rightarrow -t,$ which are different with local ones. Adding the nonlocal term into the NN, we can successfully solve the integrable nonlocal Hirota equation by PINN algorithm. The numerical results indicate the algorithm can well recover the data-driven soliton solutions of the integrable nonlocal equation. Noteworthily, the inverse problems of the integrable nonlocal equation are discussed for the first time through applying the PINN algorithm to discover the parameters of the equation in terms of its soliton solution.

preprint2021arXiv

A Robust Bayesian Copas Selection Model for Quantifying and Correcting Publication Bias

The validity of conclusions from meta-analysis is potentially threatened by publication bias. Most existing procedures for correcting publication bias assume normality of the study-specific effects that account for between-study heterogeneity. However, this assumption may not be valid, and the performance of these bias correction procedures can be highly sensitive to departures from normality. Further, there exist few measures to quantify the magnitude of publication bias based on selection models. In this paper, we address both of these issues. First, we explore the use of heavy-tailed distributions for the study-specific effects within a Bayesian hierarchical framework. The deviance information criterion (DIC) is used to determine the appropriate distribution to use for conducting the final analysis. Second, we develop a new measure to quantify the magnitude of publication bias based on Hellinger distance. Our measure is easy to interpret and takes advantage of the estimation uncertainty afforded naturally by the posterior distribution. We illustrate our proposed approach through simulation studies and meta-analyses on lung cancer and antidepressants. To assess the prevalence of publication bias, we apply our method to 1500 meta-analyses of dichotomous outcomes in the Cochrane Database of Systematic Reviews. Our methods are implemented in the publicly available R package RobustBayesianCopas.

preprint2021arXiv

A two-stage physics-informed neural network method based on conserved quantities and applications in localized wave solutions

With the advantages of fast calculating speed and high precision, the physics-informed neural network method opens up a new approach for numerically solving nonlinear partial differential equations. Based on conserved quantities, we devise a two-stage PINN method which is tailored to the nature of equations by introducing features of physical systems into neural networks. Its remarkable advantage lies in that it can impose physical constraints from a global perspective. In stage one, the original PINN is applied. In stage two, we additionally introduce the measurement of conserved quantities into mean squared error loss to train neural networks. This two-stage PINN method is utilized to simulate abundant localized wave solutions of integrable equations. We mainly study the Sawada-Kotera equation as well as the coupled equations: the classical Boussinesq-Burgers equations and acquire the data-driven soliton molecule, M-shape double-peak soliton, plateau soliton, interaction solution, etc. Numerical results illustrate that abundant dynamic behaviors of these solutions can be well reproduced and the two-stage PINN method can remarkably improve prediction accuracy and enhance the ability of generalization compared to the original PINN method.

preprint2021arXiv

Combining Cox Regressions Across a Heterogeneous Distributed Research Network Facing Small and Zero Counts

Studies of the effects of medical interventions increasingly take place in distributed research settings using data from multiple clinical data sources including electronic health records and administrative claims. In such settings, privacy concerns typically prohibit sharing of individual patient data, and instead, analyses can only utilize summary statistics from the individual databases. In the specific but very common context of the Cox proportional hazards model, we show that standard meta analysis methods then lead to substantial bias when outcome counts are small. This bias derives primarily from the normal approximations that the methods utilize. Here we propose and evaluate methods that eschew normal approximations in favor of three more flexible approximations: a skew-normal, a one-dimensional grid, and a custom parametric function that mimics the behavior of the Cox likelihood function. In extensive simulation studies we demonstrate how these approximations impact bias in the context of both fixed-effects and (Bayesian) random-effects models. We then apply these approaches to three real-world studies of the comparative safety of antidepressants, each using data from four observational healthcare databases.

preprint2021arXiv

Double and triple poles solutions for the Gerdjikov-Ivanov type of derivative nonlinear Schrödinger equation with zero/nonzero boundary conditions

In this work, the double and triple poles soliton solutions for the Gerdjikov-Ivanov(GI) type of derivative nonlinear Schrödinger equation with zero boundary conditions(ZBCs) and nonzero boundary conditions(NZBCs) are studied via Riemann-Hilbert (RH) method. Though spectral problem analysis, we first give out the Jost function and scattering matrix under ZBCs and NZBCs. Then according to the analyticity, symmetry and asymptotic behavior of Jost function and scattering matrix, the Riemann-Hilbert problem(RHP) with ZBCs and NZBCs are constructed. Further, the obtained RHP with ZBCs and NZBCs can be solved in the case that reflection coefficients have double or triple poles. Finally, we derive the general precise formulae of N-double and N-triple poles solutions corresponding to ZBCs and NZBCs, respectively. The dynamical behaviors for these solutions are further discussed by image simulation.

preprint2021arXiv

In-orbit timing calibration of the Insight-Hard X-ray Modulation Telescope

We describe the timing system and the timing calibration results of the three payloads on-board the Insight-Hard X-ray Modulation Telescope (Insight-HXMT). These three payloads are the High Energy X-ray telescope (HE, 20-250 keV), the Medium Energy X-ray telescope (ME, 5-30 keV) and the low Energy X-ray telescope (LE, 1-10 keV). We present a method to correct the temperature-dependent period response and the long-term variation of the on-board crystal oscillator, especially for ME that does not carry a temperature-compensated crystal oscillator. The time of arrivals (ToAs) of the Crab pulsar are measured to evaluate the accuracy of the timing system. As the ephemeris of the Crab pulsar given by Jodrell Bank observatory has systematic errors around 40 μs (Rots et al. 2014), we use the quasi-simultaneous observations of the X-ray Timing Instrument (XTI) on-board the Neutron star Interior Composition Explorer (NICER) to produce the Crab ephemerides and to verify the timing system of Insight-HXMT. The energy-dependent ToAs' offsets relative to the NICER measurements including physical and instrumental origins are about 24.7μs, 10.1μs and 864.7μs, and the systematic errors of the timing system are determined as 12.1μs, 8.6μs, and 15.8μs, for HE, ME and LE respectively.

preprint2021arXiv

Performance of a focal plane detector for soft X-ray imaging spectroscopy based on back-illuminated sCMOS

Spectroscopy focusing array (SFA) and Polarimetry focusing array (PFA) are the two major payloads of enhanced X-ray Timing and Polarimetry mission (eXTP). Nested Wolter-\RNum{1} X-ray mirror module is implemented in SFA and PFA to achive high effective area. When evaluating the properties of the mirror module, the alignment of the optical axis of the X-ray mirror module and a quasi-parallel X-ray beam is a prerequisite to ensure the accuracy of the results. Hence, to assist the alignment of the X-ray mirror module, an X-ray focal plane detector is designed based on the back-illuminated scientific Complementary Metal-Oxide-Semiconductor Transistor (sCMOS) sensor GSENSE6060BSI, one of the largest detection areas, is produced by \textit{Gpixel Inc}. Then the characteristics of readout noise, dark current, and split-pixel event properties of the detector are studied with the self-developed multi-target fluorescence X-ray source in a 100 m long X-ray test facility. The energy calibration is carried out with the single-pixel event and the energy non-linearity of the detector is also obtained. Eventually, the simulation of the eXTP mirror module based on the optical model is conducted and the alignment test of the Wolter-\RNum{1} X-ray mirror module designed for \textit{EP/FXT} (Einstein Probe/Follow-up X-ray Telescope) with "Burkert test" method is shown.

preprint2021arXiv

Physics-informed neural networks method in high-dimensional integrable systems

In this paper, the physics-informed neural networks (PINN) is applied to high-dimensional system to solve the (N+1)-dimensional initial boundary value problem with 2N+1 hyperplane boundaries. This method is used to solve the most classic (2+1)-dimensional integrable Kadomtsev-Petviashvili (KP) equation and (3+1)-dimensional reduced KP equation. The dynamics of (2+1)-dimensional local waves such as solitons, breathers, lump and resonance rogue are reproduced. Numerical results display that the magnitude of the error is much smaller than the wave height itself, so it is considered that the classical solutions in these integrable systems are well obtained based on the data-driven mechanism.

preprint2021arXiv

Rapid chemically selective 3D imaging in the mid-infrared with a Si-based camera

The emerging technique of mid-infrared optical coherence tomography (MIR-OCT) takes advantage of the reduced scattering of MIR light in various materials and devices, enabling tomographic imaging at deeper penetration depths. Because of challenges in MIR detection technology, the image acquisition time is however significantly longer than for tomographic imaging methods in the visible/near-infrared. Here we demonstrate an alternative approach to MIR tomography with high-speed imaging capabilities. Through femtosecond non-degenerate two-photon absorption of MIR light in a conventional Si-based CCD camera, we achieve wide-field, high-definition tomographic imaging with chemical selectivity of structured materials and biological samples in mere seconds.

preprint2021arXiv

Solving localized wave solutions of the derivative nonlinear Schrodinger equation using an improved PINN method

The solving of the derivative nonlinear Schrodinger equation (DNLS) has attracted considerable attention in theoretical analysis and physical applications. Based on the physics-informed neural network (PINN) which has been put forward to uncover dynamical behaviors of nonlinear partial different equation from spatiotemporal data directly, an improved PINN method with neuron-wise locally adaptive activation function is presented to derive localized wave solutions of the DNLS in complex space. In order to compare the performance of above two methods, we reveal the dynamical behaviors and error analysis for localized wave solutions which include one-rational soliton solution, genuine rational soliton solutions and rogue wave solution of the DNLS by employing two methods, also exhibit vivid diagrams and detailed analysis. The numerical results demonstrate the improved method has faster convergence and better simulation effect. On the bases of the improved method, the effects for different numbers of initial points sampled, residual collocation points sampled, network layers, neurons per hidden layer on the second order genuine rational soliton solution dynamics of the DNLS are considered, and the relevant analysis when the locally adaptive activation function chooses different initial values of scalable parameters are also exhibited in the simulation of the two-order rogue wave solution.

preprint2020arXiv

A regression-based method for detecting publication bias in multivariate meta-analysis

Publication bias occurs when the publication of research results depends not only on the quality of the research but also on its nature and direction. The consequence is that published studies may not be truly representative of all valid studies undertaken, and this bias may threaten the validity of systematic reviews and meta-analyses - on which evidence-based medicine increasingly relies. Multivariate meta-analysis has recently received increasing attention for its ability reducing potential bias and improving statistical efficiency by borrowing information across outcomes. However, detecting and accounting for publication bias are more challenging in multivariate meta-analysis setting because some studies may be completely unpublished whereas some studies may selectively report part of multiple outcomes. In this paper, we propose a score test for jointly testing publication bias for multiple outcomes, which is novel to the multivariate setting. The proposed test is a natural multivariate extension of the univariate Egger's test, and can handle the above mentioned scenarios simultaneously, It accounts for correlations among multivariate outcomes, while allowing different types of outcomes, and can borrow information across outcomes. The proposed test is shown to be more powerful than the Egger's test, Begg's test and Trim and Fill method through simulation studies. Two data analyses are given to illustrate the performance of the proposed test in practice.

preprint2020arXiv

Background Model for the Low-Energy Telescope of Insight-HXMT

With more than 150 blank sky observations at high Galactic latitude, we make a systematic study to the background of the Low Energy Telescope (LE) of the Hard X-ray Modulation Telescope (dubbed as Insight-HXMT). Both the on-ground simulation and the in-orbit observation indicate that the background spectrum mainly has two components. One is the particle background that dominates above 7 keV and its spectral shape is consistent in every geographical locations. Another is the diffuse X-ray background that dominates below 7 keV and has a stable spectrum less dependent of the sky region. The particle background spectral shape can be obtained from the blind detector data of all the blank sky observations, and the particle background intensity can be measured by the blind detector at 10-12.5 keV. The diffuse X-ray background in the high Galactic latitude can also be obtained from the blank sky spectra after subtracting the particle background. Based on these characteristics, we develop the background model for both the spectrum and the light curve. The systematic error for the background spectrum is investigated with different exposures (T_exp). For the spectrum with T_exp=1 ks, the average systematic errors in 1-7 keV and 1-10 keV are 4.2% and 3.7%, respectively. We also perform the systematic error analyses of the background light curves with different energy bands and time bins. The results show that the systematic errors for the light curves with different time bins are <8% in 1-10 keV.

preprint2020arXiv

Data Centers Job Scheduling with Deep Reinforcement Learning

Efficient job scheduling on data centers under heterogeneous complexity is crucial but challenging since it involves the allocation of multi-dimensional resources over time and space. To adapt the complex computing environment in data centers, we proposed an innovative Advantage Actor-Critic (A2C) deep reinforcement learning based approach called A2cScheduler for job scheduling. A2cScheduler consists of two agents, one of which, dubbed the actor, is responsible for learning the scheduling policy automatically and the other one, the critic, reduces the estimation error. Unlike previous policy gradient approaches, A2cScheduler is designed to reduce the gradient estimation variance and to update parameters efficiently. We show that the A2cScheduler can achieve competitive scheduling performance using both simulated workloads and real data collected from an academic data center.

preprint2020arXiv

Discovery of oscillations above 200 keV in a black hole X-ray binary with Insight-HXMT

Low-frequency quasi-periodic oscillations (LFQPOs) are commonly found in black hole X-ray binaries, and their origin is still under debate. The properties of LFQPOs at high energies (above 30 keV) are closely related to the nature of the accretion flow in the innermost regions, and thus play a crucial role in critically testing various theoretical models. The Hard X-ray Modulation Telescope (Insight-HXMT) is capable of detecting emissions above 30 keV, and is therefore an ideal instrument to do so. Here we report the discovery of LFQPOs above 200 keV in the new black hole MAXI J1820+070 in the X-ray hard state, which allows us to understand the behaviours of LFQPOs at hundreds of kiloelectronvolts. The phase lag of the LFQPO is constant around zero below 30 keV, and becomes a soft lag (that is, the high-energy photons arrive first) above 30 keV. The soft lag gradually increases with energy and reaches ~0.9s in the 150-200 keV band. The detection at energies above 200 keV, the large soft lag and the energy-related behaviors of the LFQPO pose a great challenge for most currently existing models, but suggest that the LFQPO probably originates from the precession of a small-scale jet.

preprint2020arXiv

Fixed-effects model: the most convincing model for meta-analysis with few studies

According to Davey et al. (2011) with a total of 22,453 meta-analyses from the January 2008 Issue of the Cochrane Database of Systematic Reviews, the median number of studies included in each meta-analysis is only three. In other words, about a half or more of meta-analyses conducted in the literature include only two or three studies. While the common-effect model (also referred to as the fixed-effect model) may lead to misleading results when the heterogeneity among studies is large, the conclusions based on the random-effects model may also be unreliable when the number of studies is small. Alternatively, the fixed-effects model avoids the restrictive assumption in the common-effect model and the need to estimate the between-study variance in the random-effects model. We note, however, that the fixed-effects model is under appreciated and rarely used in practice until recently. In this paper, we compare all three models and demonstrate the usefulness of the fixed-effects model when the number of studies is small. In addition, we propose a new estimator for the unweighted average effect in the fixed-effects model. Simulations and real examples are also used to illustrate the benefits of the fixed-effects model and the new estimator.

preprint2020arXiv

Hyperspectral Super-Resolution via Coupled Tensor Ring Factorization

Hyperspectral super-resolution (HSR) fuses a low-resolution hyperspectral image (HSI) and a high-resolution multispectral image (MSI) to obtain a high-resolution HSI (HR-HSI). In this paper, we propose a new model, named coupled tensor ring factorization (CTRF), for HSR. The proposed CTRF approach simultaneously learns high spectral resolution core tensor from the HSI and high spatial resolution core tensors from the MSI, and reconstructs the HR-HSI via tensor ring (TR) representation (Figure~\ref{fig:framework}). The CTRF model can separately exploit the low-rank property of each class (Section \ref{sec:analysis}), which has been never explored in the previous coupled tensor model. Meanwhile, it inherits the simple representation of coupled matrix/CP factorization and flexible low-rank exploration of coupled Tucker factorization. Guided by Theorem~\ref{th:1}, we further propose a spectral nuclear norm regularization to explore the global spectral low-rank property. The experiments have demonstrated the advantage of the proposed nuclear norm regularized CTRF (NCTRF) as compared to previous matrix/tensor and deep learning methods.

preprint2020arXiv

In-flight calibration of the Insight-Hard X-ray Modulation Telescope

We present the calibration of the Insight-Hard X-ray Modulation Telescope (Insight-HXMT) X-ray satellite, which can be used to perform timing and spectral studies of bright X-ray sources. Insight-HXMT carries three main payloads onboard: the High Energy X-ray telescope (HE), the Medium Energy X-ray telescope (ME) and the Low Energy X-ray telescope (LE). In orbit, the radioactive sources, activated lines, the fluorescence lines and celestial sources are used to calibrate the energy scale and energy resolution of the payloads. The Crab nebular is adopted as the primary effective area calibrator and empirical functions are constructed to modify the simulated effective areas of the three payloads respectively. The systematic errors of HE, compared to the model of the Crab nebular, are less than 2% in 28--120 keV and 2%--10% above 120 keV. The systematic errors of ME are less than 1.5% in 10--35 keV. The systematic errors of LE are less than 1% in 1--7 keV except the Si K--edge (1.839 keV, up to 1.5%) and less than 2% in 7--10 keV.

preprint2020arXiv

Missing at Random or Not: A Semiparametric Testing Approach

Practical problems with missing data are common, and statistical methods have been developed concerning the validity and/or efficiency of statistical procedures. On a central focus, there have been longstanding interests on the mechanism governing data missingness, and correctly deciding the appropriate mechanism is crucially relevant for conducting proper practical investigations. The conventional notions include the three common potential classes -- missing completely at random, missing at random, and missing not at random. In this paper, we present a new hypothesis testing approach for deciding between missing at random and missing not at random. Since the potential alternatives of missing at random are broad, we focus our investigation on a general class of models with instrumental variables for data missing not at random. Our setting is broadly applicable, thanks to that the model concerning the missing data is nonparametric, requiring no explicit model specification for the data missingness. The foundational idea is to develop appropriate discrepancy measures between estimators whose properties significantly differ only when missing at random does not hold. We show that our new hypothesis testing approach achieves an objective data oriented choice between missing at random or not. We demonstrate the feasibility, validity, and efficacy of the new test by theoretical analysis, simulation studies, and a real data analysis.

preprint2020arXiv

Parameter estimation for an Ornstein-Uhlenbeck Process driven by a general Gaussian noise

In this paper, we consider an inference problem for an Ornstein-Uhlenbeck process driven by a general one-dimensional centered Gaussian process $(G_t)_{t\ge 0}$. The second order mixed partial derivative of the covariance function $ R(t,\, s)=\mathbb{E}[G_t G_s]$ can be decomposed into two parts, one of which coincides with that of fractional Brownian motion and the other is bounded by $(ts)^{β-1}$ up to a constant factor. This condition is valid for a class of continuous Gaussian processes that fails to be self-similar or have stationary increments. Some examples include the subfractional Brownian motion and the bi-fractional Brownian motion. Under this assumption, we study the parameter estimation for drift parameter in the Ornstein-Uhlenbeck process driven by the Gaussian noise $(G_t)_{t\ge 0}$. For the least squares estimator and the second moment estimator constructed from the continuous observations, we prove the strong consistency and the asympotic normality, and obtain the Berry-Esséen bounds. The proof is based on the inner product's representation of the Hilbert space $\mathfrak{H}$ associated with the Gaussian noise $(G_t)_{t\ge 0}$, and the estimation of the inner product based on the results of the Hilbert space associated with the fractional Brownian motion.

preprint2020arXiv

Spike-and-Slab Group Lassos for Grouped Regression and Sparse Generalized Additive Models

We introduce the spike-and-slab group lasso (SSGL) for Bayesian estimation and variable selection in linear regression with grouped variables. We further extend the SSGL to sparse generalized additive models (GAMs), thereby introducing the first nonparametric variant of the spike-and-slab lasso methodology. Our model simultaneously performs group selection and estimation, while our fully Bayes treatment of the mixture proportion allows for model complexity control and automatic self-adaptivity to different levels of sparsity. We develop theory to uniquely characterize the global posterior mode under the SSGL and introduce a highly efficient block coordinate ascent algorithm for maximum a posteriori (MAP) estimation. We further employ de-biasing methods to provide uncertainty quantification of our estimates. Thus, implementation of our model avoids the computational intensiveness of Markov chain Monte Carlo (MCMC) in high dimensions. We derive posterior concentration rates for both grouped linear regression and sparse GAMs when the number of covariates grows at nearly exponential rate with sample size. Finally, we illustrate our methodology through extensive simulations and data analysis.

preprint2020arXiv

Testing for publication bias in meta-analysis under Copas selection model

In meta-analyses, publication bias is a well-known, important and challenging issue because the validity of the results from a meta-analysis is threatened if the sample of studies retrieved for review is biased. One popular method to deal with publication bias is the Copas selection model, which provides a flexible sensitivity analysis for correcting the estimates with considerable insight into the data suppression mechanism. However, rigorous testing procedures under the Copas selection model to detect bias are lacking. To fill this gap, we develop a score-based test for detecting publication bias under the Copas selection model. We reveal that the behavior of the standard score test statistic is irregular because the parameters of the Copas selection model disappear under the null hypothesis, leading to an identifiability problem. We propose a novel test statistic and derive its limiting distribution. A bootstrap procedure is provided to obtain the p-value of the test for practical applications. We conduct extensive Monte Carlo simulations to evaluate the performance of the proposed test and apply the method to several existing meta-analyses.

preprint2020arXiv

The influence of the Insight-HXMT/LE time response on timing analysis

LE is the low energy telescope of Insight-HXMT. It uses swept charge devices (SCDs) to detect soft X-ray photons. The time response of LE is caused by the structure of SCDs. With theoretical analysis and Monte Carlo simulations we discuss the influence of LE time response (LTR) on the timing analysis from three aspects: the power spectral density, the pulse profile and the time lag. After the LTR, the value of power spectral density monotonously decreases with the increasing frequency. The power spectral density of a sinusoidal signal reduces by a half at frequency 536 Hz. The corresponding frequency for QPO signals is 458 Hz. The Root mean square (RMS) of QPOs holds the similar behaviour. After the LTR, the centroid frequency and full width at half maxima (FWHM) of QPOs signals do not change. The LTR reduces the RMS of pulse profiles and shifts the pulse phase. In the time domain, the LTR only reduces the peak value of the crosscorrelation function while it does not change the peak position. Thus it will not affect the result of the time lag. When considering the time lag obtained from two instruments and one among them is LE, a 1.18 ms lag is expected caused by the LTR. The time lag calculated in the frequency domain is the same as that in the time domain.

preprint2020arXiv

Vision-Based Fall Event Detection in Complex Background Using Attention Guided Bi-directional LSTM

Fall event detection, as one of the greatest risks to the elderly, has been a hot research issue in the solitary scene in recent years. Nevertheless, there are few researches on the fall event detection in complex background. Different from most conventional background subtraction methods which depend on background modeling, Mask R-CNN method based on deep learning technique can clearly extract the moving object in noise background. We further propose an attention guided Bi-directional LSTM model for the final fall event detection. To demonstrate the efficiency, the proposed method is verified in the public dataset and self-build dataset. Evaluation of the algorithm performances in comparison with other state-of-the-art methods indicates that the proposed design is accurate and robust, which means it is suitable for the task of fall event detection in complex situation.

preprint2020arXiv

Woodpecker-DL: Accelerating Deep Neural Networks via Hardware-Aware Multifaceted Optimizations

Accelerating deep model training and inference is crucial in practice. Existing deep learning frameworks usually concentrate on optimizing training speed and pay fewer attentions to inference-specific optimizations. Actually, model inference differs from training in terms of computation, e.g. parameters are refreshed each gradient update step during training, but kept invariant during inference. These special characteristics of model inference open new opportunities for its optimization. In this paper, we propose a hardware-aware optimization framework, namely Woodpecker-DL (WPK), to accelerate inference by taking advantage of multiple joint optimizations from the perspectives of graph optimization, automated searches, domain-specific language (DSL) compiler techniques and system-level exploration. In WPK, we investigated two new automated search approaches based on genetic algorithm and reinforcement learning, respectively, to hunt the best operator code configurations targeting specific hardware. A customized DSL compiler is further attached to these search algorithms to generate efficient codes. To create an optimized inference plan, WPK systematically explores high-speed operator implementations from third-party libraries besides our automatically generated codes and singles out the best implementation per operator for use. Extensive experiments demonstrated that on a Tesla P100 GPU, we can achieve the maximum speedup of 5.40 over cuDNN and 1.63 over TVM on individual convolution operators, and run up to 1.18 times faster than TensorRT for end-to-end model inference.

preprint2019arXiv

Modulation instability, rogue waves and spectral analysis for the sixth-order nonlinear Schrodinger equation

Modulation instability, rogue wave and spectral analysis are investigated for the nonlinear Schrodinger equation with the higher-order terms. The modulation instability distribution characteristics from the sixth-order to the eighth-order nonlinear Schrodinger equations are studied. Higher-order dispersion terms are closely related to the distribution of modulation stability regime, and n-order dispersion term corresponds to $n-2$ modulation stability curves in the modulation instability band. Based on the generalized Darboux transformation method, the higher-order rational solutions are constructed. Then the compact algebraic expression of the Nth-order rogue wave is given. Dynamic phenomena of first- to third-order rogue waves are illustrated, which exhibit meaningful structures. Two arbitrary parameters play important roles in the rogue wave solution. One can control deflection of crest of rogue wave and its width, while the other can cause the change of width and amplitude of rogue wave. When it comes to the third-order rogue wave, three typical nonlinear wave constructions, namely fundamental, circular and triangular are displayed and discussed. Through the spectral analysis on first-order rogue wave, when these parameters satisfy certain conditions, it occurs a transition between W-shaped soliton and rogue wave.

preprint2019arXiv

Overview to the Hard X-ray Modulation Telescope (Insight-HXMT) Satellite

As China's first X-ray astronomical satellite, the Hard X-ray Modulation Telescope (HXMT), which was dubbed as Insight-HXMT after the launch on June 15, 2017, is a wide-band (1-250 keV) slat-collimator-based X-ray astronomy satellite with the capability of all-sky monitoring in 0.2-3 MeV. It was designed to perform pointing, scanning and gamma-ray burst (GRB) observations and, based on the Direct Demodulation Method (DDM), the image of the scanned sky region can be reconstructed. Here we give an overview of the mission and its progresses, including payload, core sciences, ground calibration/facility, ground segment, data archive, software, in-orbit performance, calibration, background model, observations and some preliminary results.

preprint2019arXiv

Photon orbits and phase transitions in Born-Infeld-dilaton black holes

Relations between photon orbits and thermodynamical phase transitions are explored in Born-Infeld-dilaton-AdS black hole. The coupling between the electromagnetic field and the dialton field is chosen such that the full phase diagram contains zeroth-order and first-order phase transitions as well as RPT. We find that there exist non-monotonic beahviors of the photon orbit radius $r_{ps}$ and the minimum impact parameter $u_{ps}$ which signal the existence of the various phase transitions. In particular, the marginal value of pressure under which RPT occur can be read off from these behaviors. Along the co-existing lines, there are changes of both $r_{ps}$ and $u_{ps}$, whose dependence on the transition temperature show characteristic behaviors signalling the existence of RPT. Moreover, the critical exponents of $Δr_{ps}$ and $Δu_{ps}$ are found to take a universal value $\frac{1}{2}$. These results imply that $r_{ps}$ and $u_{ps}$ can be used as order parameters to describe BH phase transitions.

preprint2018arXiv

Multi-dark soliton solutions for the multi-component coupled Maccari system

Based on the KP hierarchy reduction method, the general multi-dark soliton solutions in Gram type determinant forms for the multi-component coupled Maccari system are constructed. Especially, the three-component coupled Maccari system comprised of two-component short waves and single-component long waves are discussed in detail. Besides, the dynamics of one and two dark-dark solitons are analyzed. It is shown that the collisions of two dark-dark solitons are elastic by asymptotic analysis. Additionally, the two dark-dark solitons bound states are studied through two different cases (stationary and moving cases). The bound states can exist up to arbitrary order in the stationary case, however, only two-soliton bound state exists in the moving case. Besides, the oblique stationary bound state can be generated for all possible combinations of nonlinearity coefficients consisted of positive, negative and mixed cases. Nevertheless, the parallel stationary and the moving bound states are only possible when nonlinearity coefficients take opposite signs.

preprint2016arXiv

Chaoticons described by nonlocal nonlinear Schrodinger equation

It is shown that the unstable evolutions of the Hermite-Gauss-type stationary solutions for the nonlocal nonlinear Schrodinger equation with the exponential-decay response function can evolve into chaotic states. This new kind of entities are referred to as chaoticons because they exhibit not only chaotic properties (with positive Lyapunov exponents and spatial decoherence) but also soliton-like properties (with invariant statistic width and interaction of quasi-elastic collisions).

preprint2016arXiv

Luminescent and absorptive metal-coated emulsions for micro-velocimetry

Fluorescent latex beads have been widely used as tracers in microfluidics over for the last decades. They have the advantages to be density matched with water and to be easily localizable using fluorescence microscopy. We have recently synthesized silver-coated oil droplets that are both luminescent and absorptive, by first coating the oil interface with a polydopamine layer and then depositing a silver layer by a redox process. They have a mean diameter of 6 μm and a their density has been matched to the density of water by adjusting the thickness of the metallic layer. In this work we used these particles as tracers to measure the velocity profile of a water flow in a PDMS microchannel with a rectangular cross-section, which allowed us to confirm the predictions of the Stokes equation with results comparable to those of common submicronic polystyrene particles.

preprint2015arXiv

An integrable semi-discretization of the coupled Yajima--Oikawa system

In the present paper, an integrable semi-discrete analogue of the one-dimensional coupled Yajima--Oikawa system, which is comprised of multicomponent short-wave and one component long-wave, is proposed by using Hirota's bilinear method. Based on the reductions of the Bäcklund transformations of the semi-discrete BKP hierarchy, both the bright and dark soliton (for the short-wave components) solutions in terms of pfaffians are constructed.

preprint2015arXiv

Firing dynamics of an autaptic neuron

Autapses are synapses that connect a neuron to itself in the nervous system. Previously, both experimental and theoretical studies have demonstrated that autaptic connections in the nervous system have a significant physiological function. Autapses in nature provide self-delayed feedback, thus introducing an additional time scale to neuronal activities and causing many dynamic behaviors in neurons. Recently, theoretical studies have revealed that an autapse provides a control option for adjusting the response of a neuron: e.g., an autaptic connection can cause the electrical activities of the Hindmarsh-Rose neuron to switch between quiescent, periodic, and chaotic firing patterns; an autapse can enhance or suppress the mode-locking status of a neuron injected with sinusoidal current; and the firing frequency and interspike interval distributions of the response spike train can also be modified by the autapse. In this paper, we review recent studies that showed how an autapse affects the response of a single neuron.

preprint2015arXiv

General Mixed Multi-Soliton Solutions to One-Dimensional Multicomponent Yajima-Oikawa System

In this paper, we derive a general mixed (bright-dark) multi-soliton solution to a one-dimensional multicomponent Yajima-Oikawa (YO) system, i.e., the (M+1)-component YO system comprised of M-component short waves (SWs) and one-component long wave (LW) for all possible combinations of nonlinearity coefficients including positive, negative and mixed types. With the help of the KP-hierarchy reduction method, we firstly construct two types of general mixed N-soliton solution (two-bright-one-dark soliton and one-bright-two-dark one for SW components) to the (3+1)-component YO system in detail. Then by extending the corresponding analysis to the (M+1)-component YO system, a general mixed N-soliton solution in Gram determinant form is obtained. The expression of the mixed soliton solution also contains the general all bright and all dark N-soliton solution as special cases. Besides, the dynamical analysis shows that the inelastic collision can only take place among SW components when at least two SW components have bright solitons in mixed type soliton solution. Whereas, the dark solitons in SW components and the bright soliton in LW component always undergo usual elastic collision.

preprint2015arXiv

Integrable discretizations and self-adaptive moving mesh method for a coupled short pulse equation

In the present paper, integrable semi-discrete and fully discrete analogues of a coupled short pulse (CSP) equation are constructed. The key of the construction is the bilinear forms and determinant structure of solutions of the CSP equation. We also construct Nsoliton solutions for the semi-discrete and fully discrete analogues of the CSP equations in the form of Casorati determinant. In the continuous limit, we show that the fully discrete CSP equation converges to the semi-discrete CSP equation, then further to the continuous CSP equation. Moreover, the integrable semi-discretization of the CSP equation is used as a selfadaptive moving mesh method for numerical simulations. The numerical results agree with the analytical results very well.

preprint2015arXiv

Modulation instability and controllable rogue waves with multiple compression points for periodically modulated coupled Hirota equations

Based on modulation instability analysis and generalized Darboux transformation, we derive a hierarchy of rogue wave solutions for a variable-coefficients coupled Hirota equations. The explicit first-order rogue wave solution is presented, and the dark-bright and composite rogue waves with multiple compression points are shown by choosing sufficiently large periodic modulation amplitudes in the coefficients of the coupled equations. Also, the dark-bright and composite Peregrine combs are generated from the multiple compression points. For the second-order case, the darkbright three sisters, rogue wave quartets and sextets structures with one or more rogue waves involving multiple compression points are put forward, respectively. Furthermore, some wave characteristics such as the difference between light intensity and continuous wave background, and pulse energy evolution of the dark rogue wave solution features multiple compression points are discussed.

preprint2015arXiv

Multi-dark soliton solutions of the two-dimensional multi-component Yajima-Oikawa systems

We present a general form of multi-dark soliton solutions of two-dimensional multi-component soliton systems. Multi-dark soliton solutions of the two-dimensional (2D) and one-dimensional (1D) multi-component Yajima-Oikawa (YO) systems, which are often called the 2D and 1D multi-component long wave-short wave resonance interaction systems, are studied in detail. Taking the 2D coupled YO system with two short wave and one long wave components as an example, we derive the general $N$-dark-dark soliton solution in both the Gram type and Wronski type determinant forms for the 2D coupled YO system via the KP hierarchy reduction method. By imposing certain constraint conditions, the general $N$-dark-dark soliton solution of the 1D coupled YO system is further obtained. The dynamics of one dark-dark and two dark-dark solitons are analyzed in detail. In contrast with bright-bright soliton collisions, it is shown that dark-dark soliton collisions are elastic and there is no energy exchange among solitons in different components. Moreover, the dark-dark soliton bound states including the stationary and moving ones are discussed. For the stationary case, the bound states exist up to arbitrary order, whereas, for the moving case, only the two-soliton bound state is possible under the condition that the coefficients of nonlinear terms have opposite signs.

preprint2015arXiv

The Large Deviation Principle and Steady-state Fluctuation Theorem for the Entropy Production Rate of a Stochastic Process in Magnetic Fields

Fluctuation theorem is one of the major achievements in the field of nonequilibrium statistical mechanics during the past two decades. Steady-state fluctuation theorem of sample entropy production rate in terms of large deviation principle for diffusion processes have not been rigorously proved yet due to technical difficulties. Here we give a proof for the steady-state fluctuation theorem of a diffusion process in magnetic fields, with explicit expressions of the free energy function and rate function. The proof is based on the Karhunen-Loéve expansion of complex-valued Ornstein-Uhlenbeck process.

preprint2014arXiv

A digital CDS technique and the performance testing

Readout noise is a critical parameter for characterizing the performance of charge-coupled devices (CCDs), which can be greatly reduced by the correlated double sampling (CDS) circuit. However, conventional CDS circuit inevitably introduces new noises since it consists of several active analog components such as operational amplifiers. This paper proposes a digital CDS circuit technique, which transforms the pre-amplified CCD signal into a train of digital presentations by a high-speed data acquisition card directly without the noisy CDS circuit first, then implement the digital CDS algorithm through numerical method. The readout noise of 3.3 e$^{-}$ and the energy resolution of 121 eV@5.9keV can be achieved via the digital CDS technique.

preprint2014arXiv

A direct algorithm of one-dimensional optimal system for the group invariant solutions

A direct and systematic algorithm is proposed to find one-dimensional optimal system for the group invariant solutions, which is attributed to the classification of its corresponding one-dimensional Lie algebra. Since the method is based on different values of all the invariants, the process itself can both guarantee the comprehensiveness and demonstrate the inequivalence of the optimal system, with no further proof. To illustrate our method more clearly , we give a couple of well-known examples: the Korteweg-de Vries (KdV) equation and the heat equation.

preprint2014arXiv

Fluctuating Immersed Material (FIMAT) Dynamics for Fully Resolved Simulation of the Brownian Motion of Particles

Fluctuating hydrodynamics based techniques have been developed in recent years for the simulation of Brownian motion of particles. These mesoscale simulation tools are viable approaches for problems where molecular dynamics simulations may be deemed expensive. We have developed a rigid constraint-based formulation where the key idea is to assume that the entire domain is a fluctuating fluid. Rigid motion constraints are then imposed in regions that are occupied by rigid particles. The resulting solution gives the Brownian motion of the particles. This approach is shown to be viable for the simulation of long time scale diffusive behavior as well as for short time scale dynamics by using two separate solution techniques. Test cases are reported to validate the approach and to establish its efficacy.

preprint2014arXiv

Generalized Darboux transformation and higher-order rogue wave solutions of the coupled Hirota equations

This paper is dedicated to study higher-order rogue wave solutions of the coupled Hirota equations with high-order nonlinear effects like the third dispersion, self-steepening and stimulated Raman scattering terms. By using the generalized Darboux transformation, a unified representation of Nth-order rogue wave solution with 3N+1 free parameters is obtained. In particular, the first-order rogue wave solution containing polynomials of fourth order, and the second-order rogue wave solution consisting of polynomials of eighth order are explicitly presented. Through the numerical plots, we show that four or six fundamental rogue waves can coexist in the second-order rogue waves. By adjusting the values of some free parameters, different kinds of spatial-temporal distribution structures such as circular, quadrilateral, triangular, line and fundamental patterns are exhibited. Moreover, we see that nine or twelve fundamental rogue waves can synchronously emerge in the third-order rogue waves. The more intricate spatialtemporal distribution shapes are shown via adequate choices of the free parameters. Several wave characteristics such as the amplitudes and the coordinate positions of the highest peaks in the rogue waves are discussed.

preprint2014arXiv

Generalized Darboux transformation and localized waves in coupled Hirota equations

In this paper, we construct a generalized Darboux transformation to the coupled Hirota equations with high-order nonlinear effects like the third dispersion, self-steepening and inelastic Raman scattering terms. As application, an Nth-order localized wave solution on the plane backgrounds with the same spectral parameter is derived through the direct iterative rule. In particular, some semi-rational, multi-parametric localized wave solutions are obtained: (1) Vector generalization of the first- and the second-order rogue wave solution; (2) Interactional solutions between a dark-bright soliton and a rogue wave, two dark-bright solitons and a second-order rogue wave; (3) Interactional solutions between a breather and a rogue wave, two breathers and a second-order rogue wave. The results further reveal the striking dynamic structures of localized waves in complex coupled systems.

preprint2014arXiv

Proton irradiation effect on SCDs

The Low Energy X-ray Telescope is a main payload on the Hard X-ray Modulation Telescope satellite. The swept charge device is selected for the Low Energy X-ray Telescope. As swept charge devices are sensitive to proton irradiation, irradiation test was carried out on the HI-13 accelerator at the China Institute of Atomic Energy. The beam energy was measured to be 10 MeV at the SCD. The proton fluence delivered to the SCD was $3\times10^{8}\mathrm{protons}/\mathrm{cm}^{2}$ over two hours. It is concluded that the proton irradiation affects both the dark current and the charge transfer inefficiency of the SCD through comparing the performance both before and after the irradiation. The energy resolution of the proton-irradiated SCD is 212 eV@5.9 keV at $-60\,^{\circ}\mathrm{C}$, while it before irradiated is 134 eV. Moreover, better performance can be reached by lowering the operating temperature of the SCD on orbit.

preprint2014arXiv

Rational solutions to multicomponent Yajima-Oikawa systems: from two dimension to one dimension

Exact explicit rational solutions of two- and one- dimensional multicomponent Yajima-Oikawa (YO) systems, which contain multi-short-wave components and single long-wave one, are presented by using the bilinear method. For two-dimensional system, the fundamental rational solution first describes the localized lumps, which have three different patterns: bright, intermediate and dark states. Then, rogue waves can be obtained under certain parameter conditions and their behaviors are also classified to above three patterns with different definition. It is shown that the simplest (fundamental) rogue waves are line localized waves which arise from the constant background with a line profile and then disappear into the constant background again. In particular, two-dimensional intermediate and dark counterparts of rogue wave are found with the different parameter requirements. We demonstrate that multirogue waves describe the interaction of several fundamental rogue waves, in which interesting curvy wave patterns appear in the intermediate times. Different curvy wave patterns form in the interaction of different types fundamental rogue waves. Higher-order rogue waves exhibit the dynamic behaviors that the wave structures start from lump and then retreats back to it, and this transient wave possesses the patterns such as parabolas. Furthermore, different states of higher-order rogue wave result in completely distinguishing lumps and parabolas. Moreover, one-dimensional rogue wave solutions with three states are constructed through the further reduction. Specifically, higher-order rogue wave in one dimensional case is derived under the parameter constraints.

preprint2014arXiv

The construction of two-dimensional optimal systems for the invariant solutions

To search for inequivalent group invariant solutions, a general and systematic approach is established to construct two-dimensional optimal systems, which is based on commutator relations, adjoint matrix and the invariants. The details of computing all the invariants for two-dimensional subalgebras is presented and the optimality of twodimensional optimal systems is shown clearly under different values of invariants, with no further proof. Applying the algorithm to (1+1)-dimensional heat equation and (2+1)-dimensional Navier-Stokes (NS) equation, their twodimensional optimal systems are obtained, respectively. For the heat equation, eleven two-parameter elements in the optimal system are found one by one, which are discovered more comprehensive. The two-dimensional optimal system of NS equations is used to generate intrinsically different reduced ordinary differential equations and some interesting explicit solutions are provided.

preprint2014arXiv

Unbiased Correction Relations for Galaxy Cluster properties Derived from Chandra and XMM-Newton

We use a sample of 62 clusters of galaxies to investigate the discrepancies of gas temperature and total mass within r500 between XMM-Newton and Chandra data. Comparisons of the properties show that: (1) Both the de-projected and projected temperatures determined by Chandra are higher than those of XMM-Newton and there is a good linear relation for the de-projected temperature. (2) The Chandra mass is much higher than XMM-Newton mass with a bias of 0.15. To explore the reasons for the discrepancy in mass, we recalculate the Chandra mass (expressed as M_c) by modifying its temperature with the de-projected temperature relation. The results show that M_c is more close to the XMM-Newton mass with the bias reducing to 0.02. Moreover, M_c are corrected with the r500 measured by XMM-Newton and the intrinsic scatter is significantly improved with the value reducing from 0.20 to 0.12. These mean that the temperature bias may be the main factor causing the mass bias. At last, we find that M_c is consistent with the corresponding XMM-Newton mass derived directly from our mass relation at a given Chandra mass. Thus, the de-projected temperature and mass relations can provide unbiased corrections for galaxy cluster properties derived from Chandra and XMM-Newton.

preprint2013arXiv

A good mass proxy for galaxy clusters with XMM-Newton

We use a sample of 39 galaxy clusters at redshift z < 0.1 observed by XMM-Newton to investigate the relations between X-ray observables and total mass. Based on central cooling time and central temperature drop, the clusters in this sample are divided into two groups: 25 cool core clusters and 14 non-cool core clusters, respectively. We study the scaling relations of Lbol-M500, M500-T, M500-Mg and M500-YX, and also the influences of cool core on these relations. The results show that the M500-YX relation has a slope close to the standard self-similar value, has the smallest scatter and does not vary with the cluster sample. Moreover, the M500-YX relation is not affected by the cool core. Thus, the parameter of YX may be the best mass indicator.

preprint2013arXiv

A new integrable discrete generalized nonlinear Schrodinger equation and its reductions

A new integrable discrete system is constructed and studied, based on the algebraization of the difference operator. The model is named the discrete generalized nonlinear Schrodinger (GNLS) equation for which can be reduced to classical discrete nonlinear Schrodinger (NLS) equation. To show the complete integrability of the discrete GNLS equation, the recursion operator, symmetries and conservation quantities are obtained. Furthermore, all of reductions for the discrete GNLS equation are given and the discrete NLS equation is obtained by one of the reductions. At the same time, the recursion operator and symmetries of continuous GNLS equation are successfully recovered by its corresponding discrete ones.

preprint2013arXiv

Measure synchronization in a two-species bosonic Josephson junction

Measure synchronization (MS) in a two-species bosonic Josephson junction (BJJ) is studied based on semi-classical theory. Six different scenarios for MS, including two in the Josephson oscillation regime (0 phase mode) and four in the self-trapping regime ($π$ phase mode), have been clearly shown. Systematic investigations of the common features behind these different scenarios have been performed. We show that the average energies of the two species merge at the MS transition point. The scaling of the power law near the MS transition has been verified, and the critical exponent is 1/2 for all of the different scenarios for MS. We also illustrate MS in a three-dimensional phase space; from this illustration, more detailed information on the dynamical process can be obtained. Particularly, by analyzing the Poincare sections with changing interspecies interactions, we find that the two-species BJJ exhibits separatrix crossing behavior at MS transition point, and such behavior depicts the general mechanism behind the different scenarios for the MS transitions. The new critical behavior found in a two-species BJJ is expected to be found in real systems of atomic Bose gases.

preprint2013arXiv

Non-local symmetries of the Hirota-Satsuma coupled KdV system and their applications

The nonlocal symmetry is derived from the known Darboux transformation (DT) of the Hirota-Satsuma coupled KdV (HS-cKdV) system, and infinitely many nonlocal symmetries are obtained by introducing some internal parameters. By extending the HS-cKdV system to an auxiliary system with five dependent variables, the prolongation is found to localize the nonlocal symmetry related to the DT. Base on the enlarged system, the finite symmetry transformations and similarity reductions about the local symmetries are computed, which lead to some novel exact solutions of the HS-cKdV system. These solutions contain some new solutions from old ones by the finite symmetry transformation and exact interaction solutions among solitons and other complicated waves including periodic cnoidal waves and Painlevé waves through similarity reductions. Some integrable models from the nonlocal symmetry related to the DT are obtained in two aspects: the negative HS-cKdV hierarchy obtained by introducing the internal parameter and integrable models both in lower and higher dimensions given by restricting the symmetry constraints.

preprint2013arXiv

On a nonsymmetric Ornstein-Uhlenbeck semigroup and its generator

If we add a simple rotation term to both the Ornstein-Uhlenbeck semigroup and the definition of the H-derivative, then analogue to the classical Malliavin calculus on the real Wiener space [I. Shigekawa, Stochastic analysis, 2004], we get a normal but nonsymmetric Ornstein-Uhlenbeck operator $L$ on the complex Wiener space. The eigenfunctions of the operator $L$ are given. In addition, the hypercontractivity for the nonsymmetric Ornstein-Uhlenbeck semigroup is shown.

preprint2013arXiv

On the Jordan decomposition for a class of non-symmetric Ornstein-Uhlenbeck operators

In this paper, we calculate the Jordan decomposition (or say, the Jordan canonical form) for a class of non-symmetric Ornstein-Uhlenbeck operators with the drift coefficient matrix being a Jordan block and the diffusion coefficient matrix being identity multiplying a constant. For the 2-dimensional case, we present all the general eigenfunctions by the induction. For the 3-dimensional case, we divide the calculating of the Jordan decomposition into several steps (the key step is to do the canonical projection onto the homogeneous Hermite polynomials, next we use the theory of systems of linear equations). As a by-pass product, we get the geometric multiplicity of the eigenvalue of the Ornstein-Uhlenbeck operator.

preprint2013arXiv

Reducing Poisson noise and baseline drift in X-ray spectral images with bootstrap Poisson regression and robust nonparametric regression

X-ray spectral imaging provides quantitative imaging of trace elements in biological sample with high sensitivity. We propose a novel algorithm to promote the signal-to-noise ratio (SNR) of X-ray spectral images that have low photon counts. Firstly, we estimate the image data area that belongs to the homogeneous parts through confidence interval testing. Then, we apply the Poisson regression through its maximum likelihood estimation on this area to estimate the true photon counts from the Poisson noise corrupted data. Unlike other denoising methods based on regression analysis, we use the bootstrap resampling method to ensure the accuracy of regression estimation. Finally, we use a robust local nonparametric regression method to estimate the baseline and subsequently subtract it from the X-ray spectral data to further improve the SNR of the data. Experiments on several real samples show that the proposed method performs better than some state-of-the-art approaches to ensure accuracy and precision for quantitative analysis of the different trace elements in standard reference biological sample.

preprint2013arXiv

Rogue wave solutions in AB system

In this paper, the generalized Darboux transformation is established to the AB system, which mainly describes marginally unstable baroclinic wave packets in geophysical fluids and ultra-short pulses in nonlinear optics. We find a unified formula of Nth-order rogue wave solution for the AB system by the direct iterative rule. In particular, rogue waves possessing several free parameters from first to second order are calculated. The dynamics and some interesting structures of the rogue waves are illustrated through some figures.

preprint2012arXiv

AFM and Raman studies of topological insulator materials subject to argon plasma etching

Plasma etching is an important tool in nano-device fabrication. We report a study on argon plasma etching of topological insulator materials Bi2Se3, Bi2Te3, Sb2Te3 and Bi2Te2Se using exfoliated flakes (with starting thicknesses of ~100 nm) derived from bulk crystals. We present data mainly from atomic force microscopy (AFM) and Raman spectroscopy. Through AFM measurements, plasma exposure is observed to decrease the thickness of our samples and increase surface roughness (with height fluctuations reaching as large as ~20 nm). We extract an etching rate for each type of material. Plasma exposure also causes a widening (especially Eg2) of the characteristic Raman peaks, with no significant change in peak position. The overall Raman intensity is observed to initially increase, then decrease sharply after the samples are etched below ~20 nm in thickness. Our findings are valuable for understanding the effects of argon plasma etching on topological insulator materials.

preprint2012arXiv

Deprojected analysis of Abell 1835 observed with Chandra and compared with XMM-Newton

Using a deprojection technique, we study the X-ray properties of the galaxy cluster Abell 1835 observed with Chandra, including temperature, abundance, electron density, gas mass fraction, and total mass. A comparison with the results without deprojection shows that the properties do not change much. When we compare the results with those of XMM-Newton, the difference between the temperature profiles derived from Chandra and XMM-Newton data still exists, even if the point-spread function effect of XMM-Newton is corrected. To investigate the reasons for the difference, we used the spectra to cross-calibrate the temperatures. They show that the Chandra spectra can be fitted well with XMM-Newton temperatures. Furthermore, we derive the electron density profile from Chandra data with XMM-Newton temperatures and calculate the projected mass, which is consistent with the XMM-Newton mass and a little lower than the weak lensing mass at r_500. Thus, it seems that the temperature derived from XMM-Newton may be more reliable.

preprint2012arXiv

Multiscale Theory of Finite Size Bose Systems: Implications for Collective and Single-Particle Excitations

Boson droplets (i.e., dense assemblies of bosons at low temperature) are shown to mask a significant amount of single-particle behavior and to manifest collective, droplet-wide excitations. To investigate the balance between single-particle and collective behavior, solutions to the wave equation for a finite size Bose system are constructed in the limit where the ratio \varepsilon of the average nearest-neighbor boson distance to the size of the droplet or the wavelength of density disturbances is small. In this limit, the lowest order wave function varies smoothly across the system, i.e., is devoid of structure on the scale of the average nearest-neighbor distance. The amplitude of short range structure in the wave function is shown to vanish as a power of \varepsilon when the interatomic forces are relatively weak. However, there is residual short range structure that increases with the strength of interatomic forces. While the multiscale approach is applied to boson droplets, the methodology is applicable to any finite size bose system and is shown to be more direct than field theoretic methods. Conclusions for Helium-4 nanodroplets are drawn.

preprint2012arXiv

Nonlocal symmetries related to Bäcklund transformation and their applications

Starting from nonlocal symmetries related to Bäcklund transformation (BT), many interesting results can be obtained. Taking the well known potential KdV (pKdV) equation as an example, a new type of nonlocal symmetry in elegant and compact form which comes from BT is presented and used to make researches in the following three subjects: two sets of negative pKdV hierarchies and their corresponding bilinear forms are constructed; the nonlocal symmetry is localized by introduction of suitable and simple auxiliary dependent variables to generate new solutions from old ones and to consider some novel group invariant solutions; some other models both in finite dimensions and infinite dimensions are generated by comprising the original BT and evolution under new nonlocal symmetry. The finite-dimensional models are completely integrable in Liouville sense, which are shown equivalent to the results given through the nonlinearization method for Lax pair.

preprint2012arXiv

On the eigenfunctions of the complex Ornstein-Uhlenbeck operators

Starting from the 1-dimensional complex-valued Ornstein-Uhlenbeck process, we present two natural ways to imply the associated eigenfunctions of the 2-dimensional normal Ornstein-Uhlenbeck operators in the complex Hilbert space $L_{\Cnum}^2(μ)$. We call the eigenfunctions Hermite-Laguerre-Ito polynomials. In addition, the Mehler summation formula for the complex process are shown.

preprint2010arXiv

A Heuristic Algorithm for optimizing Page Selection Instructions

Page switching is a technique that increases the memory in microcontrollers without extending the address buses. This technique is widely used in the design of 8-bit MCUs. In this paper, we present an algorithm to reduce the overhead of page switching. To pursue small code size, we place the emphasis on the allocation of functions into suitable pages with a heuristic algorithm, thereby the cost-effective placement of page selection instructions. Our experimental results showed the optimization achieved a reduction in code size of 13.2 percent.

preprint2010arXiv

Layered Depth-Normal Images: a Sparse Implicit Representation of Solid Models

This paper presents a novel implicit representation of solid models. With this representation, every solid model can be effectively presented by three layered depth-normal images (LDNIs) that are perpendicular to three orthogonal axes respectively. The layered depth-normal images for a solid model, whose boundary is presented by a polygonal mesh, can be generated efficiently with help of the graphics hardware accelerated sampling. Based on this implicit representation - LDNIs, solid modeling operations including the Boolean operations and the offsetting operation have been developed. A contouring algorithm is also introduced in this paper to generate thin structure and sharp feature preserved mesh surfaces from the layered depth-normal images. Comparisons between LDNIs and other implicit representation of solid models are given at the end of the paper to demonstrate the advantages of LDNIs.

preprint2010arXiv

On the Imbedding Problem for Three-state Time Homogeneous Markov Chains with Coinciding Negative Eigenvalues

For an indecomposable $3\times 3$ stochastic matrix (i.e., 1-step transition probability matrix) with coinciding negative eigenvalues, a new necessary and sufficient condition of the imbedding problem for time homogeneous Markov chains is shown by means of an alternate parameterization of the transition rate matrix (i.e., intensity matrix, infinitesimal generator), which avoids calculating matrix logarithm or matrix square root. In addition, an implicit description of the imbedding problem for the $3\times 3$ stochastic matrix in Johansen [J. Lond. Math. Soc., 8, 345-351. (1974)] is pointed out.

preprint2009arXiv

Stability and chaotic behavior of Bose-Einstein condensates in optical lattices with two- and three-body interactions

The stability and chaotic behaviors of Bose-Einstein condensates with two- and three-atom interactions in optical lattices are discussed with analytical and numerical methods. It is found that the steady-state relative population appears tuning-fork bifurcation when the system parameters are changed to certain critical values. In particular, the existence of three-body interaction not only transforms the bifurcation point of the system but also affects greatly on the macroscopic quantum self-trapping behaviors of the system associated with the critically stable steady-state solution. In addition, we also investigated the influence of the initial conditions, three-body interaction and the energy bias on the macroscopic quantum self-trapping. Finally, by applying the periodic modulation on the energy bias, we find that the relative population oscillation exhibits a process from order to chaos, via a series of period-doubling bifurcations.

preprint2007arXiv

Transient Dynamics of Sparsely Connected Hopfield Neural Networks with Arbitrary Degree Distributions

Using probabilistic approach, the transient dynamics of sparsely connected Hopfield neural networks is studied for arbitrary degree distributions. A recursive scheme is developed to determine the time evolution of overlap parameters. As illustrative examples, the explicit calculations of dynamics for networks with binomial, power-law, and uniform degree distribution are performed. The results are good agreement with the extensive numerical simulations. It indicates that with the same average degree, there is a gradual improvement of network performance with increasing sharpness of its degree distribution, and the most efficient degree distribution for global storage of patterns is the delta function.

preprint2003arXiv

Wigner crystalization about $ν$=3

We measure a resonance in the frequency dependence of the real diagonal conductivity, Re[$σ_{xx}$], near integer filling factor, $ν=3$. This resonance depends strongly on $ν$, with peak frequency $f_{pk} \approx 1.7$ GHz at $ν=3.04$ or 2.92 close to integer $ν$, but $f_{pk} \approx$ 600 MHz at $ν=3.19$ or 2.82, the extremes of where the resonance is visible. The dependence of $f_{pk}$ upon $n^*$, the density of electrons in the partially filled level, is discussed and compared with similar measurments by Chen {\it et al.}\cite{yong} about $ν=1$ and 2. We interpret the resonance as due to a pinned Wigner crystal phase with density $n^*$ about the $ν=3$ state.

preprint2002arXiv

The attractors in sequence processing neural networks

The average length and average relaxation time of attractors in sequence processing neural networks are investigated. The simulation results show that a critical point of $α$, the loading ratio, is found. Below the turning point, the average length is equal to the number of stored patterns; conversely, the ratio of length and numbers of stored patterns, grow with an exponential dependence $\exp (Aα) $. Moreover, we find that the logarithm of average relaxation time is only linearly associated with $α$ and the turning point of coupling degree is located for examining robustness of networks.