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

Title: Curved-line search algorithm for ab initio atomic structure relaxation

; ; ;
Publication Date:
Sponsoring Org.:
OSTI Identifier:
Grant/Contract Number:
Resource Type:
Journal Article: Publisher's Accepted Manuscript
Journal Name:
Physical Review B
Additional Journal Information:
Journal Volume: 96; Journal Issue: 11; Related Information: CHORUS Timestamp: 2017-09-21 13:00:30; Journal ID: ISSN 2469-9950
American Physical Society
Country of Publication:
United States

Citation Formats

Chen, Zhanghui, Li, Jingbo, Li, Shushen, and Wang, Lin-Wang. Curved-line search algorithm for ab initio atomic structure relaxation. United States: N. p., 2017. Web. doi:10.1103/PhysRevB.96.115141.
Chen, Zhanghui, Li, Jingbo, Li, Shushen, & Wang, Lin-Wang. Curved-line search algorithm for ab initio atomic structure relaxation. United States. doi:10.1103/PhysRevB.96.115141.
Chen, Zhanghui, Li, Jingbo, Li, Shushen, and Wang, Lin-Wang. 2017. "Curved-line search algorithm for ab initio atomic structure relaxation". United States. doi:10.1103/PhysRevB.96.115141.
title = {Curved-line search algorithm for ab initio atomic structure relaxation},
author = {Chen, Zhanghui and Li, Jingbo and Li, Shushen and Wang, Lin-Wang},
abstractNote = {},
doi = {10.1103/PhysRevB.96.115141},
journal = {Physical Review B},
number = 11,
volume = 96,
place = {United States},
year = 2017,
month = 9

Journal Article:
Free Publicly Available Full Text
This content will become publicly available on September 21, 2018
Publisher's Accepted Manuscript

Save / Share:
  • The numerical integration of constitutive models in computational solid mechanics codes allows for the solution of boundary value problems involving complex material behavior. Metal plasticity models, in particular, have been instrumental in the development of these codes. Most plasticity models implemented in computational codes use an isotropic von Mises yield surface. The von Mises, or J 2, model uses a predictor–corrector algorithm–the radial return algorithm–to integrate the model. For non-quadratic yield surfaces, including anisotropic yield surfaces, no simple algorithm exists. This paper presents and analyzes a line search algorithm for the return mapping problem that shows excellent improvement over amore » Newton–Raphson model. Two non-quadratic yield surfaces–one isotropic and one anisotropic–are studied in this paper. The line search algorithm used for integrating the models is shown to be reliable and robust. The theory and implementation of the models, the details of the return mapping algorithm, and results that show the effectiveness of the method are presented. Lastly, a few simple boundary value problems verify the implementation and show the impact of the models. For the internal pressurization of a cylinder, the importance of modeling anisotropy correctly is shown.« less
  • The most important elements of {ital ab} {ital initio} calculations of x-ray-absorption fine structure (XAFS) are studied. To obtain accurate results without {ital ad} {ital hoc} adjustable parameters, we find it essential to include (i) curved-wave effects, (ii) a complex, energy-dependent self-energy, (iii) an approximate molecular potential, and (iv) a fixed energy reference for the photoelectron wave number. Based on these findings, an automated code has been developed for {ital ab} {ital initio} calculations of single-scattering XAFS, in which curved-wave effects are treated exactly in terms of effective backscattering amplitudes, inelastic losses and self-energy shifts are incorporated with use ofmore » a Hedin-Lundqvist self-energy, an automated relativistic overlapping-atom muffin-tin potential is used, and the energy threshold is estimated from electron-gas theory. The efficiency of the code is made possible by analytic expressions for the Hedin-Lundqvist self-energy. This code replaces existing tables of XAFS phases and scattering amplitudes and yields reliable theoretical XAFS standards for arbitrary pairs of atoms throughout the Periodic Table ({ital Z}{le}94). These results are comparable to those from self-consistent calculations and are valid to within about 20 eV of the absorption edge. Comparisons with experiment are presented for Cu, Ge, Pt, Br{sub 2}, and GeCl{sub 4}. The calculated XAFS amplitudes are found to be accurate to within 15%; XAFS phases are accurate to within 0.2 rad; and nearest-neighbor distances are typically accurate to within 0.02 A.« less
  • The first application of the ab initio genetic algorithm with an embedded gradient has been carried out for the elucidation of global minimum structures of a series of anionic sodium chloride clusters, NaxClx+1- (x = 1-4), produced in the gas-phase using electrospray ionization and studied by photoelectron spectroscopy. These are all superhalogen species with extremely high electron binding energies. The adiabatic electron binding energies for NaxClx+1- were measured to be 5.64, 6.22, 6.3, and 6.9 eV, for x = 1-4, respectively. Our genetic algorithm program detected the linear global minima of NaCl2- and Na2Cl3- and tree-dimensional structures for the largermore » species. Na3Cl4- was found to have C3v symmetry, which can be viewed as a Na4Cl4 cube missing a corner Na+ cation, whereas Na4Cl5- was found to have C4v symmetry, close to a 3x3 planar structure. Further accurate ab initio calculations were carried out for the elucidated global minimum structures. Excellent agreement between the theoretically calculated and the experimental spectra was observed, confirming the obtained structures and demonstrating the power of the developed genetic algorithm technique.« less
  • A novel maximum power point (MPP) search algorithm for photovoltaic (PV) array systems is introduced. The proposed algorithm determines the maximum power point of a PV array for any temperature and solar irradiation level using an on-line procedure. The method needs only the on-line values of the PV array output voltage and current, which can be obtained easily by using just current and voltage transducers. The algorithm requires neither the measurement of temperature and solar irradiation level nor a PV array model that is mostly used in look-up table based algorithms. Satisfactory results were obtained with the proposed algorithm inmore » a laboratory prototype implementation scheme consisting of a PV array computer emulation model, a chopper controlled permanent magnet dc motor, and a DT2827 data acquisition board with the ATLAB software drivers for interfacing.« less
  • We present a filter line-search algorithm that does not require inertia information of the linear system. This feature enables the use of a wide range of linear algebra strategies and libraries, which is essential to tackle large-scale problems on modern computing architectures. The proposed approach performs curvature tests along the search step to detect negative curvature and to trigger convexification. We prove that the approach is globally convergent and we implement the approach within a parallel interior-point framework to solve large-scale and highly nonlinear problems. Our numerical tests demonstrate that the inertia-free approach is as efficient as inertia detection viamore » symmetric indefinite factorizations. We also demonstrate that the inertia-free approach can lead to reductions in solution time because it reduces the amount of convexification needed.« less