BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Sihong Shao

Sihong Shao contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
math.NANumerical Analysisquant-phphysics.comp-phmath.COmath.OCmath.SPmath.STStatistics Theory

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
8works
0followers
9topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

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

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

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

Building this graph slice

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

Published work

8 published item(s)

preprint2026arXiv

Equivalent spectral theory for fundamental graph cut problems

We introduce and develop equivalent spectral graph theory for several fundamental graph cut problems including maxcut, mincut, Cheeger cut, anti-Cheeger cut, dual Cheeger problem and their useful variants. A specified strategy for achieving an equivalent eigenproblem is proposed for a general graph cut problem via the set-pair Lovász extension and the Dinkelbach scheme. For a class of 2-cut and 3-cut problems, we reveal the intrinsic difference-of-submodularity for the fractional formulations and show that their set-pair Lovász extensions yield equivalent difference-of-convex structures. Building on the Dinkelbach scheme, we finally establish a unified research roadmap for nonlinear spectral theory that provides a one-to-one correspondence between certain eigenpairs and the optimal graph cut problems. The finer structure of the eigenvectors, the Courant nodal domain theorem and the graphic feature of eigenvalues are studied systematically in the setting of these new nonlinear eigenproblems.

0 citations0 reviewsNo code linkedOpen access
preprint2023arXiv

Application of Causal Inference Techniques to the Maximum Weight Independent Set Problem

A powerful technique for solving combinatorial optimization problems is to reduce the search space without compromising the solution quality by exploring intrinsic mathematical properties of the problems. For the maximum weight independent set (MWIS) problem, using an upper bound lemma which says the weight of any independent set not contained in the MWIS is bounded from above by the weight of the intersection of its closed neighbor set and the MWIS, we give two extension theorems -- independent set extension theorem and vertex cover extension theorem. With them at our disposal, two types of causal inference techniques (CITs) are proposed on the assumption that a vertex is strongly reducible (included or not included in all MWISs) or reducible (contained or not contained in a MWIS). One is a strongly reducible state-preserving technique, which extends a strongly reducible vertex into a vertex set where all vertices have the same strong reducibility. The other, as a reducible state-preserving technique, extends a reducible vertex into a vertex set with the same reducibility as that vertex and creates some weighted packing constraints to narrow the search space. Numerical experiments show that our CITs can help reduction algorithms find much smaller remaining graphs, improve the ability of exact algorithms to find the optimal solutions and help heuristic algorithms produce approximate solutions of better quality. In particular, detailed tests on $12$ representative graphs generated from datasets in Network Data Repository demonstrate that, compared to the state-of-the-art algorithms, the size of remaining graphs is further reduced by more than 32.6%, and the number of solvable instances is increased from 1 to 5.

0 citations0 reviewsNo code linkedOpen access
preprint2023arXiv

Nonlocalization of singular potentials in quantum dynamics

Nonlocal modeling has drawn more and more attention and becomes steadily more powerful in scientific computing. In this paper, we demonstrate the superiority of a first-principle nonlocal model -- Wigner function -- in treating singular potentials which are often used to model the interaction between point charges in quantum science. The nonlocal nature of the Wigner equation is fully exploited to convert the singular potential into the Wigner kernel with weak or even no singularity, and thus highly accurate numerical approximations are achievable, which are hardly designed when the singular potential is taken into account in the local Schrödinger equation. The Dirac delta function, the logarithmic, and the inverse power potentials are considered. Numerically converged Wigner functions under all these singular potentials are obtained with an operator splitting spectral method, and display many interesting quantum behaviors as well.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

A characteristic-spectral-mixed scheme for six-dimensional Wigner-Coulomb dynamics

Numerical resolution for 6-D Wigner dynamics under the Coulomb potential faces with the combined challenges of high dimensionality, nonlocality, oscillation and singularity. In particular, the extremely huge memory storage of 6-D grids hinders the usage of all existing deterministic numerical scheme, which is well-known as the curse of dimensionality. To surmount these difficulties, we propose a massively parallel solver, termed the CHAracteristic-Spectral-Mixed (CHASM) scheme, by fully exploiting two distinct features of the Wigner equation: Locality of spatial advection and nonlocality of quantum interaction. Our scheme utilizes the local cubic B-spline basis to interpolate the local spatial advection. The key is to use a perfectly matched boundary condition to give a closure of spline coefficients, so that distributed pieces can recover the global one as accurately as possible owing to the rapid decay of wavelet basis in the dual space, and communication costs are significantly reduced. To resolve the nonlocal pseudodifferential operator with weakly singular symbol, CHASM further adopts the truncated kernel method to attain a highly efficient approximation. Several typical experiments including the quantum harmonic oscillator and Hydrogen 1s state demonstrate the accuracy and efficiency of CHASM. The non-equilibrium electron-proton couplings are also clearly displayed and reveal the uncertainty principle and quantum tunneling in phase space. Finally, the scalability of CHASM up to 16000 cores is presented.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

