Catalog footprint

What is connected

50works
25topics
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

50 published item(s)

preprint2026arXiv

Accelerated simulation of multiscale gas-radiation coupling flows via a general synthetic iterative scheme

Gas-radiation coupling critically influences hypersonic reentry flows, where extreme temperatures induce pronounced non-equilibrium gas and radiative heat transport. Accurate and efficient simulation of radiative gas dynamics is therefore indispensable for reliable design of thermal protection systems for atmospheric entry vehicles. In this study, a Boltzmann-type kinetic model for radiative gas flows is solved across a broad spectrum of flow and radiation transport regimes using the general synthetic iterative scheme (GSIS). The approach integrates an unstructured finite-volume discrete velocity method with a set of macroscopic synthetic equations. Within this framework, the kinetic model provides high-order closures for the constitutive relations in the synthetic equations. Simultaneously, the macroscopic synthetic equations drive the evolution of the mesoscopic kinetic system, significantly accelerating steady-state convergence in near-continuum regimes, as substantiated by linear Fourier stability analysis. Crucially, the algorithm is proven to be asymptotic-preserving, correctly recovering the continuum and optically thick limits, represented by the radiative Navier-Stokes-Fourier equations governing distinct translational, rotational, vibrational, and radiative temperatures, on coarse meshes independent of the mean free path. Numerical simulations of challenging benchmarks, including three-dimensional hypersonic flow over an Apollo reentry capsule, demonstrate that GSIS achieves orders-of-magnitude speedup over conventional iterative schemes in multiscale simulations of radiative gas flows while accurately capturing non-equilibrium effects and radiative heat transfer in hypersonic environments.

preprint2026arXiv

GI-Bench: A Panoramic Benchmark Revealing the Knowledge-Experience Dissociation of Multimodal Large Language Models in Gastrointestinal Endoscopy Against Clinical Standards

Multimodal Large Language Models (MLLMs) show promise in gastroenterology, yet their performance against comprehensive clinical workflows and human benchmarks remains unverified. To systematically evaluate state-of-the-art MLLMs across a panoramic gastrointestinal endoscopy workflow and determine their clinical utility compared with human endoscopists. We constructed GI-Bench, a benchmark encompassing 20 fine-grained lesion categories. Twelve MLLMs were evaluated across a five-stage clinical workflow: anatomical localization, lesion identification, diagnosis, findings description, and management. Model performance was benchmarked against three junior endoscopists and three residency trainees using Macro-F1, mean Intersection-over-Union (mIoU), and multi-dimensional Likert scale. Gemini-3-Pro achieved state-of-the-art performance. In diagnostic reasoning, top-tier models (Macro-F1 0.641) outperformed trainees (0.492) and rivaled junior endoscopists (0.727; p>0.05). However, a critical "spatial grounding bottleneck" persisted; human lesion localization (mIoU >0.506) significantly outperformed the best model (0.345; p<0.05). Furthermore, qualitative analysis revealed a "fluency-accuracy paradox": models generated reports with superior linguistic readability compared with humans (p<0.05) but exhibited significantly lower factual correctness (p<0.05) due to "over-interpretation" and hallucination of visual features. GI-Bench maintains a dynamic leaderboard that tracks the evolving performance of MLLMs in clinical endoscopy. The current rankings and benchmark results are available at https://roterdl.github.io/GIBench/.

preprint2025arXiv

Probing compressed Higgsinos at the FASER experiment

In the Minimal Supersymmetric Standard Model (MSSM), compressed Higgsinos spectrum ($Δm^0 \lesssim 1$ GeV) occurs when $|μ| \ll |M_1|, |M_2|$ and ${\rm sign}(M_1\cdot M_2)<0$, which leads to a long-lived next-to-lightest neutralino. Such a long-lived neutralino could be copiously produced at the LHC, however escape the detection at the LHC main detectors. We examine the discovery potential at the FASER experiment and find that the FASER 2 could cover the neutral Higgsino mass up to about 130 GeV with mass splitting between 4 to 30 MeV. It is complementary to both the LHC Higgsino search in the $Δm^{0,\pm} \gtrsim 1$ GeV region, and displaced vertex and disappearing track searches of charginos with $Δm^\pm \lesssim 1$ GeV.

preprint2022arXiv

An Entropy-based Measure of Intelligence Degree of System Structures

In this paper, we investigate how to measure the intelligence of systems under specific structures. Two indicators are adopted to characterize the intelligence of a given structure, namely the function diversity of the structure, and the ability to generate order under specific environments. A measure of intelligence degree is proposed, with which the intelligence degree of several basic structures is calculated. It is shown that some structures are indeed "smarter" than the others under the proposed measure. The results add a possible way of revealing the evolution mechanism of natural life and constructing life-like structures with high intelligence degree.

preprint2022arXiv

BSM Higgs Production at a Muon Collider

The potential of the non-Standard Model heavy Higgs bosons in 2HDM at a muon collider is studied. The pair production of the non-SM Higgs bosons via the universal gauge interactions is the dominant mechanism once above the kinematic threshold. On the other hand, single Higgs boson production associated with a pair of heavy fermions is also important in the parameter region with enhanced Yukawa couplings. Both $μ^+μ^-$ annihilation channels and Vector Boson Fusion processes are considered, as well as radiative return $s$-channel production. Different types of 2HDMs can also be distinguishable for moderate and large values of $\tanβ$.

preprint2022arXiv

Exotic Higgs Decays in the Type-II 2HDMs at Current and Future pp Colliders

The exotic decay modes of non-Standard Model Higgses can serve as powerful search channels to explore the parameter space of extended Higgs sectors. In this Snowmass contribution, we illustrate this using the Two-Higgs Doublet Model (2HDM) as an example. We first review the current experimental constraints on the parameter space of a Type-II 2HDM arising from existing searches for the exotic Higgs decay mode $A/H\rightarrow HZ/AZ$. We then present the sensitivity of future colliders to discover addition Higgs bosons using the exotic decay channels $A\rightarrow HZ$, $A\rightarrow H^\pm W^\mp$ and $H^\pm\rightarrow H W^\pm$. We find that a 100 TeV collider can probe almost the entire region of the Type-II 2HDM parameter space that survives current theoretical and experimental constraints and would therefore be an ideal machine to search for heavier Higgses in hierarchical scalar sectors.

preprint2022arXiv

