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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

Beamspace maximum likelihood algorithm based on sum and difference beams for elevation estimation

Beamspace super-resolution methods for elevation estimation in multipath environment has attracted significant attention, especially the beamspace maximum likelihood(BML)algorithm. However, the difference beam is rarely used in superresolution methods, especially in low elevation estimation. The target airspace information in the difference beam is different from the target airspace information in the sum beam. And the use of difference beams does not significantly increase the complexity of the system and algorithms. Thus, this paper applies the difference beam to the beamformer to improve the elevation estimation performance of BML algorithm. And the direction and number of beams can be adjusted according to the actual needs. The theoretical target elevation angle root means square error(RMSE) and the computational complexity of the proposed algorithms are analyzed. Finally, computer simulations and real data processing results demonstrate the effectiveness of the proposed algorithms.

Research on a deconvolution algorithm for laser-induced fluorescence diagnosis based on the maximum entropy principle

Laser-induced fluorescence(LIF)spectroscopy is employed for plasma diagnosis,necessitating the utilization of deconvolution algorithms to isolate the Doppler effect from the raw spectral signal.However,direct deconvolution becomes invalid in the presence of noise as it leads to infinite amplification of high-frequency noise components.To address this issue,we propose a deconvolution algorithm based on the maximum entropy principle.We validate the effectiveness of the proposed algorithm by utilizing simulated LIF spectra at various noise levels(signal-to-noise ratio,SNR=20–80 d B)and measured LIF spectra with Xe as the working fluid.In the typical measured spectrum(SNR=26.23 d B)experiment,compared with the Gaussian filter and the Richardson–Lucy(R-L)algorithm,the proposed algorithm demonstrates an increase in SNR of 1.39 d B and 4.66 d B,respectively,along with a reduction in the root-meansquare error(RMSE)of 35%and 64%,respectively.Additionally,there is a decrease in the spectral angle(SA)of 0.05 and 0.11,respectively.In the high-quality spectrum(SNR=43.96 d B)experiment,the results show that the running time of the proposed algorithm is reduced by about98%compared with the R-L iterative algorithm.Moreover,the maximum entropy algorithm avoids parameter optimization settings and is more suitable for automatic implementation.In conclusion,the proposed algorithm can accurately resolve Doppler spectrum details while effectively suppressing noise,thus highlighting its advantage in LIF spectral deconvolution applications.

Strengthening-softening transition and maximum strength in Schwarz nanocrystals

Recently,a Schwarz crystal structure with curved grain boundaries(GBs)constrained by twin-boundary(TB)networks was discovered in nanocrystalline Cu through experiments and atomistic simulations.Nanocrystalline Cu with nanosized Schwarz crystals exhibited high strength and excellent thermal stability.However,the grainsize effect and associated deformation mechanisms of Schwarz nanocrystals remain unknown.Here,we performed large-scale atomistic simulations to investigate the deformation behaviors and grain-size effect of nanocrystalline Cu with Schwarz crystals.Our simulations showed that similar to regular nanocrystals,Schwarz nanocrystals exhibit a strengthening-softening transition with decreasing grain size.The critical grain size in Schwarz nanocrystals is smaller than that in regular nanocrystals,leading to a maximum strength higher than that of regular nanocrystals.Our simulations revealed that the softening in Schwarz nanocrystals mainly originates from TB migration(or detwinning)and annihilation of GBs,rather than GB-mediated processes(including GB migration,sliding and diffusion)dominating the softening in regular nanocrystals.Quantitative analyses of simulation data further showed that compared with those in regular nanocrystals,the GB-mediated processes in Schwarz nanocrystals are suppressed,which is related to the low volume fraction of amorphous-like GBs and constraints of TB networks.The smaller critical grain size arises from the suppression of GB-mediated processes.

Stochastic Maximum Principle for Optimal Advertising Models with Delay and Non-Convex Control Spaces

In this paper we study optimal advertising problems that model the introduction of a new product into the market in the presence of carryover effects of the advertisement and with memory effects in the level of goodwill. In particular, we let the dynamics of the product goodwill to depend on the past, and also on past advertising efforts. We treat the problem by means of the stochastic Pontryagin maximum principle, that here is considered for a class of problems where in the state equation either the state or the control depend on the past. Moreover the control acts on the martingale term and the space of controls U can be chosen to be non-convex but now the space of controls U can be chosen to be non-convex. The maximum principle is thus formulated using a first-order adjoint Backward Stochastic Differential Equations (BSDEs), which can be explicitly computed due to the specific characteristics of the model, and a second-order adjoint relation.

