Three parallel algorithms for classical molecular dynamics are presented. The first assigns each processor a fixed subset of atoms; the second assigns each a fixed subset of inter-atomic forces to compute; the third assigns each a fixed spatial region. The algorithms are suitable for molecular dynamics models which can be difficult to parallelize efficiently—those with short-range forces where the neighbors of each atom change rapidly. They can be implemented on any distributed-memory parallel machine which allows for message-passing of data between independently executing processors. The algorithms are tested on a standard Lennard-Jones benchmark problem for system sizes ranging from 500 to 100,000,000 atoms on several parallel supercomputers--the nCUBE 2, Intel iPSC/860 and Paragon, and Cray T3D. Comparing the results to the fastest reported vectorized Cray Y-MP and C90 algorithm shows that the current generation of parallel machines is competitive with conventional vector supercomputers even for small problems. For large problems, the spatial algorithm achieves parallel efficiencies of 90% and a 1840-node Intel Paragon performs up to 165 faster than a single Cray C9O processor. Trade-offs between the three algorithms and guidelines for adapting them to more complex molecular dynamics simulations are also discussed.

Fast parallel algorithms for short-range molecular dynamics

Mechanical properties of nanocrystalline materials

Constitutive relations are fundamental to the solution of problems in continuum mechanics, and are required in the study of, for example, mechanically dominated clinical interventions involving soft biological tissues. Structural continuum constitutive models of arterial layers integrate information about the tissue morphology and therefore allow investigation of the interrelation between structure and function in response to mechanical loading. Collagen fibres are key ingredients in the structure of arteries. In the media (the middle layer of the artery wall) they are arranged in two helically distributed families with a small pitch and very little dispersion in their orientation (i.e. they are aligned quite close to the circumferential direction). By contrast, in the adventitial and intimal layers, the orientation of the collagen fibres is dispersed, as shown by polarized light microscopy of stained arterial tissue. As a result, continuum models that do not account for the dispersion are not able to capture accurately the stress–strain response of these layers. The purpose of this paper, therefore, is to develop a structural continuum framework that is able to represent the dispersion of the collagen fibre orientation. This then allows the development of a new hyperelastic free-energy function that is particularly suited for representing the anisotropic elastic properties of adventitial and intimal layers of arterial walls, and is a generalization of the fibre-reinforced structural model introduced by Holzapfel & Gasser (Holzapfel & Gasser 2001 Comput. Meth. Appl. Mech. Eng. 190, 4379–4403) and Holzapfel et al. (Holzapfel et al. 2000 J. Elast. 61, 1–48). The model incorporates an additional scalar structure parameter that characterizes the dispersed collagen orientation. An efficient finite element implementation of the model is then presented and numerical examples show that the dispersion of the orientation of collagen fibres in the adventitia of human iliac arteries has a significant effect on their mechanical response.

https://royalsocietypublishing.org/doi/pdf/10.1098/rsif.2005.0073

Hyperelastic modelling of arterial layers with distributed collagen fibre orientations

1: Magnetic Particles: Preparation, Properties and Applications: M. Ozaki. 2: Maghemite (gamma-Fe2O3): A Versatile Magnetic Colloidal Material C.J. Serna, M.P. Morales. 3: Dynamics of Adsorption and Oxidation of Organic Molecules on Illuminated Titanium Dioxide Particles Immersed in Water M.A. Blesa, R.J. Candal, S.A. Bilmes. 4: Colloidal Aggregation in Two-Dimensions A. Moncho-Jorda, F. Martinez-Lopez, M.A. Cabrerizo-Vilchez, R. Hidalgo Alvarez, M. Quesada-PMerez. 5: Kinetics of Particle and Protein Adsorption Z. Adamczyk.

Surface and Colloid Science

Alloying has long been used to confer desirable properties to materials. Typically, it involves the addition of relatively small amounts of secondary elements to a primary element. For the past decade and a half, however, a new alloying strategy that involves the combination of multiple principal elements in high concentrations to create new materials called high-entropy alloys has been in vogue. The multi-dimensional compositional space that can be tackled with this approach is practically limitless, and only tiny regions have been investigated so far. Nevertheless, a few high-entropy alloys have already been shown to possess exceptional properties, exceeding those of conventional alloys, and other outstanding high-entropy alloys are likely to be discovered in the future. Here, we review recent progress in understanding the salient features of high-entropy alloys. Model alloys whose behaviour has been carefully investigated are highlighted and their fundamental properties and underlying elementary mechanisms discussed. We also address the vast compositional space that remains to be explored and outline fruitful ways to identify regions within this space where high-entropy alloys with potentially interesting properties may be lurking. High-entropy alloys have greatly expanded the compositional space for alloy design. In this Review, the authors discuss model high-entropy alloys with interesting properties, the physical mechanisms responsible for their behaviour and fruitful ways to probe and discover new materials in the vast compositional space that remains to be explored.

Vlado A. Lubarda

Papers

The onset of twinning in metals: a constitutive description

Void growth by dislocation emission

Damage tensors and the crack density distribution

On the mechanics of solids with a growing mass

On the effective lattice parameter of binary alloys