A hybrid method integrating Green's function Monte Carlo and projected entangled pair states

This paper introduces a hybrid approach combining Green’s function Monte Carlo(GFMC)method with projected entangled pair state(PEPS)ansatz.This hybrid method regards PEPS as a trial state and a guiding wave function in GFMC.By leveraging PEPS’s proficiency in capturing quantum state entanglement and GFMC’s efficient parallel architecture,the hybrid method is well-suited for the accurate and efficient treatment of frustrated quantum spin systems.As a benchmark,we applied this approach to study the frustrated J_(1)–J_(2) Heisenberg model on a square lattice with periodic boundary conditions(PBCs).Compared with other numerical methods,our approach integrating PEPS and GFMC shows competitive accuracy in the performance of ground-state energy.This paper provides systematic and comprehensive discussion of the approach of our previous work[Phys.Rev.B 109235133(2024)].

Implementation of the Integrated Green’s Function Method for 3D Poisson’s Equation in a Large Aspect Ratio Computational Domain

The solution of Poisson’s Equation plays an important role in many areas, including modeling high-intensity and high-brightness beams in particle accelerators. For the computational domain with a large aspect ratio, the integrated Green’s function method has been adopted to solve the 3D Poisson equation subject to open boundary conditions. In this paper, we report on the efficient implementation of this method, which can save more than a factor of 50 computing time compared with the direct brute force implementation and its improvement under certain extreme conditions.

Analytical approach to quantum phase transitions of ultracold Bose gases in bipartite optical lattices using the generalized Green＇s function method

In order to investigate the quantum phase transitions and the time-of-flight absorption pictures analyt- ically in a systematic way for ultracold Bose gases in bipartite optical lattices, we present a generalized Green＇s function method. Utilizing this method, we study the quantum phase transitions of ultracold Bose gases in two types of bipartite optical lattices, i.e., a hexagonal lattice with normal Bose-Hubbard interaction and a d-dimensional hypercubic optical lattice with extended Bose-Hubbard interaction. Furthermore, the time-of-flight absorption pictures of ultracold Bose gases in these two types of lat- tices are also calculated analytically. In hexagonal lattice, the time-of-flight interference patterns of ultracold Bose gases obtained by our analytical method are in good qualitative agreement with the exDerimental results of Soltan-Panahi, et al. [Nat. Phys. 7, 434 （2011）]. In square optical lattice, the emergence of peaks at（±π/a,±π/a） in the time-of-flight absorption pictures, which is believed to bea sort of evidence of the existence of a supersolid phase, is clearly seen when the system enters the compressible phase from charge-density-wave phase.

Study on Characteristics of 3-D Translating-Pulsating Source Green Function of Deep-Water Havelock Form and Its Fast Integration Method

The singularities, oscillatory performances and the contributing factors to the 3-＇D translating-pulsating source Green function of deep-water Havelock form which consists of a local disturbance part and a far-field wave-like part, are analyzed systematically. Relative numerical integral methods about the two parts are presented in this paper. An improved method based on LOBATTO rule is used to eliminate singularities caused respectively by infinite discontinuity and jump discontinuous node from the local disturbance part function, which makes the improvement of calculation efficiency and accuracy possible. And variable substitution is applied to remove the singularity existing at the end of the integral interval of the far-field wave-like part function. Two auxiliary techniques such as valid interval calculation and local refinement of integral steps technique in narrow zones near false singularities are applied so as to avoid unnecessary integration of invalid interval and improve integral accordance. Numerical test results have proved the efficiency and accuracy in these integral methods that thus can be applied to calculate hydrodynamic performance of floating structures moving in waves.

Existence of Monotone Positive Solution for a Fourth-Order Three-Point BVP with Sign-Changing Green’s Function

This paper is concerned with the following fourth-order three-point boundary value problem , where , we discuss the existence of positive solutions to the above problem by applying to the fixed point theory in cones and iterative technique.

Nonequilibrium Green＇s function method for quantum thermal transport

