BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Shidong Jiang

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

ResearcherUnclaimed source record
math.NAphysics.opticscond-mat.mtrl-scicond-mat.otherNumerical Analysisphysics.atm-clusphysics.atom-phcond-mat.softastro-phcond-mat.mes-hallmath-phmath.MPphysics.chem-phphysics.comp-ph

Catalog footprint

What is connected

17works
14topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

17 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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A Quantitative Study of the Effect of Cladding Thickness on Modal Confinement Loss in Photonic Waveguides

There has been increasing interest in making the photonic devices more and more compact in the integrated photonics industry, and one of the important questions for manufacturers and design engineers is how to quantify the effect of the finite cladding thickness on the modal confinement loss of photonic waveguides. This requires at least six to seven digits accuracy for the computation of propagation constant $β$ since the modal confinement loss is proportional to the imaginary part of $β$ that is six to seven orders of magnitude smaller than its real part by the industrial standard. In this paper, we present an accurate and efficient method to compute the propagation constant of electromagnetic modes of photonic waveguides with arbitrary number of (nonsmooth) inclusions in a layered media. The method combines a well-conditioned boundary integral equation formulation for photonic waveguides which requires the discretization of the material interface only, and efficient Sommerfeld integral representations to treat the effect of the layered medium. Our scheme is capable of calculating the propagation loss of the electromagnetic modes with high fidelity, even for waveguides with corners imbedded in a cladding material of finite thickness. The numerical results, with more than $10$-digit accuracy, show quantitatively that the modal confinement loss of the rectangular waveguide increases exponentially fast as the cladding thickness decreases.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Fast evaluation of the Caputo fractional derivative and its applications to fractional diffusion equations

We present an efficient algorithm for the evaluation of the Caputo fractional derivative $_0^C\!D_t^αf(t)$ of order $α\in (0,1)$, which can be expressed as a convolution of $f'(t)$ with the kernel $t^{-α}$. The algorithm is based on an efficient sum-of-exponentials approximation for the kernel $t^{-1-α}$ on the interval $[Δt, T]$ with a uniform absolute error $\varepsilon$, where the number of exponentials $N_{\text{exp}}$ needed is of the order $O\left(\log\frac{1}{\varepsilon}\left( \log\log\frac{1}{\varepsilon}+\log\frac{T}{Δt}\right) +\log\frac{1}{Δt}\left( \log\log\frac{1}{\varepsilon}+\log\frac{1}{Δt}\right) \right)$. As compared with the direct method, the resulting algorithm reduces the storage requirement from $O(N_T)$ to $O(N_{\text{exp}})$ and the overall computational cost from $O(N_T^2)$ to $O(N_TN_{\text{exp}})$ with $N_T$ the total number of time steps. Furthermore, when the fast evaluation scheme of the Caputo derivative is applied to solve the fractional diffusion equations, the resulting algorithm requires only $O(N_SN_{\text{exp}})$ storage and $O(N_SN_TN_{\text{exp}})$ work with $N_S$ the total number of points in space; whereas the direct methods require $O(N_SN_T$) storage and $O(N_SN_T^2)$ work. The complexity of both algorithms is nearly optimal since $N_{\text{exp}}$ is of the order $O(\log N_T)$ for $T\gg 1$ or $O(\log^2N_T)$ for $T\approx 1$ for fixed accuracy $\varepsilon$. We also present a detailed stability and error analysis of the new scheme for solving linear fractional diffusion equations. The performance of the new algorithm is illustrated via several numerical examples. Finally, the algorithm can be parallelized in a straightforward manner.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Second Kind Integral Equation Formulation for the Mode Calculation of Optical Waveguides

We present a second kind integral equation (SKIE) formulation for calculating the electromagnetic modes of optical waveguides, where the unknowns are only on material interfaces. The resulting numerical algorithm can handle optical waveguides with a large number of inclusions of arbitrary irregular cross section. It is capable of finding the bound, leaky, and complex modes for optical fibers and waveguides including photonic crystal fibers (PCF), dielectric fibers and waveguides. Most importantly, the formulation is well conditioned even in the case of nonsmooth geometries. Our method is highly accurate and thus can be used to calculate the propagation loss of the electromagnetic modes accurately, which provides the photonics industry a reliable tool for the design of more compact and efficient photonic devices. We illustrate and validate the performance of our method through extensive numerical studies and by comparison with semi-analytical results and previously published results.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Computing the ground state and dynamics of the nonlinear Schrödinger equation with nonlocal interactions via the nonuniform FFT

