Source author record

Michael Engel

Michael Engel appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

32works
15topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

32 published item(s)

preprint2022arXiv

A statistical shape model for radiation-free assessment and classification of craniosynostosis

The assessment of craniofacial deformities requires patient data which is sparsely available. Statistical shape models provide realistic and synthetic data enabling comparisons of existing methods on a common dataset. We build the first publicly available statistical 3D head model of craniosynostosis patients and the first model focusing on infants younger than 1.5 years. We further present a shape-model-based classification pipeline to distinguish between three different classes of craniosynostosis and a control group on photogrammetric surface scans. To the best of our knowledge, our study uses the largest dataset of craniosynostosis patients in a classification study for craniosynostosis and statistical shape modeling to date. We demonstrate that our shape model performs similar to other statistical shape models of the human head. Craniosynostosis-specific pathologies are represented in the first eigenmodes of the model. Regarding the automatic classification of craniosynostis, our classification approach yields an accuracy of 97.8%, comparable to other state-of-the-art methods using both computed tomography scans and stereophotogrammetry. Our publicly available, craniosynostosis-specific statistical shape model enables the assessment of craniosynostosis on realistic and synthetic data. We further present a state-of-the-art shape-model-based classification approach for a radiation-free diagnosis of craniosynostosis.

preprint2020arXiv

Effects of additives on oil displacement in nanocapillaries: a mesoscale simulation study

We performed mesoscale simulations in order to investigate the effects that additives, such as surfactants and polymers, have in the oil displacement process by water and brine injection. A Many-Body Dissipative Particle Dynamics (MDPD) model was parameterized in order to reproduce physical properties obtained either by experiments or Classical Molecular Dynamics (MD). The MDPD model was then employed to simulate the displacement of n-dodecane by water and brine, with surfactant concentrations ranging from 1 wt. % to 10 wt. % SDS and polymer concentration of 0.5 wt. % HPAM. We observed that while the additives may enter the capillary without clogging, specific combinations of surfactant concentrations and injection rates may lead to poor oil displacement. This result can be understood in terms of the balance between water-wall and (surface-mediated) water-oil interactions. As a consequence, our results show an effective reversal of the wettability character, from water-wet to oil-wet, as a function of the interface speed. We conclude by summarizing the implications such results have on the fluid-fluid displacements involving additives.

preprint2020arXiv

Spontaneous Crystallization in Systems of Binary Hard Sphere Colloids

Computer simulations of the fluid-to-solid phase transition in the hard sphere system were instrumental for our understanding of crystallization processes. But while colloid experiments and theory have been predicting the stability of several binary hard sphere crystals for many years, simulations were not successful to confirm this phenomenon. Here, we report the growth of binary hard sphere crystals isostructural to Laves phases, AlB$_2$, and NaZn$_{13}$ in simulation directly from the fluid. We analyze particle kinetics during Laves phase growth using event-driven molecular dynamics simulations with and without swap moves that speed up diffusion. The crystallization process transitions from nucleation and growth to spinodal decomposition already deep within the fluid-solid coexistence regime. Finally, we present packing fraction-size ratio state diagrams in the vicinity of the stability regions of three binary crystals.

preprint2016arXiv

The Role of Short-Range Order and Hyperuniformity in the Formation of Band Gaps in Disordered Photonic Materials

We study photonic band gap formation in two-dimensional high refractive index disordered materials where the dielectric structure is derived from packing disks in real and reciprocal space. Numerical calculations of the photonic density of states demonstrate the presence of a band gap for all polarizations in both cases. We find that the band gap width is controlled by the increase in positional correlation inducing short-range order and hyperuniformity concurrently. Our findings suggest that the optimization of short-range order, in particular the tailoring of Bragg scattering at the isotropic Brillouin zone, are of key importance for designing disordered PBG materials.

preprint2015arXiv

Controlling Chirality of Entropic Crystals

