For the study of crystal formation and dynamics, we introduce a simple two-dimensional monatomic model system with a parametrized interaction potential. We find in molecular dynamics simulations that a surprising variety of crystals, a decagonal, and a dodecagonal quasicrystal are self-assembled. In the case of the quasicrystals, the particles reorder by phason flips at elevated temperatures. During annealing, the entropically stabilized decagonal quasicrystal undergoes a reversible phase transition at 65% of the melting temperature into an approximant, which is monitored by the rotation of the de Bruijn surface in hyperspace.

Using angular momentum representation a method is proposed that allows the systematic construction of a generalized Landau--de Gennes elastic free energy of liquid crystals, in powers of a symmetric and traceless tensor order parameter, polarization field, of external fields and all respective derivatives. By this method all linearly independent elastic invariants and surface terms are constructed for nematics and cholesterics up to fourth order terms. In particular it is shown that up to fourth order in the tensor order parameter there are nineteen bulk elastic constants and four surface terms in the free energy of a general, biaxial nematic. In addition, the stability of this expansion is studied in detail. Some special cases of the elastic free energy of liquid crystals, already discussed in the literature, are reexamined and discrepancies with our results are emphasized. Finally, a thermodynamically correct way of establishing contact between the generalized de Gennes elastic free energy and other theories, like those of Oseen-Frank or Meyer, is proposed by applying fluctuation theory. Thus, the degeneracy of splay and bend elastic constants is removed even when these are calculated from the standard de Gennes free energy. Restrictions on higher order elastic constants are also obtained by comparing mean field relations and stability conditions with available experimental data.

We report on implementation and performance of the program IMD, designed for short range molecular dynamics simulations on massively parallel computers. After a short explanation of the cell-based algorithm, its extension to parallel computers as well as two variants of the communication scheme are discussed. We provide performance numbers for simulations of different sizes and compare them with values found in the literature. Finally we describe two applications, namely a very large scale simulation with more than 1.23×109 atoms, to our knowledge the largest published MD simulation up to this day and a simulation of a crack propagating in a two-dimensional quasicrystal.

No abstract

No abstract