Highly Accurate Golden Section Search Algorithms and Fictitious Time Integration Method for Solving Nonlinear Eigenvalue Problems

This study sets up two new merit functions,which are minimized for the detection of real eigenvalue and complex eigenvalue to address nonlinear eigenvalue problems.For each eigen-parameter the vector variable is solved from a nonhomogeneous linear system obtained by reducing the number of eigen-equation one less,where one of the nonzero components of the eigenvector is normalized to the unit and moves the column containing that component to the right-hand side as a nonzero input vector.1D and 2D golden section search algorithms are employed to minimize the merit functions to locate real and complex eigenvalues.Simultaneously,the real and complex eigenvectors can be computed very accurately.A simpler approach to the nonlinear eigenvalue problems is proposed,which implements a normalization condition for the uniqueness of the eigenvector into the eigenequation directly.The real eigenvalues can be computed by the fictitious time integration method(FTIM),which saves computational costs compared to the one-dimensional golden section search algorithm(1D GSSA).The simpler method is also combined with the Newton iterationmethod,which is convergent very fast.All the proposed methods are easily programmed to compute the eigenvalue and eigenvector with high accuracy and efficiency.

Efficient method to calculate the eigenvalues of the Zakharov–Shabat system

A numerical method is proposed to calculate the eigenvalues of the Zakharov–Shabat system based on Chebyshev polynomials. A mapping in the form of tanh(ax) is constructed according to the asymptotic of the potential function for the Zakharov–Shabat eigenvalue problem. The mapping can distribute Chebyshev nodes very well considering the gradient for the potential function. Using Chebyshev polynomials, tanh(ax) mapping, and Chebyshev nodes, the Zakharov–Shabat eigenvalue problem is transformed into a matrix eigenvalue problem. This method has good convergence for the Satsuma–Yajima potential and the convergence rate is faster than the Fourier collocation method. This method is not only suitable for simple potential functions but also converges quickly for a complex Y-shape potential. It can also be further extended to other linear eigenvalue problems.

THE INTERIOR TRANSMISSION EIGENVALUE PROBLEM FOR AN ANISOTROPIC MEDIUM BY A PARTIALLY COATED BOUNDARY

We consider the interior transmission eigenvalue problem corresponding to the scattering for an anisotropic medium of the scalar Helmholtz equation in the case where the boundary?Ωis split into two disjoint parts and possesses different transmission conditions.Using the variational method,we obtain the well posedness of the interior transmission problem,which plays an important role in the proof of the discreteness of eigenvalues.Then we achieve the existence of an infinite discrete set of transmission eigenvalues provided that n≡1,where a fourth order differential operator is applied.In the case of n■1,we show the discreteness of the transmission eigenvalues under restrictive assumptions by the analytic Fredholm theory and the T-coercive method.

On the eigenvalues and eigendisplacement of the critical mode in horizontally layered media

Wave propagation in horizontally layered media is a classical problem in seismic-wave theory.In semi-infinite space,a nondispersive Rayleigh wave mode exists,and the eigendisplacement decays exponentially with depth.In a layered model with increasing layer velocity,the phase velocity of the Rayleigh wave varies between the S-wave velocity of the bottom half-space and that of the classical Rayleigh wave propagated in a supposed half-space formed by the parameters of the top layer.If the phase velocity is the same as the P-or S-wave velocity of the layer,which is called the critical mode or critical phase velocity of surface waves,the general solution of the wave equation is not a homogeneous(expressed by trigonometric functions)or inhomogeneous(expressed by exponential functions)plane wave,but one whose amplitude changes linearly with depth(expressed by a linear function).Theories based on a general solution containing only trigonometric or exponential functions do not apply to the critical mode,owing to the singularity at the critical phase velocity.In this study,based on the classical framework of generalized reflection and transmission coefficients,the propagation of surface waves in horizontally layered media was studied by introducing a solution for the linear function at the critical phase velocity.Therefore,the eigenvalues and eigenfunctions of the critical mode can be calculated by solving a singular problem.The eigendisplacement characteristics associated with the critical phase velocity were investigated for different layered models.In contrast to the normal mode,the eigendisplacement associated with the critical phase velocity exhibits different characteristics.If the phase velocity is equal to the S-wave velocity in the bottom half-space,the eigendisplacement remains constant with increasing depth.

