Random Green's Function Method for Large-Scale Electronic Structure Calculation

We report a linear-scaling random Green's function(rGF) method for large-scale electronic structure calculation. In this method, the rGF is defined on a set of random states and is efficiently calculated by projecting onto Krylov subspace. With the rGF method, the Fermi–Dirac operator can be obtained directly, avoiding the polynomial expansion to Fermi–Dirac function. To demonstrate the applicability, we implement the rGF method with the density-functional tight-binding method. It is shown that the Krylov subspace can maintain at small size for materials with different gaps at zero temperature, including H_(2)O and Si clusters. We find with a simple deflation technique that the rGF self-consistent calculation of H_(2)O clusters at T = 0 K can reach an error of~ 1 me V per H_(2)O molecule in total energy, compared to deterministic calculations. The rGF method provides an effective stochastic method for large-scale electronic structure simulation.

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.

Dynamic Characteristics of Functionally Graded Timoshenko Beams by Improved Differential Quadrature Method

This study proposes an effective method to enhance the accuracy of the Differential Quadrature Method(DQM)for calculating the dynamic characteristics of functionally graded beams by improving the form of discrete node distribution.Firstly,based on the first-order shear deformation theory,the governing equation of free vibration of a functionally graded beam is transformed into the eigenvalue problem of ordinary differential equations with respect to beam axial displacement,transverse displacement,and cross-sectional rotation angle by considering the effects of shear deformation and rotational inertia of the beam cross-section.Then,ignoring the shear deformation of the beam section and only considering the effect of the rotational inertia of the section,the governing equation of the beam is transformed into the eigenvalue problem of ordinary differential equations with respect to beam transverse displacement.Based on the differential quadrature method theory,the eigenvalue problem of ordinary differential equations is transformed into the eigenvalue problem of standard generalized algebraic equations.Finally,the first several natural frequencies of the beam can be calculated.The feasibility and accuracy of the improved DQM are verified using the finite element method(FEM)and combined with the results of relevant literature.

A Hybrid Level Set Optimization Design Method of Functionally Graded Cellular Structures Considering Connectivity

With the continuous advancement in topology optimization and additive manufacturing(AM)technology,the capability to fabricate functionally graded materials and intricate cellular structures with spatially varying microstructures has grown significantly.However,a critical challenge is encountered in the design of these structures–the absence of robust interface connections between adjacent microstructures,potentially resulting in diminished efficiency or macroscopic failure.A Hybrid Level Set Method(HLSM)is proposed,specifically designed to enhance connectivity among non-uniform microstructures,contributing to the design of functionally graded cellular structures.The HLSM introduces a pioneering algorithm for effectively blending heterogeneous microstructure interfaces.Initially,an interpolation algorithm is presented to construct transition microstructures seamlessly connected on both sides.Subsequently,the algorithm enables the morphing of non-uniform unit cells to seamlessly adapt to interconnected adjacent microstructures.The method,seamlessly integrated into a multi-scale topology optimization framework using the level set method,exhibits its efficacy through numerical examples,showcasing its prowess in optimizing 2D and 3D functionally graded materials(FGM)and multi-scale topology optimization.In essence,the pressing issue of interface connections in complex structure design is not only addressed but also a robust methodology is introduced,substantiated by numerical evidence,advancing optimization capabilities in the realm of functionally graded materials and cellular structures.

Transfer matrix method for free and forced vibrations of multi-level functionally graded material stepped beams with different boundary conditions

Functionally graded materials(FGMs)are a novel class of composite materials that have attracted significant attention in the field of engineering due to their unique mechanical properties.This study aims to explore the dynamic behaviors of an FGM stepped beam with different boundary conditions based on an efficient solving method.Under the assumptions of the Euler-Bernoulli beam theory,the governing differential equations of an individual FGM beam are derived with Hamilton’s principle and decoupled via the separation-of-variable approach.Then,the free and forced vibrations of the FGM stepped beam are solved with the transfer matrix method(TMM).Two models,i.e.,a three-level FGM stepped beam and a five-level FGM stepped beam,are considered,and their natural frequencies and mode shapes are presented.To demonstrate the validity of the method in this paper,the simulation results by ABAQUS are also given.On this basis,the detailed parametric analyses on the frequencies and dynamic responses of the three-level FGM stepped beam are carried out.The results show the accuracy and efficiency of the TMM.

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.

A Radial Basis Function Method with Improved Accuracy for Fourth Order Boundary Value Problems

Accurately approximating higher order derivatives is an inherently difficult problem. It is shown that a random variable shape parameter strategy can improve the accuracy of approximating higher order derivatives with Radial Basis Function methods. The method is used to solve fourth order boundary value problems. The use and location of ghost points are examined in order to enforce the extra boundary conditions that are necessary to make a fourth-order problem well posed. The use of ghost points versus solving an overdetermined linear system via least squares is studied. For a general fourth-order boundary value problem, the recommended approach is to either use one of two novel sets of ghost centers introduced here or else to use a least squares approach. When using either ghost centers or least squares, the random variable shape parameter strategy results in significantly better accuracy than when a constant shape parameter is used.