This review deals with the nonequilibrium Green＇s function （NEGF） method applied to the problems of energy transport due to atomic vibrations （phonons）, primarily for small jtmction systems. We present a pedagogical introduction to the subject, deriving some of the well-known results such as the Laudauer-like formula for heat current in ballistic systerms. The main aim of the review is to build the machinery of the method so that it can be applied to other situations, which are not directly treated here. In addition to the above, we consider a nmnber of applications of NEGF, not in routine model system calculations, but in a few new aspects showing the power and usefulness of the formalism. In partkaflar, we discuss the problems of multiple leads, coupled left-right-lead system, and system without a center. We also apply the method to the problem of full counting statisties. In the case of nonlinear svstems, we make general comments on the thermal expansion effect. phonon relaxation timv. and a certain class of mean-field approximations. Lastly, we examine the relationship between NEGF. reduced density matrix, and master equation approaches to thermal transport,

Continuum Skyrme Hartree-Fock-Bogoliubov theory with Green’s function method for neutron-rich Ca,Ni,Zr,and Sn isotopes

The possible exotic nuclear properties in the neutron-rich Ca,Ni,Zr,and Sn isotopes are examined with the continuum Skyrme Hartree-Fock-Bogoliubov theory in the framework of the Green’s function method.The pairing correlation,the couplings with the continuum,and the blocking effects for the unpaired nucleon in odd-A nuclei are properly treated.The Skyrme interaction SLy4 is adopted for the ph channel and the density-dependentinteraction is adopted for the pp chan-nel,which well reproduce the experimental two-neutron separation energies S_(2n)and one-neutron separation energies Sn.It is found that the criterion S_(n)>0 predicts a neutron drip line with neutron numbers much smaller than those for S_(2n)>0.Owing to the unpaired odd neutron,the neutron pairing energies−E_(pair)in odd-A nuclei are much lower than those in the neighbor-ing even-even nuclei.By investigating the single-particle structures,the possible halo structures in the neutron-rich Ca,Ni,and Sn isotopes are predicted,where sharp increases in the root-mean-square(rms)radii with significant deviations from the traditional rA^(1∕3)rule and diffuse spatial density distributions are observed.Analyzing the contributions of various partial waves to the total neutron densityρlj(r)∕ρ(r)reveals that the orbitals located around the Fermi surface-particularly those with small angular momenta-significantly affect the extended nuclear density and large rms radii.The number of neutrons Nλ(N_(0))occupying above the Fermi surfacen(continuum threshold)is discussed,whose evolution as a function of the mass number A in each isotope is consistent with that of the pairing energy,supporting the key role of the pairing correlation in halo phenomena.

Searching for single-particle resonances with the Green’s function method

Single-particle resonances in the continuum are crucial for studies of exotic nuclei.In this study,the Green’s function approach is employed to search for single-particle resonances based on the relativistic-mean-field model.Taking^(120)Sn as an example,we identify singleparticle resonances and determine the energies and widths directly by probing the extrema of the Green’s functions.In contrast to the results found by exploring for the extremum of the density of states proposed in our recent study[Chin.Phys.C,44:084105(2020)],which has proven to be very successful,the same resonances as well as very close energies and widths are obtained.By comparing the Green’s functions plotted in different coordinate space sizes,we also found that the results very slightly depend on the space size.These findings demonstrate that the approach by exploring for the extremum of the Green’s function is also very reliable and effective for identifying resonant states,regardless of whether they are wide or narrow.

Adomian Decomposition Method with Green’s Function for Solving Tenth-Order Boundary Value Problems

In this paper, the Adomian decomposition method with Green’s function (Standard Adomian and Modified Technique) is applied to solve linear and nonlinear tenth-order boundary value problems with boundary conditions defined at any order derivatives. The numerical results obtained with a small amount of computation are compared with the exact solutions to show the efficiency of the method. The results show that the decomposition method is of high accuracy, more convenient and efficient for solving high-order boundary value problems.

METHOD OF GREEN’S FUNCTION OF CORRUGATED SHELLS