Colloidal crystal structures with complexity and diversity rivaling atomic and molecular crystals have been predicted and obtained for hard particles by entropy maximization. However, so far homochiral colloidal crystals, which are candidates for photonic metamaterials, are absent. Using Monte Carlo simulations we show that chiral polyhedra exhibiting weak directional entropic forces self-assemble either an achiral crystal or a chiral crystal with limited control over the crystal handedness. Building blocks with stronger faceting exhibit higher selectivity and assemble a chiral crystal with handedness uniquely determined by the particle chirality. Tuning the strength of directional entropic forces by means of particle rounding or the use of depletants allows for reconfiguration between achiral and homochiral crystals. We rationalize our findings by quantifying the chirality strength of each particle, both from particle geometry and potential of mean force and torque diagrams.

preprint2015arXiv

Fast radio bursts: the observational case for a Galactic origin

There are by now ten published detections of fast radio bursts (FRBs), single bright GHz-band millisecond pulses of unknown origin. Proposed explanations cover a broad range from exotic processes at cosmological distances to atmospheric and terrestrial sources. Loeb et al. have previously suggested that FRB sources could be nearby flare stars, and pointed out the presence of a W-UMa-type contact binary within the beam of one out of three FRB fields that they examined. Using time-domain optical photometry and spectroscopy, we now find possible flare stars in additional FRB fields, with one to three such cases among eight FRB fields studied. We evaluate the chance probabilities of these possible associations to be in the range 0.1% to 9%, depending on the input assumptions. Further, we re-analyze the probability that two FRBs recently discovered 3 years apart within the same radio beam are unrelated. Contrary to other claims, we conclude with 99% confidence that the two events are from the same repeating source. The different dispersion measures between the two bursts then rule out a cosmological origin for the dispersion measure, but are consistent with the flare-star scenario with a varying plasma blanket between bursts. Finally, we review some theoretical objections that have been raised against a local flare-star FRB origin, and show that they are incorrect.

preprint2015arXiv

Power dissipation and electrical breakdown in black phosphorus

We report operating temperatures and heating coefficients measured in a multi-layer black phosphorus device as a function of injected electrical power. By combining micro-Raman spectroscopy and electrical transport measurements, we have observed a linear temperature increase up to 600K at a power dissipation rate of 0.896Kμm^3/mW. By further increasing the bias voltage, we determined the threshold power and temperature for electrical breakdown and analyzed the fracture in the black phosphorus layer that caused the device failure by means of scanning electron microscopy and atomic force microscopy. The results will benefit the research and development of electronics and optoelectronics based on novel two-dimensional materials.

preprint2014arXiv

A black phosphorus photo-detector for multispectral, high-resolution imaging

Black phosphorus is a layered semiconductor that is intensely researched in view of applications in optoelectronics. In this Letter, we investigate a multi-layer black phosphorus photo-detector that is capable of acquiring high-contrast (V>0.9) images both in the visible (λ_{VIS}=532nm) as well as in the infrared (λ_{IR}=1550nm) spectral regime. In a first step, by using photocurrent microscopy, we map the active area of the device and we characterize responsivity and gain. In a second step, by deploying the black phosphorus device as a point-like detector in a confocal microsope setup, we acquire diffraction-limited optical images with sub-micron resolution. The results demonstrate the usefulness of black phosphorus as an optoelectronic material for hyperspectral imaging applications.

preprint2014arXiv

Controlled self-assembly of periodic and aperiodic cluster crystals

Soft particles are known to overlap and form stable clusters that self-assemble into periodic crystalline phases with density-independent lattice constants. We use molecular dynamics simulations in two dimensions to demonstrate that, through a judicious design of an isotropic pair potential, one can control the ordering of the clusters and generate a variety of phases, including decagonal and dodecagonal quasicrystals. Our results confirm analytical predictions based on a mean-field approximation, providing insight into the stabilization of quasicrystals in soft macromolecular systems, and suggesting a practical approach for their controlled self-assembly in laboratory realizations using synthesized soft-matter particles.

