Source author record

Hai Zhang

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

24works
15topics
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

24 published item(s)

preprint2025arXiv

Mathematical Theory for Photonic Hall Effect in Honeycomb Photonic Crystals

In this work, we develop a mathematical theory for the photonic Hall effect and prove the existence of guided electromagnetic waves at the interface of two honeycomb photonic crystals. The guided wave resembles the edge states in electronic systems: it is induced by the topological Hall effect, and the wave propagates along the interface but not in the bulk media. Starting from a symmetric honeycomb photonic crystal that attains Dirac points at the high-symmetry points of the Brillouin zone, $K$ and $K'$, we introduce two classes of perturbations for the periodic medium. The perturbations lift the Dirac degeneracy, forming a spectral band valley at the points $K$ and $K'$ with well-defined topological phase that depends on the sign of the perturbation parameters. By employing the layer potential techniques and spectral analysis, we investigate the existence of guided wave along an interface when two honeycomb photonic crystals are glued together. In particular, we elucidate the relationship between the existence of the interface mode and the nature of perturbations imposed on the two periodic media separated by the interface.

preprint2022arXiv

A measurement decoupling based fast algorithm for super-resolving point sources with multi-cluster structure

We consider the problem of resolving closely spaced point sources in one dimension from their Fourier data in a bounded domain. Classical subspace methods (e.g., MUSIC algorithm, Matrix Pencil method, etc.) show great superiority in resolving closely spaced sources, but their computational cost is usually heavy. This is especially the case for point sources with a multi-cluster structure which requires processing large-sized data matrix resulted from highly sampled measurements. To address this issue, we propose a fast algorithm termed D-MUSIC, based on a measurement decoupling strategy. We demonstrate theoretically that for point sources with a known cluster structure, their measurement can be decoupled into local measurements of each of the clusters by solving a system of linear equations that are obtained by using a multipole basis. We further develop a subsampled MUSIC algorithm to detect the cluster structure and utilize it to decouple the global measurement. In the end, the MUSIC algorithm was applied to each local measurement to resolve point sources therein. Compared to the standard MUSIC algorithm, the proposed algorithm has comparable super-resolving capability while having a much lower computational complexity.

preprint2022arXiv

Dirac points for the honeycomb lattice with impenetrable obstacles

This work is concerned with the Dirac points for the honeycomb lattice with impenetrable obstacles arranged periodically in a homogeneous medium. We consider both the Dirichlet and Neumann eigenvalue problems and prove the existence of Dirac points for both eigenvalue problems at crossing of the lower band surfaces as well as higher band surfaces. Furthermore, we perform quantitative analysis for the eigenvalues and the slopes of two conical dispersion surfaces near each Dirac point based on a combination of the layer potential technique and asymptotic analysis. It is shown that the eigenvalues are in the neighborhood of the singular frequencies associated with the Green's function for the honeycomb lattice, and the slopes of the dispersion surfaces are reciprocal to the eigenvalues.

preprint2022arXiv

Mathematical theory for electromagnetic scattering resonances and field enhancement in a subwavelength annular gap

This work presents a mathematical theory for electromagnetic scattering resonances in a subwavelength annular hole embedded in a metallic slab, with the annulus width $h\ll1$. The model is representative among many 3D subwavelength hole structures, which are able to induce resonant scattering of electromagnetic wave and the so-called extraordinary optical transmission. We develop a multiscale framework for the underlying scattering problem based upon a combination of the integral equation in the exterior domain and the waveguide mode expansion inside the tiny hole. The matching of the electromagnetic field over the hole aperture leads to a sequence of decoupled infinite systems, which are used to set up the resonance conditions for the scattering problem. By performing rigorous analysis for the infinite systems and the resonance conditions, we characterize all the resonances in a bounded domain over the complex plane. It is shown that the resonances are associated with the TE and TEM waveguide modes in the annular hole, and they are close to the real axis with the imaginary parts of order ${\cal O}(h)$. We also investigate the resonant scattering when an incident wave is present. It is proved that the electromagnetic field is amplified with order ${\cal O}(1/h)$ at the resonant frequencies that are associated with the TE modes in the annular hole. On the other hand, one particular resonance associated with the TEM mode can not be excited by a plane wave but can be excited with a near-field electric dipole source, leading to field enhancement of order ${\cal O}(1/h)$.

preprint2022arXiv

Randomized Spectral Clustering in Large-Scale Stochastic Block Models