Global fit of 2HDM with future collider results

In this work, we summarize a global fit study of Type-II two Higgs doublet models (2HDM), and explore the impact of future SM-like Higgs and Z-pole precision measurements on the allowed parameter space. The work is based on the study results of a global fit of 2HDMs with the tool GAMBIT, utilising various current constraints including theoretical constraints (unitarity, perturbativity and vacuum stability), Higgs searches at colliders, electroweak physics and flavour constraints. We further investigate the ability of future facilities, such as the HL-LHC, CEPC, ILC and FCC-ee to explore the 2HDM parameter space.

preprint2022arXiv

Global fits of SUSY at future Higgs factories

In this work, we study the impact of electroweak and Higgs precision measurements at future electron-positron colliders on several typical supersymmetric models, including the Constrained Minimal Supersymmetric Standard Model (CMSSM), Non-Universal Higgs Mass generalisations (NUHM1, NUHM2), and the 7-dimensional Minimal Supersymmetric Standard Model (MSSM7). Using publicly-available data from the \textsf{GAMBIT} community, we post-process previous SUSY global fits with additional likelihoods to explore the discovery potential of Higgs factories, such as the Circular Electron Positron Collider (CEPC), the Future Circular Collider (FCC) and the International Linear Collider (ILC). We show that the currently allowed parameter space of these models will be further tested by future precision measurements. In particular, dark matter annihilation mechanisms may be distinguished by precise measurements of Higgs observables.

preprint2022arXiv

Heavy Higgs Bosons in 2HDM at a Muon Collider

We study the discovery potential of the non-Standard Model (SM) heavy Higgs bosons in the Two-Higgs-Doublet Models (2HDMs) at a multi-TeV muon collider and explore the discrimination power among different types of 2HDMs. We find that the pair production of the non-SM Higgs bosons via the universal gauge interactions is the dominant mechanism once above the kinematic threshold. Single Higgs boson production associated with a pair of heavy fermions could be important in the parameter region with enhanced Yukawa couplings. For both signal final states, $μ^+μ^-$ annihilation channels dominate over the vector boson fusion (VBF) processes, except at high center of mass energies where the VBF processes receive large logarithmic enhancement with the increase of energies. Single Higgs boson $s$-channel production in $μ^+μ^-$-annihilation via the radiative return can also be important for the Type-L 2HDM in the very large $\tanβ$ region, extending the kinematic reach of the heavy Higgs boson mass to the collider energy. Considering both the production and decay of non-SM Higgs bosons, signals can be identified over the Standard Model backgrounds. With the pair production channels via annihilation, 95\% C.L. exclusion reaches in the Higgs mass up to the production mass threshold of $\sqrt{s}/2$ are possible when channels with different final states are combined. Including single production modes can extended the reach further. Different types of 2HDMs can be distinguishable for moderate and large values of $\tanβ$.

preprint2022arXiv

MSSM Under Higgs Factories

The high precision measurements of the Higgs mass and couplings at the future Higgs factories are sensitive to the parameter space of the Minimal Supersymmetric Standard Model (MSSM). Focused on the dominant stop sector contributions, we study the implication of the Higgs precision measurements on MSSM using multi-variable $χ^2$ fit. The results show nice complementarity between the indirect searches at Higgs factories and the direct searches at the current LHC program.

preprint2022arXiv

Simple and statistically sound recommendations for analysing physical theories

Physical theories that depend on many parameters or are tested against data from many different experiments pose unique challenges to statistical inference. Many models in particle physics, astrophysics and cosmology fall into one or both of these categories. These issues are often sidestepped with statistically unsound ad hoc methods, involving intersection of parameter intervals estimated by multiple experiments, and random or grid sampling of model parameters. Whilst these methods are easy to apply, they exhibit pathologies even in low-dimensional parameter spaces, and quickly become problematic to use and interpret in higher dimensions. In this article we give clear guidance for going beyond these procedures, suggesting where possible simple methods for performing statistically sound inference, and recommendations of readily-available software tools and standards that can assist in doing so. Our aim is to provide any physicists lacking comprehensive statistical training with recommendations for reaching correct scientific conclusions, with only a modest increase in analysis burden. Our examples can be reproduced with the code publicly available at https://doi.org/10.5281/zenodo.4322283.

preprint2022arXiv

The physics case of a 3 TeV muon collider stage

In the path towards a muon collider with center of mass energy of 10 TeV or more, a stage at 3 TeV emerges as an appealing option. Reviewing the physics potential of such muon collider is the main purpose of this document. In order to outline the progression of the physics performances across the stages, a few sensitivity projections for higher energy are also presented. There are many opportunities for probing new physics at a 3 TeV muon collider. Some of them are in common with the extensively documented physics case of the CLIC 3 TeV energy stage, and include measuring the Higgs trilinear coupling and testing the possible composite nature of the Higgs boson and of the top quark at the 20 TeV scale. Other opportunities are unique of a 3 TeV muon collider, and stem from the fact that muons are collided rather than electrons. This is exemplified by studying the potential to explore the microscopic origin of the current $g$-2 and $B$-physics anomalies, which are both related with muons.

preprint2022arXiv

Unique quantum impurity states driven by a vortex in topological superconductors

The interplay of magnetic impurity and vortex in a topological superconductor is of fundamental interest with major implications for implementing quantum computation. There are multiple degrees of freedom interacting with the impurity in the system, including the Majorana zero mode (MZM), the Caroli de Gennes Matricon (CdGM) states, and the electron bulk states that form Cooper pairs, which makes the impurity pinned vortex state elusive to date in topological superconductors. Here, we present an accurate solution of the problem, based on a generalized mapping scheme and the density-matrix renormalization group (DMRG) method. We identify in-gap states that are distinct from the established Yu-Shiba-Rusinov (YSR) states. The newly found in-gap physics is driven by three prominent mechanisms: (i) the coupling of impurity and bulk states leading to competition between Kondo screening and Cooper pairing, resulting in a singlet-doublet quantum phase transition, (ii) the coupling of impurity and CdGM states introducing an effective Zeeman splitting to the doublet state, and (iii) the coupling of impurity and MZM turning the singlet-doublet transition into a crossover. These mechanisms cooperatively produce a unique spin-resolved local density of states. Despite the MZM, a robust nearly zero-energy peak is generated, with comparable height but opposite spin polarization as that of the MZM. These results signify novel emergent physics and offer insights to elucidate intriguing experimental phenomena.