preprint2014arXiv

Entropically Patchy Particles: Engineering Valence through Shape Entropy

Patchy particles are a popular paradigm for the design and synthesis of nanoparticles and colloids for self-assembly. In "traditional" patchy particles, anisotropic interactions arising from patterned coatings, functionalized molecules, DNA, and other enthalpic means create the possibility for directional binding of particles into higher-ordered structures. Although the anisotropic geometry of non-spherical particles contributes to the interaction patchiness through van der Waals, electrostatic, and other interactions, how particle shape contributes entropically to self-assembly is only now beginning to be understood. It has been recently demonstrated that, for hard shapes, entropic forces are directional. A newly proposed theoretical framework that defines and quantifies directional entropic forces demonstrates the anisotropic--that is, patchy--nature of these emergent, attractive forces. Here we introduce the notion of entropically patchy particles as the entropic counterpart to enthalpically patchy particles. Using three example "families" of shapes, we judiciously modify entropic patchiness by introducing geometric features to the particles so as to target specific crystal structures, which then assembled with Monte Carlo simulations. We quantify the emergent entropic valence via a potential of mean force and torque. We generalize these shape operations to shape anisotropy dimensions, in analogy with the anisotropy dimensions introduced for enthalpically patchy particles. Our findings demonstrate that entropic patchiness and emergent valence provide a way of engineering directional bonding into nanoparticle systems, whether in the presence or absence of additional, non-entropic forces.

preprint2014arXiv

Origin of photoresponse in black phosphorus photo-transistors

We study the origin of photocurrent generated in doped multilayer BP photo-transistors, and find that it is dominated by thermally driven thermoelectric and bolometric processes. The experimentally observed photocurrent polarities are consistent with photo-thermal processes. The photo-thermoelectric current can be generated up to a $μ$m away from the contacts, indicating a long thermal decay length. With an applied source-drain bias, a photo-bolometric current is generated across the whole device, overwhelming the photo-thermoelectric contribution at a moderate bias. The photo-responsivity in the multilayer BP device is two orders of magnitude larger than that observed in graphene.

preprint2014arXiv

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 has been 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. Here, 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 that they are on the order of a few kT 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. 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.

preprint2013arXiv

Complexity in surfaces of densest packings for families of polyhedra

Packings of hard polyhedra have been studied for centuries due to their mathematical aesthetic and more recently for their applications in fields such as nanoscience, granular and colloidal matter, and biology. In all these fields, particle shape is important for structure and properties, especially upon crowding. Here, we explore packing as a function of shape. By combining simulations and analytic calculations, we study three 2-parameter families of hard polyhedra and report an extensive and systematic analysis of the densest packings of more than 55,000 convex shapes. The three families have the symmetries of triangle groups (icosahedral, octahedral, tetrahedral) and interpolate between various symmetric solids (Platonic, Archimedean, Catalan). We find that optimal (maximum) packing density surfaces that reveal unexpected richness and complexity, containing as many as 130 different structures within a single family. Our results demonstrate the utility of thinking of shape not as a static property of an object in the context of packings, but rather as but one point in a higher dimensional shape space whose neighbors in that space may have identical or markedly different packings. Finally, we present and interpret our packing results in a consistent and generally applicable way by proposing a method to distinguish regions of packings and classify types of transitions between them.

preprint2013arXiv

Emergent collective phenomena in a mixture of hard shapes through active rotation

We investigate collective phenomena with rotationally driven spinners of concave shape. Each spinner experiences a constant internal torque in either a clockwise or counterclockwise direction. Although the spinners are modeled as hard, otherwise non-interacting rigid bodies, we find that their active motion induces an effective interaction that favors rotation in the same direction. With increasing density and activity, phase separation occurs via spinodal decomposition, as well as self-organization into rotating crystals. We observe the emergence of cooperative, super-diffusive motion along interfaces, which can transport inactive test particles. Our results demonstrate novel phase behavior of actively rotated particles that is not possible with linear propulsion or in non-driven, equilibrium systems of identical hard particles.

