No Arabic abstract
A comprehensive, critical study of the vibrational, thermodynamic and thermoelastic properties of bcc iron is presented, using well established semi-empirical embedded-atom method potentials available in the literature. Classical molecular dynamics simulations are used to address temperature effects, where dynamical matrices are constructed as a time average of the second moment of the atomic displacements. The $C_{11}, C_{44}, C$ elastic constants are then obtained from the sound velocities along high symmetry directions in reciprocal space. Results are compared to ultrasonic measurements and highlight the limitations of the potentials considered here in describing thermoelastic properties.
We investigate the topology of the spin-polarized charge density in bcc and fcc iron. While the total spin-density is found to possess the topology of the non-magnetic prototypical structures, in some cases the spin-polarized densities are characterized by unique topologies; for example, the spin-polarized charge densities of bcc and high-spin fcc iron are atypical of any known for non-magnetic materials. In these cases, the two spin-densities are correlated: the spin-minority electrons have directional bond paths with deep minima in the minority density, while the spin-majority electrons fill these holes, reducing bond directionality. The presence of two distinct spin topologies suggests that a well-known magnetic phase transition in iron can be fruitfully reexamined in light of these topological changes. We show that the two phase changes seen in fcc iron (paramagnetic to low-spin and low-spin to high-spin) are different. The former follows the Landau symmetry-breaking paradigm and proceeds without a topological transformation, while the latter also involves a topological catastrophe.
Large scale computer simulations are used to elucidate a longstanding controversy regarding the existence, or otherwise, of spin waves in paramagnetic BCC iron. Spin dynamics simulations of the dynamic structure factor of a Heisenberg model of Fe with first principles interactions reveal that well defined peaks persist far above Curie temperature T_c. At large wave vectors these peaks can be ascribed to propagating spin waves, at small wave vectors the peaks correspond to over-damped spin waves. Paradoxically, spin wave excitations exist despite only limited magnetic short-range order at and above T_c.
The vibrational properties of graphene fluoride and graphane are studied using ab initio calculations. We find that both sp3 bonded derivatives of graphene have different phonon dispersion relations and phonon density of states as expected from the different masses associated with the attached atoms fluorine and hydrogen, respectively. These differences manifest themselves in the predicted temperature behavior of the constant-volume specific heat of both compounds.
The vibrational density of states (VDOS) of nanoclusters and nanocrystalline materials are derived from molecular-dynamics simulations using empirical tight-binding potentials. The results show that the VDOS inside nanoclusters can be understood as that of the corresponding bulk system compressed by the capillary pressure. At the surface of the nanoparticles the VDOS exhibits a strong enhancement at low energies and shows structures similar to that found near flat crystalline surfaces. For the nanocrystalline materials an increased VDOS is found at high and low phonon energies, in agreement with experimental findings. The individual VDOS contributions from the grain centers, grain boundaries, and internal surfaces show that, in the nanocrystalline materials, the VDOS enhancements are mainly caused by the grain-boundary contributions and that surface atoms play only a minor role. Although capillary pressures are also present inside the grains of nanocrystalline materials, their effect on the VDOS is different than in the cluster case which is probably due to the inter-grain coupling of the modes via the grain-boundaries.
Large-scale atomistic simulations with classical potentials can provide valuable insights into microscopic deformation mechanisms and defect-defect interactions in materials. Unfortunately, these assets often come with the uncertainty of whether the observed mechanisms are based on realistic physical phenomena or whether they are artifacts of the employed material models. One such example is the often reported occurrence of stable planar faults (PFs) in body-centered cubic (bcc) metals subjected to high strains, e.g., at crack tips or in strained nano-objects. In this paper, we study the strain dependence of the generalized stacking fault energy (GSFE) of {110} planes in various bcc metals with material models of increasing sophistication, i.e., (modified) embedded atom method, angular-dependent, Tersoff, and bond-order potentials as well as density functional theory. We show that under applied tensile strains the GSFE curves of many classical potentials exhibit a local minimum which gives rise to the formation of stable PFs. These PFs do not appear when more sophisticated material models are used and have thus to be regarded as artifacts of the potentials. We demonstrate that the local GSFE minimum is not formed for reasons of symmetry and we recommend including the determination of the strain-dependent (110) GSFE as a benchmark for newly developed potentials.