preprint2021arXiv

Convergent Bayesian Global Fits of 4D Composite Higgs Models

Models in which the Higgs boson is a composite pseudo-Nambu-Goldstone boson offer attractive solutions to the Higgs mass naturalness problem. We consider three such models based on the minimal $SO(5) \rightarrow SO(4)$ symmetry breaking pattern, and perform convergent global fits on the models under a Bayesian framework in order to find the regions of their parameter spaces that best fit a wide range of constraints, including recent Higgs measurements. We use a novel technique to analyse the fine-tuning of the models, quantifying the tuning as the Kullback-Leibler divergence from the prior to the posterior probability on the parameter space. Each model is found to be able to satisfy all constraints at the $3σ$ level simultaneously. As a by-product of the fits, we analyse the collider phenomenology of our models in these viable regions. In two of the three models, we find that the $g g \rightarrow H \rightarrow γγ$ cross section is less than ${\sim}90$% that predicted by the SM, which is already in slight tension with experiment and could potentially be ruled out in the future high-luminosity run of the LHC. In addition, the lightest fermions $F$ arising from the new strong dynamics in these models are seen in general to lie above ${\sim}1.1$ TeV, with the $F \rightarrow tW^{+}$ and $F \rightarrow \bar{b}W^{+}$ decays offering particularly promising channels for probing these models in future collider searches.

preprint2020arXiv

An ultra-long and quite thin coronal loop without significant expansion

Context. Coronal loops are the basic building blocks of the solar corona, which are related to the mass supply and heating of solar plasmas in the corona. However, their fundamental magnetic structures are still not well understood. Most coronal loops do not expand significantly, whereas the diverging magnetic field would have an expansion factor of about 5-10 over one pressure scale height. Aims. In this study, we investigate a unique coronal loop with a roughly constant cross section, it is ultra long and quite thin. A coronal loop model with magnetic helicity is presented to explain the small expansion of the loop width. Methods. This coronal loop was predominantly detectable in the 171 A channel of the Atmospheric Imaging Assembly (AIA). Then, the local magnetic field line was extrapolated by a Potential-Field-Source-Surface model. Finally, the differential emission measure analysis made from six AIA bandpasses was applied to obtain the thermal properties of this loop. Results. This coronal loop has a projected length of roughly 130 Mm, a width of about 1.5 +(-) 0.5 Mm and a lifetime of around 90 minutes. It follows an open magnetic field line. The cross section expanded very little (i.e., 1.5-2.0) along the loop length during its whole lifetime. This loop has a nearly constant temperature at about 0.7 +(-) 0.2 MK, whereas its density exhibits the typical structure of a stratified atmosphere. Conclusions. We use a thin twisted flux tube theory to construct a model for this non-expanding loop, and find that indeed with sufficient twist a coronal loop can attain equilibrium. However, we can not rule out other possibilities such as footpoint heating by small-scale reconnection, elevated scale height by a steady flow along the loop etc.

preprint2020arXiv

Analyses of residual accelerations for TianQin based on the global MHD simulation

TianQin is a proposed space-based gravitational wave observatory. It is designed to detect the gravitational wave signals in the frequency range of 0.1 mHz -- 1 Hz. At a geocentric distance of $10^5$ km, the plasma in the earth magnetosphere will contribute as the main source of environmental noises. Here, we analyze the acceleration noises that are caused by the magnetic field of space plasma for the test mass of TianQin. The real solar wind data observed by the Advanced Composition Explorer are taken as the input of the magnetohydrodynamic simulation. The Space Weather Modeling Framework is used to simulate the global magnetosphere of the earth, from which we obtain the plasma and magnetic field parameters on the detector's orbits. We calculate the time series of the residual accelerations and the corresponding amplitude spectral densities on these orbit configurations. We find that the residual acceleration produced by the interaction between the TM's magnetic moment induced by the space magnetic field and the spacecraft magnetic field ($\bm{a}_{\rm M1}$) is the dominant term, which can approach $10^{-15}$ m/s$^2$/Hz$^{1/2}$ at $f \approx$ 0.2 mHz for the nominal values of the magnetic susceptibility ($χ_{\rm m} = 10^{-5}$) and the magnetic shielding factor ($ξ_{\rm m} = 10$) of the test mass. The ratios between the amplitude spectral density of the acceleration noise caused by the space magnetic field and the preliminary goal of the inertial sensor are 0.38 and 0.08 at 1 mHz and 10 mHz, respectively. We discuss the further reduction of this acceleration noise by decreasing $χ_{\rm m}$ and/or increasing $ξ_{\rm m}$ in the future instrumentation development for TianQin.

preprint2020arXiv

Comparative Studies of 2HDMs under the Higgs Boson Precision Measurements

We perform comparative studies for four types of the two Higgs Doublet Models (2HDMs) under the precision measurements of the Standard Model (SM) Higgs observables at the proposed Higgs factories. We explore the discovery potential based on the hypothetical deviations in the precision data for the 2HDMs up to one-loop level. We find $5σ$ observability from the $χ^2$ fitting in a significant theory parameter space at future Higgs factories. For the Type-I 2HDM, regions with $\cos(β-α)\lesssim -0.1$ or $\cos(β-α)\gtrsim 0.08$ are discoverable at more than $5σ$ level. For the other three types of 2HDMs, the $5σ$ region is even bigger: $|\cos(β-α)|\gtrsim 0.02$ for $\tanβ\sim 1$. At small and large values of $\tanβ$, the region in $\cos(β-α)$ is further tightened. We examine the extent to which the different 2HDM theories may be distinguishable from one to the other at the $95\%$ Confidence Level with four benchmark points as case studies. We show that a large part of the parameter space of the other types of 2HDMs can be distinguished from the benchmark points of the target model. The impacts of loop corrections are found to be significant in certain parameter regions.

preprint2020arXiv

Design and Practice of the Regulation in Speed with Flywheel

This paper aims at the teaching contents of "Fluctuation and Regulation in Speed of Machines" for students, explores the effect of the flywheel on the speed regulation of the mechanism system and its influencing factors, designs an experimental device with applied value, and develops a comparatively complete experimental scheme by contrasting experimental principles. Using this experimental device, we can in-depth study of the effect of the flywheel on the speed of different drive shafts, and deepen students' understanding and mastery of knowledge about regulation in speed.

