Search

Lasers are becoming a more and more important tool in cutting and shaping materials. Improving precision and effectivity is an ongoing demand in science and industry. One possibility are double pulses. Here we study laser ablation of aluminum by the two-temperature model. There the laser is modeled as a source in a continuum heat conduction equation for the electrons, whose temperature then is transferred to a molecular dynamics particle model by an electron-phonon coupling term. The melting and ablation effectivity is investigated depending on the relative intensity and the time delay between two Gaussian shaped laser pulses. It turns out that at least for aluminum the optimal pulse shapes are standard Gaussian pulses. For double pulses with delay times up to 200 ps we find a behavior as observed in experiment: the ablation depth decreases beyond a delay of 10 ps even if one does not account for the weakening at the second pulse due to laser-plasma interaction., Comment: 18 pages, 6 figures

Experiments show that shallow nitrogen implantations (<10keV) result in a negatively charged nitrogen-vacancy center (NV-) yield of 0.01-0.1%. The most succesful technique for introducing NV- centers in the carbon matrix is ion implantation followed by annealing at 1100K. We investigated the influence of channeling effects during shallow implantation and statistical diffusion using molecular dynamics (MD) and Monte Carlo (MC) approaches. Energy barriers for the diffusion process were calculated using the density functional theory (DFT). Our simulations show a significant difference in the NV yield compared to the experiment. Statistically, 25% of the implanted nitrogens form a NV center after annealing., Comment: 5 pages, 6 figures

We have studied the dependence of metal oxide properties in molecular dynamics (MD) simulations on the polarizability of oxygen ions. We present studies of both liquid and crystalline structures of silica (SiO2), magnesia (MgO) and alumina (Al2O3). For each of the three oxides, two separately optimized sets of force fields were used: (i) Long-range Coulomb interactions between oxide and metal ions combined with a short-range pair potential. (ii) Extension of force field (i) by adding polarizability to the oxygen ions. We show that while an effective potential of type (i) without polarizable oxygen ions can describe radial distributions and lattice constants reasonably well, potentials of type (ii) are required to obtain correct values for bond angles and the equation of state. The importance of polarizability for metal oxide properties decreases with increasing temperature., Comment: 8 pages, 7 figures

We present an effective atomic interaction potential for crystalline alpha-Al2O3 generated by the program potfit. The Wolf direct, pairwise summation method with spherical truncation is used for electrostatic interactions. The polarizability of oxygen atoms is included by use of the Tangney-Scandolo interatomic force field approach. The potential is optimized to reproduce the forces, energies and stresses in relaxed and strained configurations as well as {0001}, {10-10} and {11-20} surfaces of Al2O3. Details of the force field generation are given, and its validation is demonstrated. We apply the developed potential to investigate crack propagation in alpha-Al2O3 single crystals., Comment: 8 pages, 5 figures

We extend the program potfit, which generates effective atomic interaction potentials from ab initio data, to electrostatic interactions and induced dipoles. The potential parametrization algorithm uses the Wolf direct, pairwise summation method with spherical truncation. The polarizability of oxygen atoms is modeled with the Tangney-Scandolo interatomic force field approach. Due to the Wolf summation, the computational effort in simulation scales linearly in the number of particles, despite the presence of electrostatic interactions. Thus, this model allows to perform large-scale molecular dynamics simulations of metal oxides with realistic potentials. Details of the implementation are given, and the generation of potentials for SiO2 and MgO is demonstrated. The approach is validated by simulations of microstructural, thermodynamic and vibrational properties of liquid silica and magnesia., Comment: 9 pages, 4 tables, 9 figures

