Population-based local search for protein folding simulation in the MJ energy model and cubic lattices

We present experimental results on benchmark problems in 3D cubic lattice structures with the Miyazawa-Jernigan energy function for two local search procedures that utilise the pull-move set: (i) population-based local search (PLS) that traverses the energy landscape with greedy steps towards (potential) local minima followed by upward steps up to a certain level of the objective function; (ii) simulated annealing with a logarithmic cooling schedule (LSA). The parameter settings for PLS are derived from short LSA-runs executed in pre-processing and the procedure utilises tabu lists generated for each member of the population. In terms of the total number of energy function evaluations both methods perform equally well, however. PLS has the potential of being parallelised with an expected speed-up in the region of the population size. Furthermore, both methods require a significant smaller number of function evaluations when compared to Monte Carlo simulations with kink-jump moves. (C) 2009 Elsevier Ltd. All rights reserved.

How Much Multiuser Diversity is Required for Energy Limited Multiuser Systems?

Multiuser diversity (MUDiv) is one of the central concepts in multiuser (MU) systems. In particular, MUDiv allows for scheduling among users in order to eliminate the negative effects of unfavorable channel fading conditions of some users on the system performance. Scheduling, however, consumes energy (e.g., for making users' channel state information available to the scheduler). This extra usage of energy, which could potentially be used for data transmission, can be very wasteful, especially if the number of users is large. In this paper, we answer the question of how much MUDiv is required for energy limited MU systems. Focusing on uplink MU wireless systems, we develop MU scheduling algorithms which aim at maximizing the MUDiv gain. Toward this end, we introduce a new realistic energy model which accounts for scheduling energy and describes the distribution of the total energy between scheduling and data transmission stages. Using the fact that such energy distribution can be controlled by varying the number of active users, we optimize this number by either i) minimizing the overall system bit error rate (BER) for a fixed total energy of all users in the system or ii) minimizing the total energy of all users for fixed BER requirements. We find that for a fixed number of available users, the achievable MUDiv gain can be improved by activating only a subset of users. Using asymptotic analysis and numerical simulations, we show that our approach benefits from MUDiv gains higher than that achievable by generic greedy access algorithm, which is the optimal scheduling method for energy unlimited systems. © 2010 IEEE.

A Significance-Driven Programming Framework for Energy-Constrained Approximate Computing

Approximate execution is a viable technique for energy-con\-strained environments, provided that applications have the mechanisms to produce outputs of the highest possible quality within the given energy budget.
We introduce a framework for energy-constrained execution with controlled and graceful quality loss. A simple programming model allows users to express the relative importance of computations for the quality of the end result, as well as minimum quality requirements. The significance-aware runtime system uses an application-specific analytical energy model to identify the degree of concurrency and approximation that maximizes quality while meeting user-specified energy constraints. Evaluation on a dual-socket 8-core server shows that the proposed
framework predicts the optimal configuration with high accuracy, enabling energy-constrained executions that result in significantly higher quality compared to loop perforation, a compiler approximation technique.

Exploiting Significance of Computations for Energy-Constrained Approximate Computing

Approximate execution is a viable technique for environments with energy constraints, provided that applications are given the mechanisms to produce outputs of the highest possible quality within the available energy budget. This paper introduces a framework for energy-constrained execution with controlled and graceful quality loss. A simple programming model allows developers to structure the computation in different tasks, and to express the relative importance of these tasks for the quality of the end result. For non-significant tasks, the developer can also supply less costly, approximate versions. The target energy consumption for a given execution is specified when the application is launched. A significance-aware runtime system employs an application-specific analytical energy model to decide how many cores to use for the execution, the operating frequency for these cores, as well as the degree of task approximation, so as to maximize the quality of the output while meeting the user-specified energy constraints. Evaluation on a dual-socket 16-core Intel platform using 9 benchmark kernels shows that the proposed framework picks the optimal configuration with high accuracy. Also, a comparison with loop perforation (a well-known compile-time approximation technique), shows that the proposed framework results in significantly higher quality for the same energy budget.

Density functional and Monte Carlo studies of sulfur. II. Equilibrium polymerization of the liquid phase

The equilibrium polymerization of sulfur is investigated by Monte Carlo simulations. The potential energy model is based on density functional results for the cohesive energy, structural, and vibrational properties as well as reactivity of sulfur rings and chains [Part I, J. Chem. Phys. 118, 9257 (2003)]. Liquid samples of 2048 atoms are simulated at temperatures 450less than or equal toTless than or equal to850 K and P=0 starting from monodisperse S-8 molecular compositions. Thermally activated bond breaking processes lead to an equilibrium population of unsaturated atoms that can change the local pattern of covalent bonds and allow the system to approach equilibrium. The concentration of unsaturated atoms and the kinetics of bond interchanges is determined by the energy DeltaE(b) required to break a covalent bond. Equilibrium with respect to the bond distribution is achieved for 15less than or equal toDeltaE(b)less than or equal to21 kcal/mol over a wide temperature range (Tgreater than or equal to450 K), within which polymerization occurs readily, with entropy from the bond distribution overcompensating the increase in enthalpy. There is a maximum in the polymerized fraction at temperature T-max that depends on DeltaE(b). This fraction decreases at higher temperature because broken bonds and short chains proliferate and, for Tless than or equal toT(max), because entropy is less important than enthalpy. The molecular size distribution is described well by a Zimm-Schulz function, plus an isolated peak for S-8. Large molecules are almost exclusively open chains. Rings tend to have fewer than 24 atoms, and only S-8 is present in significant concentrations at all T. The T dependence of the density and the dependence of polymerization fraction and degree on DeltaE(b) give estimates of the polymerization temperature T-f=450+/-20 K. (C) 2003 American Institute of Physics.

Shell Model Study of Local and Global Energy Barriers in KDP

We perform a study of the energetics of KH2PO4 (KDP) by using a shell model (SM) which was constructed by adjusting the interaction parameters to ab initio calculations, and was fitted to reproduce phonons, polarization-inversion energies and structural properties. We calculate the energy profiles by performing global displacements and local distortions following the ferroelectric (FE) mode pattern in clusters of different sizes embedded in a paraelectric (PE) phase matrix. These properties are expected to be relevant to the PE-FE phase transition. The obtained SM results are compared to corresponding ab initio (AI) data. The global instabilities are found in good agreement for both KDP and DKDP. We also find qualitative good agreement in the KDP structure and even quantitative agreement in the expanded DKDP structure for the local distortions. The SM results reproduce well different trends like increasing instabilities as the cluster sizes grows, as the heavier atoms are included, and as the volume is increased, in accordance with the corresponding data from AI calculations.

Random neighbor sandpile model with constrained total energy

We propose as energy-constrained sandpile model with random neighbors. The critical behavior of the model is in the same universality class as the mean-field self-organized criticality sandpile. The critical energy E-c depends on the number of neighbors n of each site, but the various exponents do not. For n = 6, we got that E-c = 0.4545; and a self-similar structure of the energy distribution function with five major peaks is also observed. This is a natural result of system dynamics and the way the system is disturbed.

Universality and self-similarity of an energy-constrained sandpile model with random neighbors

We study an energy-constrained sandpile model with random neighbors. The critical behavior of the model is in the same universality class as the mean-field self-organized criticality sandpile. The critical energy E-c depends on the number of neighbors n for each site, but the various exponents are independent of n. A self-similar structure with n-1 major peaks is developed for the energy distribution p(E) when the system approaches its stationary state. The avalanche dynamics contributes to the major peaks appearing at E-Pk = 2k/(2n - 1) with k = 1,2,...,n-1, while the fine self-similar structure is a natural result of the way the system is disturbed. [S1063-651X(99)10307-6].

An application of bi-directional evolutionary structural optimisation for optimising energy absorbing structures using a material damage model

The Bi-directional Evolutionary Structural Optimisation (BESO) method is a numerical topology optimisation method developed for use in finite element analysis. This paper presents a particular application of the BESO method to optimise the energy absorbing capability of metallic structures. The optimisation objective is to evolve a structural geometry of minimum mass while ensuring that the kinetic energy of an impacting projectile is reduced to a level which prevents perforation. Individual elements in a finite element mesh are deleted when a prescribed damage criterion is exceeded. An energy absorbing structure subjected to projectile impact will fail once the level of damage results in a critical perforation size. It is therefore necessary to constrain an optimisation algorithm from producing such candidate solutions. An algorithm to detect perforation was implemented within a BESO framework which incorporated a ductile material damage model.

A Programming Model and Runtime System for Significance-Aware Energy-Efficient Computing.

We introduce a task-based programming model and runtime system that exploit the observation that not all parts of a program are equally significant for the accuracy of the end-result, in order to trade off the quality of program outputs for increased energy-efficiency. This is done in a structured and flexible way, allowing for easy exploitation of different points in the quality/energy space, without adversely affecting application performance. The runtime system can apply a number of different policies to decide whether it will execute less-significant tasks accurately or approximately.

The experimental evaluation indicates that our system can achieve an energy reduction of up to 83% compared with a fully accurate execution and up to 35% compared with an approximate version employing loop perforation. At the same time, our approach always results in graceful quality degradation.

Work Extraction and Energy Storage in the Dicke Model

We study work extraction from the Dicke model achieved using simple unitary cyclic transformations keeping into account both a non optimal unitary protocol, and the energetic cost of creating the initial state. By analyzing the role of entanglement, we find that highly entangled states can be inefficient for energy storage when considering the energetic cost of creating the state. Such surprising result holds notwithstanding the fact that the criticality of the model at hand can sensibly improve the extraction of work. While showing the advantages of using a many-body system for work extraction, our results demonstrate that entanglement is not necessarily advantageous for energy storage purposes, when non optimal processes are considered. Our work shows the importance of better understanding the complex interconnections between non-equilibrium thermodynamics of quantum systems and correlations among their subparts.

A Phase-field Model for Computer Simulation of Lamellar Structure Formation in y-TiAI

In the present paper, a phase-field model is developed to simulate the formation and evolution of lamellar microstructure in γ-TiAl alloys. The mechanism of formation of TiAl lamellae proposed by Denquin and Naka is incorporated into the model. The model describes the formation and evolution of the face-centered cubic (fcc) stacking lamellar zone followed by the subsequent appearance and growth of the γ-phase, involving both the chemical composition change by atom transfer and the ordering of the fcc lattice. The thermodynamics of the model system and the interaction between the displacive and diffusional transformations are described by a non-equilibrium free energy formulated as a function of concentration and structural order parameter fields. The long-range elastic interactions, arising from the lattice misfit between the α, fcc (A1) and the various orientation variants of the γ-phase are taken into account by incorporating of the elastic strain energy into the total free energy. Simulation studies based on the model successfully predicted some essential features of the lamellar structure. It is found that the formation and evolution of the lamellar structure are predominantly controlled by the minimization of the elastic energy of the interfaces between the different fcc stacking groups, low-symmetry product phase γ and the high-symmetry α-phase, as well as between the various orientation variants of the product phase.