Source author record

Sihong Shao

Sihong Shao appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

math.NA physics.comp-ph quant-ph Numerical Analysis math.SP math.CO nlin.PS nlin.SI hep-th math-ph math.MP math.OC math.ST physics.chem-ph Statistics Theory

Catalog footprint

What is connected

17works

15topics

4close collaborators

Actions

Connect this record

Open graph Browse works

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.

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.

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.

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.

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.

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.

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.

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.

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.

preprint2016arXiv

An advective-spectral-mixed method for time-dependent many-body Wigner simulations

As a phase space language for quantum mechanics, the Wigner function approach bears a close analogy to classical mechanics and has been drawing growing attention, especially in simulating quantum many-body systems. However, deterministic numerical solutions have been almost exclusively confined to one-dimensional one-body systems and few results are reported even for one-dimensional two-body problems. This paper serves as the first attempt to solve the time-dependent many-body Wigner equation through a grid-based advective-spectral-mixed method. The main feature of the method is to resolve the linear advection in $(\bm{x},t)$-space by an explicit three-step characteristic scheme coupled with the piecewise cubic spline interpolation, while the Chebyshev spectral element method in $\bm k$-space is adopted for accurate calculation of the nonlocal pseudo-differential term. Not only the time step of the resulting method is not restricted by the usual CFL condition and thus a large time step is allowed, but also the mass conservation can be maintained. In particular, for the system consisting of identical particles, the advective-spectral-mixed method can also rigorously preserve physical symmetry relations. The performance is validated through several typical numerical experiments, like the Gaussian barrier scattering, electron-electron interaction and a Helium-like system, where the third-order accuracy against both grid spacing and time stepping is observed.

preprint2016arXiv

Nodal Domains of Eigenvectors for $1$-Laplacian on Graphs

The eigenvectors for graph $1$-Laplacian possess some sort of localization property: On one hand, any nodal domain of an eigenvector is again an eigenvector with the same eigenvalue; on the other hand, one can pack up an eigenvector for a new graph by several fundamental eigencomponents and modules with the same eigenvalue via few special techniques. The Courant nodal domain theorem for graphs is extended to graph $1$-Laplacian for strong nodal domains, but for weak nodal domains it is false. The notion of algebraic multiplicity is introduced in order to provide a more precise estimate of the number of independent eigenvectors. A positive answer is given to a question raised in [{\sl K.~C. Chang, Spectrum of the $1$-Laplacian and Cheeger constant on graphs, J. Graph Theor., DOI: 10.1002/jgt.21871}], to confirm that the critical values obtained by the minimax principle may not cover all eigenvalues of graph $1$-Laplacian.

preprint2016arXiv

The $1$-Laplacian Cheeger Cut: Theory and Algorithms

This paper presents a detailed review of both theory and algorithms for the Cheeger cut based on the graph $1$-Laplacian. In virtue of the cell structure of the feasible set, we propose a cell descend (CD) framework for achieving the Cheeger cut. While plugging the relaxation to guarantee the decrease of the objective value in the feasible set, from which both the inverse power (IP) method and the steepest descent (SD) method can also be recovered, we are able to get two specified CD methods. Comparisons of all these methods are conducted on several typical graphs.

preprint2015arXiv

Nonlinear Dirac equation solitary waves in the presence of external driving forces

We consider the nonlinear Dirac (NLD) equation in 1+1 dimension with scalar-scalar self-interaction in the presence of external forces as well as damping of the form $ f(x,t) - i μγ^0 Ψ$, where both $f$ and $Ψ$ are two-component spinors. We develop an approximate variational approach using collective coordinates (CC) for studying the time dependent response of the solitary waves to these external forces. This approach predicts intrinsic oscillations of the solitary waves, i.e. the amplitude, width and phase all oscillate with the same frequency. The translational motion is also affected, because the soliton position oscillates around a mean trajectory. We then compare the results of the variational approximation with numerical simulations of the NLD equation, and find a good agreement, if we take into account a certain linear excitation with specific wavenumber that is excited together with the intrinsic oscillations such that the momentum in a transformed NLD equation is conserved. We also solve explicitly the CC equations of the variational approximation in the non-relativistic regime for a homogeneous external force and obtain excellent agreement with the numerical solution of the CC equations.

preprint2014arXiv

Stability of solitary waves in the nonlinear Dirac equation with arbitrary nonlinearity

We consider the nonlinear Dirac equation in 1+1 dimension with scalar-scalar self interaction $ \frac{g^2}{κ+1} ({\bar Ψ} Ψ)^{κ+1}$ and with mass $m$. Using the exact analytic form for rest frame solitary waves of the form $Ψ(x,t) = ψ(x) e^{-i ωt}$ for arbitrary $ κ$, we discuss the validity of various approaches to understanding stability that were successful for the nonlinear Schrödinger equation. In particular we study the validity of a version of Derrick's theorem, the criterion of Bogolubsky as well as the Vakhitov-Kolokolov criterion, and find that these criteria yield inconsistent results. Therefore, we study the stability by numerical simulations using a recently developed 4th-order operator splitting integration method. For different ranges of $κ$ we map out the stability regimes in $ω$. We find that all stable nonlinear Dirac solitary waves have a one-hump profile, but not all one-hump waves are stable, while all waves with two humps are unstable. We also find that the time $t_c$, it takes for the instability to set in, is an exponentially increasing function of $ω$ and $t_c$ decreases monotonically with increasing $κ$.

preprint2013arXiv

Analysis of the Time Reversible Born-Oppenheimer Molecular Dynamics