A higher-order accurate operator splitting spectral method for the Wigner-Poisson system

An accurate description of 2-D quantum transport in a double-gate metal oxide semiconductor filed effect transistor (dgMOSFET) requires a high-resolution solver to a coupled system of the 4-D Wigner equation and 2-D Poisson equation. In this paper, we propose an operator splitting spectral method to evolve such Wigner-Poisson system in 4-D phase space with high accuracy. After an operator splitting of the Wigner equation, the resulting two sub-equations can be solved analytically with spectral approximation in phase space. Meanwhile, we adopt a Chebyshev spectral method to solve the Poisson equation. Spectral convergence in phase space and a fourth-order accuracy in time are both numerically verified. Finally, we apply the proposed solver into simulating dgMOSFET, develop the steady states from long-time simulations and obtain numerically converged current-voltage (I-V) curves.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Adaptive Hermite Spectral Methods in Unbounded Domains

Recently, new adaptive techniques were developed that greatly improved the efficiency of solving PDEs using spectral methods. These adaptive spectral techniques are especially suited for accurately solving problems in unbounded domains and require the monitoring and dynamic adjustment of three key tunable parameters: the scaling factor, the displacement of the basis functions, and the spectral expansion order. There have been few analyses of numerical methods for unbounded domain problems. Specifically, there is no analysis of adaptive spectral methods to provide insight into how to increase efficiency and accuracy through dynamical adjustment of parameters. In this paper, we perform the first numerical analysis of the adaptive spectral method using generalized Hermite functions in both one- and multi-dimensional problems. Our analysis reveals why adaptive spectral methods work well when a "frequency indicator" of the numerical solution is controlled. We then investigate how the implementation of the adaptive spectral methods affects numerical results, thereby providing guidelines for the proper tuning of parameters. Finally, we further improve performance by extending the adaptive methods to allow bidirectional basis function translation, and the prospect of carrying out similar numerical analysis to solving PDEs arising from realistic difficult-to-solve unbounded models with adaptive spectral methods is also briefly discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Performance evaluations on the parallel CHAracteristic-Spectral-Mixed (CHASM) scheme

Performance evaluations on the deterministic algorithms for 6-D problems are rarely found in literatures except some recent advances in the Vlasov and Boltzmann community [Dimarco et al. (2018), Kormann et al. (2019)], due to the extremely high complexity. Thus a detailed comparison among various techniques shall be useful to the researchers in the related fields. We try to make a thorough evaluation on a parallel CHAracteristic-Spectral-Mixed (CHASM) scheme to support its usage. CHASM utilizes the cubic B-spline expansion in the spatial space and spectral expansion in the momentum space, which many potentially overcome the computational burden in solving classical and quantum kinetic equations in 6-D phase space. Our purpose is three-pronged. First, we would like show that by imposing some effective Hermite boundary conditions, the local cubic spline can approximate to the global one as accurately as possible. Second, we will illustrate the necessity of adopting the truncated kernel method in calculating the pseudodifferential operator with a singular symbol, since the widely used pseudo-spectral method [Ringhofer (1990)] might fail to properly tackle the singularity. Finally, we make a comparison among non-splitting Lawson schemes and Strang operator splitting. Our numerical results demonstrate the advantage of the one-stage Lawson predictor-corrector scheme over multi-stage ones as well as the splitting scheme in both accuracy and stability.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

SPADE: Sequential-clustering Particle Annihilation via Discrepancy Estimation

For an empirical signed measure $μ= \frac{1}{N} \left(\sum_{i=1}^P δ_{x_i} - \sum_{i=1}^M δ_{y_i}\right)$, particle annihilation (PA) removes $N_A$ particles from both $\{x_i\}_{i=1}^P$ and $\{y_i\}_{i=1}^M$ simultaneously, yielding another empirical signed measure $ν$ such that $\int f d ν$ approximates to $\int f d μ$ within an acceptable accuracy for suitable test functions $f$. Such annihilation of particles carrying opposite importance weights has been extensively utilized for alleviating the numerical sign problem in particle simulations. In this paper, we propose an algorithm for PA in high-dimensional Euclidean space based on hybrid of clustering and matching, dubbed the Sequential-clustering Particle Annihilation via Discrepancy Estimation (SPADE). It consists of two steps: Adaptive clustering of particles via controlling their number-theoretic discrepancies, and independent random matching among positive and negative particles in each cluster. Both deterministic error bounds by the Koksma-Hlawka inequality and non-asymptotic random error bounds by concentration inequalities are proved to be affected by two factors. One factor measures the irregularity of point distributions and reflects their discrete nature. The other relies on the variation of test function and is influenced by the continuity. Only the latter implicitly depends on dimensionality $d$, implying that SPADE can be immune to the curse of dimensionality for a wide class of test functions. Numerical experiments up to $d=1080$ validate our theoretical discoveries.

0 citations0 reviewsNo code linkedOpen access