By using the fundamental equations of axisymmetric shallow shells of revolution, the nonlinear bending of a shallow corrugated shell with taper under arbitrary load has been investigated. The nonlinear boundary value problem of the corrugated shell was reduced to the nonlinear integral equations by using the method of Green＇s function. To solve the integral equations, expansion method was used to obtain Green＇s function. Then the integral equations were reduced to the form with degenerate core by expanding Green＇s function as series of characteristic function. Therefore, the integral equations become nonlinear algebraic equations. Newton＇ s iterative method was utilized to solve the nonlinear algebraic equations. To guarantee the convergence of the iterative method, deflection at center was taken as control parameter. Corresponding loads were obtained by increasing deflection one by one. As a numerical example,elastic characteristic of shallow corrugated shells with spherical taper was studied.Calculation results show that characteristic of corrugated shells changes remarkably. The snapping instability which is analogous to shallow spherical shells occurs with increasing load if the taper is relatively large. The solution is close to the experimental results.

A NEW METHOD FOR SOLVING DYADIC GREEN'S FUNCTION OF ELECTROMAGNETIC FIELD

A new method for solving electromagnetic field boundary value problem is given.Byusing this method,the boundary value problem of the vector wave equation can be transformedinto the independent boundary value problem of scalar wave equations and the two additionalvector differential operations.All the dyadic Green’s functions got by eigenfunction expansionof the dyadic Green’s function can be got by this method easily and some of the dyadic Green’sfunctions for complex systems which are very difficult to get by the ordinary method have beengot by this new method.The dyadic Green’s function for a dielectric loaded cavity is one of thegiven examples.

Quasi-Green’s function method for free vibration of clamped thin plates on Winkler foundation

The quasi-Green＇s function method is used to solve the free vibration problem of clamped thin plates on the Winkler foundation. Quasi-Green＇s function is established by the fundamental solution and the boundary equation of the problem. The function satisfies the homogeneous boundary condition of tile problem. The mode-shape differential equation of the free vibration problem of clamped thin plates on the Winkler foundation is reduced to the Fredholm integral equation of the second kind by Green＇s formula. The irregularity of the kernel of the integral equation is overcome by choosing a suitable form of the normalized boundary equation. The numerical results show the high accuracy of the proposed method.

Numerical Green’s function method:Application to quantifying ground motion variations of M7 earthquakes

The concept of ＂numerical Green’s functions＂ （NGF or Green’s function database） is developed. The basic idea is： a large seismic fault is divided into subfaults of appropriate size, for which synthetic Green’s functions at the surface （NGF） are calculated and stored. Consequently, ground motions from arbitrary kinematic sources can be simulated, rapidly, for the whole fault or parts of it by superposition. The target fault is a simplified, vertical model of the Newport-Inglewood fault in the Los Angeles basin. This approach and its functionality are illustrated by investigating the variations of ground motions （e.g. peak ground velocity and synthetic seismograms） due to the source complexity. The source complexities are considered with two respects： hypocenter location and slip history. The results show a complex behavior, with dependence of absolute peak ground velocity and their variation on source process directionality, hypocenter location, local structure, and static slip asperity location. We concluded that combining effect due to 3-D structure and finite-source is necessary to quan- tify ground motion characteristics and their variations. Our results will facilitate the earthquake hazard assessment projects.

Singularity-free Green's function for EM sources embedded in a stratified medium

We present a method to unify the calculation of Green＇s functions for an electromagnetic（EM） transmitting source embedded in a homogeneous stratified medium.A virtual interface parallel to layer interfaces is introduced through the source location.The potentials for Green＇s function are derived by decomposing the partial wave solutions to Helmholtz＇s equations into upward and downward within boundaries.The amplitudes of the potentials in each stratum are obtained recursively from the initial amplitudes at the source level.The initial amplitudes are derived by coupling with the transmitting sources and following the discontinuity of the tangential electric and magnetic fields at the source interface.Only the initial terms are related to the transmitting sources and thus need to be modified for different transmitters,whereas the kernel connected with the stratified media stays unchanged.Hence,the present method can be easily applied to EM transmitting sources with little modification.The application of the proposed method to the marine controlled-source electromagnetic method（MCSEM） demonstrates its simplicity and flexibility.

The Green's function of the harmonic horizontal force applied at the interior of the saturated half-space soil