Real eigenvalues determined by recursion of eigenstates

Quantum physics is primarily concerned with real eigenvalues,stemming from the unitarity of time evolutions.With the introduction of PT symmetry,a widely accepted consensus is that,even if the Hamiltonian of the system is not Hermitian,the eigenvalues can still be purely real under specific symmetry.Hence,great enthusiasm has been devoted to exploring the eigenvalue problem of non-Hermitian systems.In this work,from a distinct perspective,we demonstrate that real eigenvalues can also emerge under the appropriate recursive condition of eigenstates.Consequently,our findings provide another path to extract the real energy spectrum of non-Hermitian systems,which guarantees the conservation of probability and stimulates future experimental observations.

A Simple Embedding Method for the Laplace-Beltrami Eigenvalue Problem onImplicit Surfaces

We propose a simple embedding method for computing the eigenvalues and eigenfunctions of the Laplace-Beltrami operator on implicit surfaces.The approach follows an embedding approach for solving the surface eikonal equation.We replace the differential operator on the interface with a typical Cartesian differential operator in the surface neighborhood.Our proposed algorithm is easy to implement and efficient.We will give some two-and three-dimensional numerical examples to demonstrate the effectiveness of our proposed approach.

Reciprocal Distance Laplacian Eigenvalue Distribution Based on Graph Parameters

Let G be a connected graph of order n and m_(RD)^(L)_(G)I denote the number of reciprocal distance Laplacian eigenvaluesof G in an interval I.For a given interval I,we mainly present several bounds on m_(RD)^(L)_(G)I in terms of various structuralparameters of the graph G,including vertex-connectivity,independence number and pendant vertices.

Petrochemical eigenvalues and diagrams for the identification of metamorphic rocks'protolith,taking the host rocks of Dashuigou tellurium deposit in China as an example

The Dashuigou tellurium deposit is the world’s only known independent tellurium deposit.By restoring metamorphic rocks’protolith,we seek to understand not only the development and evolution trajectory of the region but also the origin of the relevant deposits.While there are many ways to restore metamorphic rocks’protolith,we take the host metamorphic rocks of Dashuigou tellurium deposit and leverage various petrochemical eigenvalues and related diagrams previously proposed to reveal the deposit’s host metamorphic rocks’protolith.The petrochemical eigenvalues include molecular number,Niggli’s value,REE parity ratio,CaO/Al_(2)O_(3)ratio,Fe^(3+) /(Fe^(3+) -+Fe^(2+) )ratio,chondrite-normalized REE value,logarithmic REE value,various REE eigenvalues including scandium,Eu/Sm ratio,total REE amount,light and heavy REEs,δEu,Eu anomaly,Sm/Nd ratio,and silicon isotope δ^(30) SiNBS-29‰,etc.The petrochemical plots include ACMs,100 mg-c-(al+alk),SiO_(2)-(Na_(2)O+K_(2)O),(al+fm)-(c+alk)versus Si,FeO+Fe_(2)O^(3+) TiO)-Al_(2)O_(3)-MgO,c-mg,Al_(2)O_(3)-(Na_(2)O+K_(2)O),chondrite-normalized REE model,La/Yb-REE,and Sm/Nd ratio,etc.On the basis of these comprehensive analyses,the following conclusions are drawn,starting from the many mantle-derived types of basalt developed in the study area of different geological ages,combined with the previously published research results on the deposit s fluid inclusions and sulfur and lead isotopes.The deposit is formed by mantle degassing in the form of a mantle plume in the late Yanshanian orogeny.The degassed fluids are rich in nano-sc ale substances including Fe,Te,S,As,Bi,Au,Se,H_(2),CO_(2),N_(2),H_(2)O,and CH_(4),which are enriched by nano-effect,and then rise to a certain part of the crust in the form of mantle plume along the lithospheric fault to form the deposit.The ultimate power for tellurium mineralization was from H_(2)flow with high energy,which was produced through radiation from the melted iron of the Earth’s outer core.The H,flow results in the Earth’s degassing,as well as the mantle and crust’s uplift.