Spectral clustering has been one of the widely used methods for community detection in networks. However, large-scale networks bring computational challenges to the eigenvalue decomposition therein. In this paper, we study the spectral clustering using randomized sketching algorithms from a statistical perspective, where we typically assume the network data are generated from a stochastic block model that is not necessarily of full rank. To do this, we first use the recently developed sketching algorithms to obtain two randomized spectral clustering algorithms, namely, the random projection-based and the random sampling-based spectral clustering. Then we study the theoretical bounds of the resulting algorithms in terms of the approximation error for the population adjacency matrix, the misclassification error, and the estimation error for the link probability matrix. It turns out that, under mild conditions, the randomized spectral clustering algorithms lead to the same theoretical bounds as those of the original spectral clustering algorithm. We also extend the results to degree-corrected stochastic block models. Numerical experiments support our theoretical findings and show the efficiency of randomized methods. A new R package called Rclust is developed and made available to the public.

preprint2022arXiv

Randomized spectral co-clustering for large-scale directed networks

Directed networks are broadly used to represent asymmetric relationships among units. Co-clustering aims to cluster the senders and receivers of directed networks simultaneously. In particular, the well-known spectral clustering algorithm could be modified as the spectral co-clustering to co-cluster directed networks. However, large-scale networks pose great computational challenges to it. In this paper, we leverage sketching techniques and derive two randomized spectral co-clustering algorithms, one \emph{random-projection-based} and the other \emph{random-sampling-based}, to accelerate the co-clustering of large-scale directed networks. We theoretically analyze the resulting algorithms under two generative models -- the stochastic co-block model and the degree-corrected stochastic co-block model, and establish their approximation error rates and misclustering error rates, indicating better bounds than the state-of-the-art results of co-clustering literature. Numerically, we design and conduct simulations to support our theoretical results and test the efficiency of the algorithms on real networks with up to millions of nodes. A publicly available R package \textsf{RandClust} is developed for better usability and reproducibility of the proposed methods.

preprint2021arXiv

A mathematical theory of computational resolution limit in one dimension

Given an image generated by the convolution of point sources with a band-limited function, the deconvolution problem is to reconstruct the source number, positions, and amplitudes. This problem arises from many important applications in imaging and signal processing. It is well-known that it is impossible to resolve the sources when they are close enough in practice. Rayleigh investigated this problem and formulated a resolution limit, the so-called Rayleigh limit, for the case of two sources with identical amplitudes. On the other hand, many numerical experiments demonstrate that a stable recovery of the sources is possible even if the sources are separated below the Rayleigh limit. In this paper, a mathematical theory for the deconvolution problem in one dimension is developed. The theory addresses the issue when one can recover the source number exactly from noisy data. The key is a new concept "computational resolution limit" which is defined to be the minimum separation distance between the sources such that exact recovery of the source number is possible. This new resolution limit is determined by the signal-to-noise ratio and the sparsity of sources, in addition to the cutoff frequency of the image. Quantitative bounds for this limit is derived, which reveal the importance of the sparsity as well as the signal-to-noise ratio to the recovery problem. The stability for recovering the source positions is also analyzed under a condition on their separation distances. Moreover, a singular-value-thresholding algorithm is proposed to recover the source number of a cluster of closely spaced point sources and to verify our theoretical results on the computational resolution limit. The results are based on a multipole expansion method and a non-linear approximation theory in Vandermonde space.

preprint2021arXiv

Differentially Private SGD with Non-Smooth Losses

In this paper, we are concerned with differentially private {stochastic gradient descent (SGD)} algorithms in the setting of stochastic convex optimization (SCO). Most of the existing work requires the loss to be Lipschitz continuous and strongly smooth, and the model parameter to be uniformly bounded. However, these assumptions are restrictive as many popular losses violate these conditions including the hinge loss for SVM, the absolute loss in robust regression, and even the least square loss in an unbounded domain. We significantly relax these restrictive assumptions and establish privacy and generalization (utility) guarantees for private SGD algorithms using output and gradient perturbations associated with non-smooth convex losses. Specifically, the loss function is relaxed to have an $α$-Hölder continuous gradient (referred to as $α$-Hölder smoothness) which instantiates the Lipschitz continuity ($α=0$) and the strong smoothness ($α=1$). We prove that noisy SGD with $α$-Hölder smooth losses using gradient perturbation can guarantee $(ε,δ)$-differential privacy (DP) and attain optimal excess population risk $\mathcal{O}\Big(\frac{\sqrt{d\log(1/δ)}}{nε}+\frac{1}{\sqrt{n}}\Big)$, up to logarithmic terms, with the gradient complexity $ \mathcal{O}( n^{2-α\over 1+α}+ n).$ This shows an important trade-off between $α$-Hölder smoothness of the loss and the computational complexity for private SGD with statistically optimal performance. In particular, our results indicate that $α$-Hölder smoothness with $α\ge {1/2}$ is sufficient to guarantee $(ε,δ)$-DP of noisy SGD algorithms while achieving optimal excess risk with the linear gradient complexity $\mathcal{O}(n).$

