A Chebyshev polynomial method for computing analytic solutions to eigenvalue problems with application to the anharmonic oscillator

1978 ◽  
Vol 19 (6) ◽  
pp. 1445-1456 ◽  
Author(s):  
John P. Boyd
2008 ◽  
Vol 86 (4) ◽  
pp. 597-600 ◽  
Author(s):  
S R Valluri ◽  
M D Fried

The master equation described by Badri Krishnan et al. (Phys Rev. D, 70, 082001 (2004)) for the time-frequency pattern using the F-statistic is studied in the context of Chebyshev-polynomial modified spindown parameters for the case of gravitational wave pulsar signals. The Chebyshev-polynomial method enables an analytic and numeric evaluation of the Fourier transform (FT) for both the non-demodulated and F-statistic demodulated FT.PACS Nos.: 04.30.Tv, 95.85.sz, 02.30.Gp, 02.40.Re


Accurate eigenvalues and eigenfunctions of the anharm onic oscillator ( H = p 2 + x 2 + λx 4 , λ > 0) and the quartic oscillator ( H = p 2 + x 4 ) are obtained in all regimes of the quantum num ber n and the anharm onicity constant λ. Transition moments of comparable accuracy are obtained for the quartic oscillator. The method, applicable quite generally for eigenvalue problems, is non-perturbative and involves the use of an appropriately scaled basis for the determ ination of each eigenvalue. The appropriate scaling formula for a given regime of ( n , λ) is constructed from the oscillation properties of the eigenfunctions. More general anharm onic oscillators are also discussed.


2021 ◽  
Vol 247 ◽  
pp. 04008
Author(s):  
F. Filiciotto ◽  
A. Jinaphanh ◽  
A. Zoia

Time eigenvalues emerge in several key applications related to neutron transport problems, including reactor start-up and reactivity measurements. In this context, experimental validation and uncertainty quantification would demand to assess the variation of the dominant time eigenvalue in response to a variation of nuclear data. Recently, we proposed the use of a Generalized Iterated Fission Probability method (G-IFP) to compute adjoint-weighted tallies, such as kinetic parameters, perturbations and sensitivity coefficients, for Monte Carlo time (or alpha) eigenvalue calculations. With the massive use of parallel Monte Carlo calculations, it would be therefore useful to trade the memory burden of the G-IFP method (which is comparable to that of the standard IFP method for k-eigenvalue problems) for computation time and to rely on history-based schemes for such adjoint-weighted tallies. For this purpose, we investigate the use of the super-history method as applied to estimating adjoint-weighted tallies within the α-k power iteration, based on previous work on k-eigenvalue problems. Verification of the algorithms is performed on some simple preliminary tests where analytic solutions exist. In addition, the performances of the proposed method are assessed by comparing the super-history and the G-IFP methods for the same sets of benchmark problems.


2021 ◽  
Vol 11 (1) ◽  
Author(s):  
Xinran Zheng ◽  
Mingqi Huang ◽  
Dongqi An ◽  
Chao Zhou ◽  
Rui Li

AbstractNew analytic bending, buckling, and free vibration solutions of rectangular nanoplates with combinations of clamped and simply supported edges are obtained by an up-to-date symplectic superposition method. The problems are reformulated in the Hamiltonian system and symplectic space, where the mathematical solution framework involves the construction of symplectic eigenvalue problems and symplectic eigen expansion. The analytic symplectic solutions are derived for several elaborated fundamental subproblems, the superposition of which yields the final analytic solutions. Besides Lévy-type solutions, non-Lévy-type solutions are also obtained, which cannot be achieved by conventional analytic methods. Comprehensive numerical results can provide benchmarks for other solution methods.


Sign in / Sign up

Export Citation Format

Share Document