Physics-constrained neural network for solving discontinuous interface K-eigenvalue problem with application to reactor physics

Machine learning-based modeling of reactor physics problems has attracted increasing interest in recent years.Despite some progress in one-dimensional problems,there is still a paucity of benchmark studies that are easy to solve using traditional numerical methods albeit still challenging using neural networks for a wide range of practical problems.We present two networks,namely the Generalized Inverse Power Method Neural Network(GIPMNN)and Physics-Constrained GIPMNN(PC-GIPIMNN)to solve K-eigenvalue problems in neutron diffusion theory.GIPMNN follows the main idea of the inverse power method and determines the lowest eigenvalue using an iterative method.The PC-GIPMNN additionally enforces conservative interface conditions for the neutron flux.Meanwhile,Deep Ritz Method(DRM)directly solves the smallest eigenvalue by minimizing the eigenvalue in Rayleigh quotient form.A comprehensive study was conducted using GIPMNN,PC-GIPMNN,and DRM to solve problems of complex spatial geometry with variant material domains from the fleld of nuclear reactor physics.The methods were compared with the standard flnite element method.The applicability and accuracy of the methods are reported and indicate that PC-GIPMNN outperforms GIPMNN and DRM.

LARGE DEVIATIONS FOR TOP EIGENVALUES OFβ-JACOBI ENSEMBLES AT SCALING TEMPERATURES

Letλ=(λ_(1),…,λ_(n))beβ-Jacobi ensembles with parameters p_(1),p_(2),n andβ,withβvarying with n.Set■.Suppose that■and 0≤σγ<1.We offer the large deviation for p_(1)+p_(2)/p_(1)■λ_(i)whenγ>0 via the large deviation of the corresponding empirical measure and via a direct estimate,respectively,whenγ=0.

Estimates on the eigenvalues of complex nonlocal Sturm-Liouville problems

The present paper deals with the eigenvalues of complex nonlocal Sturm-Liouville boundary value problems.The bounds of the real and imaginary parts of eigenvalues for the nonlocal Sturm-Liouville differential equation involving complex nonlocal potential terms associated with nonlocal boundary conditions are obtained in terms of the integrable conditions of coefficients and the real part of the eigenvalues.

Minimum eigenvalue based adaptive fault compensation for hypersonic vehicles

The attitude tracking control problem is addressed for hypersonic vehicles under actuator faults that may cause an uncertain time-varying control gain matrix.An adaptive compensation scheme is developed to ensure system stability and asymptotic tracking properties,including a kinematic control signal and a dynamic control signal.To deal with the uncertainties of the control gain matrix,a new positive definite one is constructed.The minimum eigenvalue of such a new control gain matrix is estimated.Simulation results of application to an X-33 vehicle model verify the effectiveness of the proposed minimum eigenvalue based adaptive fault compensation scheme.

基于一般散射模型的Hybrid Freeman/Eigenvalue分解算法(英文)

提出了一种新的基于一般散射模型的hybrid Freeman/eigenvalue分解算法,用于分析极化合成孔径雷达(PolS AR)数据。文中,单位矩阵作为体散射模型,相干矩阵的两个较大特征值对应的特征向量作为表面散射模型和二次散射模型,并且不需要反射对称条件。新算法有三个优点:第一,表面散射和二次散射不需要反射对称条件,更符合一般散射体的建模;第二,因为散射能量是相干矩阵特征值的线性组合,所以散射能量具有旋转不变性;第三,表面散射能量和二次散射能量避免了负值现象。在San Francisco地区的AIRSAR数据上进行了实验,证明了新算法的有效性。

一种自适应的混合Freeman/Eigenvalue极化分解模型