Efficient slope reliability analysis under soil spatial variability using maximum entropy distribution with fractional moments

Spatial variability of soil properties imposes a challenge for practical analysis and design in geotechnical engineering.The latter is particularly true for slope stability assessment,where the effects of uncertainty are synthesized in the so-called probability of failure.This probability quantifies the reliability of a slope and its numerical calculation is usually quite involved from a numerical viewpoint.In view of this issue,this paper proposes an approach for failure probability assessment based on Latinized partially stratified sampling and maximum entropy distribution with fractional moments.The spatial variability of geotechnical properties is represented by means of random fields and the Karhunen-Loève expansion.Then,failure probabilities are estimated employing maximum entropy distribution with fractional moments.The application of the proposed approach is examined with two examples:a case study of an undrained slope and a case study of a slope with cross-correlated random fields of strength parameters under a drained slope.The results show that the proposed approach has excellent accuracy and high efficiency,and it can be applied straightforwardly to similar geotechnical engineering problems.

Maximum Correntropy Criterion-Based UKF for Loosely Coupling INS and UWB in Indoor Localization

Indoor positioning is a key technology in today’s intelligent environments,and it plays a crucial role in many application areas.This paper proposed an unscented Kalman filter(UKF)based on the maximum correntropy criterion(MCC)instead of the minimummean square error criterion(MMSE).This innovative approach is applied to the loose coupling of the Inertial Navigation System(INS)and Ultra-Wideband(UWB).By introducing the maximum correntropy criterion,the MCCUKF algorithm dynamically adjusts the covariance matrices of the system noise and the measurement noise,thus enhancing its adaptability to diverse environmental localization requirements.Particularly in the presence of non-Gaussian noise,especially heavy-tailed noise,the MCCUKF exhibits superior accuracy and robustness compared to the traditional UKF.The method initially generates an estimate of the predicted state and covariance matrix through the unscented transform(UT)and then recharacterizes the measurement information using a nonlinear regression method at the cost of theMCC.Subsequently,the state and covariance matrices of the filter are updated by employing the unscented transformation on the measurement equations.Moreover,to mitigate the influence of non-line-of-sight(NLOS)errors positioning accuracy,this paper proposes a k-medoid clustering algorithm based on bisection k-means(Bikmeans).This algorithm preprocesses the UWB distance measurements to yield a more precise position estimation.Simulation results demonstrate that MCCUKF is robust to the uncertainty of UWB and realizes stable integration of INS and UWB systems.

Effective approach for conformal subarray design based on maximum entropy of planar mappings

In this paper,an effective algorithm for optimizing the subarray of conformal arrays is proposed.The method first divides theconformal array into several first-level subarrays.It uses the X algorithm to solve the feasible solution of first-level subarray tiling and employs the particle swarm algorithm to optimize the conformal array subarray tiling scheme with the maximum entropy of the planar mapping as the fitness function.Subsequently,convex optimization is applied to optimize the subarray amplitude phase.Data results verify that the method can effectively find the optimal conformal array tiling scheme.

Correlation between the rock mass properties and maximum horizontal stress:A case study of overcoring stress measurements

Understanding the mechanical properties of the lithologies is crucial to accurately determine the horizontal stress magnitude.To investigate the correlation between the rock mass properties and maximum horizontal stress,the three-dimensional(3D)stress tensors at 89 measuring points determined using an improved overcoring technique in nine mines in China were adopted,a newly defined characteristic parameter C_(ERP)was proposed as an indicator for evaluating the structural properties of rock masses,and a fuzzy relation matrix was established using the information distribution method.The results indicate that both the vertical stress and horizontal stress exhibit a good linear growth relationship with depth.There is no remarkable correlation between the elastic modulus,Poisson's ratio and depth,and the distribution of data points is scattered and messy.Moreover,there is no obvious relationship between the rock quality designation(RQD)and depth.The maximum horizontal stress σ_(H) is a function of rock properties,showing a certain linear relationship with the C_(ERP)at the same depth.In addition,the overall change trend of σ_(H) determined by the established fuzzy identification method is to increase with the increase of C_(ERP).The fuzzy identification method also demonstrates a relatively detailed local relationship betweenσ_H and C_(ERP),and the predicted curve rises in a fluctuating way,which is in accord well with the measured stress data.