preprint2020arXiv

Exploring the 2HDM with Global Fits in $\texttt{GAMBIT}$

In this work, we present preliminary results of a global fit of the type-II two-Higgs-doublet model (2HDM) with the tool $\texttt{GAMBIT}$. Our study includes various constraints, including the theoretical constraints (unitarity, perturbativity and vacuum stability), Higgs searches at colliders, electroweak physics and flavour constraints. With the latest experimental results, our results not only confirm past studies but also go further in probing the model. We find, for example, that the measurements of $B\rightarrow K^*μ^+μ^-$ angular observables cannot be explained in the type-II 2HDM.

preprint2020arXiv

Fast convergence and asymptotic preserving of the General Synthetic Iterative Scheme

Recently the general synthetic iteration scheme (GSIS) is proposed to find the steady-state solution of the Boltzmann equation~\cite{SuArXiv2019}, where various numerical simulations have shown that (i) the steady-state solution can be found within dozens of iterations at any Knudsen number $K$, and (ii) the solution is accurate even when the spatial cell size in the bulk region is much larger than the molecular mean free path, i.e. Navier-Stokes solutions are recovered at coarse grids. The first property indicates that the error decay rate between two consecutive iterations decreases to zero with $K$, while the second one implies that the GSIS is asymptotically preserving the Navier-Stokes limit. This paper is dedicated to the rigorous proof of both properties.

preprint2020arXiv

General synthetic iteration scheme for non-linear gas kinetic simulation of multi-scale rarefied gas flows

The general synthetic iteration scheme (GSIS) is extended to find the steady-state solution of nonlinear gas kinetic equation, removing the long-standing problems of slow convergence and requirement of ultra-fine grids in near-continuum flows. The key ingredients of GSIS are that the gas kinetic equation and macroscopic synthetic equations are tightly coupled, and the constitutive relations in macroscopic synthetic equations explicitly contain Newton's law of shear stress and Fourier's law of heat conduction. The higher-order constitutive relations describing rarefaction effects are calculated from the velocity distribution function, however, their constructions are simpler than our previous work (Su et al. Journal of Computational Physics 407 (2020) 109245) for linearized gas kinetic equations. On the other hand, solutions of macroscopic synthetic equations are used to inform the evolution of gas kinetic equation at the next iteration step. A rigorous linear Fourier stability analysis in periodic system shows that the error decay rate of GSIS can be smaller than 0.5, which means that the deviation to steady-state solution can be reduced by 3 orders of magnitude in 10 iterations. Other important advantages of the GSIS are (i) it does not rely on the specific form of Boltzmann collision operator and (ii) it can be solved by sophisticated techniques in computational fluid dynamics, making it amenable to large scale engineering applications. In this paper, the efficiency and accuracy of GSIS is demonstrated by a number of canonical test cases in rarefied gas dynamics.

preprint2020arXiv

Minimal 4D Composite Higgs Models Under Current LHC Constraints

We present preliminary results of the first convergent global fits of several minimal composite Higgs models. Our fits are performed using the differential evolution optimisation package $\texttt{Diver}$. A variety of physical constraints are taken into account, including a wide range of exclusion bounds on heavy resonance production from Run 2 of the LHC. As a by-product of the fits, we analyse the collider phenomenology of the lightest new up-type and down-type resonances in the viable regions of our models, finding some low-mass resonances that can be probed in future collider searches.

preprint2020arXiv

Preflare very long-periodic pulsations observed in Halpha emission before the onset of a solar flare

Very long-periodic pulsations during preflare phases (preflare-VLPs) have been detected in the full-disk solar soft X-Ray (SXR) flux (see Tan et al. 2016). They may be regarded as precursors to solar flares and may help us better understand the trigger mechanism of solar flares. In this letter, we report a preflare-VLP event before the onset of an M1.1 circular-ribbon flare on 2015 October 16. It was simultaneously observed in Halpha, SXR, and extreme ultraviolet (EUV) wavelengths, which were recorded by the NVST, GOES, EVE, and AIA respectively. The preflare-VLP is identified as the repeat and quasi-periodic pulses in light curves during preflare phase, which might be modulated by LRC-circuit oscillation in the current-carrying plasma loop. The quasi-periodicity can be determined from the Fourier power spectrum with Markov chain Monte Carlo (MCMC)-based Bayesian (e.g., Liang et al. 2020), such as ~9.3 minutes. We present the first report of a preflare-VLP event in the local Halpha line and EUV wavelength, which could be considered a precursor of a solar flare. This finding should therefore prove useful for the prediction of solar flares, especially for powerful flares.

preprint2020arXiv

Quasi-synchronization of bounded confidence opinion dynamics with stochastic asynchronous rule

Recently the theory of noise-induced synchronization of Hegselmann-Krause (HK) dynamics has been well developed. As a typical opinion dynamics of bounded confidence, the HK model obeys a synchronous updating rule, i.e., \emph{all} agents check and update their opinions at each time point. However, whether asynchronous bounded confidence models, including the famous Deffuant-Weisbuch (DW) model, can be synchronized by noise have not been theoretically proved. In this paper, we propose a generalized bounded confidence model which possesses a stochastic asynchronous rule. The model takes the DW model and the HK model as special cases and can significantly generalize the bounded confidence models to practical application. We discover that the asynchronous model possesses a different noise-based synchronization behavior compared to the synchronous HK model. Generally, the HK dynamics can achieve quasi-synchronization \emph{almost surely} under the drive of noise. For the asynchronous dynamics, we prove that the model can achieve quasi-synchronization \emph{in mean}, which is a new type of quasi-synchronization weaker than the "almost surely" sense. The results unify the theory of noise-induced synchronization of bounded confidence opinion dynamics and hence proves the noise-induced synchronization of DW model theoretically for the first time. Moreover, the results provide a theoretical foundation for developing noise-based control strategy of more complex social opinion systems with stochastic asynchronous rules.

preprint2020arXiv

The TianQin project: current progress on science and technology

