skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Coordinate-space solver for superfluid many-fermion systems with the shifted conjugate-orthogonal conjugate-gradient method

Authors:
; ; ;
Publication Date:
Sponsoring Org.:
USDOE
OSTI Identifier:
1349702
Grant/Contract Number:
FG02-97ER41014; AC05-00OR22725; AC02-05CH11231
Resource Type:
Journal Article: Publisher's Accepted Manuscript
Journal Name:
Physical Review C
Additional Journal Information:
Journal Volume: 95; Journal Issue: 4; Related Information: CHORUS Timestamp: 2017-06-25 10:07:34; Journal ID: ISSN 2469-9985
Publisher:
American Physical Society
Country of Publication:
United States
Language:
English

Citation Formats

Jin, Shi, Bulgac, Aurel, Roche, Kenneth, and Wlazłowski, Gabriel. Coordinate-space solver for superfluid many-fermion systems with the shifted conjugate-orthogonal conjugate-gradient method. United States: N. p., 2017. Web. doi:10.1103/PhysRevC.95.044302.
Jin, Shi, Bulgac, Aurel, Roche, Kenneth, & Wlazłowski, Gabriel. Coordinate-space solver for superfluid many-fermion systems with the shifted conjugate-orthogonal conjugate-gradient method. United States. doi:10.1103/PhysRevC.95.044302.
Jin, Shi, Bulgac, Aurel, Roche, Kenneth, and Wlazłowski, Gabriel. Mon . "Coordinate-space solver for superfluid many-fermion systems with the shifted conjugate-orthogonal conjugate-gradient method". United States. doi:10.1103/PhysRevC.95.044302.
@article{osti_1349702,
title = {Coordinate-space solver for superfluid many-fermion systems with the shifted conjugate-orthogonal conjugate-gradient method},
author = {Jin, Shi and Bulgac, Aurel and Roche, Kenneth and Wlazłowski, Gabriel},
abstractNote = {},
doi = {10.1103/PhysRevC.95.044302},
journal = {Physical Review C},
number = 4,
volume = 95,
place = {United States},
year = {Mon Apr 03 00:00:00 EDT 2017},
month = {Mon Apr 03 00:00:00 EDT 2017}
}

Journal Article:
Free Publicly Available Full Text
Publisher's Version of Record at 10.1103/PhysRevC.95.044302

Citation Metrics:
Cited by: 3works
Citation information provided by
Web of Science

Save / Share:
  • Here, we present an iterative algorithm for computing an invariant subspace associated with the algebraically smallest eigenvalues of a large sparse or structured Hermitian matrix A. We are interested in the case in which the dimension of the invariant subspace is large (e.g., over several hundreds or thousands) even though it may still be small relative to the dimension of A. These problems arise from, for example, density functional theory (DFT) based electronic structure calculations for complex materials. The key feature of our algorithm is that it performs fewer Rayleigh–Ritz calculations compared to existing algorithms such as the locally optimalmore » block preconditioned conjugate gradient or the Davidson algorithm. It is a block algorithm, and hence can take advantage of efficient BLAS3 operations and be implemented with multiple levels of concurrency. We discuss a number of practical issues that must be addressed in order to implement the algorithm efficiently on a high performance computer.« less
  • The effects of possible superfluidity on the cranking moment of a lange many-fermion system moving under pericdic boundary conditions were investigated within the framework of the theory of superconductivity recently formulated by Bogolyubov. The Hamiltonian was initially subjected to Bogolyubov' s general unitary transformation. The coliective excitations of the fermions were then considered in the usual pnir approximation; the appropriate cranking terms were linear in the boson pair operators. On performing a unitary transformation which transformed away these linear terms, an expression was obtained for the momert of inertia of the system which included both the effects of pcssible superfluiditymore » and collective excitation. This expression, by virtue of its being stationary with respect to arbitrary variations in the amplitude associated with the latter unitary transformation, was then utilized as a variational principle for the moment of inertia. The result previousiy obtained that the moment of inertia had a rigid value, was rederived in more compact form. For the superfluid state, it was found that collective excitations effected a marked increase in the superfluid moment at intermediate coupling strengths although the resulting moment was still quite small compared to the rigid value. In contrast to the normal state case, where particle-hole pairs play a major role, this increase was almost entirely due to excitations consisting of particle-pairs or hole-pairs. The precise magnitude of the apparent resonance in the moment produced by the d- wave part of the cranking interaction was dependent to some extent on the features of the particle-particle potential which led to the superffuid state. Variational expressions for the moment were exhibited for both Yukawa and deltafunction shell potentials. These results were identical in charged and neutral Fermi systems. A calculation of the cranking moment at finite temperatures was presented in an Appendix along with an interpretation of it in terms of Bardeen's two-fluid model of superconductivity. (auth)« less
  • The technique of the initial guess calculation for the conjugate gradient method is proposed. Computational schemes of the linear system solution with symmetrical positive definite matrices are constructed on its basis. Their efficient modifications for systems with five-diagonal matrices are proposed. The investigation of the developed methods using the problem of two-dimensional numerical simulation of bipolar transistors has been carried out. Experimental evidence of the proposed method's efficiency has been obtained. 30 refs., 6 figs., 1 tab.
  • Difficulties in the straightfoward application of the complex-coordinate method to the calculation of resonance states in many-electron systems are examined. For the case of shape resonances, it is shown that many of these difficulties can be avoided by using complex coordinates only after reduction of the system to an effective one-electron problem. Further simplifications are achieved by the use of an inner-projection technique to facilitate the computation of the complex Hamiltonian matrix elements. The method is first illustrated by application to a model-potential problem. Its suitability for studying many-electron problems is demonstrated by calculation of the position and width ofmore » a low-energy /sup 2/P/sup o/ shape resonance in Be/sup -/. We discuss the modifications necessary to study core-excited (Feshbach) resonances.« less