preprint2020arXiv

A mathematical theory for Fano resonance in a periodic array of narrow slits

This work concerns resonant scattering by a perfectly conducting slab with periodically arranged subwavelength slits, with two slits per period. There are two classes of resonances, corresponding to poles of a scattering problem. A sequence of resonances has an imaginary part that is nonzero and on the order of the width $\varepsilon$ of the slits; these are associated with Fabry-Perot resonance, with field enhancement of order $1/\varepsilon$ in the slits. The focus of this study is another class of resonances which become real valued at normal incidence, when the Bloch wavenumber $κ$ is zero. These are embedded eigenvalues of the scattering operator restricted to a period cell, and the associated eigenfunctions extend to surface waves of the slab that lie within the radiation continuum. When $0<|κ|\ll 1$, the real embedded eigenvalues will be perturbed as complex-valued resonances, which induce the Fano resonance phenomenon. We derive the asymptotic expansions of embedded eigenvalues and their perturbations as resonances when the Bloch wavenumber becomes nonzero. Based on the quantitative analysis of the diffracted field, we prove that the Fano-type anomalies occurs for the transmission of energy through the slab, and show that the field enhancement is of order $1/(κ\varepsilon)$, which is stronger than Fabry-Perot resonance.

preprint2020arXiv

Fano resonance in metallic grating via strongly coupled subwavelength resonators

We investigate the Fano resonance in grating structures by using coupled resonators. The grating consists of a perfectly conducting slab with periodically arranged subwavelength slit holes, where inside each period, a pair of slits sit very close to each other. The slit holes act as resonators and are strongly coupled. It is shown rigorously that there exist two groups of resonances corresponding to poles of the scattering problem. One sequence of resonances have imaginary part on the order of $\varepsilon$, where $\varepsilon$ is the size of the slit aperture, while the other sequence have imaginary part on the order of $\varepsilon^2$. When coupled with the incident wave at resonant frequencies, the narrow-band resonant scattering induced by the latter will interfere with the broader background resonant radiation induced by the former. The interference of these two resonances generates the Fano type transmission anomaly, which persists in the whole radiation continuum of the grating structure as long as the slit aperture size is small compared to the incident wavelength.

preprint2020arXiv

Sensitivity of resonance frequency in the detection of thin layer using nano-slit structures

We derive the formulas for the resonance frequencies and their sensitivity when the nano-slit structures are used in the detection of thin layers. For a thin layer with a thickness of $H$ deposited over the nanostructure, we show quantitatively that for both single and periodic slit structures with slit aperture size $δ$, the sensitivity of resonance frequency reduces as $H$ increases. Specifically, the sensitivity is of order $O(δ/H)$ if $H >δ$ and of order $O(1+\ln H/δ)$ otherwise. The evanescent wave modes are present along the interface between the thin dielectric film and ambient medium above. From the mathematical derivations, it is observed that the sensitivity of the resonance frequency highly depends on the effect of evanescent wave modes on the tiny slit apertures.

preprint2016arXiv

A mathematical and numerical framework for bubble meta-screens

The aim of this paper is to provide a mathematical and numerical framework for the analysis and design of bubble meta-screens. An acoustic meta-screen is a thin sheet with patterned subwavelength structures, which nevertheless has a macroscopic effect on the acoustic wave propagation. In this paper, periodic subwavelength bubbles mounted on a reflective surface (with Dirichlet boundary condition) is considered. It is shown that the structure behaves as an equivalent surface with Neumann boundary condition at the Minnaert resonant frequency which corresponds to a wavelength much greater than the size of the bubbles. Analytical formula for this resonance is derived. Numerical simulations confirm its accuracy and show how it depends on the ratio between the periodicity of the lattice, the size of the bubble, and the distance from the reflective surface. The results of this paper formally explain the super-absorption behavior observed in [V. Leroy et al., Phys. Rev. B, 2015].

