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 solve...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.展开更多
Direct numerical simulations are carried out with different disturbance forms introduced into the inlet of a flat plate boundary layer with the Mach number 4.5. According to the biorthogonal eigenfunction system of th...Direct numerical simulations are carried out with different disturbance forms introduced into the inlet of a flat plate boundary layer with the Mach number 4.5. According to the biorthogonal eigenfunction system of the linearized Navier-Stokes equations and the adjoint equations, the decomposition of the direct numerical simulation results into the discrete normal mode is easily realized. The decomposition coefficients can be solved by doing the inner product between the numerical results and the eigenfunctions of the adjoint equations. For the quadratic polynomial eigenvalue problem, the inner product operator is given in a simple form, and it is extended to an Nth-degree polynomial eigenvalue problem. The examples illustrate that the simplified mode decomposition is available to analyze direct numerical simulation results.展开更多
The dynamic analysis of damped structural system by using finite element method leads to nonlinear eigenvalue problem(NEP)(particularly,quadratic eigenvalue problem).In general,the parameters of NEP are considered as ...The dynamic analysis of damped structural system by using finite element method leads to nonlinear eigenvalue problem(NEP)(particularly,quadratic eigenvalue problem).In general,the parameters of NEP are considered as exact values.But in actual practice because of different errors and incomplete information,the parameters may have uncertain or vague values and such uncertain values may be considered in terms of fuzzy numbers.This article proposes an efficient fuzzy-affine approach to solve fully fuzzy nonlinear eigenvalue problems(FNEPs)where involved parameters are fuzzy numbers viz.triangular and trapezoidal.Based on the parametric form,fuzzy numbers have been transformed into family of standard intervals.Further due to the presence of interval overestimation problem in standard interval arithmetic,affine arithmetic based approach has been implemented.In the proposed method,the FNEP has been linearized into a generalized eigenvalue problem and further solved by using the fuzzy-affine approach.Several application problems of structures and also general NEPs with fuzzy parameters are investigated based on the proposed procedure.Lastly,fuzzy eigenvalue bounds are illustrated with fuzzy plots with respect to its membership function.Few comparisons are also demonstrated to show the reliability and efficacy of the present approach.展开更多
We introduce a type of full multigrid method for the nonlinear eigenvalue problem. The main idea is to transform the solution of the nonlinear eigenvalue problem into a series of solutions of the corresponding linear ...We introduce a type of full multigrid method for the nonlinear eigenvalue problem. The main idea is to transform the solution of the nonlinear eigenvalue problem into a series of solutions of the corresponding linear boundary value problems on the sequence of finite element spaces and nonlinear eigenvalue problems on the coarsest finite element space. The linearized boundary value problems are solved by some multigrid iterations.Besides the multigrid iteration, all other efficient iteration methods for solving boundary value problems can serve as the linear problem solver. We prove that the computational work of this new scheme is truly optimal,the same as solving the linear corresponding boundary value problem. In this case, this type of iteration scheme certainly improves the overfull efficiency of solving nonlinear eigenvalue problems. Some numerical experiments are presented to validate the efficiency of the new method.展开更多
We study the smooth LU decomposition of a given analytic functional A-matrix A(A) and its block-analogue. Sufficient conditions for the existence of such matrix decompositions are given, some differentiability about...We study the smooth LU decomposition of a given analytic functional A-matrix A(A) and its block-analogue. Sufficient conditions for the existence of such matrix decompositions are given, some differentiability about certain elements arising from them are proved, and several explicit expressions for derivatives of the specified elements are provided. By using these smooth LU decompositions, we propose two numerical methods for computing multiple nonlinear eigenvalues of A(A), and establish their locally quadratic convergence properties. Several numerical examples are provided to show the feasibility and effectiveness of these new methods.展开更多
In this paper we consider the existence of infinitely many positive solutions for second order nonlinear eigenvalue problem with singular coefficient function. By the use of Krasnosel'skii fixed point theorem of c...In this paper we consider the existence of infinitely many positive solutions for second order nonlinear eigenvalue problem with singular coefficient function. By the use of Krasnosel'skii fixed point theorem of cone expansion-compression type we give several sufficient conditions.展开更多
Nonlinear rank-one modification of the symmetric eigenvalue problem arises from eigenvibrations of mechanical structures with elastically attached loads and calculation of the propagation modes in optical fiber. In th...Nonlinear rank-one modification of the symmetric eigenvalue problem arises from eigenvibrations of mechanical structures with elastically attached loads and calculation of the propagation modes in optical fiber. In this paper, we first study the existence and uniqueness of eigenvalues, and then investigate three numerical algorithms, namely Picard iteration, nonlinear Rayleigh quotient iteration and successive linear approximation method (SLAM). The global convergence of the SLAM is proven under some mild assumptions. Numerical examples illustrate that the SLAM is the most robust method.展开更多
The transmission eigenvalue problem is an eigenvalue problem that arises in the scatter- ing of time-harmonic waves by an inhomogeneous medium of compact support. Based on a fourth order formulation, the transmission ...The transmission eigenvalue problem is an eigenvalue problem that arises in the scatter- ing of time-harmonic waves by an inhomogeneous medium of compact support. Based on a fourth order formulation, the transmission eigenvalue problem is discretized by the Mor- ley element. For the resulting quadratic eigenvalue problem, a recursive integral method is used to compute real and complex eigenvalues in prescribed regions in the complex plane. Numerical examples are presented to demonstrate the effectiveness of the proposed method.展开更多
In this paper we propose a uniformly convergent numerical method for discretizing singularly perturbed nonlinear eigenvalue problems under constraints with applications in Bose-Einstein condensation and quantum chemis...In this paper we propose a uniformly convergent numerical method for discretizing singularly perturbed nonlinear eigenvalue problems under constraints with applications in Bose-Einstein condensation and quantum chemistry.We begin with the time-independent Gross-Pitaevskii equation and show how to reformulate it into a singularly perturbed nonlinear eigenvalue problem under a constraint.Matched asymptotic approximations for the problem are presented to locate the positions and characterize the widths of boundary layers and/or interior layers in the solution.A uniformly convergent numerical method is proposed by using the normalized gradient flow and piecewise uniform mesh techniques based on the asymptotic approximations for the problem.Extensive numerical results are reported to demonstrate the effectiveness of our numerical method for the problems.Finally,the method is applied to compute ground and excited states of Bose-Einstein condensation in the semiclassical regime and some conclusive findings are reported.展开更多
In this paper, we investigate a fourth-order three-point nonlinear eigenvalue problem. By the Krasnosel’skii’s fixed-point theorem in cone, some new results on the existence of positive solutions are obtained.
In this paper we study the existence of nontrivial solutions of a class of asymptotically linear elliptic resonant problems at higher eigenvalues with the nonlinear term which may be un- bounded by making use of the M...In this paper we study the existence of nontrivial solutions of a class of asymptotically linear elliptic resonant problems at higher eigenvalues with the nonlinear term which may be un- bounded by making use of the Morse theory for a C^2-function at both isolated critical point and infinity.展开更多
The focus of this paper is on determination of the dynamic parameters of structural systems with viscoelastic (VE) dampers described by Maxwell rheological models. Such parameters could be obtained after solving the a...The focus of this paper is on determination of the dynamic parameters of structural systems with viscoelastic (VE) dampers described by Maxwell rheological models. Such parameters could be obtained after solving the appropriately defined nonlinear eigenvalue problem for frames with VE dampers. The solution to the nonlinear eigenvalue problem is obtained by equating to zero the determinant of the considered system of equations. Apart from complex conjugate eigenvalues, the real ones occurred when dampers that are described by the classic Maxwell model, are also determined.展开更多
We propose a finite element method to compute the band structures of dispersive photonic crystals in 3D.The nonlinear Maxwell’s eigenvalue problem is formulated as the eigenvalue problem of a holomorphic operator fun...We propose a finite element method to compute the band structures of dispersive photonic crystals in 3D.The nonlinear Maxwell’s eigenvalue problem is formulated as the eigenvalue problem of a holomorphic operator function.The N´ed´elec edge elements are employed to discretize the operators,where the divergence free condition for the electric field is realized by a mixed form using a Lagrange multiplier.The convergence of the eigenvalues is proved using the abstract approximation theory for holomorphic operator functions with the regular approximation of the edge elements.The spectral indicator method is then applied to compute the discrete eigenvalues.Numerical examples are presented demonstrating the effectiveness of the proposed method.展开更多
The necessary condition established in Part I of this paper for the global maximizers of the maximization problem max V tr(VTAV)/tr(VTBV)+tr(VTCV)over the Stiefel manifold{V∈Rm×l |VTV=Il}(l〈m),natural...The necessary condition established in Part I of this paper for the global maximizers of the maximization problem max V tr(VTAV)/tr(VTBV)+tr(VTCV)over the Stiefel manifold{V∈Rm×l |VTV=Il}(l〈m),naturally leads to a self-consistent-field(SCF)iteration for computing a maximizer.In this part,we analyze the global and local convergence of the SCF iteration,and show that the necessary condition for the global maximizers is fulfilled at any convergent point of the sequences of approximations generated by the SCF iteration.This is one of the advantages of the SCF iteration over optimization-based methods.Preliminary numerical tests are reported and show that the SCF iteration is very efficient by comparing with some manifold-based optimization methods.展开更多
We are concerned with the maximization of tr(V T AV)/tr(V T BV)+tr(V T CV) over the Stiefel manifold {V ∈ R m×l | V T V = Il} (l 〈 m), where B is a given symmetric and positive definite matrix, A and...We are concerned with the maximization of tr(V T AV)/tr(V T BV)+tr(V T CV) over the Stiefel manifold {V ∈ R m×l | V T V = Il} (l 〈 m), where B is a given symmetric and positive definite matrix, A and C are symmetric matrices, and tr(. ) is the trace of a square matrix. This is a subspace version of the maximization problem studied in Zhang (2013), which arises from real-world applications in, for example, the downlink of a multi-user MIMO system and the sparse Fisher discriminant analysis in pattern recognition. We establish necessary conditions for both the local and global maximizers and connect the problem with a nonlinear extreme eigenvalue problem. The necessary condition for the global maximizers offers deep insights into the problem, on the one hand, and, on the other hand, naturally leads to a self-consistent-field (SCF) iteration to be presented and analyzed in detail in Part II of this paper.展开更多
The(continuous) finite element approximations of different orders for the computation of the solution to electronic structures were proposed in some papers and the performance of these approaches is becoming appreciab...The(continuous) finite element approximations of different orders for the computation of the solution to electronic structures were proposed in some papers and the performance of these approaches is becoming appreciable and is now well understood.In this publication,the author proposes to extend this discretization for full-potential electronic structure calculations by combining the refinement of the finite element mesh,where the solution is most singular with the increase of the degree of the polynomial approximations in the regions where the solution is mostly regular.This combination of increase of approximation properties,done in an a priori or a posteriori manner,is well-known to generally produce an optimal exponential type convergence rate with respect to the number of degrees of freedom even when the solution is singular.The analysis performed here sustains this property in the case of Hartree-Fock and Kohn-Sham problems.展开更多
In this work the one-band effective Hamiltonian governing the electronic states of a quantum dot/ring in a homogenous magnetic field is used to derive a pair/quadruple of nonlinear eigenvalue problems corresponding to...In this work the one-band effective Hamiltonian governing the electronic states of a quantum dot/ring in a homogenous magnetic field is used to derive a pair/quadruple of nonlinear eigenvalue problems corresponding to different spin orientations and in case of rotational symmetry additionally to quantum number±ℓ.We show,that each of those pair/quadruple of nonlinear problems allows for the minmax characterization of its eigenvalues under certain conditions,which are satisfied for our examples and the common InAs/GaAs heterojunction.Exploiting the minmax property we devise efficient iterative projection methods simultaneously handling the pair/quadruple of nonlinear problems and thereby saving up to 40%of the computational time as compared to the nonlinear Arnoldi method applied to each of the problems separately.展开更多
In this paper we apply a cone theoretic fixed-point theorem proved by Krasnosel'skii to obtain sufficient conditions for the existence of positive solutions to some boundary value problems for a class of functional d...In this paper we apply a cone theoretic fixed-point theorem proved by Krasnosel'skii to obtain sufficient conditions for the existence of positive solutions to some boundary value problems for a class of functional difference equations. We consider the case that the nonlinear term satisfies asvrnntoticallv linear growth.展开更多
基金the National Science and Tech-nology Council,Taiwan for their financial support(Grant Number NSTC 111-2221-E-019-048).
文摘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.
基金supported by the National Natural Science Foundation of China(Nos.1133200711202147+2 种基金and 9216111)the Specialized Research Fund for the Doctoral Program of Higher Education(No.20120032120007)the Open Fund from State Key Laboratory of Aerodynamics(Nos.SKLA201201 and SKLA201301)
文摘Direct numerical simulations are carried out with different disturbance forms introduced into the inlet of a flat plate boundary layer with the Mach number 4.5. According to the biorthogonal eigenfunction system of the linearized Navier-Stokes equations and the adjoint equations, the decomposition of the direct numerical simulation results into the discrete normal mode is easily realized. The decomposition coefficients can be solved by doing the inner product between the numerical results and the eigenfunctions of the adjoint equations. For the quadratic polynomial eigenvalue problem, the inner product operator is given in a simple form, and it is extended to an Nth-degree polynomial eigenvalue problem. The examples illustrate that the simplified mode decomposition is available to analyze direct numerical simulation results.
文摘The dynamic analysis of damped structural system by using finite element method leads to nonlinear eigenvalue problem(NEP)(particularly,quadratic eigenvalue problem).In general,the parameters of NEP are considered as exact values.But in actual practice because of different errors and incomplete information,the parameters may have uncertain or vague values and such uncertain values may be considered in terms of fuzzy numbers.This article proposes an efficient fuzzy-affine approach to solve fully fuzzy nonlinear eigenvalue problems(FNEPs)where involved parameters are fuzzy numbers viz.triangular and trapezoidal.Based on the parametric form,fuzzy numbers have been transformed into family of standard intervals.Further due to the presence of interval overestimation problem in standard interval arithmetic,affine arithmetic based approach has been implemented.In the proposed method,the FNEP has been linearized into a generalized eigenvalue problem and further solved by using the fuzzy-affine approach.Several application problems of structures and also general NEPs with fuzzy parameters are investigated based on the proposed procedure.Lastly,fuzzy eigenvalue bounds are illustrated with fuzzy plots with respect to its membership function.Few comparisons are also demonstrated to show the reliability and efficacy of the present approach.
基金supported by National Natural Science Foundation of China (Grant Nos. 91330202, 11371026, 11201501, 11571389, 11001259 and 11031006)National Basic Research Program of China (Grant No. 2011CB309703)the National Center for Mathematics and Interdisciplinary Science, Chinese Academy of Sciences, the President Foundation of Academy of Mathematics and Systems Science, Chinese Academy of Sciences and the Program for Innovation Research in Central University of Finance and Economics
文摘We introduce a type of full multigrid method for the nonlinear eigenvalue problem. The main idea is to transform the solution of the nonlinear eigenvalue problem into a series of solutions of the corresponding linear boundary value problems on the sequence of finite element spaces and nonlinear eigenvalue problems on the coarsest finite element space. The linearized boundary value problems are solved by some multigrid iterations.Besides the multigrid iteration, all other efficient iteration methods for solving boundary value problems can serve as the linear problem solver. We prove that the computational work of this new scheme is truly optimal,the same as solving the linear corresponding boundary value problem. In this case, this type of iteration scheme certainly improves the overfull efficiency of solving nonlinear eigenvalue problems. Some numerical experiments are presented to validate the efficiency of the new method.
基金supported by the National Basic Research Program(No.2005CB321702)the China Outstanding Young Scientist F0undation(No.10525102)the National Natural Science Foundation (No.10471146),P.R.China
文摘We study the smooth LU decomposition of a given analytic functional A-matrix A(A) and its block-analogue. Sufficient conditions for the existence of such matrix decompositions are given, some differentiability about certain elements arising from them are proved, and several explicit expressions for derivatives of the specified elements are provided. By using these smooth LU decompositions, we propose two numerical methods for computing multiple nonlinear eigenvalues of A(A), and establish their locally quadratic convergence properties. Several numerical examples are provided to show the feasibility and effectiveness of these new methods.
文摘In this paper we consider the existence of infinitely many positive solutions for second order nonlinear eigenvalue problem with singular coefficient function. By the use of Krasnosel'skii fixed point theorem of cone expansion-compression type we give several sufficient conditions.
基金supported in part by NSF grants DMS-0611548 and OCI-0749217 and DOE grant DE-FC02-06ER25794supported in part by NSF of China under the contract number 10871049 and Shanghai Down project 200601.
文摘Nonlinear rank-one modification of the symmetric eigenvalue problem arises from eigenvibrations of mechanical structures with elastically attached loads and calculation of the propagation modes in optical fiber. In this paper, we first study the existence and uniqueness of eigenvalues, and then investigate three numerical algorithms, namely Picard iteration, nonlinear Rayleigh quotient iteration and successive linear approximation method (SLAM). The global convergence of the SLAM is proven under some mild assumptions. Numerical examples illustrate that the SLAM is the most robust method.
文摘The transmission eigenvalue problem is an eigenvalue problem that arises in the scatter- ing of time-harmonic waves by an inhomogeneous medium of compact support. Based on a fourth order formulation, the transmission eigenvalue problem is discretized by the Mor- ley element. For the resulting quadratic eigenvalue problem, a recursive integral method is used to compute real and complex eigenvalues in prescribed regions in the complex plane. Numerical examples are presented to demonstrate the effectiveness of the proposed method.
基金Singapore Ministry of Education grant No.R-146-000-083-112 and would like to thank Professor Tao Tang for very helpful discussion on the subject.
文摘In this paper we propose a uniformly convergent numerical method for discretizing singularly perturbed nonlinear eigenvalue problems under constraints with applications in Bose-Einstein condensation and quantum chemistry.We begin with the time-independent Gross-Pitaevskii equation and show how to reformulate it into a singularly perturbed nonlinear eigenvalue problem under a constraint.Matched asymptotic approximations for the problem are presented to locate the positions and characterize the widths of boundary layers and/or interior layers in the solution.A uniformly convergent numerical method is proposed by using the normalized gradient flow and piecewise uniform mesh techniques based on the asymptotic approximations for the problem.Extensive numerical results are reported to demonstrate the effectiveness of our numerical method for the problems.Finally,the method is applied to compute ground and excited states of Bose-Einstein condensation in the semiclassical regime and some conclusive findings are reported.
基金the National Natural Science Foundation of China(11261053)The Natural Science Foundation of Gansu Province(1308RJZA125)
文摘In this paper, we investigate a fourth-order three-point nonlinear eigenvalue problem. By the Krasnosel’skii’s fixed-point theorem in cone, some new results on the existence of positive solutions are obtained.
文摘In this paper we study the existence of nontrivial solutions of a class of asymptotically linear elliptic resonant problems at higher eigenvalues with the nonlinear term which may be un- bounded by making use of the Morse theory for a C^2-function at both isolated critical point and infinity.
基金the financial support received from the Poznan University of Technology(Grant No.DS 11-088/12)in connection with this work.
文摘The focus of this paper is on determination of the dynamic parameters of structural systems with viscoelastic (VE) dampers described by Maxwell rheological models. Such parameters could be obtained after solving the appropriately defined nonlinear eigenvalue problem for frames with VE dampers. The solution to the nonlinear eigenvalue problem is obtained by equating to zero the determinant of the considered system of equations. Apart from complex conjugate eigenvalues, the real ones occurred when dampers that are described by the classic Maxwell model, are also determined.
基金China Postdoctoral Science Foundation Grant 2019M650460the NSF grant DMS-2011148.The research of J.Sun is supported partially by the Simons Foundation Grant 711922.
文摘We propose a finite element method to compute the band structures of dispersive photonic crystals in 3D.The nonlinear Maxwell’s eigenvalue problem is formulated as the eigenvalue problem of a holomorphic operator function.The N´ed´elec edge elements are employed to discretize the operators,where the divergence free condition for the electric field is realized by a mixed form using a Lagrange multiplier.The convergence of the eigenvalues is proved using the abstract approximation theory for holomorphic operator functions with the regular approximation of the edge elements.The spectral indicator method is then applied to compute the discrete eigenvalues.Numerical examples are presented demonstrating the effectiveness of the proposed method.
基金Acknowledgements The first author was supported by National Natural Science Foundation of China(Grant Nos.11101257 and 11371102)the Basic Academic Discipline Program,the 11th Five Year Plan of 211 Project for Shanghai University of Finance and Economics+1 种基金supported by National Science Foundation of USA(Grant Nos.1115834and 1317330)a Research Gift Grant from Intel Corporation
文摘The necessary condition established in Part I of this paper for the global maximizers of the maximization problem max V tr(VTAV)/tr(VTBV)+tr(VTCV)over the Stiefel manifold{V∈Rm×l |VTV=Il}(l〈m),naturally leads to a self-consistent-field(SCF)iteration for computing a maximizer.In this part,we analyze the global and local convergence of the SCF iteration,and show that the necessary condition for the global maximizers is fulfilled at any convergent point of the sequences of approximations generated by the SCF iteration.This is one of the advantages of the SCF iteration over optimization-based methods.Preliminary numerical tests are reported and show that the SCF iteration is very efficient by comparing with some manifold-based optimization methods.
基金supported by National Natural Science Foundation of China(Grant Nos.11101257 and 11371102)the Basic Academic Discipline Program+3 种基金the 11th Five Year Plan of 211 Project for Shanghai University of Finance and Economicsa visiting scholar at the Department of Mathematics,University of Texas at Arlington from February 2013 toJanuary 2014supported by National Science Foundation of USA(Grant Nos.1115834and 1317330)a Research Gift Grant from Intel Corporation
文摘We are concerned with the maximization of tr(V T AV)/tr(V T BV)+tr(V T CV) over the Stiefel manifold {V ∈ R m×l | V T V = Il} (l 〈 m), where B is a given symmetric and positive definite matrix, A and C are symmetric matrices, and tr(. ) is the trace of a square matrix. This is a subspace version of the maximization problem studied in Zhang (2013), which arises from real-world applications in, for example, the downlink of a multi-user MIMO system and the sparse Fisher discriminant analysis in pattern recognition. We establish necessary conditions for both the local and global maximizers and connect the problem with a nonlinear extreme eigenvalue problem. The necessary condition for the global maximizers offers deep insights into the problem, on the one hand, and, on the other hand, naturally leads to a self-consistent-field (SCF) iteration to be presented and analyzed in detail in Part II of this paper.
文摘The(continuous) finite element approximations of different orders for the computation of the solution to electronic structures were proposed in some papers and the performance of these approaches is becoming appreciable and is now well understood.In this publication,the author proposes to extend this discretization for full-potential electronic structure calculations by combining the refinement of the finite element mesh,where the solution is most singular with the increase of the degree of the polynomial approximations in the regions where the solution is mostly regular.This combination of increase of approximation properties,done in an a priori or a posteriori manner,is well-known to generally produce an optimal exponential type convergence rate with respect to the number of degrees of freedom even when the solution is singular.The analysis performed here sustains this property in the case of Hartree-Fock and Kohn-Sham problems.
基金We would like to thank Oleksandr Voskoboynikov for his comments on the physical relevance of the model under consideration.We also thank the anonymous referees for their comments helping us to improve this manuscript.
文摘In this work the one-band effective Hamiltonian governing the electronic states of a quantum dot/ring in a homogenous magnetic field is used to derive a pair/quadruple of nonlinear eigenvalue problems corresponding to different spin orientations and in case of rotational symmetry additionally to quantum number±ℓ.We show,that each of those pair/quadruple of nonlinear problems allows for the minmax characterization of its eigenvalues under certain conditions,which are satisfied for our examples and the common InAs/GaAs heterojunction.Exploiting the minmax property we devise efficient iterative projection methods simultaneously handling the pair/quadruple of nonlinear problems and thereby saving up to 40%of the computational time as compared to the nonlinear Arnoldi method applied to each of the problems separately.
文摘In this paper we apply a cone theoretic fixed-point theorem proved by Krasnosel'skii to obtain sufficient conditions for the existence of positive solutions to some boundary value problems for a class of functional difference equations. We consider the case that the nonlinear term satisfies asvrnntoticallv linear growth.