-
Unraveling the interlayer and intralayer coupling in two-dimensional layered MoS$_2$ by X-ray absorption spectroscopy and ab initio molecular dynamics simulations
Authors:
Inga Pudza,
Dmitry Bocharov,
Andris Anspoks,
Matthias Krack,
Aleksandr Kalinko,
Edmund Welter,
Alexei Kuzmin
Abstract:
Understanding interlayer and intralayer coupling in two-dimensional layered materials (2DLMs) has fundamental and technological importance for their large-scale production, engineering heterostructures, and development of flexible and transparent electronics. At the same time, the quantification of weak interlayer interactions in 2DMLs is a challenging task, especially, from the experimental point…
▽ More
Understanding interlayer and intralayer coupling in two-dimensional layered materials (2DLMs) has fundamental and technological importance for their large-scale production, engineering heterostructures, and development of flexible and transparent electronics. At the same time, the quantification of weak interlayer interactions in 2DMLs is a challenging task, especially, from the experimental point of view. Herein, we demonstrate that the use of X-ray absorption spectroscopy in combination with reverse Monte Carlo (RMC) and ab initio molecular dynamics (AIMD) simulations can provide useful information on both interlayer and intralayer coupling in 2DLM 2H$_c$-MoS$_2$. The analysis of the low-temperature (10-300 K) Mo K-edge extended X-ray absorption fine structure (EXAFS) using RMC simulations allows for obtaining information on the means-squared relative displacements $σ^2$ for nearest and distant Mo-S and Mo-Mo atom pairs. This information allowed us further to determine the strength of the interlayer and intralayer interactions in terms of the characteristic Einstein frequencies $ω_E$ and the effective force constants $κ$ for the nearest ten coordination shells around molybdenum. The studied temperature range was extended up to 1200 K employing AIMD simulations which were validated at 300 K using the EXAFS data. Both RMC and AIMD results provide evidence of the reduction of correlation in thermal motion between distant atoms and suggest strong anisotropy of atom thermal vibrations within the plane of the layers and in the orthogonal direction.
△ Less
Submitted 2 June, 2023;
originally announced June 2023.
-
How to verify the precision of density-functional-theory implementations via reproducible and universal workflows
Authors:
Emanuele Bosoni,
Louis Beal,
Marnik Bercx,
Peter Blaha,
Stefan Blügel,
Jens Bröder,
Martin Callsen,
Stefaan Cottenier,
Augustin Degomme,
Vladimir Dikan,
Kristjan Eimre,
Espen Flage-Larsen,
Marco Fornari,
Alberto Garcia,
Luigi Genovese,
Matteo Giantomassi,
Sebastiaan P. Huber,
Henning Janssen,
Georg Kastlunger,
Matthias Krack,
Georg Kresse,
Thomas D. Kühne,
Kurt Lejaeghere,
Georg K. H. Madsen,
Martijn Marsman
, et al. (20 additional authors not shown)
Abstract:
In the past decades many density-functional theory methods and codes adopting periodic boundary conditions have been developed and are now extensively used in condensed matter physics and materials science research. Only in 2016, however, their precision (i.e., to which extent properties computed with different codes agree among each other) was systematically assessed on elemental crystals: a firs…
▽ More
In the past decades many density-functional theory methods and codes adopting periodic boundary conditions have been developed and are now extensively used in condensed matter physics and materials science research. Only in 2016, however, their precision (i.e., to which extent properties computed with different codes agree among each other) was systematically assessed on elemental crystals: a first crucial step to evaluate the reliability of such computations. We discuss here general recommendations for verification studies aiming at further testing precision and transferability of density-functional-theory computational approaches and codes. We illustrate such recommendations using a greatly expanded protocol covering the whole periodic table from Z=1 to 96 and characterizing 10 prototypical cubic compounds for each element: 4 unaries and 6 oxides, spanning a wide range of coordination numbers and oxidation states. The primary outcome is a reference dataset of 960 equations of state cross-checked between two all-electron codes, then used to verify and improve nine pseudopotential-based approaches. Such effort is facilitated by deploying AiiDA common workflows that perform automatic input parameter selection, provide identical input/output interfaces across codes, and ensure full reproducibility. Finally, we discuss the extent to which the current results for total energies can be reused for different goals (e.g., obtaining formation energies).
△ Less
Submitted 26 May, 2023;
originally announced May 2023.
-
CP2K: An Electronic Structure and Molecular Dynamics Software Package -- Quickstep: Efficient and Accurate Electronic Structure Calculations
Authors:
Thomas D. Kühne,
Marcella Iannuzzi,
Mauro Del Ben,
Vladimir V. Rybkin,
Patrick Seewald,
Frederick Stein,
Teodoro Laino,
Rustam Z. Khaliullin,
Ole Schütt,
Florian Schiffmann,
Dorothea Golze,
Jan Wilhelm,
Sergey Chulkov,
Mohammad Hossein Bani-Hashemian,
Valéry Weber,
Urban Borstnik,
Mathieu Taillefumier,
Alice Shoshana Jakobovits,
Alfio Lazzaro,
Hans Pabst,
Tiziano Müller,
Robert Schade,
Manuel Guidon,
Samuel Andermatt,
Nico Holmberg
, et al. (14 additional authors not shown)
Abstract:
CP2K is an open source electronic structure and molecular dynamics software package to perform atomistic simulations of solid-state, liquid, molecular and biological systems. It is especially aimed at massively-parallel and linear-scaling electronic structure methods and state-of-the-art ab-initio molecular dynamics simulations. Excellent performance for electronic structure calculations is achiev…
▽ More
CP2K is an open source electronic structure and molecular dynamics software package to perform atomistic simulations of solid-state, liquid, molecular and biological systems. It is especially aimed at massively-parallel and linear-scaling electronic structure methods and state-of-the-art ab-initio molecular dynamics simulations. Excellent performance for electronic structure calculations is achieved using novel algorithms implemented for modern high-performance computing systems. This review revisits the main capabilities of CP2k to perform efficient and accurate electronic structure simulations. The emphasis is put on density functional theory and multiple post-Hartree-Fock methods using the Gaussian and plane wave approach and its augmented all-electron extension.
△ Less
Submitted 11 March, 2020; v1 submitted 8 March, 2020;
originally announced March 2020.
-
Ab initio molecular dynamics simulations of negative thermal expansion in ScF3: the effect of the supercell size
Authors:
D. Bocharov,
M. Krack,
Yu. Rafalskij,
A. Kuzmin,
J. Purans
Abstract:
Scandium fluoride (ScF3) belongs to a class of negative thermal expansion (NTE) materials. It shows a strong lattice contraction up to about 1000 K switching to expansion at higher temperatures. Here the NTE effect in ScF3 is studied in the temperature range from 300 K to 1600 K using ab initio molecular dynamics (AIMD) simulations in the isothermal-isobaric (NpT) ensemble. The temperature depende…
▽ More
Scandium fluoride (ScF3) belongs to a class of negative thermal expansion (NTE) materials. It shows a strong lattice contraction up to about 1000 K switching to expansion at higher temperatures. Here the NTE effect in ScF3 is studied in the temperature range from 300 K to 1600 K using ab initio molecular dynamics (AIMD) simulations in the isothermal-isobaric (NpT) ensemble. The temperature dependence of the lattice constant, inter-atomic Sc-F-Sc bond angle distributions and the Sc-F and Sc-Sc radial distribution functions is obtained as a function of supercell size from 2a x 2a x 2a to 5a x 5a x 5a where a is the lattice parameter of ScF3. A comparison with the experimental Sc K-edge EXAFS data at 600 K is used to validate the accuracy of the AIMD simulations. Our results suggest that the AIMD calculations are able to reproduce qualitatively the NTE effect in ScF3, however a supercell size larger than 2a x 2a x 2a should be used to account accurately for dynamic disorder. The origin of the NTE in ScF3 is explained by the interplay between expansion and rotation of ScF6 octahedra.
△ Less
Submitted 24 February, 2020;
originally announced February 2020.
-
Comparative classical and ab initio Molecular Dynamics study of molten and glassy germanium dioxide
Authors:
Michael Hawlitzky,
Juergen Horbach,
Simona Ispas,
Matthias Krack,
Kurt Binder
Abstract:
A Molecular Dynamics (MD) study of static and dynamic properties of molten and glassy germanium dioxide is presented. The interactions between the atoms are modelled by the classical pair potential proposed by Oeffner and Elliott (OE) [Oeffner R D and Elliott S R 1998, Phys. Rev. B, 58, 14791]. We compare our results to experiments and previous simulations. In addition, an ab initio method, the…
▽ More
A Molecular Dynamics (MD) study of static and dynamic properties of molten and glassy germanium dioxide is presented. The interactions between the atoms are modelled by the classical pair potential proposed by Oeffner and Elliott (OE) [Oeffner R D and Elliott S R 1998, Phys. Rev. B, 58, 14791]. We compare our results to experiments and previous simulations. In addition, an ab initio method, the so-called Car-Parrinello Molecular Dynamics (CPMD), is applied to check the accuracy of the structural properties, as obtained by the classical MD simulations with the OE potential. As in a similar study for SiO2, the structure predicted by CPMD is only slightly softer than that resulting from the classical MD. In contrast to earlier simulations, both the static structure and dynamic properties are in very good agreement with pertinent experimental data. MD simulations with the OE potential are also used to study the relaxation dynamics. As previously found for SiO2, for high temperatures the dynamics of molten GeO2 is compatible with a description in terms of mode coupling theory.
△ Less
Submitted 16 April, 2008;
originally announced April 2008.
-
Coexistence of tetrahedral and octahedral-like sites in amorphous phase change materials
Authors:
S. Caravati,
M. Bernasconi,
T. D. Kuehne,
M. Krack,
M. Parrinello
Abstract:
Chalcogenide alloys are materials of interest for optical recording and non-volatile memories. We perform ab-initio molecular dynamics simulations aiming at shading light onto the structure of amorphous Ge2Sb2Te5 (GST), the prototypical material in this class. First principles simulations show that amorphous GST obtained by quenching from the liquid phase displays two types of short range order.…
▽ More
Chalcogenide alloys are materials of interest for optical recording and non-volatile memories. We perform ab-initio molecular dynamics simulations aiming at shading light onto the structure of amorphous Ge2Sb2Te5 (GST), the prototypical material in this class. First principles simulations show that amorphous GST obtained by quenching from the liquid phase displays two types of short range order. One third of Ge atoms are in a tetrahedral environment while the remaining Ge, Sb and Te atoms display a defective octahedral environment, reminiscent of cubic crystalline GST.
△ Less
Submitted 1 April, 2008; v1 submitted 9 August, 2007;
originally announced August 2007.
-
An Efficient and Accurate Car-Parrinello-like Approach to Born-Oppenheimer Molecular Dynamics
Authors:
Thomas D. Kühne,
Matthias Krack,
Fawzi R. Mohamed,
Michele Parrinello
Abstract:
We present a new method which combines Car-Parrinello and Born-Oppenheimer molecular dynamics in order to accelerate density functional theory based ab-initio simulations. Depending on the system a gain in efficiency of one to two orders of magnitude has been observed, which allows ab-initio molecular dynamics of much larger time and length scales than previously thought feasible. It will be dem…
▽ More
We present a new method which combines Car-Parrinello and Born-Oppenheimer molecular dynamics in order to accelerate density functional theory based ab-initio simulations. Depending on the system a gain in efficiency of one to two orders of magnitude has been observed, which allows ab-initio molecular dynamics of much larger time and length scales than previously thought feasible. It will be demonstrated that the dynamics is correctly reproduced and that high accuracy can be maintained throughout for systems ranging from insulators to semiconductors and even to metals in condensed phases. This development considerably extends the scope of ab-initio simulations.
△ Less
Submitted 20 December, 2006; v1 submitted 19 October, 2006;
originally announced October 2006.