Search NASA⌕ Search

Engineering topics

Engel, Michael

Publications and source records attributed to Engel, Michael.

23 records · Page 2

Computational self-assembly of a one-component icosahedral quasicrystal

Icosahedral quasicrystals (IQCs) are a form of matter that is ordered but not periodic in any direction. All reported IQCs are intermetallic compounds and either of face-centred-icosahedral or primitive-icosahedral type, and the positions of their atoms have been resolved from diffraction data. However, unlike axially symmetric quasicrystals, IQCs have not been observed in non-atomic (that is, micellar or nanoparticle) systems, where real-space information would be directly available. Here, we show that an IQC can be assembled by means of molecular dynamics simulations from a one-component system of particles interacting via a tunable, isotropic pair potential extending only to the third-neighbour shell. The IQC is body-centred, self-assembles from a fluid phase, and in parameter space neighbours clathrates and other tetrahedrally bonded crystals. Our findings elucidate the structure and dynamics of the IQC, and suggest routes to search for it and design it in soft matter and nanoscale systems.

36 MATERIALS SCIENCE↗

Understanding shape entropy through local dense packing

Entropy drives the phase behavior of colloids ranging from dense suspensions of hard spheres or rods to dilute suspensions of hard spheres and depletants. Entropic ordering of anisotropic shapes into complex crystals, liquid crystals, and even quasicrystals was demonstrated recently in computer simulations and experiments. The ordering of shapes appears to arise from the emergence of directional entropic forces (DEFs) that align neighboring particles, but these forces have been neither rigorously defined nor quantified in generic systems. In this paper, we show quantitatively that shape drives the phase behavior of systems of anisotropic particles upon crowding through DEFs. We define DEFs in generic systems and compute them for several hard particle systems. We show they are on the order of a few times the thermal energy (k B T) at the onset of ordering, placing DEFs on par with traditional depletion, van der Waals, and other intrinsic interactions. In experimental systems with these other interactions, we provide direct quantitative evidence that entropic effects of shape also contribute to self-assembly. We use DEFs to draw a distinction between self-assembly and packing behavior. We show that the mechanism that generates directional entropic forces is the maximization of entropy by optimizing local particle packing. Finally, we show that this mechanism occurs in a wide class of systems and we treat, in a unified way, the entropy-driven phase behavior of arbitrary shapes, incorporating the well-known works of Kirkwood, Onsager, and Asakura and Oosawa.

36 MATERIALS SCIENCE↗

Massively parallel Monte Carlo for many-particle simulations on GPUs

Current trends in parallel processors call for the design of efficient massively parallel algorithms for scientific computing. Parallel algorithms for Monte Carlo simulations of thermodynamic ensembles of particles have received little attention because of the inherent serial nature of the statistical sampling. In this paper, we present a massively parallel method that obeys detailed balance and implement it for a system of hard disks on the GPU. We reproduce results of serial high-precision Monte Carlo runs to verify the method. This is a good test case because the hard disk equation of state over the range where the liquid transforms into the solid is particularly sensitive to small deviations away from the balance conditions. On a Tesla K20, our GPU implementation executes over one billion trial moves per second, which is 148 times faster than on a single Intel Xeon E5540 CPU core, enables 27 times better performance per dollar, and cuts energy usage by a factor of 13. With this improved performance we are able to calculate the equation of state for systems of up to one million hard disks. These large system sizes are required in order to probe the nature of the melting transition, which has been debated for the last forty years. In this paper we present the details of our computational method, and discuss the thermodynamics of hard disks separately in a companion paper.

97 MATHEMATICS AND COMPUTING↗