We analyze the time reversible Born-Oppenheimer molecular dynamics (TRBOMD) scheme, which preserves the time reversibility of the Born-Oppenheimer molecular dynamics even with non-convergent self-consistent field iteration. In the linear response regime, we derive the stability condition as well as the accuracy of TRBOMD for computing physical properties such as the phonon frequency obtained from the molecular dynamic simulation. We connect and compare TRBOMD with the Car-Parrinello molecular dynamics in terms of accuracy and stability. We further discuss the accuracy of TRBOMD beyond the linear response regime for non-equilibrium dynamics of nuclei. Our results are demonstrated through numerical experiments using a simplified one dimensional model for Kohn-Sham density functional theory.

preprint2013arXiv

Efficient iterative method for solving the Dirac-Kohn-Sham density functional theory

We present for the first time an efficient iterative method to directly solve the four-component Dirac-Kohn-Sham (DKS) density functional theory. Due to the existence of the negative energy continuum in the DKS operator, the existing iterative techniques for solving the Kohn-Sham systems cannot be efficiently applied to solve the DKS systems. The key component of our method is a novel filtering step (F) which acts as a preconditioner in the framework of the locally optimal block preconditioned conjugate gradient (LOBPCG) method. The resulting method, dubbed the LOBPCG-F method, is able to compute the desired eigenvalues and eigenvectors in the positive energy band without computing any state in the negative energy band. The LOBPCG-F method introduces mild extra cost compared to the standard LOBPCG method and can be easily implemented. We demonstrate our method in the pseudopotential framework with a planewave basis set which naturally satisfies the kinetic balance prescription. Numerical results for Pt$_{2}$, Au$_{2}$, TlF, and Bi$_{2}$Se$_{3}$ indicate that the LOBPCG-F method is a robust and efficient method for investigating the relativistic effect in systems containing heavy elements.

preprint2013arXiv

Multi-hump solitary waves of nonlinear Dirac equation

This paper concentrates on a (1+1)-dimensional nonlinear Dirac (NLD) equation with a general self-interaction, being a linear combination of the scalar, pseudoscalar, vector and axial vector self-interactions to the power of the integer $k+1$. The solitary wave solutions to the NLD equation are analytically derived, and the upper bounds of the hump number in the charge, energy and momentum densities for the solitary waves are proved in theory. The results show that: (1) for a given integer $k$, the hump number in the charge density is not bigger than $4$, while that in the energy density is not bigger than $3$; (2) those upper bounds can only be achieved in the situation of higher nonlinearity, namely, $k\in\{5,6,7,\cdots \}$ for the charge density and $k\in\{3,5,7,\cdots\}$ for the energy density; (3) the momentum density has the same multi-hump structure as the energy density; (4) more than two humps (resp. one hump) in the charge (resp. energy) density can only happen under the linear combination of the pseudoscalar self-interaction and at least one of the scalar and vector (or axial vector) self-interactions. Our results on the multi-hump structure will be interesting in the interaction dynamics for the NLD solitary waves.

preprint2013arXiv

Numerical methods for nonlinear Dirac equation

This paper presents a review of the current state-of-the-art of numerical methods for nonlinear Dirac (NLD) equation. Several methods are extendedly proposed for the (1+1)-dimensional NLD equation with the scalar and vector self-interaction and analyzed in the way of the accuracy and the time reversibility as well as the conservation of the discrete charge, energy and linear momentum. Those methods are the Crank-Nicolson (CN) schemes, the linearized CN schemes, the odd-even hopscotch scheme, the leapfrog scheme, a semi-implicit finite difference scheme, and the exponential operator splitting (OS) schemes. The nonlinear subproblems resulted from the OS schemes are analytically solved by fully exploiting the local conservation laws of the NLD equation. The effectiveness of the various numerical methods, with special focus on the error growth and the computational cost, is illustrated on two numerical experiments, compared to two high-order accurate Runge-Kutta discontinuous Galerkin methods. Theoretical and numerical comparisons show that the high-order accurate OS schemes may compete well with other numerical schemes discussed here in terms of the accuracy and the efficiency. A fourth-order accurate OS scheme is further applied to investigating the interaction dynamics of the NLD solitary waves under the scalar and vector self-interaction. The results show that the interaction dynamics of two NLD solitary waves depend on the exponent power of the self-interaction in the NLD equation; collapse happens after collision of two equal one-humped NLD solitary waves under the cubic vector self-interaction in contrast to no collapse scattering for corresponding quadric case.

Sihong Shao

What is connected

Connect this record

See the researcher in context

Building this map preview

17 published item(s)

Equivalent spectral theory for fundamental graph cut problems

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

Nonlocalization of singular potentials in quantum dynamics

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

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

Adaptive Hermite Spectral Methods in Unbounded Domains

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

SPADE: Sequential-clustering Particle Annihilation via Discrepancy Estimation

An advective-spectral-mixed method for time-dependent many-body Wigner simulations

Nodal Domains of Eigenvectors for $1$-Laplacian on Graphs

The $1$-Laplacian Cheeger Cut: Theory and Algorithms

Nonlinear Dirac equation solitary waves in the presence of external driving forces

Stability of solitary waves in the nonlinear Dirac equation with arbitrary nonlinearity

Analysis of the Time Reversible Born-Oppenheimer Molecular Dynamics

Efficient iterative method for solving the Dirac-Kohn-Sham density functional theory

Multi-hump solitary waves of nonlinear Dirac equation

Numerical methods for nonlinear Dirac equation