Source author record

Tao Yin

Tao Yin 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

13works
9topics
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

13 published item(s)

preprint2022arXiv

On the hyper-singular boundary integral equation methods for dynamic poroelasticity: three dimensional case

In our previous work [SIAM J. Sci. Comput. 43(3) (2021) B784-B810], an accurate hyper-singular boundary integral equation method for dynamic poroelasticity in two dimensions has been developed. This work is devoted to studying the more complex and difficult three-dimensional problems with Neumann boundary condition and both the direct and indirect methods are adopted to construct combined boundary integral equations. The strongly-singular and hyper-singular integral operators are reformulated into compositions of weakly-singular integral operators and tangential-derivative operators, which allow us to prove the jump relations associated with the poroelastic layer potentials and boundary integral operators in a simple manner. Relying on both the investigated spectral properties of the strongly-singular operators, which indicate that the corresponding eigenvalues accumulate at three points whose values are only dependent on two Lamé constants, and the spectral properties of the Calderón relations of the poroelasticity, we propose low-GMRES-iteration regularized integral equations. Numerical examples are presented to demonstrate the accuracy and efficiency of the proposed methodology by means of a Chebyshev-based rectangular-polar solver.

preprint2022arXiv

Spectral Galerkin method for solving elastic wave scattering problems with multiple open arcs

We study the elastic time-harmonic wave scattering problems on unbounded domains with boundaries composed of finite collections of disjoints finite open arcs (or cracks) in two dimensions. Specifically, we present a fast spectral Galerkin method for solving the associated weakly- and hyper-singular boundary integral equations (BIEs) arising from Dirichlet and Neumann boundary conditions, respectively. Discretization bases of the resulting BIEs employ weighted Chebyshev polynomials that capture the solutions' edge behavior. We show that these bases guarantee exponential convergence in the polynomial degree when assuming analyticity of sources and arcs geometries. Numerical examples demonstrate the accuracy and robustness of the proposed method with respect to number of arcs and wavenumber.

preprint2020arXiv

A Windowed Green Function method for elastic scattering problems on a half-space

This paper presents a windowed Green function (WGF) method for the numerical solution of problems of elastic scattering by "locally-rough surfaces" (i.e., local perturbations of a half space), under either Dirichlet or Neumann boundary conditions, and in both two and three spatial dimensions. The proposed WGF method relies on an integral-equation formulation based on the free-space Green function, together with smooth operator windowing (based on a "slow-rise" windowing function) and efficient high-order singular-integration methods. The approach avoids the evaluation of the expensive layer Green function for elastic problems on a half-space, and it yields uniformly fast convergence for all incident angles. Numerical experiments for both two and three dimensional problems are presented, demonstrating the accuracy and super-algebraically fast convergence of the proposed method as the window-size grows.

preprint2020arXiv

An accurate hyper-singular boundary integral equation method for dynamic poroelasticity in two dimensions

This paper is concerned with the boundary integral equation method for solving the exterior Neumann boundary value problem of dynamic poroelasticity in two dimensions. The main contribution of this work consists of two aspescts: the proposal of a novel regularized boundary integral equation, and the presentation of new regularized formulations of the strongly-singular and hyper-singular boundary integral operators. Firstly, turning to the spectral properties of the double-layer operator and the corresponding Calderón relation of the poroelasticity, we propose the novel low-GMRES-iteration integral equation whose eigenvalues are bounded away from zero and infinity. Secondly, with the help of the Günter derivatives, we reformulate the strongly-singular and hyper-singular integral operators into combinations of the weakly-singular operators and the tangential derivatives. The accuracy and efficiency of the proposed methodology are demonstrated through several numerical examples.

preprint2020arXiv

Hybrid quantum-classical algorithms for solving quantum chemistry in Hamiltonian-wavefunction space