By using integral transform methods, the Green(s functions of horizontal harmonic force applied at the interior of the saturated half-space soil are obtained in the paper. The general solutions of the Biot dynamic equations in frequency domain are established through the use of Hankel integral transforms technique. Utilizing the above- mentioned general solutions, and the boundary conditions of the surface of the half-space and the continuous conditions at the plane of the horizontal force, the solutions of the boundary value problem can be determined. By the numerical inverse Hankel transforms method, the Green(s functions of the harmonic horizontal force are obtainable. The degenerate case of the results deduced from this paper agrees well with the known results. Two numerical examples are given in the paper.

纳米分子器件量子传导中的Green's Function模型

将纳米分子器件的电极和分子本体作为一个系统,考虑到入射电子在原子电极和分子本体中的传导,在电极与分子之间界面上的透射的量子状态,在满足Fisher-Lee关系式的基础上,一个新型的Green’sFunction被设计出来。这一理论模型把三个区域上的不同量子状态完整地表达在一个方程式里,符合物理学基本原理,可以系统地用来研究纳米分子器件的量子传导特性。

Green's functions for uniformly distributed loads acting on an inclined line in a poroelastic layered site

Based on one type of practical Biot＇s equation and the dynamic-stiffness matrices ofa poroelastic soil layer and half-space, Green＇s functions were derived for unitformly distributed loads acting on an inclined line in a poroelastie layered site. This analysis overcomes significant problems in wave scattering due to local soil conditions and dynamic soil-structure interaction. The Green＇s functions can be reduced to the case of an elastic layered site developed by Wolf in 1985. Parametric studies are then carried out through two example problems.

In-plane dynamic Green's functions for inclined and uniformly distributed loads in a multi-layered transversely isotropic half-space

The dynamic stiffness method combined with the Fourier transform is utilized to derive the in-plane Green’s functions for inclined and uniformly distributed loads in a multi-layered transversely isotropic（TI）half-space.The loaded layer is fixed to obtain solutions restricted in it and the corresponding reactions forces,which are then applied to the total system with the opposite sign.By adding solutions restricted in the loaded layer to solutions from the reaction forces,the global solutions in the wavenumber domain are obtained,and the dynamic Green’s functions in the space domain are recovered by the inverse Fourier transform.The presented formulations can be reduced to the isotropic case developed by Wolf（1985）,and are further verified by comparisons with existing solutions in a uniform isotropic as well as a layered TI halfspace subjected to horizontally distributed loads which are special cases of the more general problem addressed.The deduced Green’s functions,in conjunction with boundary element methods,will lead to significant advances in the investigation of a variety of wave scattering,wave radiation and soil-structure interaction problems in a layered TI site.Selected numerical results are given to investigate the influence of material anisotropy,frequency of excitation,inclination angle and layered on the responses of displacement and stress,and some conclusions are drawn.

Green's function of multi-layered poroelastic half-space for models of ground vibration due to railway traffic

This study proposes a Green＇s function, an essential representation of water-saturated ground under moving excitation, to simulate ground borne vibration from trains. First, general solutions to the governing equations of poroelastic medium are derived by means of integral transform. Secondly, the transmission and reflection matrix approach is used to formulate the relationship between displacement and stress of the stratified ground, which results in the matrix of the Green＇s function. Then the Green＇s function is combined into a train-track-ground model, and is verified by typical examples and a field test. Additional simulations show that the computed ground vibration attenuates faster in the immediate vicinity of the track than in the surrounding area. The wavelength of wheel-rail unevenness has a notable effect on computed displacement and pore pressure. The variation of vibration intensity with the depth of ground is significantly influenced by the layering of the strata soil. When the train speed is equal to the velocity of the Rayleigh wave, the Mach cone appears in the simulated wave field. The proposed Green＇s function is an appropriate representation for a layered ground with shallow ground water table, and will be helpful to understand the dynamic responses of the ground to complicated moving excitation.

NONLINEAR EVOLUTION SYSTEMS AND GREEN’S FUNCTION

In this paper, we will introduce how to apply Green＇s function method to get the pointwise estimates for the solutions of Cauchy problem of nonlinear evolution equations with dissipative structure. First of all, we introduce the pointwise estimates of the time-asymptotic shape of the solutions of the isentropic Navier-Stokes equations and show to exhibit the generalized Huygen＇s principle. Then, for other nonlinear dissipative evolution equations, we will only introduce the result and give some brief explanations. Our approach is based on the detailed analysis of the Green＇s function of the linearized system and micro-local analysis, such as frequency decomposition and so on.