PACKMOL: A Package for Building Initial Configurations for Molecular Dynamics Simulations

Adequate initial configurations for molecular dynamics simulations consist of arrangements of molecules distributed in space in such a way to approximately represent the system`s overall structure. In order that the simulations are not disrupted by large van der Waals repulsive interactions, atoms from different molecules Must keep safe pairwise distances. Obtaining Such a molecular arrangement can be considered it packing problem: Each type molecule must satisfy spatial constraints related to the geometry of the system, and the distance between atoms of different molecules Must be greater than some specified tolerance. We have developed a code able to pack millions of atoms. grouped in arbitrarily complex molecules, inside a variety of three-dimensional regions. The regions may be intersections of spheres, ellipses, cylinders, planes, or boxes. The user must provide only the structure of one molecule of each type and the geometrical constraints that each type of molecule must satisfy. Building complex mixtures, interfaces, solvating biomolecules in water, other solvents, or mixtures of solvents, is straight forward. In addition. different atoms belonging to the same molecule may also be restricted to different spatial regions, in Such a way that more ordered molecular arrangements call be built, as micelles. lipid double-layers, etc. The packing time for state-of-the-art molecular dynamics systems varies front a few seconds to a few Minutes in a personal Computer. The input files are simple and Currently compatible with PDB, Tinker, Molden, or Moldy coordinate files. The package is distributed as free software and call be downloaded front http://www.ime.unicamp.br/similar to martinez/packmol/. (C) 2009 Wiley Periodicals. Inc. J Comput Chem 30: 2157-2164, 2009

Evolutionary dynamics on rugged fitness landscapes: Exact dynamics and information theoretical aspects

The parallel mutation-selection evolutionary dynamics, in which mutation and replication are independent events, is solved exactly in the case that the Malthusian fitnesses associated to the genomes are described by the random energy model (REM) and by a ferromagnetic version of the REM. The solution method uses the mapping of the evolutionary dynamics into a quantum Ising chain in a transverse field and the Suzuki-Trotter formalism to calculate the transition probabilities between configurations at different times. We find that in the case of the REM landscape the dynamics can exhibit three distinct regimes: pure diffusion or stasis for short times, depending on the fitness of the initial configuration, and a spin-glass regime for large times. The dynamic transition between these dynamical regimes is marked by discontinuities in the mean-fitness as well as in the overlap with the initial reference sequence. The relaxation to equilibrium is described by an inverse time decay. In the ferromagnetic REM, we find in addition to these three regimes, a ferromagnetic regime where the overlap and the mean-fitness are frozen. In this case, the system relaxes to equilibrium in a finite time. The relevance of our results to information processing aspects of evolution is discussed.

Multiclass Hammersley-Aldous-Diaconis process and multiclass-customer queues

In the Hammersley-Aldous-Diaconis process, infinitely many particles sit in R and at most one particle is allowed at each position. A particle at x, whose nearest neighbor to the right is at y, jumps at rate y - x to a position uniformly distributed in the interval (x, y). The basic coupling between trajectories with different initial configuration induces a process with different classes of particles. We show that the invariant measures for the two-class process can be obtained as follows. First, a stationary M/M/1 queue is constructed as a function of two homogeneous Poisson processes, the arrivals with rate, and the (attempted) services with rate rho > lambda Then put first class particles at the instants of departures (effective services) and second class particles at the instants of unused services. The procedure is generalized for the n-class case by using n - 1 queues in tandem with n - 1 priority types of customers. A multi-line process is introduced; it consists of a coupling (different from Liggett's basic coupling), having as invariant measure the product of Poisson processes. The definition of the multi-line process involves the dual points of the space-time Poisson process used in the graphical construction of the reversed process. The coupled process is a transformation of the multi-line process and its invariant measure is the transformation described above of the product measure.

The natural force density method for the shape finding of taut structures

This work presents, with the aid of the natural approach, an extension of the force density method for the initial shape finding of cable and membrane structures, which leads to the solution of a system of linear equations. This method, here called the natural force density method, preserves the linearity which characterizes the original force density method. At the same time, it overcomes the difficulties that the original procedure presents to cope with irregular triangular finite element meshes. Furthermore, if this method is applied iteratively in the lines prescribed herewith, it leads to a viable initial configuration with a uniform, isotropic plane Cauchy stress state. This means that a minimal surface for the membrane can be achieved through a succession of equilibrated configurations. Several numerical examples illustrate the simplicity and robustness of the method. (C) 2008 Elsevier B.V. All rights reserved.

Spatiotemporal solitons in multidimensional optical media with a quadratic nonlinearity

We consider solutions to the second-harmonic generation equations in two-and three-dimensional dispersive media in the form of solitons localized in space and time. As is known, collapse does not take place in these models, which is why the solitons may be stable. The general solution is obtained in an approximate analytical form by means of a variational approach, which also allows the stability of the solutions to be predicted. Then, we directly simulate the two-dimensional case, taking the initial configuration as suggested by the variational approximation. We thus demonstrate that spatiotemporal solitons indeed exist and are stable. Furthermore, they are not, in the general case, equivalent to the previously known cylindrical spatial solitons. Direct simulations generate solitons with some internal oscillations. However, these oscillations neither grow nor do they exhibit any significant radiative damping. Numerical solutions of the stationary version of the equations produce the same solitons in their unperturbed form, i.e., without internal oscillations. Strictly stable solitons exist only if the system has anomalous dispersion at both the fundamental harmonic and second harmonic (SH), including the case of zero dispersion at SH. Quasistationary solitons, decaying extremely slowly into radiation, are found in the presence of weak normal dispersion at the second-harmonic frequency.

An evolutionary approach for the motion planning of redundant and hyper-redundant manipulators

The trajectory planning of redundant robots is an important area of research and efficient optimization algorithms are needed. The pseudoinverse control is not repeatable, causing drift in joint space which is undesirable for physical control. This paper presents a new technique that combines the closed-loop pseudoinverse method with genetic algorithms, leading to an optimization criterion for repeatable control of redundant manipulators, and avoiding the joint angle drift problem. Computer simulations performed based on redundant and hyper-redundant planar manipulators show that, when the end-effector traces a closed path in the workspace, the robot returns to its initial configuration. The solution is repeatable for a workspace with and without obstacles in the sense that, after executing several cycles, the initial and final states of the manipulator are very close.

Attitude change in arbitrarily large organizations

The alignment of collective goals and individual behavior has been extensively studied by economists under a principal-agent framework. Two main solutions have been presented: explicit incentive contracts and monitoring. These solutions correspond to changes in the objective situation faced by individuals. However, an extensive literature in social psychology provides evidence that behavior is influenced, not only by situational constraints, but also by attitudes. Therefore, an important aspect of organization is to choose the structures and procedures that best contribute to the dissemination of the desired attitudes throughout the organization. This paper studies how the initial configuration of attitudes and the size of the organization affect the optimal organizational structure and the timing of information flows when the objective is to align the members' attitudes. We identify and characterize three factors that affect the optimal organizational structures and procedures and the degree of alignment of attitudes: (1) clustering effects; (2) member cross-influence effects; and (3) leader cross-influence effects.

Stable isotope composition of smectite in suevites at the Ries crater, Germany: Implications for hydrous alteration of impactites

The 24-km diameter Ries crater, Germany, exhibits well-preserved crater filling and surficial melt-rich breccia deposits that are believed to have been altered by post-impact hydrothermal fluids. The alteration mineralogy of the crater filling breccias is characterized by clay (smectite, chlorite) and a zeolite assemblage, and secondary clay phases (smectite, minor halloysite) in surficial melt-bearing breccia deposits. Using stable isotope analysis of secondary smectitic clay fractions, evidence of significant hydrous alteration of impactites at large water/rock ratios was found. The estimated fluid temperatures, using data derived by delta(18)O and delta D fractionation, suggest smectite precipitation in surficial breccias in equilibrium with meteoric fluids at temperatures 16 +/- 5 degrees C in agreement with the long-term variation of modern precipitation in the area. The stable isotope composition of smectite in crater-fill breccia, however, suggests a trend of monotonously increasing temperatures from 43 to 112 degrees C. with increasing depth through the breccia sequence. This demonstrates a different origin of alteration and temperature distribution for the surficial and crater filling melt-bearing impact breccias in the Ries crater. Our results suggest that the inverted structure of hydrothermal systems observed in some terrestrial impact craters, including the Ries crater, could indicate the initial configuration of a thermal anomaly in the crater filling sequence, but which is replaced with a normal hydrothermal convection in crater proper, during the course of post-impact cooling. (C) 2010 Elsevier B.V. All rights reserved.

RADIC: un middleware de tolerancia a fallos que preserva el rendimiento

La tolerancia a fallos es una línea de investigación que ha adquirido una importancia relevante con el aumento de la capacidad de cómputo de los súper-computadores actuales. Esto es debido a que con el aumento del poder de procesamiento viene un aumento en la cantidad de componentes que trae consigo una mayor cantidad de fallos. Las estrategias de tolerancia a fallos actuales en su mayoría son centralizadas y estas no escalan cuando se utiliza una gran cantidad de procesos, dado que se requiere sincronización entre todos ellos para realizar las tareas de tolerancia a fallos. Además la necesidad de mantener las prestaciones en programas paralelos es crucial, tanto en presencia como en ausencia de fallos. Teniendo en cuenta lo citado, este trabajo se ha centrado en una arquitectura tolerante a fallos descentralizada (RADIC – Redundant Array of Distributed and Independant Controllers) que busca mantener las prestaciones iniciales y garantizar la menor sobrecarga posible para reconfigurar el sistema en caso de fallos. La implementación de esta arquitectura se ha llevado a cabo en la librería de paso de mensajes denominada Open MPI, la misma es actualmente una de las más utilizadas en el mundo científico para la ejecución de programas paralelos que utilizan una plataforma de paso de mensajes. Las pruebas iniciales demuestran que el sistema introduce mínima sobrecarga para llevar a cabo las tareas correspondientes a la tolerancia a fallos. MPI es un estándar por defecto fail-stop, y en determinadas implementaciones que añaden cierto nivel de tolerancia, las estrategias más utilizadas son coordinadas. En RADIC cuando ocurre un fallo el proceso se recupera en otro nodo volviendo a un estado anterior que ha sido almacenado previamente mediante la utilización de checkpoints no coordinados y la relectura de mensajes desde el log de eventos. Durante la recuperación, las comunicaciones con el proceso en cuestión deben ser retrasadas y redirigidas hacia la nueva ubicación del proceso. Restaurar procesos en un lugar donde ya existen procesos sobrecarga la ejecución disminuyendo las prestaciones, por lo cual en este trabajo se propone la utilización de nodos spare para la recuperar en ellos a los procesos que fallan, evitando de esta forma la sobrecarga en nodos que ya tienen trabajo. En este trabajo se muestra un diseño propuesto para gestionar de un modo automático y descentralizado la recuperación en nodos spare en un entorno Open MPI y se presenta un análisis del impacto en las prestaciones que tiene este diseño. Resultados iniciales muestran una degradación significativa cuando a lo largo de la ejecución ocurren varios fallos y no se utilizan spares y sin embargo utilizándolos se restablece la configuración inicial y se mantienen las prestaciones.

Action publique concertée et gestion des sports de nature : une ethnographie participative au Conseil général du Bas-Rhin

En France, la décentralisation et la territorialisation de l'action publique ont fait des sports de nature un objet d'action publique légitime en donnant naissance à de nouveaux outils de management public dédiés à la concertation et à la planification des usages de la nature. Nés de l'article 52 de la Loi sur le sport modifiée en 2000, la Commission Départementale des Espaces, Sites et Itinéraires relatifs aux sports de nature (CDESI) et le Plan Départemental des Espaces Sites et Itinéraires relatifs aux sports de nature (PDESI) sont des outils de concertation territoriale dédiés à la gestion publique des sports de nature au niveau départemental. Un enjeu de ce travail tient à l'appréhension des transformations de l'action publique en s'attachant à l'étude des dispositifs de concertation sur les sports de nature. Un deuxième enjeu de ce travail s'attache à mettre en évidence les effets de la concertation en analysant les interactions et les différents modes d'engagements des acteurs au cours de la « chose publique en train de se faire » (Cefaï, 2002). Les acteurs s'engagent non seulement dans la concertation comprise comme une activité sociale faite d'interactions, mais ils s'engagent également dans la concertation en tant que processus d'action publique. Aussi, un autre enjeu de ce travail est d'appréhender les effets de la concertation par une analyse processuelle des engagements (Fillieule, 2004) des acteurs et des organisations. En mobilisant les outils conceptuels de la sociologie interactionniste, de la sociologie pragmatique, ainsi que de la sociologie structuraliste, l'analyse des situations interactionnelles a notamment permis d'identifier les procédures de cadrage et les techniques dramaturgiques mises en oeuvre par les interactants, ainsi que les répertoires argumentatifs mobilisés par ces acteurs pendant l « 'épreuve » de la concertation. Les confrontations des points de vue et les justifications des prises de positions des acteurs peuvent faire évoluer la configuration initiale des jeux d'acteurs même si, pour certains, ces changements ne restent parfois qu'éphémères. Les organisations s'engagent dans la concertation en fonction de la revendication d'une légitimité qui est à comprendre comme une forme militantisme institutionnel s'articulant autour de la valorisation d'une expertise militante, environnementale, institutionnelle, ou encore de leur statut de partenaire institutionnel. In France, decentralization and territorialization of public action have made outdoor sports become an object of public policies justifiable by giving birth to new tools of public management dedicated to the public consultation, the dialogue, and the planning of the uses of the landscapes. Indeed, born of article 52 of the Law on sport modified in 2000, the Departmental committee for Spaces, Sites and Routes relative to natural sports ( CDESI) and the Departmental Plan of Spaces Sites and Routes relative to natural sports ( PDESI) are governance tools dedicated to the public management of outdoor sports for counties. A challenge of this work is to understand the changes of public policy by focusing on the study of mechanisms for consultation on outdoor sports. A second item of this work is to highlight the effects of cooperation by focusing on the analysis of interactions and actor's commitments during the "public thing in the making" (Cefaï, 2002). Actors commit themselves not only in the dialogue included as a social activity made by interactions, but they also take part to the dialogue included as a process of public action. Also, another issue of this work is to understand the effects of consultation by a processual approach of individual commitments (Fillieule, 2004) of actors and organizations. Using the conceptual tools of symbolic interactionism, pragmatic sociology, and structuralist sociology, the analysis of interactional situations has highlighted the framing work and procedures implemented by the interactants, as well as the dramaturgical techniques and argumentative directories which, they mobilize during the "test" of the consultation. Confrontation of viewpoints and justifications of interactants' positions can evolve from their initial configuration sets, even if for some of them these changes are sometimes ephemeral. Organizations involve themselves according to demands of legitimacy which, are to understand as a shape institutional militancy articulating around the valuation of a militant, environmental, institutional expertise, or still around their status of institutional partner.

Integration of Bridge Damage Detection Concepts and Components (3 volumes), TR-636, 2013

This work is divided into three volumes: Volume I: Strain-Based Damage Detection; Volume II: Acceleration-Based Damage Detection; Volume III: Wireless Bridge Monitoring Hardware. Volume I: In this work, a previously-developed structural health monitoring (SHM) system was advanced toward a ready-for-implementation system. Improvements were made with respect to automated data reduction/analysis, data acquisition hardware, sensor types, and communication network architecture. The statistical damage-detection tool, control-chart-based damage-detection methodologies, were further investigated and advanced. For the validation of the damage-detection approaches, strain data were obtained from a sacrificial specimen attached to the previously-utilized US 30 Bridge over the South Skunk River (in Ames, Iowa), which had simulated damage,. To provide for an enhanced ability to detect changes in the behavior of the structural system, various control chart rules were evaluated. False indications and true indications were studied to compare the damage detection ability in regard to each methodology and each control chart rule. An autonomous software program called Bridge Engineering Center Assessment Software (BECAS) was developed to control all aspects of the damage detection processes. BECAS requires no user intervention after initial configuration and training. Volume II: In this work, a previously developed structural health monitoring (SHM) system was advanced toward a ready-for-implementation system. Improvements were made with respect to automated data reduction/analysis, data acquisition hardware, sensor types, and communication network architecture. The objective of this part of the project was to validate/integrate a vibration-based damage-detection algorithm with the strain-based methodology formulated by the Iowa State University Bridge Engineering Center. This report volume (Volume II) presents the use of vibration-based damage-detection approaches as local methods to quantify damage at critical areas in structures. Acceleration data were collected and analyzed to evaluate the relationships between sensors and with changes in environmental conditions. A sacrificial specimen was investigated to verify the damage-detection capabilities and this volume presents a transmissibility concept and damage-detection algorithm that show potential to sense local changes in the dynamic stiffness between points across a joint of a real structure. The validation and integration of the vibration-based and strain-based damage-detection methodologies will add significant value to Iowa’s current and future bridge maintenance, planning, and management Volume III: In this work, a previously developed structural health monitoring (SHM) system was advanced toward a ready-for-implementation system. Improvements were made with respect to automated data reduction/analysis, data acquisition hardware, sensor types, and communication network architecture. This report volume (Volume III) summarizes the energy harvesting techniques and prototype development for a bridge monitoring system that uses wireless sensors. The wireless sensor nodes are used to collect strain measurements at critical locations on a bridge. The bridge monitoring hardware system consists of a base station and multiple self-powered wireless sensor nodes. The base station is responsible for the synchronization of data sampling on all nodes and data aggregation. Each wireless sensor node include a sensing element, a processing and wireless communication module, and an energy harvesting module. The hardware prototype for a wireless bridge monitoring system was developed and tested on the US 30 Bridge over the South Skunk River in Ames, Iowa. The functions and performance of the developed system, including strain data, energy harvesting capacity, and wireless transmission quality, were studied and are covered in this volume.

Coisotropic regularization of singular Lagrangians

We present an alternative approach to the usual treatments of singular Lagrangians. It is based on a Hamiltonian regularization scheme inspired on the coisotropic embedding of presymplectic systems. A Lagrangian regularization of a singular Lagrangian is a regular Lagrangian defined on an extended velocity phase space that reproduces the original theory when restricted to the initial configuration space. A Lagrangian regularization does not always exists, but a family of singular Lagrangians is studied for which such a regularization can be described explicitly. These regularizations turn out to be essentially unique and provide an alternative setting to quantize the corresponding physical systems. These ideas can be applied both in classical mechanics and field theories. Several examples are discussed in detail. 1995 American Institute of Physics.

Merkel cells and the individuality of friction ridge skin

There is no definite theory yet for the mechanism by which the pattern of epidermal ridges on fingers, palms and soles forming friction ridge skin (FRS) patterns is created. For a long time growth forces in the embryonal epidermis have been believed to be involved in FRS formation. More recent evidence suggests that Merkel cells play an important part in this process as well. Here we suggest a model for the formation of FRS patterns that links Merkel cells to the epidermal stress distribution. The Merkel cells are modeled as agents in an agent based model that move anisotropically where the anisotropy is created by the epidermal stress tensor. As a result ridge patterns are created with pattern defects as they occur in real FRS patterns. As a consequence we suggest why the topology of FRS patterns is indeed unique as the arrangement of pattern defects is sensitive to the initial configuration of Merkel cells.

Influence of the feeding mechanism on deposits of square particles

In a previous paper [Hidalgo et al., Phys. Rev. Lett. 103, 118001 (2009)] it was shown that square particles deposited in a silo tend to align with a diagonal parallel to the gravity, giving rise to a deposit with very particular properties. Here we explore, both experimentally and numerically, the effect on these properties of the filling mechanism. In particular, we modify the volume fraction of the initial configuration from which the grains are deposited. Starting from a very dilute case, increasing the volume fraction results in an enhancement of the disorder in the final deposit characterized by a decrease of the final packing fraction and a reduction of the number of particles oriented with their diagonal in the direction of gravity. However, for very high initial volume fractions, the final packing fraction increases again. This result implies that two deposits with the same final packing fraction can be obtained from very different initial conditions. The structural properties of such deposits are analyzed, revealing that, although the final volume fraction is the same, their micromechanical properties notably differ.

Structural modelling of liquid ZrO2 and glassy B2O3 by using AB INITIO Molecular Dynamics method

Atomic structure of ZrO2 and B2O3 was investigated in this work. New data under extreme conditions (T = 3100 K) was obtained for the liquid ZrO2 structure. A fractional number of boron was investigated for glassy structure of B2O3. It was shown that it is possible to obtain an agreement for the fractional number between NMR and DFT techniques using a suitable initial configuration.