We present efficient and accurate numerical methods for computing the ground state and dynamics of the nonlinear Schrödinger equation (NLSE) with nonlocal interactions based on a fast and accurate evaluation of the long-range interactions via the nonuniform fast Fourier transform (NUFFT). We begin with a review of the fast and accurate NUFFT based method in \cite{JGB} for nonlocal interactions where the singularity of the Fourier symbol of the interaction kernel at the origin can be canceled by switching to spherical or polar coordinates. We then extend the method to compute other nonlocal interactions whose Fourier symbols have stronger singularity at the origin that cannot be canceled by the coordinate transform. Many of these interactions do not decay at infinity in the physical space, which adds another layer of complexity since it is more difficult to impose the correct artificial boundary conditions for the truncated bounded computational domain. The performance of our method against other existing methods is illustrated numerically, with particular attention on the effect of the size of the computational domain in the physical space. Finally, to study the ground state and dynamics of the NLSE, we propose efficient and accurate numerical methods by combining the NUFFT method for potential evaluation with the normalized gradient flow using backward Euler Fourier pseudospectral discretization and time-splitting Fourier pseudospectral method, respectively. Extensive numerical comparisons are carried out between these methods and other existing methods for computing the ground state and dynamics of the NLSE with various nonlocal interactions. Numerical results show that our scheme performs much better than those existing methods in terms of both accuracy and efficiency.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Efficient Brownian Dynamics Simulation of Single DNA with Hydrodynamic Interactions in Linear Flows

The coarse-grained molecular dynamics (MD) or Brownian dynamics (BD) simulation is a particle-based approach that has been applied to a wide range of biological problems that involve interactions with surrounding fluid molecules or the so-called hydrodynamic interactions (HIs). In this paper, an efficient algorithm is proposed to simulate the motion of a single DNA molecule in linear flows. The algorithm utilizes the integraing factor to cope with the effect of the linear flow of the surrounding fluid and applies the Metropolis method (MM) in [N. Bou-Rabee, A. Donev, and E. Vanden-Eijnden, Multiscale Model. Simul. 12, 781 (2014)] to achieve more efficient BD simulation. Thus our method permits much larger time step size than previous methods while still maintaining the stability of the BD simulation, which is advantageous for long-time BD simulation. Our numerical results on $λ$-DNA agree very well with both experimental data and previous simulation results. Finally, when combined with fast algorithms such as the fast multipole method which has nearly optimal complexity in the total number of beads, the resulting method is parallelizable, scalable to large systems, and stable for large time step size, thus making the long-time large-scale BD simulation within practical reach. This will be useful for the study of membranes, long-chain molecules, and a large collection of molecules in the fluids.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Fast and Accurate Evaluation of Nonlocal Coulomb and Dipole-Dipole Interactions via the Nonuniform FFT

We present a fast and accurate algorithm for the evaluation of nonlocal (long-range) Coulomb and dipole-dipole interactions in free space. The governing potential is simply the convolution of an interaction kernel $U(\bx)$ and a density function $ρ(\bx)=|ψ(\bx)|^2$, for some complex-valued wave function $ψ(\bx)$, permitting the formal use of Fourier methods. These are hampered by the fact that the Fourier transform of the interaction kernel $\widehat{U}(\bk)$ has a singularity at the origin $\bk={\bf 0}$ in Fourier (phase) space. Thus, accuracy is lost when using a uniform Cartesian grid in $\bk$ which would otherwise permit the use of the FFT for evaluating the convolution. Here, we make use of a high-order discretization of the Fourier integral, accelerated by the nonuniform fast Fourier transform (NUFFT). By adopting spherical and polar phase-space discretizations in three and two dimensions, respectively, the singularity in $\hat{U}(\bk)$ at the origin is canceled, so that only a modest number of degrees of freedom are required to evaluate the Fourier integral, assuming that the density function $ρ(\bx)$ is smooth and decays sufficiently fast as $\bx \rightarrow \infty$. More precisely, the calculation requires $O(N\log N)$ operations, where $N$ is the total number of discretization points in the computational domain. Numerical examples are presented to demonstrate the performance of the algorithm.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Efficient sum-of-exponentials approximations for the heat kernel and their applications

In this paper, we show that efficient separated sum-of-exponentials approximations can be constructed for the heat kernel in any dimension. In one space dimension, the heat kernel admits an approximation involving a number of terms that is of the order $O(\log(\frac{T}δ) (\log(\frac{1}ε)+\log\log(\frac{T}δ)))$ for any $x\in\bbR$ and $δ\leq t \leq T$, where $ε$ is the desired precision. In all higher dimensions, the corresponding heat kernel admits an approximation involving only $O(\log^2(\frac{T}δ))$ terms for fixed accuracy $ε$. These approximations can be used to accelerate integral equation-based methods for boundary value problems governed by the heat equation in complex geometry. The resulting algorithms are nearly optimal. For $N_S$ points in the spatial discretization and $N_T$ time steps, the cost is $O(N_S N_T \log^2 \frac{T}δ)$ in terms of both memory and CPU time for fixed accuracy $ε$. The algorithms can be parallelized in a straightforward manner. Several numerical examples are presented to illustrate the accuracy and stability of these approximations.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The solution of the scalar wave equation in the exterior of a sphere