preprint2013arXiv

Hard-disk equation of state: First-order liquid-hexatic transition in two dimensions with three simulation methods

We report large-scale computer simulations of the hard-disk system at high densities in the region of the melting transition. Our simulations reproduce the equation of state, previously obtained using the event-chain Monte Carlo algorithm, with a massively parallel implementation of the local Monte Carlo method and with event-driven molecular dynamics. We analyze the relative performance of these simulation methods to sample configuration space and approach equilibrium. Our results confirm the first-order nature of the melting phase transition in hard disks. Phase coexistence is visualized for individual configurations via the orientational order parameter field. The analysis of positional order confirms the existence of the hexatic phase.

preprint2013arXiv

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.

preprint2013arXiv

Packing and Self-assembly of Truncated Triangular Bipyramids

Motivated by breakthroughs in the synthesis of faceted nano- and colloidal particles, as well as theoretical and computational studies of their packings, we investigate a family of truncated triangular bipyramids. We report dense periodic packings with small unit cells that were obtained via numerical and analytical optimization. The maximal packing fraction $ϕ_{\max}$ changes continuously with the truncation parameter $t$. Eight distinct packings are identified based on discontinuities in the first and second derivatives of $ϕ_{\max}(t)$. These packings differ in the number of particles in the fundamental domain (unit cell) and the type of contacts between the particles. In particular, we report two packings with four particles in the unit cell for which both $ϕ_{\max}(t)$ and $ϕ'_{\max}(t)$ are continuous and the discontinuity occurs in the second derivative only. In the self-assembly simulations that we perform for larger boxes with 2048 particles, only one out of eight packings is found to assemble. In addition, the degenerate quasicrystal reported previously for triangular bipyramids without truncation [{http://prl.aps.org/abstract/PRL/v107/i21/e215702}{Haji-Akbari \emph{et al.}, \emph{Phys. Rev. Lett.} \textbf{107}: 215702 (2011)}] assembles for truncations as high as 0.45. The self-assembly propensities for the structures formed in the thermodynamic limit are explained using the isoperimetric quotient of the particles and the coordination number in the disordered fluid and in the assembled structure.

preprint2012arXiv

Confirmation of the random tiling hypothesis for a decagonal quasicrystal

Mechanisms that stabilize quasicrystals are much discussed but not finally resolved. We confirm the random tiling hypothesis and its predictions in a fully atomistic decagonal quasicrystal model by calculating the free energy and the phason elastic constants over a wide range of temperatures. The Frenkel-Ladd method is applied for the phonon part and an approach of uncorrelated phason flips for the configurational part. When lowering the temperature, a phase transition to an approximant occurs. Close to the transition temperature one of the phason elastic constants becomes soft.

preprint2012arXiv

High-frequency performance of scaled carbon nanotube array field-effect transistors

We report the radio-frequency performance of carbon nanotube array transistors that have been realized through the aligned assembly of highly separated, semiconducting carbon nanotubes on a fully scalable device platform. At a gate length of 100 nm, we observe output current saturation and obtain as-measured, extrinsic current gain and power gain cut-off frequencies, respectively, of 7 GHz and 15 GHz. While the extrinsic current gain is comparable to the state-of-the-art the extrinsic power gain is improved. The de-embedded, intrinsic current gain and power gain cut-off frequencies of 153 GHz and 30 GHz are the highest values experimentally achieved to date. We analyze the consistency of DC and AC performance parameters and discuss the requirements for future applications of carbon nanotube array transistors in high-frequency electronics.

preprint2012arXiv

Light-matter interaction in a microcavity-controlled graphene transistor

