skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Efficient optimization of natural resonance theory weightings and bond orders by gram‐based convex programming
We describe the formal algorithm and numerical applications of a novel convex quadratic programming (QP) strategy for performing the variational minimization that underlies natural resonance theory (NRT). The QP algorithm vastly improves the numerical efficiency, thoroughness, and accuracy of variational NRT description, which now allows uniform treatment ofallreference structures at the high level of detail previously reserved only for leading “reference” structures, with little or no user guidance. We illustrate overall QPNRT search strategy, program I/O, and numerical results for a specific application to adenine, and we summarize more extended results for a data set of 338 species from throughout the organic, bioorganic, and inorganic domain. The improved QP‐based implementation of NRT is a principal feature of the newly releasedNBO 7.0program version. © 2019 Wiley Periodicals, Inc.  more » « less
Award ID(s):
1740707
PAR ID:
10460459
Author(s) / Creator(s):
 ;  ;  
Publisher / Repository:
Wiley Blackwell (John Wiley & Sons)
Date Published:
Journal Name:
Journal of Computational Chemistry
Volume:
40
Issue:
23
ISSN:
0192-8651
Page Range / eLocation ID:
p. 2028-2035
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract The interplay of correlated spatial modulation and symmetry breaking leads to quantum critical phenomena intermediate between those of the clean and randomly disordered cases. By performing a detailed analytic and numerical case study of the quasi-periodically (QP) modulated transverse field Ising chain, we provide evidence for the conjectures of reference (Crowleyet al2018Phys. Rev. Lett.120175702) regarding the QP-Ising universality class. In the generic case, we confirm that the logarithmic wandering coefficientwgoverns both the macroscopic critical exponents and the energy-dependent localisation length of the critical excitations. However, for special values of the phase difference Δ between the exchange and transverse field couplings, the QP-Ising transition has different properties. For Δ = 0, a generalised Aubry–André duality prevents the finite energy excitations from localising despite the presence of logarithmic wandering. For Δ such that the fields and couplings are related by a lattice shift, the wandering coefficientwvanishes. Nonetheless, the presence of small couplings leads to non-trivial exponents and localised excitations. Our results add to the rich menagerie of quantum Ising transitions in the presence of spatial modulation. 
    more » « less
  2. Intermediate data structures are a common cause of inefficiency in functional programming.Fusionattempts to eliminate intermediate data structures by combining adjacent data traversals into one; existing fusion techniques, however, are based on predefined rewrite rules and hence are limited in expressiveness. In this work we explore a different approach to eliminating intermediate data structures, based on inductive program synthesis. We dub this approachsuperfusion(by analogy withsuperoptimization, which uses inductive synthesis for program optimization). Starting from a reference program annotated with data structures to be eliminated, superfusion first generates asketchwhere program fragments operating on those data structures are replaced with holes; it then fills the holes with constant-time expressions such that the resulting program is equivalent to the reference. The main technical challenge here is scalability because optimized programs are often complex, making the search space intractably large for naive enumeration. To address this challenge, our key insight is to first synthesize aghost functionthat describes the relationship between the original intermediate data structure and its compressed version; this function, although not used in the final program, serves to decompose the joint sketch filling problem into independent simpler problems for each hole. We implement superfusion in a tool calledSuFuand evaluate it on a dataset of 290 tasks collected from prior work on deductive fusion and program restructuring. The results show that SuFu solves 264 out of 290 tasks, exceeding the capabilities of rewriting-based fusion systems and achieving comparable performance with specialized approaches to program restructuring on their respective domains. 
    more » « less
  3. We study one-dimensional hybrid quantum circuits perturbed by quenched quasiperiodic (QP) modulations across the measurement-induced phase transition (MIPT). Considering non-Pisot QP structures, characterized by unbounded fluctuations, allows us to tune the wandering exponent β to exceed the Luck bound ν ≥ 1/(1−β) for the stability of the MIPT, where ν = 1.28(2). Via robust numerical simulations of random Clifford circuits interleaved with local projective measurements, we find that sufficiently large QP structural fluctuations destabilize the MIPT and induce a flow to a broad family of critical dynamical phase transitions of the infinite QP type that is governed by the wandering exponent β. We numerically determine the associated critical properties, including the correlation length exponent consistent with saturating the Luck bound, and a universal activated dynamical scaling with activation exponent ψ ≅ β, finding excellent agreement with the conclusions of real-space renormalization group calculations. 
    more » « less
  4. Abstract Observations of 30‐MHz coherent backscatter from sporadic‐Eionization layers were obtained with a VHF imaging radar located in Ithaca, New York. The volume probed by the radar lies at relatively high magnetic latitudes, on the northern edge of the mid‐latitude region and underneath the ionospheric trough. Banded, quasi‐periodic (QP) echoes observed from Ithaca are similar to those found in lower midlatitude regions. The Doppler shifts observed are smaller and, so far, do not appear to reach the threshold for Farley‐Buneman instability. However, many of the echoes exhibit fine‐scale structure, with secondary bands or braids oriented obliquely to the primary bands. Secondary bands have been seen only rarely at lower middle latitudes. In previous observations, the QP scattering has been linked to unstable neutral wind shears. Neutral wind shear commonly found in the lower thermosphere could play a key role in the formation of these irregularities and explain some morphological features of the resulting plasma density irregularities and the radar echoes. We consider whether neutral instability and turbulence in the lower thermosphere is the likely cause for some of the structuring in the sporadic‐Elayers. Results of 3D numerical simulations of atmospheric dynamics in the mesosphere to lower thermosphere support the proposition. In particular, we focus on Ekman‐type instabilities that, like the more common Kelvin‐Helmholtz instabilities, are inflection point instabilities, although specifically associated with turning shears, and result in convective rolls aligned close to the mean wind direction, with smaller‐scale secondary waves aligned normal to the primary structures. 
    more » « less
  5. We model a three-link fully actuated biped as a hybrid system and propose a prediction-based control algorithm for global tracking of reference trajectories. The proposed control strategy consists of a reference system that generates the desired periodic gait, a virtual system that generates a suitable reference trajectory using prediction, and a tracking control law that steers the biped to the virtual trajectory. The proposed algorithms achieves, in finite time, tracking in two steps. We present mathematical properties that define the main elements in the hybrid predictive controller for achieving convergence to the reference within the first two steps. The results are validated through numerical simulations. 
    more » « less