We derive new, explicit representations for the solution to the scalar wave equation in the exterior of a sphere, subject to either Dirichlet or Robin boundary conditions. Our formula leads to a stable and high-order numerical scheme that permits the evaluation of the solution at an arbitrary target, without the use of a spatial grid and without numerical dispersion error. In the process, we correct some errors in the analytic literature concerning the asymptotic behavior of the logarithmic derivative of the spherical modified Hankel function. We illustrate the performance of the method with several numerical examples.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Breaking of the overall permutation symmetry in nonlinear optical susceptibilities of one-dimensional periodic dimerized Huckel model

Based on infinite one-dimensional single-electron periodic models of trans-polyacetylene, we show analytically that the overall permutation symmetry of nonlinear optical susceptibilities is, albeit preserved in the molecular systems with only bound states, no longer generally held for the periodic systems. The overall permutation symmetry breakdown provides a fairly natural explanation to the widely observed large deviations of Kleinman symmetry for periodic systems in off-resonant regions. Physical conditions to experimentally test the overall permutation symmetry break are discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Size dependence of second-order hyperpolarizability of finite periodic chain under Su-Schrieffer-Heeger model

The second hyperpolarizability $γ_N(-3ωω,ω,ω)$ of $N$ double-bond finite chain of trans-polyactylene is analyzed using the Su-Schrieffer-Heeger model to explain qualitative features of the size-dependence behavior of $γ_N$. Our study shows that $γ_N/N$ is {\it nonmonotonic} with $N$ and that the nonmonotonicity is caused by the dominant contribution of the intraband transition to $γ_N$ in polyenes. Several important physical effects are discussed to reduce quantitative discrepancies between experimental and our results

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

Hyperpolarizabilities for the one-dimensional infinite single-electron periodic systems: I. Analytical solutions under dipole-dipole correlations

The analytical solutions for the general-four-wave-mixing hyperpolarizabilities $χ^{(3)}(-(w_1+w_2+w_3);w_1,w_2,w_3)$ on infinite chains under both Su-Shrieffer-Heeger and Takayama-Lin-Liu-Maki models of trans-polyacetylene are obtained through the scheme of dipole-dipole correlation. Analytical expressions of DC Kerr effect $χ^{(3)}(-w;0,0,w)$, DC-induced second harmonic generation $χ^{(3)}(-2w;0,w,w)$, optical Kerr effect $χ^{(3)}(-w;w,-w,w)$ and DC-electric-field-induced optical rectification $χ^{(3)}(0;w,-w,0)$ are derived. By including or excluding ${\bf \nabla_k}$ terms in the calculations, comparisons show that the intraband contributions dominate the hyperpolarizabilities if they are included. $\nabla_k$ term or intraband transition leads to the break of the overall permutation symmetry in $χ^{(3)}$ even for the low frequency and non-resonant regions. Hence it breaks the Kleinman symmetry that is directly based on the overall permutation symmetry. Our calculations provide a clear understanding of the Kleinman symmetry breaks that are widely observed in many experiments. We also suggest a feasible experiment on $χ^{(3)}$ to test the validity of overall permutation symmetry and our theoretical prediction. Finally, our calculations show the following trends for the various third-order nonlinear optical processes in the low frequency and non-resonant region: $χ^{(3)}(-3w;w,w,w)> χ^{(3)}(-2w;0,w,w)> χ^{(3)}(-w;w,-w,w)>χ^{(3)}(-w; 0,0,w)>= χ^{(3)}(0;w,-w,0)$, and in the resonant region: $χ^{(3)}(-w;0,0,w)> χ^{(3)}(-w;w,-w,w)> χ^{(3)}(-2w;0,w,w)>χ^{(3)}(0;w,-w,0)>χ^{(3)}(-3w;w,w,w)$. (w=ω)

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

Hyperpolarizabilities for the one-dimensional infinite single-electron periodic systems: II. Dipole-dipole versus current-current correlations

Based on Takayama-Lin-Liu-Maki model, analytical expressions for the third-harmonic generation, DC Kerr effect, DC-induced second harmonic optical Kerr effect, optical Kerr effect or intensity-dependent index of refraction and DC-electric-field-induced optical rectification are derived under the static current-current($J_0J_0$) correlation for one-dimensional infinite chains. The results of hyperpolarizabilities under $J_0J_0$ correlation are then compared with those obtained using the dipole-dipole ($DD$) correlation. The comparison shows that the conventional $J_0J_0$ correlation, albeit quite successful for the linear case, is incorrect for studying the nonlinear optical properties of periodic systems.

0 citations0 reviewsNo code linkedOpen access