Graphene has extraordinary electronic and optical properties and holds great promise for applications in photonics and optoelectronics. Demonstrations including high-speed photodetectors, optical modulators, plasmonic devices, and ultrafast lasers have now been reported. More advanced device concepts would involve photonic elements such as cavities to control light-matter interaction in graphene. Here we report the first monolithic integration of a graphene transistor and a planar, optical microcavity. We find that the microcavity-induced optical confinement controls the efficiency and spectral selection of photocurrent generation in the integrated graphene device. A twenty-fold enhancement of photocurrent is demonstrated. The optical cavity also determines the spectral properties of the electrically excited thermal radiation of graphene. Most interestingly, we find that the cavity confinement modifies the electrical transport characteristics of the integrated graphene transistor. Our experimental approach opens up a route towards cavity-quantum electrodynamics on the nanometre scale with graphene as a current-carrying intra-cavity medium of atomic thickness.

preprint2012arXiv

Structural diversity and the role of particle shape and dense fluid behavior in assemblies of hard polyhedra

A fundamental characteristic of matter is its ability to form ordered structures under the right thermodynamic conditions. Predicting these structures - and their properties - from the attributes of a material's building blocks is the holy grail of materials science. Here, we investigate the self-assembly of 145 hard convex polyhedra whose thermodynamic behavior arises solely from their anisotropic shape. Our results extend previous works on entropy-driven crystallization by demonstrating a remarkably high propensity for self-assembly and an unprecedented structural diversity, including some of the most complex crystalline phases yet observed in a non-atomic system. In addition to 22 Bravais and non-Bravais crystals, we report 66 plastic crystals (both Bravais and topologically close-packed), 21 liquid crystals (nematic, smectic, and columnar), and 44 glasses. We show that from simple measures of particle shape and local order in the disordered fluid, the class of ordered structure can be predicted.

preprint2011arXiv

Controlled fabrication of single-walled carbon nanotube electrodes by electron-beam-induced oxidation

The fabrication of metallic single-walled carbon nanotube electrodes separated by gaps of typically 20nm width by electron-beam-induced oxidation is studied within an active device configuration. The tube conductance is measured continuously during the process. The experiment provides a statistical evaluation of gap sizes as well as the electron dose needed for gap formation. Also, the ability to precisely cut many carbon nanotubes in parallel is demonstrated.

preprint2011arXiv

Crystalline Assemblies and Densest Packings of a Family of Truncated Tetrahedra and the Role of Directional Entropic Forces

Polyhedra and their arrangements have intrigued humankind since the ancient Greeks and are today important motifs in condensed matter, with application to many classes of liquids and solids. Yet, little is known about the thermodynamically stable phases of polyhedrally-shaped building blocks, such as faceted nanoparticles and colloids. Although hard particles are known to organize due to entropy alone, and some unusual phases are reported in the literature, the role of entropic forces in connection with polyhedral shape is not well understood. Here, we study thermodynamic self-assembly of a family of truncated tetrahedra and report several atomic crystal isostructures, including diamond, β-tin, and high- pressure lithium, as the polyhedron shape varies from tetrahedral to octahedral. We compare our findings with the densest packings of the truncated tetrahedron family obtained by numerical compression and report a new space filling polyhedron, which has been overlooked in previous searches. Interestingly, the self-assembled structures differ from the densest packings. We show that the self-assembled crystal structures can be understood as a tendency for polyhedra to maximize face-to-face alignment, which can be generalized as directional entropic forces.

preprint2011arXiv

Degenerate Quasicrystal of Hard Triangular Bipyramids

We report a degenerate quasicrystal in Monte Carlo simulations of hard triangular bipyramids each composed of two regular tetrahedra sharing a single face. The dodecagonal quasicrystal is similar to that recently reported for hard tetrahedra [Haji-Akbari et al., Nature (London) 462, 773 (2009)] but degenerate in the pairing of tetrahedra, and self-assembles at packing fractions above 54%. This notion of degeneracy differs from the degeneracy of a quasiperiodic random tiling arising through phason flips. Free energy calculations show that a triclinic crystal is preferred at high packing fractions.