水禽细小病毒SYBR Green Ⅰ荧光定量PCR检测方法的建立与应用

为了建立水禽细小病毒(WPV)快速检测方法,根据序列比对结果在水禽细小病毒NS基因SF3保守区域内设计特异性引物,建立SYBR Green Ⅰ荧光定量PCR通用检测方法。该方法的扩增效率(E)为90.0%,相关系数(R~2)=0.99,标准曲线方程为y=-3.607x+38.77;除WPV出现S形扩增曲线外,新城疫病毒(NDV)、H9亚型禽流感病毒(H9 AIV)、鸭坦布苏病毒(DTMUV)、鸭肝炎病毒(DHAV)、鸭肠炎病毒(DEV)、鸭呼肠孤病毒(DRV)样品均未出现S形阳性扩增曲线;批内变异系数(CV)为0.15%~0.23%,批间变异系数为0.09%~0.28%。结果表明,SYBR Green Ⅰ荧光定量PCR检测方法重复性好、灵敏度高和特异性强。临床样品检测结果表明,SYBR Green Ⅰ荧光定量PCR与普通PCR的符合率达98.4%,灵敏度是普通PCR的1 000倍。SYBR Green Ⅰ荧光定量PCR检测方法不仅能定性检测WPV,还可以进行定量检测,可用于种鸭场、种鹅场的WPV净化检测,也可用于WPV临床大量样品的快速检测。

One-step green method to prepare progressive burning gun propellant through gradient denitration strategy

Gradiently denitrated gun propellant(GDGP)prepared by a“gradient denitration”strategy is obviously superior in progressive burning performance to the traditional deterred gun propellant.Currently,the preparation of GDGP employed a tedious two-step method involving organic solvents,which hinders the large-scale preparation of GDGP.In this paper,GDGP was successfully prepared via a novelty and environmentally friendly one-step method.The obtained samples were characterized by FT-IR,Raman,SEM and XPS.The results showed that the content of nitrate groups gradiently increased from the surface to the core in the surface layer of GDGP and the surface layer of GDGP exhibited a higher compaction than that of raw gun propellant,with a well-preserved nitrocellulose structure.The denitration process enabled the propellant surface with regressive energy density and good progressive burning performance,as confirmed by oxygen bomb and closed bomb test.At the same time,the effects of different solvents on the component loss of propellant were compared.The result showed that water caused the least component loss.Finally,the stability of GDGP was confirmed by methyl-violet test.This work not only provided environmentally friendly,simple and economic preparation of GDGP,but also confirmed the stability of GDGP prepared by this method.

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.

约束Green函数与非线性地基梁模态分析

在生物和医学领域,微机电系统(MEMS)中的微梁结构在植入人体的使用时,由于体内的细胞环境类似于水凝胶,在这种环境下工作,设备和仪器的精度和稳定性很大程度上受到细胞弹性的影响.为了分析此类地基梁的动力学问题,本文建立了非线性基础上的梁振动模型,研究了任意位置弹簧和非线性弹簧基础上的梁模态.通过Laplace变换和线性叠加原理,得到了一种约束Green函数,利用数值计算验证方案的有效性,并研究了各种重要物理参数的影响,发现弹簧位置向跨中移动时,模态对称性被打破,弹簧刚度增加,模态阶数改变.

GREEN'S FUNCTION AND THE POINTWISE BEHAVIORS OF THE VLASOV-POISSON-FOKKER-PLANCK SYSTEM

The pointwise space-time behaviors of the Green’s function and the global solution to the Vlasov-Poisson-Fokker-Planck(VPFP)system in three dimensional space are studied in this paper.It is shown that the Green’s function consists of the diffusion waves decaying exponentially in time but algebraically in space,and the singular kinetic waves which become smooth for all(t,x,v)when t>0.Furthermore,we establish the pointwise space-time behaviors of the global solution to the nonlinear VPFP system when the initial data is not necessarily smooth in terms of the Green’s function.

An improved empirical green function method to extract microseismic signals

In mining, reservoir impoundment and hydraulic fracturing exploitation, the induced microseismic clusters may present complex waveforms at receivers since the individual events generally arrive very close. And the routine methods for arrival picking are insufficiently efficient due to these strong influences. Here, we modified the empirical green function method and applied it to extract the green functions and the radiation coefficients of microseismic events with the largest energy in clusters.Multiple-channel records were used to estimate an "average" source spectrum and then the "average" source wavelet was removed from the records by deconvolution. We applied this method to the real data,and the result indicated clear improvement in extracting the dominant event of the clusters.

Characteristic analysis of scattering field in two-layer media by Green's function