preprint2016arXiv

Effective medium theory for acoustic waves in bubbly fluids near Minnaert resonant frequency

We derive an effective medium theory for acoustic wave propagation in bubbly fluid near Minnaert resonant frequency. We start with a multiple scattering formulation of the scattering problem of an incident wave by a large number of identical small bubbles in a homogeneous fluid. Under certain conditions on the configuration of the bubbles, we justify the point interaction approximation and establish an effective medium theory for the bubbly fluid as the number of bubbles tends to infinity. The convergence rate is also derived. As a consequence, we show that near and below the Minnaert resonant frequency, the obtained effective media can have a high refractive index, which is the reason for the super-focusing experiment observed in "Subwavelength focusing in bubbly media using broadband time reversal", Physical Review, B (2015), by M. Lanoy et al. Moreover, our results indicate that the obtained effective medium can be dissipative above the Minnaert resonant frequency, while at that frequency, effective medium theory does not hold. Our theory sheds light on the mechanism of the extraordinary wave properties of metamaterials, which include bubbly fluid as an example, near resonant frequencies.

preprint2016arXiv

Minnaert resonances for acoustic waves in bubbly media

Through the application of layer potential techniques and Gohberg-Sigal theory we derive an original formula for the Minnaert resonance frequencies of arbitrarily shaped bubbles. We also provide a mathematical justification for the monopole approximation of scattering of acoustic waves by bubbles at their Minnaert resonant frequency. Our results are complemented by several numerical examples which serve to validate our formula in two dimensions.

preprint2016arXiv

Shape reconstruction of nanoparticles from their associated plasmonic resonances

We prove by means of a couple of examples that plasmonic resonances can be used on one hand to classify shapes of nanoparticles with real algebraic boundaries and on the other hand to reconstruct the separation distance between two nanoparticles from measurements of their first collective plasmonic resonances. To this end, we explicitly compute the spectral decompositions of the Neumann-Poincaré operators associated with a class of quadrature domains and two nearly touching disks. Numerical results are included in support of our main findings.

preprint2016arXiv

Theory of plasmonic metasurfaces

In this paper we derive an impedance boundary condition to approximate the optical scattering effect of an array of plasmonic nanoparticles mounted on a perfectly conducting plate. We show that at some resonant frequencies the impedance blows up, allowing for a significant reduction of the scattering from the plate. Using the spectral properties of a Neumann-Poincare type operator, we investigate the dependency of the impedance with respect to changes in the nanoparticle geometry and configuration.

preprint2015arXiv

Mathematical analysis of plasmonic nanoparticles: the scalar case

Localized surface plasmons are charge density oscillations confined to metallic nanoparticles. Excitation of localized surface plasmons by an electromagnetic field at an incident wavelength where resonance occurs results in a strong light scattering and an enhancement of the local electromagnetic fields. This paper is devoted to the mathematical modeling of plasmonic nanoparticles. Its aim is threefold: (i) to mathematically define the notion of plasmonic resonance and to analyze the shift and broadening of the plasmon resonance with changes in size and shape of the nanoparticles; (ii) to study the scattering and absorption enhancements by plasmon resonant nanoparticles and express them in terms of the polarization tensor of the nanoparticle. Optimal bounds on the enhancement factors are also derived; (iii) to show, by analyzing the imaginary part of the Green function, that one can achieve super-resolution and super-focusing using plasmonic nanoparticles. For simplicity, the Helmholtz equation is used to model electromagnetic wave propagation.

preprint2015arXiv

Mathematical analysis of plasmonic resonances for nanoparticles: the full Maxwell equations

In this paper we use the full Maxwell equations for light propagation in order to analyze plasmonic resonances for nanoparticles. We mathematically define the notion of plasmonic resonance and analyze its shift and broadening with respect to changes in size, shape, and arrangement of the nanoparticles, using the layer potential techniques associated with the full Maxwell equations. We present an effective medium theory for resonant plasmonic systems and derive a condition on the volume fraction under which the Maxwell-Garnett theory is valid at plasmonic resonances.

preprint2014arXiv

A mathematical theory of super-resolution by using a system of sub-wavelength Helmholtz resonators