Two-dimensional colloidal suspensions subject to laser interference patterns with decagonal symmetry can form an Archimedean-like tiling phase where rows of squares and triangles order aperiodically along one direction [J. Mikhael et al., Nature 454, 501 (2008)]. In experiments as well as in Monte-Carlo and Brownian dynamics simulations, we identify a similar phase when the laser field possesses tetradecagonal symmetry. We characterize the structure of both Archimedean-like tilings in detail and point out how the tilings differ from each other. Furthermore, we also estimate specific particle densities where the Archimedean-like tiling phases occur. Finally, using Brownian dynamics simulations we demonstrate how phasonic distortions of the decagonal laser field influence the Archimedean-like tiling. In particular, the domain size of the tiling can be enlarged by phasonic drifts and constant gradients in the phasonic displacement. We demonstrate that the latter occurs when the interfering laser beams are not adjusted properly.

We extend the Wolf direct, pairwise r^(-1) summation method with spherical truncation to dipolar interactions in silica. The Tangney-Scandolo interatomic force field for silica takes regard of polarizable oxygen atoms whose dipole moments are determined by iteration to a self-consistent solution. With Wolf summation, the computational effort scales linearly in the system size and can easily be distributed among many processors, thus making large-scale simulations of dipoles possible. The details of the implementation are explained. The approach is validated by estimations of the error term and simulations of microstructural and thermodynamic properties of silica., Comment: See http://link.aip.org/link/?JCP/132/194109 - 8 pages, 6 figures. Changes in v3: Copyright notice added, minor typographical changes. Changes in v2: 1. Inserted Paragraph in Sec. IV B describing the limitations of the TS potential. 2. We corrected transcription errors in Tab. II, and adjusted the deviation percentages mentioned in Sec. IV B, first paragraph, accordingly

The coincidence problem for planar patterns with $N$-fold symmetry is considered. For the N-fold symmetric module with $N<46$, all isometries of the plane are classified that result in coincidences of finite index. This is done by reformulating the problem in terms of algebraic number fields and using prime factorization. The more complicated case $N \ge 46$ is briefly discussed and N=46 is described explicitly. The results of the coincidence problem also solve the problem of colour lattices in two dimensions and its natural generalization to colour modules., Comment: 38 pages, 4 figures; updated and slightly corrected version of a paper published in J. Math. Phys. (see below)

Motivated by recent experimental findings, we investigate the possible occurrence and characteristics of quasicrystalline order in two-dimensional mixtures of point dipoles with two sorts of dipole moments. Despite the fact that the dipolar interaction potential does not exhibit an intrinsic length scale and cannot be tuned a priori to support the formation of quasicrystalline order, we find that configurations with long--range quasicrystallinity yield minima in the potential energy surface of the many particle system. These configurations emanate from an ideal or perturbed ideal decoration of a binary tiling by steepest descent relaxation. Ground state energy calculations of alternative ordered states and parallel tempering Monte-Carlo simulations reveal that the quasicrystalline configurations do not correspond to a thermodynamically stable state. On the other hand, steepest descent relaxations and conventional Monte-Carlo simulations suggest that they are rather robust against fluctuations. Local quasicrystalline order in the disordered equilibrium states can be strong., Comment: 10 pages, 7 figures

The coincidence problem for three-dimensional discrete structures with icosahedral symmetry is reinvestigated. We present a parametric description of the coincidence rotations based on special quaternions, called icosian numbers. In particular, we give a characterization of the possible coincidence indices (Sigma-factors) and present a complete list of the possible rotations with Sigma < 100., Comment: A version with abbreviated tables will occur in Z.Krist

We present a molecular dynamics study on atomic self-diffusion in Frank-Kasper type dodecagonal quasicrystals. It is found that the quasicrystal-specific flip mechanism for atomic diffusion as predicted by Kalugin and Katz, indeed occurs in this system. However, in order to be effective, this mechanism needs to be catalyzed by other defects such as half-vacancies if the structure is truly three-dimensional. For this reason, flip diffusion is difficult to distinguish from standard vacancy diffusion. In a quasi-two-dimensional setup, however, the flips may occur without other defects. Activation energies and flip frequencies are also determined., Comment: 6 pages, LaTeX, with 4 postscript figures