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

Title: Local-order metric for condensed-phase environments

; ; ;
Publication Date:
Sponsoring Org.:
OSTI Identifier:
Grant/Contract Number:
SC0008626; AC02-05CH11231
Resource Type:
Journal Article: Publisher's Accepted Manuscript
Journal Name:
Physical Review B
Additional Journal Information:
Journal Volume: 97; Journal Issue: 6; Related Information: CHORUS Timestamp: 2018-02-12 10:17:22; Journal ID: ISSN 2469-9950
American Physical Society
Country of Publication:
United States

Citation Formats

Martelli, Fausto, Ko, Hsin-Yu, Oğuz, Erdal C., and Car, Roberto. Local-order metric for condensed-phase environments. United States: N. p., 2018. Web. doi:10.1103/PhysRevB.97.064105.
Martelli, Fausto, Ko, Hsin-Yu, Oğuz, Erdal C., & Car, Roberto. Local-order metric for condensed-phase environments. United States. doi:10.1103/PhysRevB.97.064105.
Martelli, Fausto, Ko, Hsin-Yu, Oğuz, Erdal C., and Car, Roberto. 2018. "Local-order metric for condensed-phase environments". United States. doi:10.1103/PhysRevB.97.064105.
title = {Local-order metric for condensed-phase environments},
author = {Martelli, Fausto and Ko, Hsin-Yu and Oğuz, Erdal C. and Car, Roberto},
abstractNote = {},
doi = {10.1103/PhysRevB.97.064105},
journal = {Physical Review B},
number = 6,
volume = 97,
place = {United States},
year = 2018,
month = 2

Journal Article:
Free Publicly Available Full Text
This content will become publicly available on February 12, 2019
Publisher's Accepted Manuscript

Save / Share:
  • The forces acting on the atoms as well as the stress tensor are crucial ingredients for calculating the structural and dynamical properties of systems in the condensed phase. Here, these derivatives of the total energy are evaluated for the second-order Møller-Plesset perturbation energy (MP2) in the framework of the resolution of identity Gaussian and plane waves method, in a way that is fully consistent with how the total energy is computed. This consistency is non-trivial, given the different ways employed to compute Coulomb, exchange, and canonical four center integrals, and allows, for example, for energy conserving dynamics in various ensembles.more » Based on this formalism, a massively parallel algorithm has been developed for finite and extended system. The designed parallel algorithm displays, with respect to the system size, cubic, quartic, and quintic requirements, respectively, for the memory, communication, and computation. All these requirements are reduced with an increasing number of processes, and the measured performance shows excellent parallel scalability and efficiency up to thousands of nodes. Additionally, the computationally more demanding quintic scaling steps can be accelerated by employing graphics processing units (GPU’s) showing, for large systems, a gain of almost a factor two compared to the standard central processing unit-only case. In this way, the evaluation of the derivatives of the RI-MP2 energy can be performed within a few minutes for systems containing hundreds of atoms and thousands of basis functions. With good time to solution, the implementation thus opens the possibility to perform molecular dynamics (MD) simulations in various ensembles (microcanonical ensemble and isobaric-isothermal ensemble) at the MP2 level of theory. Geometry optimization, full cell relaxation, and energy conserving MD simulations have been performed for a variety of molecular crystals including NH{sub 3}, CO{sub 2}, formic acid, and benzene.« less
  • Ge{sub 0.15}Sb{sub 0.85} is a phase change material that undergoes a transition from a semi-conducting glass to a crystalline metal when being densified. In this work, we investigate some parameters controlling the crystallization pressure P{sub c} in related compounds by conducting high-pressure experiments on initially amorphous, germanium- and silicon-doped antimony. We find that the amorphous phase is stabilized to higher pressures when the dopant size is decreased. This result can be easily rationalized under the assumption that in the glass, Ge and Si atoms occupy tetrahedral or related small coordination shells. When pressure increases, the larger Ge atoms move soonermore » into the larger octahedral shells, which are characteristic for the crystal, than the smaller Si atoms do. We also find that P{sub c} increases quickly with Ge and Si concentrations. This observation implies that the pressure-induced change of coordination cannot be a local, elementary event, but that the four-coordination of the group-14 atoms is stabilized by the presence of other four-coordinated atoms.« less
  • Using deep narrow-band (NB) imaging and optical spectroscopy from the Keck telescope and MMT, we identify a sample of 20 emission-line galaxies at z = 0.065-0.90 where the weak auroral emission line, [O III] λ4363, is detected at ≥3σ. These detections allow us to determine the gas-phase metallicity using the 'direct' method. With electron temperature measurements, and dust attenuation corrections from Balmer decrements, we find that four of these low-mass galaxies are extremely metal-poor with 12 + log (O/H) ≤ 7.65 or one-tenth solar. Our most metal-deficient galaxy has 12 + log (O/H) = 7.24{sub −0.30}{sup +0.45} (95% confidence), similarmore » to some of the lowest metallicity galaxies identified in the local universe. We find that our galaxies are all undergoing significant star formation with average specific star formation rate (SFR) of (100 Myr){sup –1}, and that they have high central SFR surface densities (average of 0.5 M {sub ☉} yr{sup –1} kpc{sup –2}). In addition, more than two-thirds of our galaxies have between one and four nearby companions within a projected radius of 100 kpc, which we find is an excess among star-forming galaxies at z = 0.4-0.85. We also find that the gas-phase metallicities for a given stellar mass and SFR lie systematically lower than the local M {sub *}-Z-(SFR) relation by ≈0.2 dex (2σ significance). These results are partly due to selection effects, since galaxies with strong star formation and low metallicity are more likely to yield [O III] λ4363 detections. Finally, the observed higher ionization parameter and high electron density suggest that they are lower redshift analogs to typical z ≳ 1 galaxies.« less
  • We prove the existence of a Foldy-Wouthuysen (FW) transformation which decouples the mass states of all the Bhabha Poincare generators and diagonalizes the Hamiltonian. Since the Bhabha fields operate in an indefinite-metric space, such an existence is not a priori guaranteed. The FW-transformed generators are given and satisfy the Poincare Lie algebra. We observe that although the FW transformation expressed as a power series in c/sup -1/ more clearly exhibits the physics of the situation, it is only in the Dirac and Duffin-Kemmer-Petiau special cases of the Bhabha fields that the FW-transformation power series can easily be summed to yieldmore » closed-form expression. The general closed-form expression surrenders more readily to another technique. Therefore, as a first calculation, in this paper we present a method of generating the FW transformation as a power series in c/sup -1/. Our discussion concentrates on the indefinite metric, the physics which is evident in the power-series form (such as size and types of Zitterbewegung), and on a detailed examination of special cases up to spin 3/2. In all the above, a special handling of the built-in subsidiary components of the integer-spin fields is once again necessary. We also comment on what the indefinite metric may be implying about the possibility of finding a totally consistent high-spin field theory. In a later paper we present a derivation of the exact, closed-form transformation. (AIP)« less
  • A local, second-order (truncated) approximation is applied to the Hubbard model in three dimensions. Lowering the temperature, at half-filling, the paramagnetic ground state becomes unstable towards the formation of a commensurate spin-density-wave (SDW) state (antiferromagnetism) and sufficiently far away from half-filling towards the formation of incommensurate SDW states. The incommensurate-ordering wave vector does not deviate much from the commensurate one, which is in accordance with the experimental data for the SDW in chromium alloys. {copyright} {ital 1997} {ital The American Physical Society}