Researcher profile

Shidong Jiang

Shidong Jiang contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 19 - UnverifiedVerification L1Unclaimed author
5works
0followers
6topics
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

5 published item(s)

preprint2026arXiv

An $O(\log N)$ Monte Carlo method for periodic Coulomb systems

Efficient Monte Carlo (MC) sampling of many-body systems with long-range electrostatics is often limited by the cost of per-move energy-difference evaluation under periodic boundary conditions. We present DMK-MC, an accelerated MC method that adapts the dual-space multilevel kernel-splitting (DMK) framework to single-particle Metropolis updates. DMK-MC computes the energy change and, upon acceptance, updates the stored incoming plane-wave fields with $O(1)$ work per tree level, yielding an overall $O(\log N)$ expected work per trial move for fixed accuracy. The method decomposes the Coulomb kernel into three components: a global, periodized smooth part; a multilevel sequence of smooth difference kernels whose interactions are restricted to same-level colleague boxes; and a singular residual kernel whose short-range interactions are evaluated directly. Benchmarks on uniform, highly nonuniform, and implicit-solvent electrolyte and colloidal configurations show that DMK-MC consistently outperforms a recent FMM-based $O(\log N)$ Monte Carlo method, delivering several-fold speedups at comparable tolerances.

preprint2026arXiv

Fast Ewald Summation with Prolates for Charged Systems in the NPT Ensemble

We present an NPT extension of Ewald summation with prolates (ESP), a spectrally accurate and scalable particle-mesh method for molecular dynamics simulations of periodic, charged systems. Building on the recently introduced ESP framework, this work focuses on rigorous and thermodynamically consistent pressure/stress evaluation in the isothermal--isobaric ensemble. ESP employs prolate spheroidal wave functions as both splitting and spreading kernels, reducing the Fourier grid size needed to reach a prescribed pressure accuracy compared with current widely used mesh-Ewald methods based on Gaussian splitting and B-spline spreading. We derive a unified pressure-tensor formulation applicable to isotropic, semi-isotropic, anisotropic, and fully flexible cells, and show that the long-range pressure can be evaluated with a single forward FFT followed by diagonal scaling, whereas force evaluation requires both forward and inverse transforms. We provide production implementations in LAMMPS and GROMACS and validate pressure and force accuracy on bulk water, LiTFSI ionic liquids, and a transmembrane system. Benchmarks on up to $3\times 10^3$ CPU cores demonstrate strong scaling and reduced communication cost at matched accuracy, particularly for NPT pressure evaluation.

preprint2020arXiv

An Integral Equation Method for the Cahn-Hilliard Equation in the Wetting Problem

We present an integral equation approach to solving the Cahn-Hilliard equation equipped with boundary conditions that model solid surfaces with prescribed Young's angles. The discretization of the system in time using convex splitting leads to a modified biharmonic equation at each time step. To solve it, we split the solution into a volume potential computed with free space kernels, plus the solution to a second kind integral equation (SKIE). The volume potential is evaluated with the help of a box-based volume-FMM method. For non-box domains, source density is extended by solving a biharmonic Dirichlet problem. The near-singular boundary integrals are computed using quadrature by expansion (QBX) with FMM acceleration. Our method has linear complexity in the number of surface/volume degrees of freedom and can achieve high order convergence with adaptive refinement to manage error from function extension.

preprint2020arXiv

On time-domain NRBC for Maxwell's equations and its application in accurate simulation of electromagnetic invisibility cloaks

In this paper, we present analytic formulas of the temporal convolution kernel functions involved in the time-domain non-reflecting boundary condition (NRBC) for the electromagnetic scattering problems. Such exact formulas themselves lead to accurate and efficient algorithms for computing the NRBC for domain reduction of the time-domain Maxwell's system in $\mathbb R^3$. A second purpose of this paper is to derive a new time-domain model for the electromagnetic invisibility cloak. Different from the existing models, it contains only one unknown field and the seemingly complicated convolutions can be computed as efficiently as the temporal convolutions in the NRBC. The governing equation in the cloaking layer is valid for general geometry, e.g., a spherical or polygonal layer. Here, we aim at simulating the spherical invisibility cloak. We take the advantage of radially stratified dispersive media and special geometry, and develop an efficient vector spherical harmonic (VSH)-spectral-element method for its accurate simulation. Compared with limited results on FDTD simulation, the proposed method is optimal in both accuracy and computational cost. Indeed, the saving in computational time is significant.

preprint2019arXiv

Evaluation of Abramowitz functions in the right half of the complex plane

A numerical scheme is developed for the evaluation of Abramowitz functions $J_n$ in the right half of the complex plane. For $n=-1,\, \ldots,\, 2$, the scheme utilizes series expansions for $|z|<1$ and asymptotic expansions for $|z|>R$ with $R$ determined by the required precision, and modified Laurent series expansions which are precomputed via a least squares procedure to approximate $J_n$ accurately and efficiently on each sub-region in the intermediate region $1\le |z| \le R$. For $n>2$, $J_n$ is evaluated via a recurrence relation. The scheme achieves nearly machine precision for $n=-1, \ldots, 2$, with the cost about four times of evaluating a complex exponential per function evaluation.