The problem of three-dimensional(3D) acoustic scattering in a complex medium has aroused considerable interest of researchers for many years. An ultrasonic scattered field calculating technique is proposed to study the scattering echo from strongly scattered materials in a two-layer medium in this work. Firstly, with the high frequency stationary phase method,the Green's function of two-layer fluid media is derived. And then based on the idea of integral equation discretization,the Green's function method is extended to two-layer fluid media to derive the scattering field expression of defects in a complex medium. With this method, the scattering field of 3D defect in a two-layer medium is calculated and the characteristics of received echoes are studied. The results show that this method is able to solve the scattering P wave field of 3D defect with arbitrary shape at any scattering intensity in two-layer media. Considering the circumstance of waterimmersion ultrasonic non-destructive test(NDT), the scattering sound field characteristics of different types of defects are analyzed by simulation, which will help to optimize the detection scheme and corresponding imaging method in practice so as to improve the detection quality.

Green's functions of two-dimensional piezoelectric quasicrystal in half-space and bimaterials

In this paper,we obtain Green’s functions of two-dimensional(2D)piezoelectric quasicrystal(PQC)in half-space and bimaterials.Based on the elastic theory of QCs,the Stroh formalism is used to derive the general solutions of displacements and stresses.Then,we obtain the analytical solutions of half-space and bimaterial Green’s functions.Besides,the interfacial Green’s function for bimaterials is also obtained in the analytical form.Before numerical studies,a comparative study is carried out to validate the present solutions.Typical numerical examples are performed to investigate the effects of multi-physics loadings such as the line force,the line dislocation,the line charge,and the phason line force.As a result,the coupling effect among the phonon field,the phason field,and the electric field is prominent,and the butterfly-shaped contours are characteristic in 2D PQCs.In addition,the changes of material parameters cause variations in physical quantities to a certain degree.

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.

Finite Volume Time Domain with the Green Function Method for Electromagnetic Scattering in Schwarzschild Spacetime

The finite volume time domain(FVTD) algorithm and Green function algorithm are extended to Schwarzschild spacetime for numerical simulation of electromagnetic scattering. The FVTD method in Schwarzschild spacetime is developed by filling the flat spacetime with an equivalent medium. The Green function in Schwarzschild spacetime is acquired by solving initial value problems. Both the FVTD code and the Green function code are validated by numerical results. Scattering in Schwarzschild spacetime is simulated with these methods.

An Improved Preisach Distribution Function Identification Method Considering the Reversible Magnetization

This paper presents an identification method of the scalar Preisach model to consider the effect of reversible magnetization in the process of distribution function identification.By reconsidering the identification process by stripping the influence of reversible components from the measurement data,the Preisach distribution function is identified by the pure irreversible components.In this way,the simulation accuracy of both limiting hysteresis loops and the inner internal symmetrical small hysteresis loop is ensured.Furthermore,through a discrete Preisach plane with a hybrid discretization method,the irreversible magnetic flux density components are computed more efficiently through the improved Preisach model.Finally,the proposed method results are compared with the traditional method and the traditional method considering reversible magnetization and validated by the laboratory test for the B30P105 electrical steel by Epstein frame.

Generation of endurance time excitation functions using spectral representation method

In this study,application of the spectral representation method for generation of endurance time excitation functions is introduced.Using this method,the intensifying acceleration time series is generated so that its acceleration response spectrum in any desired time duration is compatible with a time-scaled predefined acceleration response spectrum.For this purpose,simulated stationary acceleration time series is multiplied by the time dependent linear modulation function,then using a simple iterative scheme,it is forced to match a target acceleration response spectrum.It is shown that the generated samples have excellent conformity in low frequency,which is useful for nonlinear endurance time analysis.In the second part of this study,it is shown that this procedure can be extended to generate a set of spatially correlated endurance time excitation functions.This makes it possible to assess the performance of long structures under multi-support seismic excitation using endurance time analysis.

Comprehensive Evaluation of New Maize Varieties for Grain and Fodder Based on Membership Function Method

To screen new maize(Zea mays L.)varieties suitable for food and fodder dual-purpose in Du'an Yao Autonomous County of Guangxi,the agronomic characters,yield and quality indexes of 12 new maize varieties were measured,and the correlation between various indexes were analyzed,and the comprehensive performance of tested varieties was evaluated by membership function method.The results showed that Guidan 671 had the highest grain yield and whole-plant biomass at 10908 and 49965 kg/hm^(2),respectively,and the second was Zhaoyu 215 with a grain yield and whole-plant biomass of 10086 and 47175 kg/hm^(2),respectively.Grain yield was highly significantly positively correlated with ear diameter and 100-grain weight(P<0.01),and significantly correlated with whole-plant biomass,starch content,ear length and grain number per row(P<0.05);and the whole-plant biomass was highly significantly correlated with the number of grains per row(P<0.01),and significantly correlated with starch content,panicle length,plant height and panicle height(P<0.05).The comprehensive performance scores of the tested varieties from high to low were Guidan 671,Zhaoyu 215,Guidan 669,Guidan 6208,Guidan 666,Guidan 6205,Guidan 660,Guidan 6203,Guidan 6206,Guidan 162,Guidan 668 and Guidan 673.According to the values of membership function and combined with various indexes,Guidan 671 and Zhaoyu 215 had good comprehensive performance,and could be used as the first choice for food and fodder dual-purpose maize varieties in Du'an Yao Autonomous County.