TianQin is a planned space-based gravitational wave (GW) observatory consisting of three earth orbiting satellites with an orbital radius of about $10^5~{\rm km}$. The satellites will form a equilateral triangle constellation the plane of which is nearly perpendicular to the ecliptic plane. TianQin aims to detect GWs between $10^{-4}~{\rm Hz}$ and $1~{\rm Hz}$ that can be generated by a wide variety of important astrophysical and cosmological sources, including the inspiral of Galactic ultra-compact binaries, the inspiral of stellar-mass black hole binaries, extreme mass ratio inspirals, the merger of massive black hole binaries, and possibly the energetic processes in the very early universe or exotic sources such as cosmic strings. In order to start science operations around 2035, a roadmap called the 0123 plan is being used to bring the key technologies of TianQin to maturity, supported by the construction of a series of research facilities on the ground. Two major projects of the 0123 plan are being carried out. In this process, the team has created a new generation $17~{\rm cm}$ single-body hollow corner-cube retro-reflector which has been launched with the QueQiao satellite on 21 May 2018; a new laser ranging station equipped with a $1.2~{\rm m}$ telescope has been constructed and the station has successfully ranged to all the five retro-reflectors on the Moon; and the TianQin-1 experimental satellite has been launched on 20 December 2019 and the first round result shows that the satellite has exceeded all of its mission requirements.

preprint2019arXiv

Can we find steady-state solutions to multiscale rarefied gas flows within dozens of iterations?

One of the central problems in the study of rarefied gas dynamics is to find the steady-state solution of the Boltzmann equation quickly. When the Knudsen number is large, i.e. the system is highly rarefied, the conventional iteration scheme can lead to convergence within a few iterations. However, when the Knudsen number is small, i.e. the flow falls in the near-continuum regime, hundreds of thousands iterations are needed, and yet the "converged" solutions are prone to be contaminated by accumulated error and large numerical dissipation. Recently, based on the gas kinetic models, the implicit unified gas kinetic scheme (UGKS) and its variants have significantly reduced the iterations in the near-continuum flow regime, but still much higher than that of the highly rarefied gas flows. In this paper, we put forward a general synthetic iteration scheme (GSIS) to find the steady-state solutions of general rarefied gas flows within dozens of iterations at any Knudsen number. As the GSIS does not rely on the specific kinetic model/collision operator, it can be naturally extended to quickly find converged solutions for mixture flows and even flows involving chemical reactions. These two superior advantages are also expected to accelerate the slow convergence in simulation of near-continuum flows via the direct simulation Monte Carlo method and its low-variance version.

preprint2018arXiv

High-Order Implicit Hybridizable Discontinuous Galerkin Method for the Boltzmann Equation

The high-order hybridizable discontinuous Galerkin (HDG) method combining with an implicit iterative scheme is used to find the steady-state solution of the Boltzmann equation with full collision integral on two-dimensional triangular meshes. The velocity distribution function and its trace are approximated in the piecewise polynomial space of degree up to 4. The fast spectral method (FSM) is incorporated into the DG discretization to evaluate the collision operator. Specific polynomial approximation is proposed for the collision term to reduce the computational cost. The proposed scheme is proved to be accurate and efficient.

preprint2016arXiv

Detection of Single vs Multiple Antenna Transmission Systems Using Pilot Data

In this paper, we consider the problem of classifying the transmission system when it is not known a priori whether the transmission is via a single antenna or multiple antennas. The receiver is assumed to be employed with a known number of antennas. In a data frame transmitted by most multiple input multiple output (MIMO) systems, some pilot or training data is inserted for symbol timing synchronization and estimation of the channel. Our goal is to perform MIMO transmit antenna classification using this pilot data. More specifically, the problem of determining the transmission system is cast as a multiple hypothesis testing problem where the number of hypotheses is equal to the maximum number of transmit antennas. Under the assumption of receiver having the exact knowledge of pilot data used for timing synchronization and channel estimation, we consider maximum likelihood (ML) and correlation test statistics to classify the MIMO transmit system. When only probabilistic knowledge of pilot data is available at the receiver, a hybrid-maximum likelihood (HML) based test statistic is constructed using the expectation-maximization (EM) algorithm. The performance of the proposed algorithms is illustrated via simulations and comparative merits of different techniques in terms of the computational complexity and performance are discussed.

preprint2016arXiv

Interpreting the 750 GeV diphoton excess by the singlet extension of the Manohar-Wise Model

The evidence of a new scalar particle $X$ from the 750 GeV diphoton excess, and the absence of any other signal of new physics at the LHC so far suggest the existence of new colored scalars, which may be moderately light and thus can induce sizable $X g g$ and $X γγ$ couplings without resorting to very strong interactions. Motivated by this speculation, we extend the Manohar-Wise model by adding one gauge singlet scalar field. The resulting theory then predicts one singlet dominated scalar $ϕ$ as well as three kinds of color-octet scalars, which can mediate through loops the $ϕgg$ and $ϕγγ$ interactions. After fitting the model to the diphoton data at the LHC, we find that in reasonable parameter regions the excess can be explained at $1σ$ level by the process $ g g \to ϕ\to γγ$, and the best points predict the central value of the excess rate with $χ_{min}^2=2.32$, which corresponds to a $p$-value of $0.68$. We also consider the constraints from various LHC Run I signals, and we conclude that, although these constraints are powerful in excluding the parameter space of the model, the best points are still experimentally allowed.

preprint2016arXiv

Interpreting the 750 GeV diphoton excess in the Minimal Dilaton Model

We try to interpret the 750 GeV diphoton excess in the Minimal Dilaton Model, which extends the SM by adding one linearized dilaton field and vector-like fermions. We first show by analytic formulae in this framework that the production rates of the $γγ$, $gg$, $Zγ$, $ZZ$, $WW^\ast$, $t\bar{t}$ and $hh$ signals at the $750 {\rm GeV}$ resonance are only sensitive to the dilaton-Higgs mixing angle $θ_S$ and the parameter $η\equiv v N_X/f$, where $f$ is the dilaton decay constant and $N_X$ denotes the number of the fermions. Then we scan the two parameters by considering various theoretical and experimental constraints to find the solutions to the diphoton excess. We conclude that the model can predict the central value of the diphoton rate without conflicting with any constraints. The signatures of our explanation at the LHC Run II and the vacuum stability at high energy scale are also discussed.

preprint2016arXiv

Interpreting The 750 GeV Diphoton Excess Within Topflavor Seesaw Model

