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.

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.

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.

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.

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.

可对角化矩阵特征值分解扰动问题的快速求解方法

针对特征值扰动计算的传统方法收敛速度慢的问题,提出了一种求解特征值扰动问题的快速迭代算法.首先,通过矩阵变换将初始矩阵的特征值扰动问题转化为对角矩阵的特征值扰动问题.然后,提出了一种快速迭代算法求解扰动参数,同时对算法的收敛性进行分析,并将其与基于摄动级数展开法导出的方法进行对比.再次,采用逐一求解特征值并进行矩阵降阶的策略,有效降低运算量.最后,通过2个算例分别展示算法的计算过程及其在结构模态参数追踪方面的应用效果.

基于ANSYS Workbench稳定性分析的大型真空容器优化设计

真空容器是石油、化工等行业中常见的设备之一,与内压容器不同,真空容器往往发生的是稳定性破坏。文中首先采用有限元软件ANSYS Workbench对DN22000 mm的真空容器进行了特征值屈曲分析,结果表明:其外压承载系数为17.427,远大于外压安全系数2.5。进而采用Response Surface Optimization模块对该真空容器进行优化设计,分析了外压承载系数、容器质量随穹顶壁厚和筒体壁厚的变化规律,以及外压承载系数、容器质量、应力随穹顶辐射筋数量的变化规律。最终确定穹顶壁厚取20 mm、筒体壁厚取22 mm、穹顶辐射筋取16个,优化后结构较最初结构质量减小了83140 kg,减小了17.02%,并且稳定性满足ASME规范要求,实现了容器轻量化的设计目标,从而显著降低制造成本。

Givens矩阵的性质及其在迭代法中的应用

研究了2阶Givens矩阵的一些性质,该矩阵的特征值为复数,谱半径恰好为1;在此基础上通过多个例子展示了Givens矩阵在迭代法中的应用,并从理论上给出了分析.

基于导数张量的多元函数极值点判定

针对多元函数极值点判定问题,定义了多元函数的导数张量,给出了张量形式的多元函数带皮亚诺余项的泰勒公式,提出了基于张量的最大和最小Z特征值对多元函数极值点进行判定的方法.

基于特征值的动态数字信道化子带检测算法

动态数字信道化接收结构通过检测所有子带,判断信号是否存在,为综合滤波器组处理提供依据,因此子带检测在接收结构中起着关键作用。针对传统检测算法存在低信噪比(signal-to-noise ratio,SNR)下检测性能不高的问题,依据随机矩阵理论,提出了基于最大最小特征值之差与平均特征值之比的检测算法,利用平均特征值和最小特征值的极限分布规律来推导算法的检测门限。其次,根据所有子带数据获取的特征值信息对所提算法进行了优化。最后,在动态数字信道化接收结构中,分析不同因素下算法的性能,表明了所提算法能够克服低SNR的影响,子带检测的性能更好。

计算Householder变换之行列式的六种解法

本文探讨Householder变换之行列式的几种解法.目的是提高学生综合运用线性代数知识的能力,培养学生灵活掌握并融会贯通整个线性代数知识点的能力.

基于ANSYS的波纹管波形参数对柱失稳影响分析

为了研究波纹管波形参数对波纹管柱失稳的影响,使用ANSYS Workbench对波纹管进行了模态分析与特征值屈曲分析。采用控制变量的方法,改变波纹管波形参数,分别进行有限元分析;并对相同波数下的复式波纹管进行有限元分析,且与单式波纹管稳定性做分析比较。有限元结果表明:在实际工作中,可通过增加波纹管壁厚、波距,或减小波高、减少波数的方法,减少波纹管柱失稳的发生;当波纹管长度较长、波数较多时,使用复式波纹管是更优的选择。

Kinematic analysis of flexible bipedal robotic systems

