No Arabic abstract
Driven by the unprecedented computational power available to scientific research, the use of computers in solid-state physics, chemistry and materials science has been on a continuous rise. This review focuses on the software used for the simulation of matter at the atomic scale. We provide a comprehensive overview of major codes in the field, and analyze how citations to these codes in the academic literature have evolved since 2010. An interactive version of the underlying data set is available at https://atomistic.software .
Fidimag is an open-source scientific code for the study of magnetic materials at the nano- or micro-scale using either atomistic or finite difference micromagnetic simulations, which are based on solving the Landau-Lifshitz-Gilbert equation. In addition, it implements simple procedures for calculating energy barriers in the magnetisation through variants of the nudged elastic band method. This computer software has been developed with the aim of creating a simple code structure that can be readily installed, tested, and extended. An agile development approach was adopted, with a strong emphasis on automated builds and tests, and reproducibility of results. The main code and interface to specify simulations are written in Python, which allows simple and readable simulation and analysis configuration scripts. Computationally costly calculations are written in C and exposed to the Python interface as Cython extensions. Docker containers are shipped for a convenient setup experience. The code is freely available on GitHub and includes documentation and examples in the form of Jupyter notebooks.
3C-SiC epitaxially grown on Si displays a large wealth of extended defects. In particular, single, double and triple stacking faults (SFs) are observed in several experiments to coexist. Overabundance of defects has so far limited the exploitation of 3C-SiC/Si for power electronics, in spite of its several ideal properties (mainly in terms of wide gap, high breakdown fields and thermal properties) and the possibility of a direct integration in the Si technology. Here we use a multiscale approach, based on both classical molecular dynamics (MD) simulations and first-principle calculations, to investigate in-depth the origin, nature and properties of most common 3C-SiC/Si(001) extended defects. Our MD simulations reveal a natural path for the formation of partial dislocation complexes terminating both double and triple SFs. MD results are used as input for superior DFT calculations, allowing us to better determine the core structure and to investigate electronic properties. It turns out that the partial dislocation complexes terminating double and triple SFs are responsible for the introduction of electronic states significantly filling the gap. On the other hand, individual partial dislocations terminating single SFs only induce states very close to the gap edge. We conclude that partial dislocation complexes, in particular the most abundant triple ones, are killer defects in terms of favoring leakage currents. Suggestions coming from theory/simulations for devising a strategy to lower their occurrence are discussed.
We study the spin dynamics of a Heisenberg model at finite temperature in the presence of an external field or a uniaxial anisotropy. For the case of the uniaxial anisotropy our simulations show that the macro moment picture breaks down. An effect which we refer to as a spin-wave instability (SWI) results in a non-dissipative Bloch-Bloembergen type relaxation of the macro moment where the size of the macro moment changes, and can even be made to disappear. This relaxation mechanism is studied in detail by means of atomistic spin dynamics simulations.
First-principles electronic structure calculations are very widely used thanks to the many successful software packages available. Their traditional coding paradigm is monolithic, i.e., regardless of how modular its internal structure may be, the code is built independently from others, from the compiler up, with the exception of linear-algebra and message-passing libraries. This model has been quite successful for decades. The rapid progress in methodology, however, has resulted in an ever increasing complexity of those programs, which implies a growing amount of replication in coding and in the recurrent re-engineering needed to adapt to evolving hardware architecture. The Electronic Structure Library (esl) was initiated by CECAM (European Centre for Atomic and Molecular Calculations) to catalyze a paradigm shift away from the monolithic model and promote modularization, with the ambition to extract common tasks from electronic structure programs and redesign them as free, open-source libraries. They include heavy-duty ones with a high degree of parallelisation, and potential for adaptation to novel hardware within them, thereby separating the sophisticated computer science aspects of performance optimization and re-engineering from the computational science done by scientists when implementing new ideas. It is a community effort, undertaken by developers of various successful codes, now facing the challenges arising in the new model. This modular paradigm will improve overall coding efficiency and enable specialists (computer scientists or computational scientists) to use their skills more effectively. It will lead to a more sustainable and dynamic evolution of software as well as lower barriers to entry for new developers.
Thermal management is extremely important for designing high-performance devices. The lattice thermal conductivity of materials is strongly dependent on the structural defects at different length scales, particularly point defects like vacancies, line defects like dislocations, and planar defects such as grain boundaries. Traditionally, the McKelvey-Shockley phonon Boltzmanns transport equation (BTE) method combined with molecular dynamics simulations has been widely used to evaluate the phonon mean free paths (MFPs) in defective systems. However, this method can only provide the aggregate MFPs of the whole sample. It is, therefore, challenging to extract the MFPs in the different regions with different thermal properties. In this study, the 1D McKelvey-Shockley phonon BTE method was extended to model inhomogeneous materials, where the effect of defects on the phonon MFPs is explicitly obtained. Then, the method was used to study the phonon interactions with the core structure of an edge dislocation. The phonon MFPs in the dislocation core were obtained and consistent with the analytical model such that high frequency phonons are likely to be scattered in this area. This method not only advances the knowledge of phonon-dislocation scattering but also shows the potential to investigate phonon transport behaviors in more complicated materials.