preprint2011arXiv

Electroluminescence from chirality-sorted (9,7)-semiconducting carbon nanotube devices

We have measured the electroluminescence and photoluminescence of (9,7) semiconducting carbon nanotube devices and demonstrate that the electroluminescence wavelength is determined by the nanotube's chiral index (n,m). The devices were fabricated on Si3N4 membranes by dielectrophoretic assembly of tubes from monochiral dispersion. Electrically driven (9,7) devices exhibit a single Lorentzian shaped emission peak at 825 nm in the visible part of the spectrum. The emission could be assigned to the excitonic E22 interband transition by comparison of the electroluminescence spectra with corresponding photoluminescence excitation maps. We show a linear dependence of the EL peak width on the electrical current, and provide evidence for the inertness of Si3N4 surfaces with respect to the nanotubes optical properties.

preprint2011arXiv

Entropic Stabilization of Tunable Planar Modulated Superstructures

Self-assembling novel ordered structures with nanoparticles has recently received much attention. Here we use computer simulations to study a two-dimensional model system characterized by a simple isotropic interaction that could be realized with building blocks on the nanoscale. We find that the particles arrange themselves into hexagonal superstructures of twin boundaries whose superlattice vector can be tuned reversibly by changing the temperature. Thermodynamic stability is confirmed by calculating the free energy with a combination of thermodynamic integration and the Frenkel-Ladd method. Different contributions to the free energy difference are discussed.

preprint2011arXiv

Low-temperature structure of xi'-Al-Pd-Mn optimized by ab-initio methods

We have studied and resolved occupancy correlations in the existing average structure model of the complex metallic alloy xi'-Al-Pd-Mn [Boudard et al., Phil. Mag. A, 74, 939 (1996)], which has approximately 320 atoms in the unit cell and many fractionally occupied sites. Model variants were constructed systematically in a tiling-decoration approach and subjected to simulated annealing by use of both density functional theory and molecular dynamics with empirical potentials. To obtain a measure for thermodynamic stability, we reproduce the Al-Pd-Mn phase diagram at T=0K, and derive an enthalpy of formation for each structure. Our optimal structure resolves a cloud of fractionally occupied sites in pseudo-Mackay clusters. In particular, we demonstrate the presence of rotational degrees of freedom of an Al9 inner shell, which is caged within two icosahedrally symmetric outer shells Al30 and Pd12. Outside these clusters, the chemical ordering on a chain of three nearby sites surprisingly breaks the inversion symmetry of the surrounding structure, and couples to an Al/vacancy site nearby. Our refined tiling-decoration model applies to any structure within the epsilon-phases family, including the metastable decagonal quasicrystalline phase.

preprint2011arXiv

Phase Diagram of Hard Tetrahedra

Advancements in the synthesis of faceted nanoparticles and colloids have spurred interest in the phase behavior of polyhedral shapes. Regular tetrahedra have attracted particular attention because they prefer local symmetries that are incompatible with periodicity. Two dense phases of regular tetrahedra have been reported recently. The densest known tetrahedron packing is achieved in a crystal of triangular bipyramids (dimers) with packing density 4000/4671=85.63%. In simulation a dodecagonal quasicrystal is observed; its approximant, with periodic tiling (3.4.3^2.4), can be compressed to a packing fraction of 85.03%. Here, we show that the quasicrystal approximant is more stable than the dimer crystal for packing densities below 84% using Monte Carlo computer simulations and free energy calculations. To carry out the free energy calculations, we use a variation of the Frenkel-Ladd method for anisotropic shapes and thermodynamic integration. The enhanced stability of the approximant can be attributed to a network substructure, which maximizes the free volume (and hence the 'wiggle room') available to the particles and facilitates correlated motion of particles, which further contributes to entropy and leads to diffusion for packing densities below 65%. The existence of a solid-solid transition between structurally distinct phases not related by symmetry breaking -- the approximant and the dimer crystal-- is unusual for hard particle systems.