全极化SAR数据的极化分解在土地利用分类、目标检测与识别以及地表参数反演等领域得到了广泛应用。目前,主要有基于特征值分解和基于模型分解2类极化分解方法。混合Freeman/Eigenvalue极化分解结合了两者的优势,避免了基于模型的极化分解中负功率问题并且能够利用已知的散射机制解释分解后的散射分量。为了进一步拓展该分解在不同地表类型中的应用,通过引入参数Neumann一般化体散射模型,提出了一种自适应的极化分解模型。利用德国Black Forest地区的L波段AirSAR(airborne synthetic aperture Radar)全极化数据进行实验,并与现有的Yamaguchi三分量模型和自适应非负分解(adaptive nonnegative eigenvalue decomposition,ANNED)对比分析,以验证模型的有效性。研究表明,自适应的混合Freeman/Eigenvalue极化分解模型保证了分解能量的非负性及完全分解,适应于不同类型的地表,能有效地区分不同地类。

Derivatives of repeated eigenvalues and corresponding eigenvectors of damped systems

A procedure is presented for computing the derivatives of repeated eigenvalues and the corresponding eigenvectors of damped systems. The derivatives are calculated in terms of the eigenvalues and eigenvectors of the second-order system, and the use of rather undesirable state space representation is avoided. Hence the cost of computation is greatly reduced. The efficiency of the proposed procedure is illustrated by considering a 5-DOF non-proportionally damped system.

A NEW METHOD FOR ESTIMATING BOUNDS OF EIGENVALUES

A new method for estimating the bounds of eigenvalues ispresented. In order to show that the method proposed is as effectiveas Qiu's an undamping spring-mass system with 5 nodes and 5 degrees ofreedom is given. To illustrate that the present method can beapplied to structures which cannot be treated by non-negativedecomposition, a plane frame with 202 nodes and 357 beam elements isgiven. The results show that the present method is effective forestimating the bounds of eigenvalues and is more common than Qiu's.

Joint eigenvalue estimation by balanced simultaneous Schur decomposition

The problem of joint eigenvalue estimation for the non-defective commuting set of matrices A is addressed. A procedure revealing the joint eigenstructure by simultaneous diagonalization of. A with simultaneous Schur decomposition （SSD） and balance procedure alternately is proposed for performance considerations and also for overcoming the convergence difficulties of previous methods based only on simultaneous Schur form and unitary transformations, it is shown that the SSD procedure can be well incorporated with the balancing algorithm in a pingpong manner, i. e., each optimizes a cost function and at the same time serves as an acceleration procedure for the other. Under mild assumptions, the convergence of the two cost functions alternately optimized, i. e., the norm of A and the norm of the left-lower part of A is proved. Numerical experiments are conducted in a multi-dimensional harmonic retrieval application and suggest that the presented method converges considerably faster than the methods based on only unitary transformation for matrices which are not near to normality.

Some Trace Formulas of a Eigenvalue Problem

In present paper, using some methods of approximation theory, the trace formulas for eigenvalues of a eigenvalue problem are calculated under the periodic condition and the decaying condition at x∞.

There Is a Way to Comprise Half-Integer Eigenvalues for Photon Spin

In this article, an attempt based on Spin Topological Space, STS, to give areasonable detailed account of the cause of photonic fermionization phenomena of light photon is made. STS is an unconventional spin space in quantum mechanics, which can be used to account for where the unconventional half-integer spin eigenvalues phenomenon of light photon comes from. We suggest to dectect the possible existence of photonic one-third-spinization phenomenon of light photon, by using three beams of light photon in interference experiment.

One Computational Method of the Eigenvalues of the Horizontal Vibration Problem of Beam

This paper considers one computational method of the eigenvalues approximate value of the horizontal vibration problem of beam. The proof of our main result is based on the variational formula. First of all, Cauchy inequality is used to obtain a basic inequality. Secondly, the functions of basis are made by Galerkin method, and the error estimates of eignevalues are obtained by Cauchy inequality. At last, the computational method of the approximate value of the eigenvalues turns out immediately, and acc...