Variational quantum eigensolver~(VQE) typically optimizes variational parameters in a quantum circuit to prepare eigenstates for a quantum system. Its applications to many problems may involve a group of Hamiltonians, e.g., Hamiltonian of a molecule is a function of nuclear configurations. In this paper, we incorporate derivatives of Hamiltonian into VQE and develop some hybrid quantum-classical algorithms, which explores both Hamiltonian and wavefunction spaces for optimization. Aiming for solving quantum chemistry problems more efficiently, we first propose mutual gradient descent algorithm for geometry optimization by updating parameters of Hamiltonian and wavefunction alternatively, which shows a rapid convergence towards equilibrium structures of molecules. We then establish differential equations that governs how optimized variational parameters of wavefunction change with intrinsic parameters of the Hamiltonian, which can speed up calculation of energy potential surface. Our studies suggest a direction of hybrid quantum-classical algorithm for solving quantum systems more efficiently by considering spaces of both Hamiltonian and wavefunction.

preprint2020arXiv

Inverse conductivity problem with internal data

This paper concerns the reconstruction of a scalar coefficient of a second-order elliptic equation in divergence form posed on a bounded domain from internal data. This theory finds applications in multi-wave imaging, greedy methods to approximate parameter-dependent elliptic problems, and image treatment with partial differential equations. We first show that the inverse problem for smooth coefficients can be rewritten as a linear transport equation. Assuming that the coefficient is known near the boundary, we study the well-posedness of associated transport equation as well as its numerical resolution using discontinuous Galerkin method. We propose a regularized transport equation that allow us to derive rigorous convergence rates of the numerical method in terms of the order of the polynomial approximation as well as the regularization parameter. We finally provide numerical examples for the inversion assuming a lower regularity of the coefficient, and using synthetic data.

preprint2020arXiv

Variational quantum eigensolvers by variance minimization

Variational quantum eigensolver(VQE) typically minimizes energy with hybrid quantum-classical optimization, which aims to find the ground state. Here, we propose a VQE by minimizing energy variance, which is called as variance-VQE(VVQE). The VVQE can be viewed as an self-verifying eigensolver for arbitrary eigenstate by designing, since an eigenstate for a Hamiltonian should have zero energy variance. We demonstrate properties and advantages of VVQE for solving a set of excited states with quantum chemistry problems. Remarkably, we show that optimization of a combination of energy and variance may be more efficient to find low-energy excited states than those of minimizing energy or variance alone. We further reveal that the optimization can be boosted with stochastic gradient descent by Hamiltonian sampling, which uses only a few terms of the Hamiltonian and thus significantly reduces the quantum resource for evaluating variance and its gradients.

preprint2019arXiv

Collective optimization for variational quantum eigensolvers

Variational quantum eigensolver (VQE) optimizes parameterized eigenstates of a Hamiltonian on a quantum processor by updating parameters with a classical computer. Such a hybrid quantum-classical optimization serves as a practical way to leverage up classical algorithms to exploit the power of near-term quantum computing. Here, we develop a hybrid algorithm for VQE, emphasizing the classical side, that can solve a group of related Hamiltonians simultaneously. The algorithm incorporates a snake algorithm into many VQE tasks to collectively optimize variational parameters of different Hamiltonians. Such so-called collective VQEs~(cVQEs) is applied for solving molecules with varied bond length, which is a standard problem in quantum chemistry. Numeral simulations show that cVQE is not only more efficient in convergence, but also trends to avoid single VQE task to be trapped in local minimums. The collective optimization utilizes intrinsic relations between related tasks and may inspire advanced hybrid quantum-classical algorithms for solving practical problems.

preprint2019arXiv

Regularized integral equation methods for elastic scattering problems in three dimensions