A rigorous mathematical theory is developed to explain the super-resolution phenomenon observed in the experiment by F.Lemoult, M.Fink and G.Lerosey (Acoustic resonators for far-field control of sound on a subwavelength scale, Phys. Rev. Lett., 107 (2011)). A key ingredient is the calculation of the resonances and the Green function in the half space with the presence of a system of Helmholtz resonators in the quasi-stationary regime. By using boundary integral equations and generalized Rouche's theorem, the existence and the leading asymptotic of the resonances are rigorously derived. The integral equation formulation also yields the leading order terms in the asymptotics of the Green function. The methodology developed in the paper provides an elegant and systematic way for calculating resonant frequencies for Helmholtz resonators in assorted space settings as well as in various frequency regimes. By using the asymptotics of the Green function, the analysis of the imaging functional of the time-reversal wave fields becomes possible, which clearly demonstrates the super-resolution property. The result provides the first mathematical theory of super-resolution in the context of a deterministic medium and sheds light to the mechanism of super-resolution and super-focusing for waves in deterministic complex media.

preprint2014arXiv

Sensitivity Analysis of an Inverse Problem for the Wave Equation with Caustics

The paper investigates the sensitivity of the inverse problem of recovering the velocity field in a bounded domain from the boundary dynamic Dirichlet-to-Neumann map (DDtN) for the wave equation. Three main results are obtained: (1) assuming that two velocity fields are non-trapping and are equal to a constant near the boundary, it is shown that the two induced scattering relations must be identical if their corresponding DDtN maps are sufficiently close; (2) a geodesic X-ray transform operator with matrix-valued weight is introduced by linearizing the operator which associates each velocity field with its induced Hamiltonian flow. A selected set of geodesics whose conormal bundle can cover the cotangent space at an interior point is used to recover the singularity of the X-ray transformed function at the point; a local stability estimate is established for this case. Although fold caustics are allowed along these geodesics, it is required that these caustics contribute to a smoother term in the transform than the point itself. The existence of such a set of geodesics is guaranteed under some natural assumptions in dimension greater than or equal to three by the classification result on caustics and regularity theory of Fourier Integral Operators. The interior point with the above required set of geodesics is called "fold-regular"; (3) assuming that a background velocity field with every interior point fold-regular is fixed and another velocity field is sufficiently close to it and satisfies a certain orthogonality condition, it is shown that if the two corresponding DDtN maps are sufficiently close then they must be equal.

preprint2014arXiv

Stability Analysis for Magnetic Resonance Elastography

We consider the inverse problem of finding unknown elastic parameters from internal measurements of displacement fields for tissues. The measurements are made on the entirety of a smooth domain. Since tissues can be modeled as quasi-incompressible fluids, we examine the Stokes system and consider only the recovery of shear modulus distributions. Our main result is to establish Lipschitz stable estimates on the shear modulus distributions from internal measurements of displacement fields. These estimates imply convergence of a numerical scheme known as the Landweber iteration scheme for reconstructing the shear modulus distributions.

preprint2013arXiv

Architecture of Network and Client-Server model

With the development of Internet technology, the Web is becoming more and more important in our lives so that it has even become an essential element. The application of the Web has never been limited to computers; it has been opened to all kinds of intelligent digital devices like mobile ones. In this paper, we discuss the main and substantial difference between UDP and TCP, and how to implement Client-Server model. We also discuss the efficiency of multi-thread server and its relationship with internet.

preprint2013arXiv

Asymptotic normality of maximum likelihood and its variational approximation for stochastic blockmodels

Variational methods for parameter estimation are an active research area, potentially offering computationally tractable heuristics with theoretical performance bounds. We build on recent work that applies such methods to network data, and establish asymptotic normality rates for parameter estimates of stochastic blockmodel data, by either maximum likelihood or variational estimation. The result also applies to various sub-models of the stochastic blockmodel found in the literature.

preprint2013arXiv

The Economic Trend of Video Game Industry

In recent years the game industry has had a huge growth. We've seen new game consoles, great looking games and an increase in the number of people playing them. We are presently in the seventh generation of video games which focuses on consoles released since 2004. For home consoles,the seventh generation began on November 22, 2005 with the release of Xbox 360 and continued with the release of PlayStation 3 on November 11, 2006, and Wii on November 19, 2006. The current generation is having a console battle between Nintendo's Wii, Microsoft's Xbox 360, and Sony's PlayStation 3.The appearance of the three new consoles not only offers various purchase choices, but also greatly affects economy and culture.