We propose to interpret the 750 GeV diphoton excess in a typical topflavor seesaw model. The new resonance X can be identified as a CP-even scalar emerging from a certain bi-doublet Higgs field. Such a scalar can couple to charged scalars, fermions as well as heavy gauge bosons predicted by the model, and consequently all of the particles contribute to the diphoton decay mode of the X. Numerical analysis indicates that the model can predict the central value of the diphoton excess without contradicting any constraints from 8 TeV LHC, and among the constraints, the tightest one comes from the Z γchannel, σ_{8 {\rm TeV}}^{Z γ} \lesssim 3.6 {\rm fb}, which requires σ_{13 {\rm TeV}}^{γγ} \lesssim 6 {\rm fb} in most of the favored parameter space.

preprint2016arXiv

Modulation Classification for MIMO-OFDM Signals via Approximate Bayesian Inference

The problem of modulation classification for a multiple-antenna (MIMO) system employing orthogonal frequency division multiplexing (OFDM) is investigated under the assumption of unknown frequency-selective fading channels and signal-to-noise ratio (SNR). The classification problem is formulated as a Bayesian inference task, and solutions are proposed based on Gibbs sampling and mean field variational inference. The proposed methods rely on a selection of the prior distributions that adopts a latent Dirichlet model for the modulation type and on the Bayesian network formalism. The Gibbs sampling method converges to the optimal Bayesian solution and, using numerical results, its accuracy is seen to improve for small sample sizes when switching to the mean field variational inference technique after a number of iterations. The speed of convergence is shown to improve via annealing and random restarts. While most of the literature on modulation classification assume that the channels are flat fading, that the number of receive antennas is no less than that of transmit antennas, and that a large number of observed data symbols are available, the proposed methods perform well under more general conditions. Finally, the proposed Bayesian methods are demonstrated to improve over existing non-Bayesian approaches based on independent component analysis and on prior Bayesian methods based on the `superconstellation' method.

preprint2016arXiv

Natural NMSSM after LHC Run I and the Higgsino dominated dark matter scenario

We investigate the impact of the direct searches for SUSY at LHC Run I on the naturalness of the Next-to-Minimal Supersymmetric Standard Model (NMSSM). For this end, we first scan the vast parameter space of the NMSSM to get the region where the fine tuning measures $Δ_Z$ and $Δ_h$ at the electroweak scale are less than about 50, then we implement by simulations the constraints of the direct searches on the parameter points in the region. Our results indicate that although the direct search experiments are effective in excluding the points, the parameter intervals for the region and also the minimum reaches of $Δ_Z$ and $Δ_h$ are scarcely changed by the constraints, which implies that the fine tuning of the NMSSM does not get worse after LHC Run I. Moreover, based on the results we propose a natural NMSSM scenario where the lightest neutralino $\tildeχ_1^0$ as the dark matter (DM) candidate is Higgsino-dominated. In this scenario, $Δ_Z$ and $Δ_h$ may be as low as 2 without conflicting with any experimental constraints, and intriguingly $\tildeχ_1^0$ can easily reach the measured DM relic density due to its significant Singlino component. We exhibit the features of the scenario which distinguish it from the other natural SUSY scenario, including the properties of its neutralino-chargino sector and scalar top quark sector. We emphasize that the scenario can be tested either through searching for $3 l + E_T^{miss}$ signal at 14 TeV LHC or through future DM direct detection experiments.

preprint2016arXiv

Noise Leads to Quasi-Consensus of Hegselmann-Krause Opinion Dynamics

This paper aims at providing rigorous theoretical analysis to investigate the consensus behavior of opinion dynamics in noisy environments. It is known that the well-known Hegselmann-Krause (HK) opinion dynamics demonstrates various agreement or disagreement behaviors in the deterministic case. Here we strictly show how noises provide great help to "synchronize" the opinions of the HK model. In fact, we prove a "critical phenomena" of the noisy HK dynamics, that is, the opinions merge as a quasi-consensus in finite time in noisy environment when the noise strength is below a critical value, which implies the fragmentation phenomenon of the HK dynamics could eventually vanish in the presence of noise. On the other hand, the opinions almost surely diverge when the noise strength exceeds the critical value.

preprint2016arXiv

SUSY effects in $R_b$: revisited under current experimental constraints

In this note we revisit the SUSY effects in $R_b$ under current experimental constraints including the LHC Higgs data, the $B$-physics measurements, the dark matter relic density and direct detection limits, as well as the precision electroweak data. We first perform a scan to figure out the currently allowed parameter space and then display the SUSY effects in $R_b$. We find that although the SUSY parameter space has been severely restrained by current experimental data, both the general MSSM and the natural-SUSY scenario can still alter $R_b$ with a magnitude sizable enough to be observed at future $Z$-factories (ILC, CEPC, FCC-ee, Super $Z$-factory) which produce $10^9-10^{12}$ $Z$-bosons. To be specific, assuming a precise measurement $δR_b = 2.0 \times 10^{-5}$ at FCC-ee, we can probe a right-handed stop up to 530 GeV through chargino-stop loops, probe a sbottom to 850 GeV through neutralino-sbottom loops and a charged Higgs to 770 GeV through the Higgs-top quark loops for a large $\text{tan}β$. The full one-loop SUSY correction to $R_b$ can reach $1 \times 10^{-4}$ in natural SUSY and $2 \times 10^{-4}$ in the general MSSM.

preprint2016arXiv

The Formation and Early Evolution of a Coronal Mass Ejection and its Associated Shock Wave on 2014 January 8

In this paper, we study the formation and early evolution of a limb coronal mass ejection (CME) and its associated shock wave that occurred on 2014 January 8. The extreme ultraviolet (EUV) images provided by the Atmospheric Imaging Assembly (AIA) on board \textit{Solar Dynamics Observatory} disclose that the CME first appears as a bubble-like structure. Subsequently, its expansion forms the CME and causes a quasi-circular EUV wave. Interestingly, both the CME and the wave front are clearly visible at all of the AIA EUV passbands. Through a detailed kinematical analysis, it is found that the expansion of the CME undergoes two phases: a first phase with a strong but transient lateral over-expansion followed by a second phase with a self-similar expansion. The temporal evolution of the expansion velocity coincides very well with the variation of the 25--50 keV hard X-ray flux of the associated flare, which indicates that magnetic reconnection most likely plays an important role in driving the expansion. Moreover, we find that, when the velocity of the CME reaches $\sim$600 km s$^{-1}$, the EUV wave starts to evolve into a shock wave, which is evidenced by the appearance of a type II radio burst. The shock's formation height is estimated to be $\sim$0.2$R_{sun}$, which is much lower than the height derived previously. Finally, we also study the thermal properties of the CME and the EUV wave. We find that the plasma in the CME leading front and the wave front has a temperature of $\sim$2 MK, while that in the CME core region and the flare region has a much higher temperature of $\ge$8 MK.

preprint2015arXiv

A Metastability Result for the Contact Process on a Random Regular Graph

In this paper we study the metastability of the contact process on a random regular graph. We show that the extinction time of the contact process, when initialized so that all vertices are infected at time 0, grows exponentially with the vertex number. Moreover, we show that the extinction time divided by its mean converges to a unit exponential distribution in law.

preprint2015arXiv

Contact Processes on Random Regular Graphs

We show that the contact process on a random $d$-regular graph initiated by a single infected vertex obeys the "cutoff phenomenon" in its supercritical phase. In particular, we prove that when the infection rate is larger than the critical value of the contact process on the infinite $d$-regular tree there are positive constants $C, p$ depending on the infection rate such that for sufficiently small $ε> 0$, when the number $n$ of vertices is large then (a) at times $t< (C - ε) \log n$ the fraction of infected vertices is vanishingly small, but (b) at time $(C + ε) \log n$ the fraction of infected vertices is within $ε$ of $p$, with probability $p$.

preprint2015arXiv

Testing the light dark matter scenario of the MSSM at the LHC

In the light dark matter (DM) scenario of the MSSM, the DM relic density puts non-trivial requirements on the spectrum of supersymmetric particles. As a result, the direct search for multi-lepton signals at the LHC has great impact on the scenario. In this work, we concentrate on the searches for sleptons and electroweak-inos at the LHC, and investigate the constraints on the light DM scenario from the 8 TeV LHC data as well as the capability of the 14 TeV LHC to test the scenario. We first get the samples of the scenario by scanning the vast parameter space of the MSSM and considering some easily available constraints, such as those from the DM relic density, the LUX experiment and the Higgs searches at colliders. Then for the surviving samples, we simulate the $2l+E_T^{miss}$ signal from slepton pair production process and the $2l+E_T^{miss}$ and $3l+E_T^{miss}$ signals from chargino and neutralino associated production processes at both the 8 TeV LHC and the 14 TeV LHC. Our simulations indicate that the 8 TeV LHC has excluded a sizable portion of the samples, and the 14 TeV LHC can be even much more powerful in testing the scenario. For example, in case that no multi-lepton signals are observed at the 14 TeV LHC, most samples of the light DM scenario will be excluded, especially a lower limit on the lightest neutralino mass will be set at 42 GeV and 44 GeV with 30 $fb^{-1}$ and 100 $fb^{-1}$ data respectively, and this limit can be further pushed up to $55\ {\rm GeV}$ with 300 $fb^{-1}$ data.

preprint2014arXiv

Branching random walks and contact processes on Galton-Watson trees

We consider branching random walks and contact processes on infinite, connected, locally finite graphs whose reproduction and infectivity rates across edges are inversely proportional to vertex degree. We show that when the ambient graph is a Galton-Watson tree then, in certain circumstances, the branching random walks and contact processes will have weak survival phases. We also provide bounds on critical values.

preprint2014arXiv

Dual-polarization light emission from InAs quantum dots in a annular photonic crystal cavity

The annular photonic crystals have been regarded as a satisfactory candidate to realize dual-polarization photonic device. In this letter, we focus our attention on the study of annular photonic crystal cavity to verify its application in light emission. We proposed a two-dimensional photonic crystal model with annular air units and a point-line defect to construct a cavity for the enhancement of light emission of InAs quantum dots. With the help of global optimization method, we have obtained an annular photonic crystal cavity design which can show a high in-plane quality factor of about 1.3*105 and 2.8*106 for transverse electric and transverse magnetic polarizations, respectively. Based on the Electron Beam Lithography and Reactive Ion Etching techniques, such cavity pattern was transferred into the top of InAs/GaAs active layer. The photoluminescence spectra of sample demonstrated clear light emission at around 1.3 um for both polarizations. Such dual-polarization light emitter has potential applications in optical communications, quantum information processing and other relative areas.

preprint2014arXiv

Modulation Classification via Gibbs Sampling Based on a Latent Dirichlet Bayesian Network

A novel Bayesian modulation classification scheme is proposed for a single-antenna system over frequency-selective fading channels. The method is based on Gibbs sampling as applied to a latent Dirichlet Bayesian network (BN). The use of the proposed latent Dirichlet BN provides a systematic solution to the convergence problem encountered by the conventional Gibbs sampling approach for modulation classification. The method generalizes, and is shown to improve upon, the state of the art.

preprint2014arXiv

Waiting time distribution of solar energetic particle events modeled with a non-stationary Poisson process

We present a study of the waiting time distributions (WTDs) of solar energetic particle (SEP) events observed with the spacecraft $WIND$ and $GOES$. Both the WTDs of solar electron events (SEEs) and solar proton events (SPEs) display a power-law tail $\sim Δt^{-γ}$. The SEEs display a broken power-law WTD. The power-law index is $γ_{1} =$ 0.99 for the short waiting times ($<$70 hours) and $γ_{2} =$ 1.92 for large waiting times ($>$100 hours). The break of the WTD of SEEs is probably due to the modulation of the corotating interaction regions (CIRs). The power-law index $γ\sim$ 1.82 is derived for the WTD of SPEs that is consistent with the WTD of type II radio bursts, indicating a close relationship between the shock wave and the production of energetic protons. The WTDs of SEP events can be modeled with a non-stationary Poisson process which was proposed to understand the waiting time statistics of solar flares (Wheatland 2000; Aschwanden $\&$ McTiernan 2010). We generalize the method and find that, if the SEP event rate $λ= 1/Δt$ varies as the time distribution of event rate $f(λ) = A λ^{-α}exp(-βλ)$, the time-dependent Poisson distribution can produce a power-law tail WTD $\sim Δt^{α- 3}$, where $0 \leq α< 2$.

preprint2013arXiv

Modulation Classification of MIMO-OFDM Signals by Independent Component Analysis and Support Vector Machines

A modulation classification (MC) scheme based on Independent Component Analysis (ICA) in conjunction with either maximum likelihood (ML) or Support Vector Machines (SVM) is proposed for MIMO-OFDM signals over frequency selective, time varying channels. The method is blind in the sense that it is assumed that the receiver has no information about the channel and transmitted signals other than that the spatial streams of signals are statistically independent. The processing consists of separation of the MIMO streams followed by modulation classification of the separated signals. While in general, blind separation of signals over frequency selective channels is a difficult problem, the non-frequency selective nature of the channel experienced by individual symbols in a MIMO-OFDM system enables the application of well-known ICA algorithms. Modulation classification is implemented by maximum likelihood and by an SVM-based modulation classification method relying on pre-selected modulation-dependent features. To improve performance in time varying channels, the invariance of the channel is exploited across the coherence bandwidth and the time coherence. The proposed method is shown to perform with high probability of correct classification over realistic ITU pedestrian and vehicular channels. An upper bound on the probability of correct classification is developed based on the Cramer Rao bound of channel estimation.

preprint2012arXiv

A Note on a Conjecture for Balanced Elementary Symmetric Boolean Functions

In 2008, Cusick {\it et al.} conjectured that certain elementary symmetric Boolean functions of the form $σ_{2^{t+1}l-1, 2^t}$ are the only nonlinear balanced ones, where $t$, $l$ are any positive integers, and $σ_{n,d}=\bigoplus_{1\le i_1<...<i_d\le n}x_{i_1}x_{i_2}...x_{i_d}$ for positive integers $n$, $1\le d\le n$. In this note, by analyzing the weight of $σ_{n, 2^t}$ and $σ_{n, d}$, we prove that ${\rm wt}(σ_{n, d})<2^{n-1}$ holds in most cases, and so does the conjecture. According to the remainder of modulo 4, we also consider the weight of $σ_{n, d}$ from two aspects: $n\equiv 3({\rm mod\}4)$ and $n\not\equiv 3({\rm mod\}4)$. Thus, we can simplify the conjecture. In particular, our results cover the most known results. In order to fully solve the conjecture, we also consider the weight of $σ_{n, 2^t+2^s}$ and give some experiment results on it.

preprint2012arXiv

Asymptotic Properties of Likelihood Based Linear Modulation Classification Systems

The problem of linear modulation classification using likelihood based methods is considered. Asymptotic properties of most commonly used classifiers in the literature are derived. These classifiers are based on hybrid likelihood ratio test (HLRT) and average likelihood ratio test (ALRT), respectively. Both a single-sensor setting and a multi-sensor setting that uses a distributed decision fusion approach are analyzed. For a modulation classification system using a single sensor, it is shown that HLRT achieves asymptotically vanishing probability of error (Pe) whereas the same result cannot be proven for ALRT. In a multi-sensor setting using soft decision fusion, conditions are derived under which Pe vanishes asymptotically. Furthermore, the asymptotic analysis of the fusion rule that assumes independent sensor decisions is carried out.

preprint2012arXiv

Characterization of Negabent Functions and Construction of Bent-Negabent Functions with Maximum Algebraic Degree

We present necessary and sufficient conditions for a Boolean function to be a negabent function for both even and odd number of variables, which demonstrate the relationship between negabent functions and bent functions. By using these necessary and sufficient conditions for Boolean functions to be negabent, we obtain that the nega spectrum of a negabent function has at most 4 values. We determine the nega spectrum distribution of negabent functions. Further, we provide a method to construct bent-negabent functions in $n$ variables ($n$ even) of algebraic degree ranging from 2 to $\frac{n}{2}$, which implies that the maximum algebraic degree of an $n$-variable bent-negabent function is equal to $\frac{n}{2}$. Thus, we answer two open problems proposed by Parker and Pott and by Stǎnicǎ \textit{et al.} respectively.

preprint2010arXiv

The ROSAT Bright Source 1RXS J201607.0+251645: An Active Algol-type Binary

1RXS J201607.0251645 is identified to be an eclipsing binary for the first time. We present the preliminary observations in V band with the 0.6-m telescope for three years and the extensive observations in V and R band with the 0.8-m telescope for six nights respectively. The light curve of the system is EB type. Five light minimum times were obtained and the orbital period of 0.388058d(0.00044d) is determined. The photometric solution given by the 2003-version Wilson-Devinney program suggests the binary is a semi-detached system with the photometric mass ratio 0.895(0.006), which probably comprise a G5 primary and an oversize K5 secondary. The less massive component has completely filled up its Roche lobe, while the other almost fills its Roche lobe with the filling factor of 93.4%. The system shows a varying O'Connell effect in its phase folded diagrams from 2005 to 2007, and is X-ray luminous with log(LXLbol)3.27. Possible mechanisms to account for these two phenomena are discussed. Finally, we infer the binary may be in thermal oscillation predicted by thermal relaxation oscillation (TRO) theory or may evolve into a contact binary.

preprint2007arXiv

Resolving the nature of the Rosette HH1 jet facing strong UV dissipation

The Rosette HH1 jet is a collimated flow immersed in the strong UV radiation field of the Rosette Nebula. We investigate the physical properties of the Rosette HH1 jet using high-quality narrow-band images and high-dispersion spectroscopy. The new images show that the axis of the jet is not precisely aligned with the star near the base of the jet. The high resolution of the spectra allows us to accurately determine the contributions from the \ion{H}{2} region, jet, and star. The appoaching and receding sides of the expanding shell of the Rosette Nebula are at heliocentric velocities of 13 and 40 km s$^{-1}$, while the jet reaches a maximum velocity offset at a heliocentric velocity of $-$30 km s$^{-1}$. The [\ion{S}{2}] doublet ratios indicate an electron density of $\sim$1000 cm$^{-3}$ in the jet and $\le$100 cm$^{-3}$ in the \ion{H}{2} region. With a careful subtraction of the nebular and jet components, we find the stellar H$α$ line is dominated by a broad absorption profile with little or no emission component, indicating a lack of substantial circumstellar material. The circumstellar material has most likely been photo-evaporated by the strong UV radiation field in the Rosette Nebula. The evaporation time scale is 10$^3$ -- 10$^4$ yr. The Rosette HH1 jet source provides evidence for an accelerated evolution from a CTTS to a WTTS due to the strong UV radiation field; therefore, both CTTSs and WTTSs can be spatially mixed in regions with massive star formation.