This paper presents novel methodologies for the numerical simulation of scattering of elastic waves by both closed and open surfaces in three-dimensional space. The proposed approach utilizes new integral formulations as well as an extension to the elastic context of the efficient high-order singular-integration methods~\cite{BG18} introduced recently for the acoustic case. In order to obtain formulations leading to iterative solvers (GMRES) which converge in small numbers of iterations we investigate, theoretically and computationally, the character of the spectra of various operators associated with the elastic-wave Calderón relation---including some of their possible compositions and combinations. In particular, by relying on the fact that the eigenvalues of the composite operator $NS$ are bounded away from zero and infinity, new uniquely-solvable, low-GMRES-iteration integral formulation for the closed-surface case are presented. The introduction of corresponding low-GMRES-iteration equations for the open-surface equations additionally requires, for both spectral quality as well as accuracy and efficiency, use of weighted versions of the classical integral operators to match the singularity of the unknown density at edges. Several numerical examples demonstrate the accuracy and efficiency of the proposed methodology.

preprint2016arXiv

Boundary integral equation methods for the two dimensional fluid-solid interaction problem

This paper is concerned with boundary integral equation methods for solving the two-dimensional fluid-solid interaction problem. We reduce the problem to three differential systems of boundary integral equations via direct and indirect approaches. Existence and uniqueness results for variational solutions of boundary integral equations equations are established. Since in all these boundary variational formulations, the hypersingular boundary integral operator associated with the time-harmonic Navier equation is a dominated integral operator, we also include a new regularization formulation for this hypersingular operator, which allows us to treat the hypersingular kernel by a wealkly singular kernel. Numerical examples are presented to verify and validate the theoretical results.

preprint2016arXiv

Direct and Inverse Elastic Scattering From Anisotropic Media

Assume a time-harmonic elastic wave is incident onto a penetrable anisotropic body embedded into a homogeneous isotropic background medium. We propose an equivalent variational formulation in a truncated bounded domain and show the uniqueness and existence of weak solutions by applying the Fredholm alternative and using properties of the Dirichlet-to-Neumann map in both two and three dimensions. The Fréchet derivative of the near-field solution operator with respect to the boundary of the scatterer is derived. As an application, we design a descent algorithm for recovering the interface from the near-field data of one or several incident directions and frequencies. Numerical examples in 2D are presented to show the validity and accuracy of our methods.

preprint2015arXiv

Polaronic effects in one- and two-band quantum systems

In this work we study the formation and dynamics of polarons in a system with a few impurities in a lattice immersed in a Bose-Einstein condensate (BEC). This system has been experimentally realized using ultracold atoms and optical lattices. Here we consider a two-band model for the impurity atoms, along with a Bogoliubov approximation for the BEC, with phonons coupled to impurities via both intra- and inter-band transitions. We decouple this Fröhlich-like term by an extended two-band Lang-Firsov polaron transformation using a variational method. The new effective Hamiltonian with two (polaron) bands differs from the original Hamiltonian by modified coherent transport, polaron energy shifts and induced long-range interaction. A Lindblad master equation approach is used to take into account residual incoherent coupling between polaron and bath. This polaronic treatment yields a renormalized inter-band relaxation rate compared to Fermi's Golden Rule. For a strongly coupled two-band Fröhlich Hamiltonian, the polaron is tightly dressed in each band and can not tunnel between them, leading to an inter-band self-trapping effect.

preprint2011arXiv

Stable Heteronuclear Few-Atom Bound States in Mixed Dimensions

We study few-body problems in mixed dimensions with $N \ge 2$ heavy atoms trapped individually in parallel one-dimensional tubes or two-dimensional disks, and a single light atom travels freely in three dimensions. By using the Born-Oppenheimer approximation, we find three- and four-body bound states for a broad region of heavy-light atom scattering length combinations. Specifically, the existence of trimer and tetramer states persist to negative scattering lengths regime, where no two-body bound state is present. These few-body bound states are analogous to the Efimov states in three dimensions, but are stable against three-body recombination due to geometric separation. In addition, we find that the binding energy of the ground trimer and tetramer state reaches its maximum value when the scattering lengths are comparable to the separation between the low-dimensional traps. This resonant behavior is a unique feature for the few-body bound states in mixed dimensions.