In spite of its intrinsic complexities,the passive gait of bipedal robots on a sloping ramp is a subject of interest for numerous researchers.What distinguishes the present research from similar works is the consideration of flexibility in the constituent links of this type of robotic systems.This is not a far-fetched assumption because in the transient(impact)phase,due to the impulsive forces which are applied to the system,the likelihood of exciting the vibration modes increases considerably.Moreover,the human leg bones that are involved in walking are supported by viscoelastic muscles and ligaments.Therefore,for achieving more exact results,it is essential to model the robot links with viscoelastic properties.To this end,the Gibbs-Appell formulation and Newton's kinematic impact law are used to derive the most general form of the system's dynamic equations in the swing and transient phases of motion.The most important issue in the passive walking motion of bipedal robots is the determination of the initial robot configuration with which the system could accomplish a periodic and stable gait solely under the effect of gravitational force.The extremely unstable nature of the system studied in this paper and the vibrations caused by the impulsive forces induced by the impact of robot feet with the inclined surface are some of the very serious challenges encountered for achieving the above-mentioned goal.To overcome such challenges,an innovative method that uses a combination of the linearized equations of motion in the swing phase and the algebraic motion equations in the transition phase is presented in this paper to obtain an eigenvalue problem.By solving this problem,the suitable initial conditions that are necessary for the passive gait of this bipedal robot on a sloping surface are determined.The effects of the characteristic parameters of elastic links including the modulus of elasticity and the Kelvin-Voigt coefficient on the walking stability of this type of robotic systems are also studied.The findings of this parametric study reveal that the increase in the Kelvin-Voigt coefficient enhances the stability of the robotic system,while the increase in the modulus of elasticity has an opposite effect.

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.

Review of Computational Approaches to Optimization Problems in Inhomogeneous Rods and Plates

In this paper,we review computational approaches to optimization problems of inhomogeneous rods and plates.We consider both the optimization of eigenvalues and the localization of eigenfunctions.These problems are motivated by physical problems including the determination of the extremum of the fundamental vibration frequency and the localization of the vibration displacement.We demonstrate how an iterative rearrangement approach and a gradient descent approach with projection can successfully solve these optimization problems under different boundary conditions with different densities given.

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.

Trace formula of the integro-differential operator on a quantum graph

In this paper we study the eigenvalue problem for integro-differential operators on a lasso graph.The trace formula of the operator is established by applying the residual technique in complex analysis.

Stability Analysis of Inverse Lax-Wendroff Procedure for a High order Compact Finite Difference Schemes

This paper considers the finite difference(FD)approximations of diffusion operators and the boundary treatments for different boundary conditions.The proposed schemes have the compact form and could achieve arbitrary even order of accuracy.The main idea is to make use of the lower order compact schemes recursively,so as to obtain the high order compact schemes formally.Moreover,the schemes can be implemented efficiently by solving a series of tridiagonal systems recursively or the fast Fourier transform(FFT).With mathematical induction,the eigenvalues of the proposed differencing operators are shown to be bounded away from zero,which indicates the positive definiteness of the operators.To obtain numerical boundary conditions for the high order schemes,the simplified inverse Lax-Wendroff(SILW)procedure is adopted and the stability analysis is performed by the Godunov-Ryabenkii method and the eigenvalue spectrum visualization method.Various numerical experiments are provided to demonstrate the effectiveness and robustness of our algorithms.

用于抵抗扭转采取边缘弹性约束下的单轴压缩下板材极限强度:一种新的综合方法

The ultimate strength of platings under compression is one of the most important factors to be addressed in the ship design.Current Rules for ship structural design generally provide explicit strength check criteria against buckling for simply supported and clamped platings.Nevertheless,ship platings generally exhibit an intermediate behaviour between the simple support and the clamped conditions,which implies that the torsional stiffness of supporting members should be duly considered.Hence,the main aim of this study is the development of new design formulas for the ultimate strength of platings under uniaxial compression,with short and/or long edges elastically restrained against torsion.In this respect,two benchmark studies are performed.The former is devoted to the development of new equations for the elastic buckling coefficients of platings with edges elastically restrained against torsion,based on the results of the eigenvalue buckling analysis,performed by Ansys Mechanical APDL.The latter investigates the ultimate strength of platings with elastically restrained edges,by systematically varying the plate slenderness ratio and the torsional stiffness of supporting members.Finally,the effectiveness of the new formulation is checked against a wide number of finite element(FE)simulations,to cover the entire design space of ship platings.