preprint2010arXiv

Dense crystalline dimer packings of regular tetrahedra

We present the densest known packing of regular tetrahedra with density phi = 4000/4671 = 0.856347... Like the recently discovered packings of Kallus et al. [arXiv:0910.5226] and Torquato-Jiao [arXiv:0912.4210], our packing is crystalline with a unit cell of four tetrahedra forming two triangular dipyramids (dimer clusters). We show that our packing has maximal density within a three-parameter family of dimer packings. Numerical compressions starting from random configurations suggest that the packing may be optimal at least for small cells with up to 16 tetrahedra and periodic boundaries.

preprint2010arXiv

Disordered, Quasicrystalline and Crystalline Phases of Densely Packed Tetrahedra

All hard, convex shapes are conjectured by Ulam to pack more densely than spheres, which have a maximum packing fraction of ϕ = π/\sqrt18 ~ 0.7405. For many shapes, simple lattice packings easily surpass this packing fraction. For regular tetrahedra, this conjecture was shown to be true only very recently; an ordered arrangement was obtained via geometric construction with ϕ = 0.7786, which was subsequently compressed numerically to ϕ = 0.7820. Here we show that tetrahedra pack much better than this, and in a completely unexpected way. Following a conceptually different approach, using thermodynamic computer simulations that allow the system to evolve naturally towards high-density states, we observe that a fluid of hard tetrahedra undergoes a first-order phase transition to a dodecagonal quasicrystal, which can be compressed to a packing fraction of ϕ = 0.8324. By compressing a crystalline approximant of the quasicrystal, the highest packing fraction we obtain is ϕ = 0.8503. If quasicrystal formation is suppressed, the system remains disordered, jams, and compresses to ϕ = 0.7858. Jamming and crystallization are both preceded by an entropy-driven transition from a simple fluid of independent tetrahedra to a complex fluid characterized by tetrahedra arranged in densely packed local motifs that form a percolating network at the transition. The quasicrystal that we report represents the first example of a quasicrystal formed from hard or non-spherical particles. Our results demonstrate that particle shape and entropy can produce highly complex, ordered structures.

preprint2010arXiv

Dynamics of particle flips in two-dimensional quasicrystals

The dynamics of quasicrystals is more complicated than the dynamics of periodic solids and difficult to study in experiments. Here, we investigate a decagonal and a dodecagonal quasicrystal using molecular dynamics simulations of the Lennard-Jones-Gauss interaction system. We observe that the short time dynamics is dominated by stochastic particle motion, so-called phason flips, which can be either single-particle jumps or correlated ring-like multi-particle moves. Over long times, the flip mechanism is efficient in reordering the quasicrystals and can generate diffusion. The temperature dependence of diffusion is described by an Arrhenius law. We also study the spatial distribution and correlation of mobile particles by analyzing the dynamic propensity.

preprint2010arXiv

Phason Dynamics in One-Dimensional Lattices

In quasicrystals, the phason degree of freedom and the inherent anharmonic potentials lead to complex dynamics which cannot be described by the usual phonon modes of motion. We have constructed simple one-dimensional model systems, the dynamic Fibonacci chain (DFC) and approximants thereof. They allow us to study the dynamics of periodic and quasiperiodic structures with anharmonic double well interactions both by analytical calculations and by molecular dynamics simulations. We found soliton modes like breathers and kink solitons and we were able to obtain closed analytical solutions for special cases, which are in good agreement with our simulations. Calculation of the specific heat shows an increase above the Dulong-Petit value, which is due to the anharmonicity of the potential and not caused by the phason degree of freedom.