en:
We propose a variational algorithm to study the real time dynamics of quantum systems as a ground-state problem on quantum devices. The method is based on the original proposal of Feynman and Kitaev to encode time into a register of auxiliary qubits.
We prepare the Feynman-Kitaev Hamiltonian acting on the composed system as a qubit operator and find an approximate ground state using the Variational Quantum Eigensolver. We apply the algorithm to the study of the dynamics of a transverse field Ising chain with an increasing number of spins and time steps, proving a favorable scaling in terms of the number of two qubit gates.
Through numerical experiments, we investigate its robustness against hardware noise, showing that the method can be use to evaluate dynamical properties of quantum systems and detect the presence of dynamical quantum phase transitions by measuring Loschmidt echoes.
The scripts provided implement the algorithm both in Python, using the Qiskit library, and in Julia, using the Yao.jl package.
Moreover, we included the data to reproduce the figures of the manuscript without running the simulations.

en: Variational dynamics as a ground-state problem on a quantum computer

quantum-computing

display_name: quantum-computing

quantum-simulation

display_name: quantum-simulation

real-time-evolution

display_name: real-time-evolution


$oid: OBSOLETE

en: The piecewise linearity condition on the total energy with respect to the total magnetization of finite quantum systems is derived using the infinite-separation-limit technique. This generalizes the well-known constancy condition, related to static correlation error, in approximate density functional theory. The magnetic analog of Koopmans' theorem in density functional theory is also derived. Moving to fractional electron count, the tilted-plane condition is derived, lifting certain assumptions in previous works. This generalization of the flat-plane condition characterizes the total energy surface of a finite system for all values of electron count N and magnetization M. This result is used in combination with tabulated spectroscopic data to show the flat-plane structure of the oxygen atom, among others. We find that derivative discontinuities with respect to electron count sometimes occur at noninteger values. A diverse set of tilted-plane structures is shown to occur in d-orbital subspaces, depending on chemical coordination. General occupancy-based total-energy expressions are demonstrated thereby to be necessarily dependent on the symmetry-imposed degeneracies.

en: Tilted-plane structure of the energy of finite quantum systems

materialscloud

approval_status: approved

en:
Materials Cloud is a platform built to enable the seamless sharing and dissemination of resources in computational materials science, offering educational, research, and archiving tools; simulation software and services; and curated and raw data. These underpin published results and empower data-based discovery, compliant with data management plans and the FAIR principles.

In the Materials Cloud you can browse, explore, download, or deposit raw and curated data (in the Discover, Explore, and Archive sections); you can access cloud or redeployable simulation services (from AiiDAlab turnkey simulations to the Quantum Mobile virtual machine, the AiiDA registry of plugins and workflows, and the OSSCAR notebooks, in the Work section), and you can use all our educational material, tutorials and lectures (in the Learn section).

en: Materials Cloud

de: Materials Cloud

fr: Materials Cloud

it: Materials Cloud

political_level: other

state: active

type: organization

psi

capacity: public

type: dataset

key: guid

key: uri

value: DCAT Harvester Materials Cloud

envi

de: Umwelt

en: Environment

fr: Environnement

it: Ambiente

density-functional-theory

display_name: density-functional-theory

static-correlation-error

display_name: static-correlation-error


en: Recent theory proposes that edge dislocations in random body-centered cubic (BCC) high entropy alloys have high barriers for motion, conveying high strengths up to high temperatures. Here, the energy barriers for edge motion are computed for two model alloys, NbTaV and MoNbTaW as represented by interatomic potentials, using the Nudged Elastic Band method and compared to theoretical predictions. The average magnitude of the barriers and the average spacing of the barriers along the glide direction agree well with the analytical theory, with no adjustable parameters. The evolution of the barriers versus applied stress is modeled, and the mean strength is in reasonable agreement with the predicted zero-temperature strength. These findings validate the analytic theory. A reduced analytic model based on solute misfit volumes is then applied to Hf-Mo-Nb-Ta-Ti-Zr and Mo-Nb-Ta-Ti-V-W alloys, rationalizing the observed significant strength increases at room temperature and 1000 ∘C upon addition of solutes with large misfit into a base alloy. The analytic theory for edge motion is thus a powerful validated tool for guiding alloy selection.

en: High energy barriers for edge dislocation motion in body-centered cubic high entropy alloys

atomistic-models

display_name: atomistic-models

high-entropy-alloys

display_name: high-entropy-alloys

snsf

display_name: snsf

solute-strengthening

display_name: solute-strengthening


en: The electronic density of states (DOS) provides information regarding the distribution of electronic energy levels in a material, and can be used to approximate its optical and electronic properties and therefore guide computational material design. Given its usefulness and relative simplicity, it has been one of the first electronic properties used as target for machine-learning approaches going beyond interatomic potentials. A subtle but important point, well-appreciated in the condensed matter community but usually overlooked in the construction of data-driven models, is that for bulk configurations the absolute energy reference of single-particle energy levels is ill-defined. Only energy differences matter, and quantities derived from the DOS are typically independent on the absolute alignment. We introduce an adaptive scheme that optimizes the energy reference of each structure as part of the training process, and show that it consistently improves the quality of ML models compared to traditional choices of energy reference, for different classes of materials and different model architectures. On a practical level, we trace the improved performance to the ability of this self-aligning scheme to match the most prominent features in the DOS. More broadly, we believe that this work highlights the importance of incorporating insights on the nature of the physical target into the definition of the architecture and of the appropriate figures of merit for machine-learning models, that translate in better transferability and overall performance.

en: Adaptive energy reference for machine-learning models of the electronic density of states

density-of-states

display_name: density-of-states

electronic-structure

display_name: electronic-structure

epfl

display_name: epfl

erc

display_name: erc

machine-learning

display_name: machine-learning

marvel

display_name: marvel


en: Atomistic simulations provide insights into structure-property relations on an atomic size and length scale that are complementary to the macroscopic observables that can be obtained from experiments.  Quantitative predictions, however, are usually hindered by the need to strike a balance between the accuracy of the calculation of the interatomic potential and the modelling of realistic thermodynamic conditions.  Machine-learning techniques make it possible to efficiently approximate the outcome of accurate electronic-structure calculations that can, therefore, be combined with extensive thermodynamic sampling.  We take elemental nickel as a prototypical material, whose alloys have applications from cryogenic temperatures up to close to their melting point and use it to demonstrate how a combination of machine-learning models of electronic properties and statistical sampling methods makes it possible to compute accurate finite-temperature properties at an affordable cost.  We demonstrate the calculation of a broad array of bulk, interfacial, and defect properties over a temperature range from 100 to 2500 K, modelling also, when needed, the impact of nuclear quantum fluctuations and electronic entropy.  The framework we demonstrate here can be easily generalized to more complex alloys and different classes of materials.

en: Finite-temperature materials modeling from the quantum nuclei to the hot electrons regime

cosmo

display_name: cosmo

electronic-density-of-states

display_name: electronic-density-of-states

materials-science

display_name: materials-science

metal-alloys

display_name: metal-alloys

molecular-dynamics-simulation

display_name: molecular-dynamics-simulation

neural-network-potentials

display_name: neural-network-potentials


en: With the recent advances in machine learning for quantum chemistry, it is now possible to predict the chemical properties of compounds and to generate novel molecules. Existing generative models mostly use a string- or graph-based representation, but the precise three-dimensional coordinates of the atoms are usually not encoded. First attempts in this direction have been proposed, where autoregressive or GAN-based models generate atom coordinates. Those either lack a latent space in the autoregressive setting, such that a smooth exploration of the compound space is not possible, or cannot generalize to varying chemical compositions. We propose a new approach to efficiently generate molecular structures that are not restricted to a fixed size or composition. Our model is based on the variational autoencoder which learns a translation-, rotation-, and permutation-invariant low-dimensional representation of molecules. Our experiments yield a mean reconstruction error below 0.05 Angstrom, outperforming the current state-of-the-art methods by a factor of four, and which is even lower than the spatial quantization error of most chemical descriptors. The compositional and structural validity of newly generated molecules has been confirmed by quantum chemical methods in a set of experiments.

en: 3DMolNet: a generative network for molecular structures

generative-model

display_name: generative-model

molecules-discovery

display_name: molecules-discovery

three-dimensional-structures

display_name: three-dimensional-structures


en: Time-dependent density functional theory has become state-of-the-art for describing photophysical and photochemical processes in extended materials due to its affordable cost. The inclusion of exact exchange was shown to be essential for the correct description of the long-range asymptotics of electronic interactions and thus a well-balanced description of valence, Rydberg and charge-transfer excitations. Several approaches for an efficient treatment of exact exchange have been established for the ground state, while implementations for excited-state properties are rare. Furthermore, the high computational costs required for excited-state properties in comparison to ground-state computations often hinder large-scale applications on periodic systems with hybrid functional accuracy. We therefore propose two approximate schemes for improving computational efficiency for the treatment of exact exchange. Within the auxiliary density matrix method (ADMM), exact exchange is estimated using a relatively small auxiliary basis and the introduced basis-set incompleteness error is compensated by an exchange density functional correction term. Benchmark results for a test set of 35 molecules demonstrate that the mean absolute error introduced by ADMM is smaller than 0.3 pm for excited-state bond lengths and in the range of 0.02 - 0.07 eV for vertical excitation, adiabatic excitation and fluorescence energies. Computational timings for a series of covalent-organic frameworks demonstrate that a speed-up of at least one order of magnitude can be achieved for ES geometry optimizations in comparison to conventional hybrid functionals. The second method is to use a semi-empirical tight binding approximation for both Coulomb and exchange contributions to the excited-state kernel. This simplified Tamm-Dancoff approximation (sTDA) achieves an accuracy comparable to approximated hybrid density functional theory when referring to highly accurate coupled-cluster reference data. We find that excited-state bond lengths deviate by 1.1 pm on average and mean absolute errors in vertical excitation, adiabatic excitation and fluorescence energies are in the range of 0.2 - 0.5 eV. In comparison to ADMM-approximated hybrid functional theory, sTDA accelerates the computation of broad-band excitation spectra by one order of magnitude, suggesting its potential use for large-scale screening purposes.

en: Excited-state properties for extended systems: efficient hybrid density functional methods

excited-state-properties

display_name: excited-state-properties

hybrid-density-functional-theory

display_name: hybrid-density-functional-theory

marie-curie-fellowship

display_name: marie-curie-fellowship

theoretical-spectroscopy

display_name: theoretical-spectroscopy


en: Doping and compositing are two universal design strategies used to engineer the electronic state of a material and mitigate its disadvantages. These two strategies have been extensively applied to the design of efficient electrocatalysts for water splitting. Using cobalt oxide (CoO) as a model catalyst, we prove that the oxygen evolution reaction (OER) performance could be progressively improved, first by Fe-doping to form Fe-CoO solid solution, and further by the addition of CeO2 to produce a Fe-CoO/CeO2 composite. X-ray adsorption spectroscopy (XAS) reveals that distinct electronic interactions are induced by the processes of doping and compositing. Fe-doping of CoO can break down the structural symmetry in the pristine material, changing the electronic structure of both Co and O species at the surface and decreasing the flat-band potential (Vfb). In comparison, subsequent compositing of Fe-CoO with CeO2 induces negligible electronic changes in the as-synthesized Fe-CoO (as seen in ex-situ characterizations), but significantly modifies the oxidative transformations of both Co and Fe under OER conditions. Our spectroscopic investigations reveal that Fe-doping and CeO2 compositing play different roles in modifying the electronic properties of CoO in its pristine state and during OER catalysis, respectively, in return, providing useful guidance for the design of more efficient electrocatalysts using these two strategies.

en: Spectroscopic investigations of complex electronic interactions by elemental doping and material compositing of cobalt oxide for enhanced oxygen evolution reaction activity

compositing

display_name: compositing

doping

display_name: doping

electronic-interaction

display_name: electronic-interaction

oxygen-evolution-reaction

display_name: oxygen-evolution-reaction

x-ray-adsorption-spectroscopy

display_name: x-ray-adsorption-spectroscopy


en: This database contains phosphorus and sulfur K-edge X-ray absorption near-edge structure (XANES) of delithiated Lithium-Phosphorus-Sulfur compounds. The structures were generated by a computational delithiation procedure of β-Li3PS4 with density functional theory calculations, where we systematically enumerated distinct lithium/vacancy orderings in the super cells with the general composition Li12-xP4S16 (i.e., four Li3PS4 formula units and x from 0 to 12). The database contains a total of 2227 P K-edge and 8885 S K-edge XANES spectra of symmetrically inequivalent absorbing P and S sites. The XANES spectra were calculated using the excited electron and core hole method as implemented in The Vienna Ab initio Simulation Package (VASP) 6.2.1. Details of the structure generation procedure as well as the VASP simulations can be found in the associated manuscript (see reference below).

en: Atomic insights into the oxidative degradation mechanisms of sulfide solid electrolytes

lithium-ion-batteries

display_name: lithium-ion-batteries

lithium-phosphorus-sulfide

display_name: lithium-phosphorus-sulfide

solid-electrolyte

display_name: solid-electrolyte

vasp

display_name: vasp

x-ray-absorption-spectroscopy

display_name: x-ray-absorption-spectroscopy


en: As exemplified by the growing interest in the quantum anomalous Hall effect, the research on topology as an organizing principle of quantum matter is greatly enriched from the interplay with magnetism. In this vein, we present a combined electrical and thermoelectrical transport study on the magnetic Weyl semimetal EuCd₂As₂. Unconventional contribution to the anomalous Hall and anomalous Nernst effects were observed both above and below the magnetic transition temperature of EuCd₂As₂, indicating the existence of significant Berry curvature. EuCd₂As₂ represents a rare case in which this unconventional transverse transport emerges both above and below the magnetic transition temperature in the same material. The transport properties evolve with temperature and field in the antiferromagnetic phase in a different manner than in the paramagnetic phase, suggesting different mechanisms to their origin. Our results indicate EuCd₂As₂ is a fertile playground for investigating the interplay between magnetism and topology, and potentially a plethora of topologically nontrivial phases rooted in this interplay. This record contains the raw data for this manuscript [Physical Review Letters 126, 076602 (2021)].

en: Unconventional transverse transport above and below the magnetic transition temperature in Weyl semimetal EuCd₂As₂

materialscloud

en:
Materials Cloud is a platform built to enable the seamless sharing and dissemination of resources in computational materials science, offering educational, research, and archiving tools; simulation software and services; and curated and raw data. These underpin published results and empower data-based discovery, compliant with data management plans and the FAIR principles.

In the Materials Cloud you can browse, explore, download, or deposit raw and curated data (in the Discover, Explore, and Archive sections); you can access cloud or redeployable simulation services (from AiiDAlab turnkey simulations to the Quantum Mobile virtual machine, the AiiDA registry of plugins and workflows, and the OSSCAR notebooks, in the Work section), and you can use all our educational material, tutorials and lectures (in the Learn section).

anomalous-hall-effect

display_name: anomalous-hall-effect

anomalous-nernst-effect

display_name: anomalous-nernst-effect

experimental

display_name: experimental

magnetic-weyl-semimetal

display_name: magnetic-weyl-semimetal

topological-spin-texture

display_name: topological-spin-texture


en: The development of enhanced sampling methods has greatly extended the scope of atomistic simulations, allowing long-time phenomena to be studied with accessible computational resources. Many such methods rely on the identification of an appropriate set of collective variables. These are meant to describe the system's modes that most slowly approach equilibrium. Once identified, the equilibration of these modes is accelerated by the enhanced sampling method of choice. An attractive way of determining the collective variables is to relate them to the eigenfunctions and eigenvalues of the transfer operator. Unfortunately, this requires knowing the long-term dynamics of the system beforehand, which is generally not available. However, we have recently shown that it is indeed possible to determine efficient collective variables starting from biased simulations. In this paper, we bring the power of machine learning and the efficiency of the recently developed on-the-fly probability enhanced sampling method to bear on this approach. The result is a powerful and robust algorithm that, given an initial enhanced sampling simulation performed with trial collective variables or generalized ensembles, extracts transfer operator eigenfunctions using a neural network ansatz and then accelerates them to promote sampling of rare events. To illustrate the generality of this approach we apply it to several systems, ranging from the conformational transition of a small molecule to the folding of a mini-protein and the study of materials crystallization.

en: Deep learning the slow modes for rare events sampling

materialscloud

en: Materials Cloud

envi

en: Environment

collective-variables

display_name: collective-variables

enhanced-sampling

display_name: enhanced-sampling

machine-learning

display_name: machine-learning

molecular-dynamics

display_name: molecular-dynamics

protein-folding

display_name: protein-folding

silicon-crystallization

display_name: silicon-crystallization


en:
Despite their incredible performance, it is well reported that deep neural networks tend to be overoptimistic about their prediction confidence. Finding effective and efficient calibration methods for neural networks is therefore an important endeavour towards better uncertainty quantification in deep learning. In this manuscript, we introduce a novel calibration technique named expectation consistency (EC), consisting of a post-training rescaling of the last layer weights by enforcing that the average validation confidence coincides with the average proportion of correct labels. First, we show that the EC method achieves similar calibration performance to temperature scaling (TS) across different neural network architectures and data sets, all while requiring similar validation samples and computational resources. However, we argue that EC provides a principled method grounded on a Bayesian optimality principle known as the Nishimori identity. Next, we provide an asymptotic characterization of both TS and EC in a synthetic setting and show that their performance crucially depends on the target function. In particular, we discuss examples where EC significantly outperforms TS.
This record provides the code for the paper "Expectation consistency for calibration of neural networks".

en: Expectation consistency for calibration of neural networks (code)

neural-networks

display_name: neural-networks

temperature-scaling

display_name: temperature-scaling

uncertainty-quantification

display_name: uncertainty-quantification


en: We propose Landau levels as a probe for the topological character of electronic bands in two-dimensional moiré superlattices. We consider two configurations of twisted double bilayer graphene (TDBG) that have very similar band structures, but show different valley Chern numbers of the flat bands. These differences between the AB-AB and AB-BA configurations of TDBG clearly manifest as different Landau level sequences in the Hofstadter butterfly spectra calculated using the tight-binding model. The Landau level sequences are explained from the point of view of the distribution of orbital magnetization in momentum space that is governed by the rotational C2 and time-reversal T symmetries. Our results can be readily extended to other twisted graphene multilayers and h-BN/graphene heterostructures thus establishing the Hofstadter butterfly spectra as a powerful tool for detecting the nontrivial valley band topology.

en: Landau levels as a probe for band topology in graphene moiré superlattices

band-topology

display_name: band-topology

landau-level

display_name: landau-level

moire-superlattice

display_name: moire-superlattice

twisted-double-bilayer-graphene

display_name: twisted-double-bilayer-graphene


en: CeO₂ greatly enhances the electrocatalytic oxygen evolution reaction (OER) activity of CoOx, though the enhancement mechanism beyond this synergy is yet to be understood. Here, operando hard X-ray adsorption spectroscopy (hXAS) is applied to monitor the Co K edge and Ce L₃ edge in CoOx/CeO₂ to shed light on the evolution of the Co and Ce oxidation states during OER. In addition, ex situ soft XAS (sXAS) characterizations provide information on the irreversible surface-specific transformations of the Co L₃ edge as well as of the O K edge. Combining the operando and ex situ spectroscopic characterizations with comprehensive electrochemical analyses, we confirm CeO₂ is not the active center for the OER. However, coupling CeO₂ with CoOx introduces significant modifications in the Co and O species at the CoOx surface and alters the flat band potential (Efb), leading to more favorable Co oxidation state transformations during OER and possibly modifying the preferential reaction pathway. This work establishes the connections between electronic structures, Co oxidation state and the OER reaction mechanism for CoOx/CeO₂ composites electrodes.

en: Operando tracking the interactions between CoOx and CeO₂ during oxygen evolution reaction

co-oxidation

display_name: co-oxidation

co-redox

display_name: co-redox

operando-xas

display_name: operando-xas

oxygen-evolution-reaction

display_name: oxygen-evolution-reaction


en: An ab initio  study of the coupled electrons and phonon transport properties of MoS2-hBN devices. A comparison of the device characteristics when hBN is treated as a perfectly insulating, non vibrating layer and one where it is included in the DFT together with MoS2.

en: Influence of the hBN dielectric layers on the quantum transport properties of MoS2 transistors.

ab-initio

display_name: ab-initio

quantum-transport

display_name: quantum-transport

thermal-properties

display_name: thermal-properties


en: A catalyst possessing a broad substrate scope, in terms of both turnover and enantioselectivity, is sometimes called "general". Despite their great utility in asymmetric synthesis, truly general catalysts are difficult or expensive to discover via traditional high-throughput screening and are, therefore, rare. Existing computational tools accelerate the evaluation of reaction conditions from a pre-defined set of experiments to identify the most general ones, but cannot generate entirely new catalysts with enhanced substrate breadth. For these reasons, we report an inverse design strategy based on the open-source genetic algorithm NaviCatGA and on the OSCAR database of organocatalysts to simultaneously probe the catalyst and substrate scope and optimize generality as primary target. We apply this strategy to the Pictet–Spengler condensation, for which we curate a database of 820 reactions, used to train statistical models of selectivity and activity. Starting from OSCAR, we define a combinatorial space of millions of catalyst possibilities, and perform evolutionary experiments on a diverse substrate scope that is representative of the whole chemical space of tetrahydro-β-carboline products. While privileged catalysts emerge, we show how genetic optimization can address the broader question of generality in asymmetric synthesis, extracting structure–performance relationships from the challenging areas of chemical space.

en: A genetic optimization strategy with generality in asymmetric organocatalysis as primary target

catalysis

display_name: catalysis

enantioselectivity

display_name: enantioselectivity

generality

display_name: generality

genetic-optimization

display_name: genetic-optimization

machine-learning

display_name: machine-learning

nccr-catalysis

display_name: nccr-catalysis

organic-molecules

display_name: organic-molecules

organocatalysis

display_name: organocatalysis


en: Raman spectroscopy is a widely-used non-destructive material characterization method, which provides information about the vibrational modes of the material and therefore of its atomic structure and chemical composition. Raman spectra can be simulated using atomistic first-principles methods but these are computationally demanding and thus the existing databases of computational Raman spectra are fairly small. We developed an optimized workflow to efficiently calculate the Raman tensors, from which the Raman spectra can be straightforwardly simulated. The workflow was benchmarked and validated by comparison to experiments and previous computational methods for select technologically relevant material systems. Using the workflow, we performed high-throughput calculations for a large set of materials (5099) belonging to many different material classes, and collected the results to a database.

en: High-throughput computation of Raman spectra from first principles

density-functional-theory

display_name: density-functional-theory

high-throughput

display_name: high-throughput

raman-tensors

display_name: raman-tensors


en: We present a density-functional theory (DFT) study of the structural, electronic, and chemical bonding behavior in germanium (Ge)-doped vanadium dioxide (VO2). Our motivation is to explain the reported increase of the metal-insulator transition temperature under Ge doping and to understand how much of the fundamental physics and chemistry behind it can be captured at the conventional DFT level. We model doping using a supercell approach, with various concentrations and different spatial distributions of Ge atoms in VO2. Our results suggest that the addition of Ge atoms strongly perturbs the high-symmetry metallic rutile phase and induces structural distortions that partially resemble the dimerization of the experimental insulating structure. Our work, therefore, hints at a possible explanation of the observed increase in transition temperature under Ge doping, motivating further studies into understanding the interplay of structural and electronic transitions in VO2.

en: Influence of germanium substitution on the structural and electronic stability of the competing vanadium dioxide phases

materialscloud

en:
Materials Cloud is a platform built to enable the seamless sharing and dissemination of resources in computational materials science, offering educational, research, and archiving tools; simulation software and services; and curated and raw data. These underpin published results and empower data-based discovery, compliant with data management plans and the FAIR principles.

In the Materials Cloud you can browse, explore, download, or deposit raw and curated data (in the Discover, Explore, and Archive sections); you can access cloud or redeployable simulation services (from AiiDAlab turnkey simulations to the Quantum Mobile virtual machine, the AiiDA registry of plugins and workflows, and the OSSCAR notebooks, in the Work section), and you can use all our educational material, tutorials and lectures (in the Learn section).

dft

display_name: dft

doping

display_name: doping


many-body-self-interaction-and-polarons

en: We address the many-body self-interaction in relation to polarons in density functional theory. Our study provides (i) a unified theoretical framework encompassing many-body and one-body forms of self-interaction and (ii) an efficient semilocal scheme for charge localization. Our theoretical formulation establishes a quantitative connection between the many-body and one-body forms of self-interaction in terms of electron screening, thereby conferring superiority to the concept of many-body self-interaction. Our semilocal methodology involves the use of a weak localized potential and applies equally to electron and hole polarons. We find that polarons free from many-body self-interaction have formation energies that are robust with respect to the functional adopted.

en: Many-body self-interaction and polarons

materialscloud

en: Materials Cloud

title_for_slug: many-body-self-interaction-and-polarons

electron-screening

display_name: electron-screening

many-body-self-interaction

display_name: many-body-self-interaction

mgo

display_name: mgo

one-body-self-interaction

display_name: one-body-self-interaction

polaron-localization

display_name: polaron-localization

polarons

display_name: polarons

self-interaction

display_name: self-interaction


en: Hydroxyapatite (HA, Ca10PO4(OH)2) is a widely explored material in the experimental domain of biomaterials science, because of its resemblance with natural bone minerals. Specifically, in the bioceramic community, HA doped with multivalent cations (e.g., Mg+2, Fe+2, Sr+2, etc.) has been extensively investigated in the last few decades. Experimental research largely established the critical role of dopant content on the changes in mechanical and biocompatibility properties. The plethora of experimental measurements of mechanical response on doped HA is based on compression or indentation testing of polycrystalline materials. Such measurements, as well as computational predictions of me,  on single crystalline (doped) HA are scarce. On that premise, the present study aims to build atomistic models of Fe2+-doped HA, a model system, with varying Fe content (10, 20, 30, and 40 mol%) and to explore their uniaxial tensile response by means of molecular dynamics (MD) simulation, together with the calculation of IR spectrum. In the equilibrated unit cell structures, Ca(1) sites were found to be energetically favourable for Fe2+ substitution. The local distribution of Fe2+ ions significantly affect the atomic partial charge distribution and chemical symmetry surrounding the functional groups. These signatures are reflected in the significant decrease in the intensity of IR peaks and are found in the Fe-doped HA, together with peak splitting because of the symmetry change in the crystal structure. Another important objective of this work is to computationally predict the mechanical response of doped HA in their single crystal format. An interesting observation is that the elastic anisotropy of undoped HA was not compromised with Fe-doping. Tensile strength (TS) is systematically reduced in doped HA with Fe+2 dopant content and a decrease in TS with temperature can be attributed to the increased thermal agitation of atoms at elevated temperatures. The physics of the tensile response was rationalized in terms of the strain dependent changes in covalent/ionic framework (Ca-P distance, P-O bond strain, O-P-O angular strain, O-H bond distance). Further, the dynamic changes in covalent bond network were energetically analyzed by calculating the changes in O-H and P-O bond vibrational energy. Summarizing, the current work develops our foundational understanding of the atomistic phenomena involved in the phase stability and tensile response of Fe-doped HA single crystals.

en: In silico study on probing atomistic insight into structural stability and tensile properties of Fe-doped hydroxyapatite single crystals

envi

en: Environment

hydroxyapatite

display_name: hydroxyapatite

iron-doping

display_name: iron-doping

md-simulation

display_name: md-simulation

mechanical-response

display_name: mechanical-response


en: We present a first-principles method to calculate electron-phonon coupling elements in atomic systems, and showcase its application to the evaluation of the phonon-limited mobility of n-type single-layer MoS₂. The method combines a density functional theory (DFT) plane-wave supercell approach with a real-space maximally localized Wannier basis. It enables the calculation of electronic structure, phonon displacements with their corresponding frequencies, and real-space electron-phonon coupling elements on the same footing, without the need for density functional perturbation theory (DFPT) or Wannier interpolation. We report a low-field, intrinsic mobility of 274 cm²/Vs at room temperature for MoS₂, and highlight its dependence on carrier density and temperature. In addition, we compare our findings to the latest available modeling data and put them in perspective with the experimentally measured values. Based on these observations, the mobilities presented in this work appear to be compatible with experimental results, when taking into account other scattering sources. Hence, the proposed approach provides a reliable framework for mobility calculations that can be extended towards large-scale device simulations.

en: Electron-phonon calculations using a wannier-based supercell approach: applications to the monolayer MoS₂ mobility

ab-initio

display_name: ab-initio

carrier-transport

display_name: carrier-transport

electron-mobility

display_name: electron-mobility

electron-phonon-coupling

display_name: electron-phonon-coupling

snsf

display_name: snsf


en:
In a recent preprint, entitled: "Development of a machine learning interatomic potential for exploring pressure-dependent kinetics of phase transitions in Germanium", we presented a novel Neural-Network (NN)  interatomic potential for Ge. We recall that Ge phases different from the cubic-diamond one are of particular interest for applications. Hexagonal Ge, for instance, displays superior optical properties. It is therefore important to investigate how, exploiting pressure, Ge can be transformed into different allotropes. In order to build a potential tackling kinetics of pressure-induced phase transformations, several kinetic paths (mainly sampled using the solid-state Nudged Elastic Band method) were added to the database, following a suitable active-learning procedure. Energies, forces, and stressed relative to the various configurations were computed ab initio using VASP with the PBE functional. The NN potential was trained using the Deep Potential Molecular Dynamic package (DeePMDkit). The potential greatly reproduces the relative stability of several Ge phases and yields at least a semi-quantitative description of the energetics along complex phase-transformation paths.
In the present archive, we provide the full potential for use in LAMMPS and ASE, together with the full database produced using VASP.

en: A NN-Potential for phase transformations in Ge

deepmd

display_name: deepmd

germanium

display_name: germanium

neural-network

display_name: neural-network

phase-transitions

display_name: phase-transitions


en: The functionality of atomic quantum emitters is intrinsically linked to their host lattice coordination. Structural distortions that spontaneously break the lattice symmetry strongly impact their optical emission properties and spin-photon interface. In a recent manuscript, we report on the direct imaging of charge state-dependent symmetry breaking of two prototypical atomic quantum emitters in mono- and bilayer MoS₂ by scanning tunneling microscopy (STM) and non-contact atomic force microscopy (nc-AFM). By changing the built-in substrate chemical potential, different charge states of sulfur vacancies (VacS) and substitutional rhenium dopants (ReMo) can be stabilized. VacS⁻¹ as well as ReMo⁰ and ReMo⁻¹ exhibit local lattice distortions and symmetry-broken defect orbitals attributed to a Jahn-Teller effect (JTE) and pseudo-JTE, respectively. By mapping the electronic and geometric structure of single point defects, we disentangle the effects of spatial averaging, charge multistability, configurational dynamics, and external perturbations that often mask the presence of local symmetry breaking. This record contains data to support the results discussed in our manuscript.

en: Charge state-dependent symmetry breaking of atomic defects in transition metal dichalcogenides

dfg

display_name: dfg

erc

display_name: erc

jahn-teller-effect

display_name: jahn-teller-effect

stm

display_name: stm


en: Field-effect transistors (FETs) based on two-dimensional materials (2DMs) with atomically thin channels have emerged as a promising platform for beyond-silicon electronics. However, low carrier mobility in 2DM transistors driven by phonon scattering remains a critical challenge. To address this issue, we propose the controlled introduction of localized tensile strain as an effective mean to inhibit electron-phonon scattering in 2DM. Strain is achieved by conformally adhering the 2DM via van-der-Waals forces to a dielectric layer previously nanoengineered with a gray-tone topography. Our results show that monolayer MoS₂ FETs under tensile strain achieve an 8-fold increase in on-state current, reaching mobilities of 185 cm²/Vs at room temperature, in good agreement with theoretical calculations. The present work on nanotopographic grayscale surface engineering and the use of high-quality dielectric materials has the potential to find application in the nanofabrication of photonic and nanoelectronic devices. Here we report the main data and codes needed to reproduce the theoretical findings of this work.

en: Deterministic grayscale nanotopography to engineer mobilities in strained MoS₂ FETs

density-functional-theory

display_name: density-functional-theory

first-principles

display_name: first-principles

transport

display_name: transport


en: Mn₅Si₃ displays an unusual and complex magnetic ground state, which is considered to be the origin of the anomalous transport and thermodynamic properties that it exhibits. We report the magnetic exchange couplings of the noncollinear antiferromagnetic phase of Mn₅Si₃ using inelastic neutron scattering measurements and density functional theory calculations. We determine the ground-state spin configuration and compute its magnon dispersion relations which are in good agreement with the ones obtained experimentally. Furthermore, we investigate the evolution of the spin texture under the application of an external magnetic field to demonstrate theoretically the multiple field-induced phase transitions that Mn₅Si₃ undergoes. Finally, we model the stability of some of the material's magnetic moments under a magnetic field and we find that very susceptible magnetic moments in a frustrated arrangement can be tuned by the field. This data set contains the data relevant to perform the DFT calculations with the JuKKR code, the spin dynamics simulations with the Spirit code, the spin-wave calculations with the SWIS code, and the neutron scattering experimental data.

en: Complex magnetic structure and spin waves of the noncollinear antiferromagnet Mn₅Si₃

antiferromagnetism

display_name: antiferromagnetism

frustration

display_name: frustration

inelastic-neutron-scattering

display_name: inelastic-neutron-scattering

intersect

display_name: intersect

marvel

display_name: marvel

noncollinear-magnetism

display_name: noncollinear-magnetism

spin-waves

display_name: spin-waves


fatigue-database-of-high-entropy-alloys

en: Fatigue failure of metallic structures is of great concern to industrial applications. A material will not be able to practically useful if it is prone to fatigue failure. To take the advantage of lately emerged high entropy alloys (HEAs) for designing novel fatigue-resistant alloys, we compiled a fatigue database of HEAs from the literature reported till the yearend of 2021. The database is subdivided into three categories, i.e., low-cycle fatigue (LCF), high-cycle fatigue (HCF), and fatigue crack growth rate (FCGR), which contains 15, 23, and 28 distinct data records, respectively. Each data record in any of three categories is characteristic of a summary, which is comprised of alloy composition, key fatigue properties, and additional information influential to or interrelated with fatigue (e.g., material processing history, phase constitution, grain size, uniaxial tensile properties, and fatigue testing conditions), and an individual dataset, which makes up the original fatigue testing curve.

en: Fatigue database of high entropy alloys

materialscloud

en:
Materials Cloud is a platform built to enable the seamless sharing and dissemination of resources in computational materials science, offering educational, research, and archiving tools; simulation software and services; and curated and raw data. These underpin published results and empower data-based discovery, compliant with data management plans and the FAIR principles.

In the Materials Cloud you can browse, explore, download, or deposit raw and curated data (in the Discover, Explore, and Archive sections); you can access cloud or redeployable simulation services (from AiiDAlab turnkey simulations to the Quantum Mobile virtual machine, the AiiDA registry of plugins and workflows, and the OSSCAR notebooks, in the Work section), and you can use all our educational material, tutorials and lectures (in the Learn section).

title_for_slug: fatigue-database-of-high-entropy-alloys

fatigue-crack-growth-rate

display_name: fatigue-crack-growth-rate

high-cycle-fatigue

display_name: high-cycle-fatigue

high-entropy-alloy

display_name: high-entropy-alloy

low-cycle-fatigue

display_name: low-cycle-fatigue


en:
We use a combination of symmetry analysis, phenomenological modeling, and first-principles density functional theory to explore the interplay between the magnetic ground state and the detailed atomic structure in the hexagonal rare-earth manganites. We find that the magnetic ordering is sensitive to a breathing mode distortion of the Mn and O ions in the ab plane, which is described by the K1 mode of the high-symmetry structure. Our density functional calculations of the magnetic interactions indicate that this mode particularly affects the single-ion anisotropy and the interplanar symmetric exchanges. By extracting the parameters of a magnetic model Hamiltonian from our first-principles results, we develop a phase diagram to describe the magnetic structure as a function of the anisotropy and exchange interactions. This in turn allows us to explain the dependence of the magnetic ground state on the identity of the rare-earth ion and on the K1 mode.
The attached files contain VASP input files and bash scripts for running all the DFT calculations in the manuscript.

en: Influence of the triangular Mn-O breathing mode on magnetic ordering in multiferroic hexagonal manganites

materialscloud

en: Materials Cloud

ab-initio

display_name: ab-initio

dftu

display_name: dftu

erc

display_name: erc


en: Dimensionality provides a clear fingerprint on the dispersion of infrared-active, polar-optical phonons. For these phonons, the local dipoles parametrized by the Born effective charges drive the LO-TO splitting of bulk materials; this splitting actually breaks down in two-dimensional materials. Here, we develop the theory for one-dimensional (1D) systems -nanowires, nanotubes, and atomic and polymeric chains. Combining an analytical model with the implementation of density-functional perturbation theory in 1D boundary conditions, we show that the dielectric splitting in the dispersion relations collapses as x²log(x) at the zone center. The dielectric properties and the radius of the 1D materials are linked by the present work to these red shifts, opening infrared and Raman characterization avenues.

en: Infrared-active phonons in one-dimensional materials and their spectroscopic signatures

envi

en: Environment

density-functional-theory

display_name: density-functional-theory

one-dimensional-materials

display_name: one-dimensional-materials

phonons

display_name: phonons

raman

display_name: raman


en:
Dynamical potentials appear in many advanced electronic-structure methods, including self-energies from many-body perturbation theory, dynamical mean-field theory, electronic-transport formulations, and many embedding approaches. Here, we propose a novel treatment for the frequency dependence, introducing an algorithmic inversion method that can be applied to dynamical potentials expanded as sum-over-poles. This approach allows for an exact solution of Dyson-like equations at all frequencies via a mapping to a matrix diagonalization, and provides simultaneously frequency-dependent (spectral) and frequency-integrated (thermodynamic) properties of the Dyson-inverted propagators.
The transformation to a sum-over-poles is performed introducing n-th order generalized Lorentzians as an improved basis set to represent the spectral function of a propagator.
Numerical results for the homogeneous electron gas at the G0W0 level are provided to argue for the accuracy and efficiency of such unified approach.
In this record, we provide the input and outputs of the AGWX suite (see article) used in the work. Also, we add the plots of the obtained Green's function for all densities, k-points, and frequencies.

en: A unified Green's function approach for spectral and thermodynamic properties from algorithmic inversion of dynamical potentials

computational-materials-science

display_name: computational-materials-science

homogeneous-electron-gas

display_name: homogeneous-electron-gas

many-body-perturbation-theory-gw

display_name: many-body-perturbation-theory-gw

max

display_name: max

snsf

display_name: snsf


en: The microscopic mechanism of heavy band formation, relevant for unconventional superconductivity in CeCoIn₅ and other Ce-based heavy fermion materials, depends strongly on the efficiency with which f electrons are delocalized from the rare earth sites and participate in a Kondo lattice. Replacing Ce³⁺ (4f 1, J = 5/2) with Sm³⁺ (4f 5, J = 5/2), we show that a combination of crystal field and on-site Coulomb repulsion causes SmCoIn₅ to exhibit a Γ7 ground state similar to CeCoIn5 with multiple f electrons. Remarkably, we also find that with this ground state, SmCoIn₅ exhibits a temperature-induced valence crossover consistent with a Kondo scenario, leading to increased delocalization of f holes below a temperature scale set by the crystal field, Tv ≈ 60 K. Our result provides evidence that in the case of many f electrons, the crystal field remains the most important tuning knob in controlling the efficiency of delocalization near a heavy fermion quantum critical point, and additionally clarifies that charge fluctuations are present and may play a generally important role in the ground state of "115" materials.

en: Charge fluctuations in the intermediate-valence ground state of SmCoIn₅

heavy-fermions

display_name: heavy-fermions

marvel

display_name: marvel

mixed-valence

display_name: mixed-valence


en: Linear scaling relations have led to an understanding of trends in catalytic activity and selectivity of many reactions in heterogeneous and electro-catalysis. Yet, linear scaling between the chemisorption energies of any two small molecule adsorbates is not guaranteed. A prominent example is the lack of scaling between the chemisorption energies of carbon and oxygen on transition metal surfaces. In this work, we show that this lack of scaling originates from different re-normalised adsorbate valence energies of lower-lying oxygen versus higher-lying carbon. We develop a model for chemisorption of small molecule adsorbates within the d-band model by combining a modified form of the Newns-Anderson hybridisation energy with an effective orthogonalization term. We develop a general descriptor to a priori determine if two adsorbates are likely to scale with each other. This record contains the AiiDA archive required to reproduce all calculations in the manuscript.

en: Limits to scaling relations between adsorption energies?

catalysis

display_name: catalysis

d-band-model

display_name: d-band-model

scaling-relations

display_name: scaling-relations


en: Existing, classical interatomic potentials for bcc iron predict contradicting crack-tip mechanisms (i.e. cleavage, dislocation emission, phase transition) for the same crack systems, thus leaving the crack propagation mechanism in bcc iron unclear. In this work, we develop a Gaussian approximation potential (GAP) by extending a DFT database for ferromagnetic bcc iron to include highly distorted primitive bcc cells and surface separation, along with small crack-tip configurations that are identified by means of a fully automated active learning workflow. Our GAP (referred to as Fe-GAP22) predicts crack propagation within 8 meV/atom accuracy. The fully automated, active learning workflow is made publicly available on GitHub. With the newly developed Fe-GAP22, we find that in absence of other defects around the crack tip (e.g. nanovoids, dislocations), the static (T=0K) crack-tip mechanism is cleavage, thus settling the contradictions in the literature. Our work also highlights the need for multi-scale modelling to predict fracture at finite temperatures and finite strain rates.

en: Atomistic fracture in bcc iron revealed by active learning of Gaussian approximation potential

bcc-iron

display_name: bcc-iron

fracture

display_name: fracture

gaussian-approximation-potential

display_name: gaussian-approximation-potential


en: The temperature-dependent properties of the aqueous electron have been extensively studied using mixed quantum-classical simulations in a wide range of thermodynamic conditions based on one-electron pseudopotentials. While the cavity model appears to explain most of the physical properties of the aqueous electron, only a non-cavity model has so far been successful in accounting for the temperature dependence of the absorption spectrum. Here, we present an accurate and efficient description of the aqueous electron under various thermodynamic conditions by combining hybrid functional-based molecular dynamics, machine learning techniques, and multiple time-step methods. Our advanced simulations accurately describe the temperature dependence of the absorption maximum in the presence of cavity formation. Specifically, our work reveals that the red shift of the absorption maximum results from an increasing gyration radius with temperature, rather than from global density variations as previously suggested.

en: Temperature dependent properties of the aqueous electron

ab-initio-molecular-dynamics

display_name: ab-initio-molecular-dynamics

aqueous-electron

display_name: aqueous-electron

hybrid-functional

display_name: hybrid-functional

machine-learning

display_name: machine-learning


en:
Machine-learned interatomic potentials are fast becoming an indispensable tool in computational materials science. One approach is the ephemeral data-derived potential (EDDP), which was designed to accelerate atomistic structure prediction. The EDDP is simple and cost-efficient. It relies on training data generated in small unit cells and is fit using a lightweight neural network, leading to smooth interactions which exhibit the robust transferability essential for structure prediction. Here, we present a variety of applications of EDDPs, enabled by recent developments of the open-source EDDP software. New features include interfaces to phonon and molecular dynamics codes, as well as deployment of the ensemble deviation for estimating the confidence in EDDP predictions. Through case studies ranging from elemental carbon and lead to the binary scandium hydride and the ternary zinc cyanide, we demonstrate that EDDPs can be trained to cover wide ranges of pressures and stoichiometries, and used to evaluate phonons, phase diagrams, superionicity, and thermal expansion. These developments complement continued success in accelerated structure prediction.
This record provides the data associated with the different case studies demonstrating the uses of ephemeral data-derived potentials in the article "Developments and Further Applications of Ephemeral Data Derived Potentials".

en: Developments and further applications of ephemeral data derived potentials

materialscloud

en:
Materials Cloud is a platform built to enable the seamless sharing and dissemination of resources in computational materials science, offering educational, research, and archiving tools; simulation software and services; and curated and raw data. These underpin published results and empower data-based discovery, compliant with data management plans and the FAIR principles.

In the Materials Cloud you can browse, explore, download, or deposit raw and curated data (in the Discover, Explore, and Archive sections); you can access cloud or redeployable simulation services (from AiiDAlab turnkey simulations to the Quantum Mobile virtual machine, the AiiDA registry of plugins and workflows, and the OSSCAR notebooks, in the Work section), and you can use all our educational material, tutorials and lectures (in the Learn section).

hydrides

display_name: hydrides

machine-learned-interatomic-potentials

display_name: machine-learned-interatomic-potentials

metal-organic-frameworks

display_name: metal-organic-frameworks

molecular-dynamics

display_name: molecular-dynamics

phase-diagram

display_name: phase-diagram

phonons

display_name: phonons

superionicity

display_name: superionicity


en: Motivated by the recent experimental discovery of strongly surface-plane-dependent superconductivity at surfaces of KTaO3 single crystals, we calculate the electron-phonon coupling strength, λ, of doped KTaO3 along the reciprocal-space high-symmetry directions. Using the Wannier-function approach implemented in the EPW package, we calculate λ across the experimentally covered doping range and compare its mode-resolved distribution along the [001], [110] and [111] reciprocal-space directions. We find that the electron-phonon coupling is strongest in the optical modes around the Γ point, with some distribution to higher k values in the [001] direction. The electron-phonon coupling strength as a function of doping has a dome-like shape in all three directions and its integrated total is largest in the [001] direction and smallest in the [111] direction, in contrast to the experimentally measured trends in critical temperatures. This disagreement points to a non-BCS character of the superconductivity. Instead, the strong localization of λ in the soft optical modes around Γ suggests an importance of ferroelectric soft-mode fluctuations, which is supported by our findings that the mode-resolved λ values are strongly enhanced in polar structures. The inclusion of spin-orbit coupling has negligible influence on our calculated mode-resolved λ values.

en: First-principles calculation of electron-phonon coupling in doped KTaO₃

materialscloud

en: Materials Cloud

ab-initio

display_name: ab-initio

electron-phonon-coupling

display_name: electron-phonon-coupling


en: This dataset includes all data presented in the figures of the main manuscript: XPS work function changes, ab-initio calculated work functions, O2p band center shifts, Densities of state for LSC and PCO, i-PLD results for LSC, LSF, STF and PCO, Densities of state for adsorbed O2 on differently decorated PCO surfaces. In this work, we investigated surface dipole changes which were induced by modification of mixed conducting surfaces with binary oxides via pulsed laser deposition.

en: Engineering surface dipoles on mixed conducting oxides with ultra-thin oxide decoration layers

envi

en: Environment

oxygen-exchange

display_name: oxygen-exchange

surface-chemistry

display_name: surface-chemistry

work-function

display_name: work-function


en: We show that hematite, 𝛼-Fe₂O₃, below its Morin transition, has a ferroic ordering of rank-5 magnetic triakontadipoles on the Fe ions. In the absence of spin-orbit coupling, these are the lowest-order ferroically aligned magnetic multipoles, and they give rise to the g-wave non-relativistic spin splitting in hematite. We find that the ferroically ordered magnetic triakontadipoles result from the simultaneous antiferroic ordering of the charge hexadecapoles and the magnetic dipoles, providing a route to manipulating the magnitude and the sign of the magnetic triakontadipoles as well as the spin splitting. Furthermore, we find that both the ferroic ordering of the magnetic triakontadipoles and many of the spin-split features persist in the weak ferromagnetic phase above the Morin transition temperature.

en: Non-relativistic ferromagnetotriakontadipolar order and spin splitting in hematite

erc

display_name: erc

multipoles

display_name: multipoles


en: Efficient O₂ reduction to H₂O₂, vital for energy conversion and environmental cleanup, relies on precise control of heterogeneous catalysts interacting with reaction species. Through high-throughput density functional theory calculations, consisting of 369 single atom catalysts, we identified the polarized descriptor (electric dipole strength) on two-dimensional carbon materials, revealing insights into the catalytic effect of support polarization. Surprisingly, this descriptor exhibits advanced scaling relationships towards H₂O₂ synthesis, incorporating factors such as active metals, coordination environments, and surface curvatures, highlighting its widespread significance. Furthermore, it demonstrates reliable predictability for O₂ adsorption in dynamic water environments, with optimal reactivity observed within the range of -1.40 to -1.00 e·Å, as confirmed by dynamic and static simulations of the 2e- pathway of O₂ reduction. In essence, these findings offer valuable insights for the rational design of electrocatalysts tailored for selective O₂ reduction.

en: Scaling relations and dynamical predictiveness of electric dipole strength on 2e- ORR catalytic property

density-functional-theory-calculations

display_name: density-functional-theory-calculations

electric-dipole-strength

display_name: electric-dipole-strength

oxygen-reduction-reaction

display_name: oxygen-reduction-reaction

support-polarization

display_name: support-polarization


en: The nature of the gap observed at the zone border in the spin excitation spectrum of CrI₃ quasi-two-dimensional single crystals is still controversial. We perform first-principles calculations based on time-dependent density functional perturbation theory, which indicate that the observed gap results from a combination of spin-orbit and interlayer interaction effects. The former give rise to the anisotropic spin-spin interactions that are responsible for its very existence, while the latter determine both its displacement from the K point of the Brillouin zone, due to the in-plane lattice distortions induced by them, and an enhancement of its magnitude, in agreement with experiments and previous theoretical work based on a lattice model.

en: First-principles study of the gap in the spin excitation spectrum of the CrI₃ honeycomb ferromagnet

dirac-cone

display_name: dirac-cone

magnons

display_name: magnons

marvel

display_name: marvel

quantum-espresso

display_name: quantum-espresso

tddft

display_name: tddft


en: The reactivity of silicates in aqueous solution is relevant to various chemistries ranging from silicate minerals in geology, to the C-S-H phase in cement, nanoporous zeolite catalysts, or highly porous precipitated silica. While simulations of chemical reactions can provide insight at the molecular level, balancing accuracy and scale in reactive simulations in the condensed phase is a challenge. Here, we demonstrate how a machine-learning reactive interatomic potential can accurately capture silicate-water reactivity. The model was trained on a new dataset comprising 400,000 energies and forces of molecular clusters at the 𝜔-B97XD\def2-TVZP level. To ensure the robustness of the model, we introduce a new and general active learning strategy based on the attribution of the model uncertainty, that automatically isolates uncertain regions of bulk simulations to be calculated as small-sized clusters. Our trained potential is found to reproduce static and dynamic properties of liquid water and solid crystalline silicates, despite having been trained exclusively on cluster data. Furthermore, we utilize enhanced sampling simulations to recover the self-ionization reactivity of water accurately, and the acidity of silicate oligomers, and lastly study the silicate dimerization reaction in a water solution at neutral conditions and find that the reaction occurs through a flanking mechanism.

en: Learning a reactive potential for silica-water through uncertainty attribution

active-learning

display_name: active-learning

chemical-reactions

display_name: chemical-reactions

free-energy

display_name: free-energy

neural-network-potential

display_name: neural-network-potential

silicates

display_name: silicates


en: Materials demonstrating positive thermal expansion (PTE) or negative thermal expansion (NTE) are quite common, whereas those exhibiting zero thermal expansion (ZTE) are notably scarce. In this work, we identify the mechanical descriptors, namely in-plane tensile stiffness and out-of-plane bending stiffness, that can effectively classify PTE and NTE 2D crystals. By utilizing high throughput calculations and the state-of-the-art symbolic regression method, these descriptors aid in the discovery of ZTE or 2D Invar monolayers with the linear thermal expansion coefficient (LTEC) within ±2×10⁻⁶ K⁻¹ in the middle range of temperatures. Additionally, the descriptors assist the discovery of large PTE and NTE 2D monolayers with the LTEC larger than ±15×10⁻⁶ K⁻¹, which are so-called 2D anti-Invar monolayers. Advancing our understanding of materials with exceptionally low or high thermal expansion is of substantial scientific and technological interest, particularly in developing next-generation electronics at the nanometer even Ångstrom scale.

en: Machine learning enables the discovery of 2D Invar and anti-Invar monolayers

bending-stiffness

display_name: bending-stiffness

linear-coefficient-of-thermal-expansion

display_name: linear-coefficient-of-thermal-expansion

machine-learning

display_name: machine-learning

symbolic-regression

display_name: symbolic-regression

thermal-expansion

display_name: thermal-expansion

zero-thermal-expansion

display_name: zero-thermal-expansion


en:
We performed a computational study using Density Functional Theory calculations on a copper-graphene system. A global minima search was performed using the Minima Hopping algorithm to find multiple graphene isomers that can appear on different substrates. As a model of copper substrate, we used a 112 copper atom cluster and as a model of graphene flake, we used 40 carbon atoms on top of the copper surface. The system was placed in vacuum to guarantee no interaction between images.

The dataset contains coordinates of graphene flake isomers on crystalline and amorphous Cu substrates.

en: Comparative study of defects in graphene flake grown on amorphous and crystalline copper surfaces

materialscloud

en:
Materials Cloud is a platform built to enable the seamless sharing and dissemination of resources in computational materials science, offering educational, research, and archiving tools; simulation software and services; and curated and raw data. These underpin published results and empower data-based discovery, compliant with data management plans and the FAIR principles.

In the Materials Cloud you can browse, explore, download, or deposit raw and curated data (in the Discover, Explore, and Archive sections); you can access cloud or redeployable simulation services (from AiiDAlab turnkey simulations to the Quantum Mobile virtual machine, the AiiDA registry of plugins and workflows, and the OSSCAR notebooks, in the Work section), and you can use all our educational material, tutorials and lectures (in the Learn section).

copper

display_name: copper

density-functional-theory

display_name: density-functional-theory

graphene

display_name: graphene

graphene-growth

display_name: graphene-growth

metal-substrate

display_name: metal-substrate


en: Fluctuations of electric fields can change the position of a gate-defined quantum dot in a semiconductor heterostructure. In the presence of magnetic field gradient, these stochastic shifts of electron's wavefunction lead to fluctuations of electron's spin splitting. The resulting spin dephasing due to charge noise limits the coherence times of spin qubits in isotopically purified Si/SiGe quantum dots. We investigate the spin splitting noise caused by such process caused by microscopic motion of charges at the semiconductor-oxide interface. We compare effects of isotropic and planar displacement of the charges, and estimate their densities and typical displacement magnitudes that can reproduce experimentally observed spin splitting noise spectra. We predict that for defect density of 10¹⁰ cm⁻², visible correlations between noises in spin splitting and in energy of electron's ground state in the quantum dot, are expected.

en: Correlations of spin splitting and orbital fluctuations due to 1/f charge noise in the Si/SiGe Quantum Dot

materialscloud

en: Materials Cloud

charge-noise

display_name: charge-noise

quantum-dots

display_name: quantum-dots

simulations

display_name: simulations


en: Recently a highly ordered Moiré dislocation lattice was identified at the interface between a SrTiO₃ (STO) thin film and the (LaAlO₃)₀.₃(Sr₂TaAlO₆)₀.₇ (LSAT) substrate. A fundamental understanding of the local ionic and electronic structure around the dislocation cores is crucial to further engineer the properties of these complex multifunctional heterostructures. Here we combine experimental characterization via analytical scanning transmission electron microscopy with results of molecular dynamics and density functional theory calculations to gain insights into the structure and defect chemistry of these dislocation arrays. Our results show that these dislocations lead to undercoordinated Ta/Al cations at the dislocation core, where oxygen vacancies can easily be formed, further facilitated by the presence of cation vacancies. The reduced Ti³⁺ observed experimentally at the dislocations by electron energy-loss spectroscopy are a consequence of both the structure of the dislocation itself, as well as of the electron-doping due to oxygen vacancy formation. Finally, the experimentally observed Ti diffusion into LSAT around the dislocation core occurs only together with cation-vacancy formation in LSAT or Ta diffusion into STO.

en: Combined theoretical and experimental study of the Moiré dislocation network at the SrTiO₃-(La,Sr)(Al,Ta)O₃ interface

envi

en: Environment

complex-oxide-interface

display_name: complex-oxide-interface

dislocation

display_name: dislocation

point-defects

display_name: point-defects


en: The electronic and vibrational properties and electron-phonon couplings of one-dimensional materials will be key to many prospective applications in nanotechnology. Dimensionality strongly affects these properties and has to be correctly accounted for in first-principles calculations. Here we develop and implement a formulation of density-functional and density-functional perturbation theory that is tailored for one-dimensional systems. A key ingredient is the inclusion of a Coulomb cutoff, a reciprocal-space technique designed to correct for the spurious interactions between periodic images in periodic-boundary conditions. This restores the proper one-dimensional open-boundary conditions, letting the true response of the isolated one-dimensional system emerge. In addition to total energies, forces and stress tensors, phonons and electron-phonon interactions are also properly accounted for. We demonstrate the relevance of the present method on a portfolio of realistic systems: BN atomic chains, BN armchair nanotubes, and GaAs nanowires. Notably, we highlight the critical role of the Coulomb cutoff by studying previously inaccessible polar-optical phonons and Fröhlich electron-phonon couplings. We also develop and apply analytical models to support the physical insights derived from the calculations and we discuss their consequences on electronic lifetimes. The present work unlocks the possibility to accurately simulate the linear response properties of one-dimensional systems, sheds light on the transition between dimensionalities and paves the way for further studies in several fields, including charge transport, optical coupling and polaritronics.

en: Density functional perturbation theory for one-dimensional systems: implementation and relevance for phonons and electron-phonon interactions

ab-initio

display_name: ab-initio

charge-transport

display_name: charge-transport

electron-phonon-coupling

display_name: electron-phonon-coupling

marvel

display_name: marvel

one-dimensional-materials

display_name: one-dimensional-materials

phonons

display_name: phonons


neural-network-potential-for-zr-h

en: The introduction of Hydrogen (H) into Zirconium (Zr) influences many mechanical properties, especially due to low H solubility and easy formation of Zirconium hydride phases. Understanding the various effects of H requires studies with atomistic resolution but at scales that incorporate defects such as cracks, interfaces, and dislocations. Such studies thus demand accurate interatomic potentials. Here, a neural network potential (NNP) for the Zr-H system is developed within the  Behler-Parrinello framework. The Zr-H NNP retains the accuracy of a recent NNP for hcp Zr and exhibits excellent agreement with first-principles density functional theory (DFT) for (i) H interstitials and their diffusion in hcp Zr, (ii) formation energies, elastic constants, and surface energies of relevant Zr hydrides, and (iii) energetics of a common Zr/Zr-H interface. The Zr-H NNP shows physical behavior for many different crack orientations in the most-stable ε-hydride and structures and reasonable relative energetics for the ⟨a⟩ screw dislocation in pure Zr. This Zr-H NNP should thus be very powerful for future study of many phenomena driving H degradation in Zr that require atomistic detail at scales far above those accessible by first-principles

en: Neural network potential for Zr-H

title_for_slug: neural-network-potential-for-zr-h

molecular-dynamics-simulation

display_name: molecular-dynamics-simulation

neural-network-potentials

display_name: neural-network-potentials

zirconium-hydrides

display_name: zirconium-hydrides


delta-project-archive-of-old-website

en: The widespread popularity of density functional theory has given rise to an extensive range of dedicated codes for predicting molecular and crystalline properties. However, each code implements the formalism in a different way, raising questions about the reproducibility of such predictions. In K. Lejaeghere et al., Science 351 (6280), aad3000 (2016) (https://doi.org/10.1126/science.aad3000), we reported the results of a community-wide effort that compared 15 solid-state codes, using 40 different potentials or basis set types, to assess the quality of the Perdew-Burke-Ernzerhof equations of state for 71 elemental crystals. We conclude that predictions from recent codes and pseudopotentials agree very well, with pairwise differences that are comparable to those between different high-precision experiments. Older methods, however, have less precise agreement. Our benchmark provides a framework for users and developers to document the precision of new applications and methodological improvements. The underlying data of this paper were accessible during 2016-2023 at https://molmod.ugent.be/deltacodesdft. Several additions that made use of the same protocol were added later. As this website is obsolete and to ensure long-term continued access, these data are stored now under the present record of the Materials Cloud Archive.

en: Delta project — archive of old website

title_for_slug: delta-project-archive-of-old-website

delta-project

display_name: delta-project

dft-benchmarking

display_name: dft-benchmarking

precision

display_name: precision

unary-crystals

display_name: unary-crystals

verification

display_name: verification


en: Spin excitations play a fundamental role in understanding magnetic properties of materials, and have significant technological implications for magnonic devices. However, accurately modeling these in transition-metal and rare-earth compounds remains a formidable challenge. Here, we present a fully first-principles approach for calculating spin-wave spectra based on time-dependent (TD) density-functional perturbation theory (DFPT), using nonempirical Hubbard functionals. This approach is implemented in a general noncollinear formulation, enabling the study of magnons in both collinear and noncollinear magnetic systems. Unlike methods that rely on empirical Hubbard U parameters to describe the ground state, and Heisenberg Hamiltonians for describing magnetic excitations, the methodology developed here probes directly the dynamical spin susceptibility (efficiently evaluated with TDDFPT throught the Liouville-Lanczos approach), and treats the linear variation of the Hubbard augmentation (in itself calculated non-empirically) in full at a self-consistent level. Furthermore, the method satisfies the Goldstone condition without requiring empirical rescaling of the exchange-correlation kernel or explicit enforcement of sum rules, in contrast to existing state-of-the-art techniques. We benchmark the novel computational scheme on prototypical transition-metal monoxides NiO and MnO, showing remarkable agreement with experiments and highlighting the fundamental role of these newly implemented Hubbard corrections. The method holds great promise for describing collective spin excitations in complex materials containing localized electronic states.

en: Magnons from time-dependent density-functional perturbation theory and nonempirical Hubbard functionals

first-principles

display_name: first-principles

hubbard-functionals

display_name: hubbard-functionals

magnetism

display_name: magnetism


en: Cement hydration is crucial for the strength development of cement-based materials; however, the mechanism that underlies this complex reaction remains poorly understood at the molecular level. An in-depth understanding of cement hydration is required for the development of environmentally friendly cement and consequently the reduction of carbon emissions in the cement industry. Here, we use molecular dynamics simulations with a reactive force field to investigate the initial hydration processes of tricalcium silicate (C₃S) and dicalcium silicate (C₂S) up to 40 ns. Our simulations provide theoretical support for the rapid initial hydration of C₃S compared to C₂S at the molecular level. The dissolution pathways of calcium ions in C₃S and C₂S are revealed, showing that, two dissolution processes are required for the complete dissolution of calcium ions in C₃S. Our findings promote the understanding of the calcium dissolution stage and serve as a valuable reference for the investigation of the initial cement hydration.

en: The initial stages of cement hydration at the molecular level

cement-hydration

display_name: cement-hydration

reaxff-reactive-force-field

display_name: reaxff-reactive-force-field


en: Covalent organic frameworks (COFs) are materials of growing interest for electronic applications due to their tunable structures, chemical stability, and layered architectures that support extended π-systems and directional charge transport. While their electronic properties are strongly influenced by the choice of molecular building blocks and the stacking arrangement, experimental control over these features remains limited, and the number of well-characterized COFs is still relatively small. Here, we explore two alternative strategies, hydrostatic pressure and metal intercalation, to tune the electronic structure of COFs. Using periodic density functional theory (DFT) calculations, we show that the band gap of pristine COF-1 decreases by ∼1 eV under compression up to 10 GPa. Metal intercalation induces an even greater reduction, in some cases leading to metallic behavior. We demonstrate that pressure and intercalation offer effective, continuous control over COF electronic properties, providing powerful means to complement and extend conventional design approaches.

en: Extensive band gap tunability in covalent organic frameworks via metal intercalation and high pressure

materialscloud

en:
Materials Cloud is a platform built to enable the seamless sharing and dissemination of resources in computational materials science, offering educational, research, and archiving tools; simulation software and services; and curated and raw data. These underpin published results and empower data-based discovery, compliant with data management plans and the FAIR principles.

In the Materials Cloud you can browse, explore, download, or deposit raw and curated data (in the Discover, Explore, and Archive sections); you can access cloud or redeployable simulation services (from AiiDAlab turnkey simulations to the Quantum Mobile virtual machine, the AiiDA registry of plugins and workflows, and the OSSCAR notebooks, in the Work section), and you can use all our educational material, tutorials and lectures (in the Learn section).

band-gap

display_name: band-gap

covalent-organic-frameworks

display_name: covalent-organic-frameworks

dft

display_name: dft

electronic-bands

display_name: electronic-bands

high-pressure

display_name: high-pressure

intercalation

display_name: intercalation


en: We report a first-principles study of the electronic and optical properties of recently identified one-dimensional semiconducting materials exfoliable from van der Waals-bonded bulk crystals. Specifically, we investigate four chalcogenide-based atomic chains, the covalently bonded S3 and Te3 chains, and polar-bonded As2S3 and Bi2Te3 chains, using a fully first-principles approach that combines density functional theory (DFT), density functional perturbation theory (DFPT), and many-body perturbation theory within the GW approximation and Bethe−Salpeter equation (BSE). Our vibrational analysis shows that the isolated, freestanding wires remain dynamically stable, with the zone-center optical phonon modes leading to infrared activity. The main finding of this study is the presence of very strong exciton binding energies (1−3 eV), which make these exfoliable 1D materials suitable platforms for room-temperature excitonic applications. Interestingly, the exciton character remains Wannier−Mott-like, as indicated by average electron−hole separations greater than the lattice constant. Notably, the optical gaps of these materials span a wide range - from infrared (0.8 eV, Bi2Te3), through the visible spectrum (yellow: 2.17 eV, Te3; blue: 2.71 eV, As2S3), up to ultraviolet (4.07 eV, S3) - highlighting their versatility for broadband optoelectronic applications. Our results offer a detailed, many-body perspective on the optoelectronic behavior of these low-dimensional materials and underscore their potential for applications in nanoscale optoelectronic devices.

en: Quasiparticle effects and strong excitonic features in exfoliable 1D semiconducting materials

materialscloud

en: Materials Cloud

exciton

display_name: exciton

gw

display_name: gw

one-dimensional-materials

display_name: one-dimensional-materials

phonons

display_name: phonons

van-der-waals-wires

display_name: van-der-waals-wires


en:
Noncollinear magnets are notoriously difficult to describe within first-principles approaches based on density-functional theory (DFT) because of the presence of low-lying spin excitations. At the level of ground-state calculations, several methods exist to constrain the magnetic moments to a predetermined configuration, and thereby accelerate convergence towards self-consistency. Their use in a perturbative context, however, remains very limited. Here we present a general methodological framework to achieve parametric control over the local spin moments at the linear-response level. Our strategy builds on the concept of Legendre transform to switch between various flavors of magnetic functionals, and to relate their second derivatives via simple linear-algebra operations. Thereby, we can address an arbitrary response function at the time-dependent DFT level of theory with optimal accuracy and minimal computational effort. In the low frequency limit, we identify the leading correction to the existing adiabatic formulation of the problem [S. Ren , Phys. Rev. X 14, 011041 (2024)], consisting in a renormalization of the phonon and magnon masses due to electron inertia. As a demonstration, we apply our methodology to the THz optical response of bulk CrI3 and Cr2O3, where we identify contributions from hybrid (electro)magnons with mixed spin-lattice character.

This dataset contains the input and output files used to produce the results presented in sections V and VI of the manuscript.

en: Dynamical response of noncollinear spin systems at constrained magnetic moments

envi

en: Environment

density-functional-perturbation-theory

display_name: density-functional-perturbation-theory

magnetic-insulators

display_name: magnetic-insulators

magnon-phonon-coupling

display_name: magnon-phonon-coupling

nonadiabatic-regime

display_name: nonadiabatic-regime


en:
A first-principles approach for calculating ion separation in solution through 2D membranes is proposed. Ionic energy profiles across the membrane are obtained first, where solvation effects are explicitly simulated by machine-learning molecular dynamics, electrostatic corrections are applied to remove finite-size capacitive effects, and a mean-field treatment of the electrochemical double layer charging is used. Entropic contributions are assessed analytically and through a thermodynamic integration scheme. Ionic separations are then inferred through a microkinetic model of the filtration process, accounting for steady-state charge separation effects across the membrane. The approach is applied to Li+, Na+, K+ sieving through a crown-ether functionalized graphene membrane, with a case study of the mechanisms for a highly selective and efficient extraction of lithium from aqueous solutions.
This record contains the MD trajectories used to generate the energy and free energy profiles of Fig. 4.

en: Ion sieving in 2D membranes from first principles

ion-sieving

display_name: ion-sieving

machine-learning

display_name: machine-learning

molecular-dynamics

display_name: molecular-dynamics


en:
Reproducibility in running scientific simulations on high-performance computing (HPC) environments is a persistent challenge due to variations in software and hardware stacks. Differences in software versions or hardware-specific optimizations often lead to discrepancies in simulation outputs. While Linux containers are commonly used to standardize software environments, tools like Docker lack reproducibility in image creation, requiring archiving of binary image blobs for future use. This method turns containers into black boxes, preventing verification of how the contained software was built.

In the linked paper, we demonstrate how we use GNU Guix to create our software stack bit-by-bit reproducible from a source bootstrap. Our approach incorporates a portable OpenMPI implementation, optimized software builds, and deployment via Apptainer images across three HPC environments. We show that our reproducible software stack facilitates consistent multi-physics simulations and complex workflows on diverse HPC platforms, exemplified by the OpenGeoSys software project. To ensure provenance of our findings, we utilized the AiiDA workflow manager.

This dataset includes the complete AiiDA provenance database underlying the results presented in the paper. The AiiDA workflow itself is defined in and can be reproduced with this repository: https://gitlab.opengeosys.org/bilke/hpc-container-study.

Version 1 of this record did not use optimized packages (Guix `--tune` option was not passed properly). This has been fixed in version 2 of this record.

en: Reproducible HPC software deployments, simulations and workflows

aiida

display_name: aiida

gnu-guix

display_name: gnu-guix

hpc

display_name: hpc

opengeosys

display_name: opengeosys

reproducibility

display_name: reproducibility

workflows

display_name: workflows


en: Precise modeling and understanding of heat transport in the superionic phase are of great interest. Although simulations combining Green-Kubo (GK) molecular dynamics with machine-learned potentials (MLPs) stand as a promising approach, substantial challenges remain due to the crucial impact of atomic diffusion. Here, we first show that the thermal conductivity (κ) of superionic materials calculated via conventional GK integral of the energy flux varies notably with the MLP model. Subsequently, we highlight that reliable, model-independent κ values can be obtained by applying Onsager’s reciprocal relations to correctly capture the coupled heat and mass transport. Remarkably, an anomalously invariant κ is observed over a wide temperature range, distinct from the characteristic trends in traditional crystals and glasses. Finally, we illustrate that conventional κ decompositions into kinetic, potential, and cross terms suffer from ambiguities in the physical interpretation, despite their mathematical rigor.

en: Heat transport in superionic materials via machine-learned molecular dynamics


en: The substantial consumption of lithium (Li) ions and sluggish reaction kinetics at the anode detrimentally impact the deliverable energy and fast-charging capability of lithium-ion batteries (LIBs) with silicon-based anodes. The prevailing contact prelithiation method using an electrolyte medium can replenish the active Li, but it may cause materials/electrode instability and bring barrier for Li-ion transport due to the nonuniform reaction. Here we explore a contact prelithiation methodology employing cyclic carbonate mediums that can enable spatially and temporally uniform prelithiation reaction. These mediums enable a delicate equilibrium between a Li-ion diffusion and the intrinsic prelithiation reaction rate of the material throughout the electrode depth. To gain an in-depth insight into the mechanism for producing Li-ions under the contact between Li metal and cyclic carbonate solvents, we undertook molecular dynamics simulations. Our calculation results of the self-produced Li alkyl carbonates reveal that while the long-chain Li alkyl carbonates are capable of forming an extensive polymer network in FEC and participating in the formation of a stable SEI on the Li metal surface, short-chain Li alkyl carbonates can dissolve into the solvent to produce Li-ions.

en: A salt-free medium facilitating electrode prelithiation towards fast-charging, high-energy lithium-ion batteries

battery

display_name: battery

lithium

display_name: lithium


en: Two-dimensional topological insulators are characterized by an insulating bulk and conductive edge states protected by the nontrivial topology of the bulk electronic structure. They remain robust against moderate disorder until Anderson localization occurs and destroys the topological phase. Interestingly, disorder can also induce a topological phase—known as a topological Anderson insulator—starting from an otherwise pristine trivial phase. While topological invariants are generally regarded as global quantities, we argue that space-resolved topological markers can act as local order parameters, revealing the role of fluctuations and correlations in the local topology under Anderson disorder and vacancies. With this perspective, we perform numerical simulations of disorder-driven topological phase transitions in the Haldane and Kane-Mele models, using supercells with both open and periodic boundary conditions. We find that short-scale fluctuations of topological markers vanish upon coarse-graining, except at the topological phase transition, where their correlation length peaks and large-scale fluctuations remain. Notably, such a topological correlation function is characterized by critical exponents that appear universal across disorder types, yet they can resolve different topological phase transitions.

en: Fluctuations and correlations of local topological order parameters in disordered two-dimensional topological insulators

anderson-disorder

display_name: anderson-disorder

critical-exponents

display_name: critical-exponents

fluctuations-and-correlations

display_name: fluctuations-and-correlations

haldane-model

display_name: haldane-model

kane-mele-model

display_name: kane-mele-model

local-topology

display_name: local-topology

topological-order-parameters

display_name: topological-order-parameters

vacancies

display_name: vacancies


en: Fullerides forms a big familiy of molecular crystals exhibiting various useful electro- and magnetochemical properties. Therefore, an efficient in silico method is desirable to look into such chemical system. To alleviate the burden of computationally-heavy first-principles calculations, here we present the successful attempts of using Machine Learning Force Field (MLFF) in predicting dynamical properties of the alkali-doped fulleride K₃C₆₀. Two on-the-fly Gaussian-Process-MLFF schemes based on different atomistic descriptors, Smooth Overlap of Atomic Position (SOAP) and Atomic Cluster Expansion (ACE), have been experimented. The performance of generated K₃C₆₀ MLFFs are validated by accurate prediction on energy and forces of 1,000 randomly disturbed K₃C₆₀ structures compared to respective DFT results. Several other dynamical properties including phonon dispersion, elastic moduli and heat capacity obtained by MLFFs have also shown good agreement with results from DFT calculations, whose computational cost on such post-processing is several times more expensive than training a single MLFF that performs effortless post-processing. This shows the potential of applying MLFF to complex molecular crystal systems and pave the way to the investigation of much more intricate fulleride properties.

en: Efficient modeling of dynamic properties in K₃C₆₀ using machine learning force fields

materialscloud

en:
Materials Cloud is a platform built to enable the seamless sharing and dissemination of resources in computational materials science, offering educational, research, and archiving tools; simulation software and services; and curated and raw data. These underpin published results and empower data-based discovery, compliant with data management plans and the FAIR principles.

In the Materials Cloud you can browse, explore, download, or deposit raw and curated data (in the Discover, Explore, and Archive sections); you can access cloud or redeployable simulation services (from AiiDAlab turnkey simulations to the Quantum Mobile virtual machine, the AiiDA registry of plugins and workflows, and the OSSCAR notebooks, in the Work section), and you can use all our educational material, tutorials and lectures (in the Learn section).

dynamical-properties

display_name: dynamical-properties

machine-learning-force-field

display_name: machine-learning-force-field

molecular-crystal

display_name: molecular-crystal


en: The electron-phonon Wannier interpolation (EPWI) method is an efficient way to compute the properties of electron-phonon interactions (EPIs) accurately. This study presents a GPU-accelerated implementation of the EPWI method for computing transport properties, followed by a performance analysis. The implementation is tested on common systems such as aluminum andsilicon. The results show complete consistency with those obtained through CPU computations. The proposed algorithm has the capability of computing the conductivity of aluminum in 20 minutes on a single NVIDIA Tesla V100 GPU, adopting a 2003 electron and phonon sampling grid. This speed is 54 times higher than the CPU-based algorithm, running on two nodes of the Intel Xeon Platinum 8260 CPU. Such impressive acceleration is achieved by carefully designing the algorithm to exploit the GPU’s specific features. The idea of this approach can also be applied to the GPU code for other EPI-related properties. This record includes the input/output files of this work. The code snippets of the core GPU implementation can be found in https://github.com/zheliu137/EPWI_BTE_GPU

en: A high-performance GPU implementation of the electron-phonon Wannier interpolation and the related transport properties

materialscloud

en: Materials Cloud

boltzmann-transport-equation

display_name: boltzmann-transport-equation

electron-phonon-coupling

display_name: electron-phonon-coupling

gpu-acceleration

display_name: gpu-acceleration


en: The existence of the hydroxyl ligand on a single- atom catalyst embedded in graphene, especially Fe(OH)−N4−C is suggested to improve the oxygen reduction reaction (ORR) from the pristine Fe−N4−C. However, the theoretical electrocatalytic activity of the ORR process is shown to be highly dependent on the electrolyte solution and the applied electrode potential. Herein, we performed the constant-potential simulation and microkinetic modeling to investigate the mechanism of ORR on well studied Fe(OH)−N4−C and Co(OH)−N4−C catalysts and compared to the pristine ones, using density functional theory (DFT) calculation combined with the effective screening medium method and the reference interaction site model (ESM-RISM). It was found that the Fe(OH)−N4−C and Co(OH)−N4−C have comparable ORR activities to Fe−N4−C and Co−N4−C, and the calculated limiting potentials and half-wave potentials agree with the experimental values, in contrast to the results obtained using the constant-charge method (i.e., under the neutral condition) regardless of consideration of the solvation effect.

en: Effect of the OH functionalization on the catalytic activity of single-atom catalyst

envi

en: Environment

density-functional-theory

display_name: density-functional-theory

implicit-solvation-model

display_name: implicit-solvation-model

oxygen-reduction-reaction

display_name: oxygen-reduction-reaction

single-atom-catalyst

display_name: single-atom-catalyst


en: Polyacene analogues, consisting of short acene segments separated by nonbenzenoid rings, offer intriguing electronic properties and magnetic interactions. Pentalene-bridged polyacenes (PPs), in particular, hold promise for enhancing the electrical conductivity and potential open-shell ground states. However, PPs have remained elusive in solution chemistry due to poor solubility and limited synthetic protocols. In a recent publication,  we report the on-surface synthesis of PPs through the annulation between ortho-xylene groups. Scanning tunneling microscopy and atomic force microscopy reveal that the reconstructed Au(110) surface significantly enhances the chemoselectivity of the annulation process. Scanning tunneling spectroscopy combined with density functional theory suggests that PP exhibits a narrow direct band gap, similar to long acenes. This work demonstrates the potential for band structure engineering in polyacene analogues by incorporating nonbenzenoid rings, paving the way to advancements in organic electronics and spintronics. This record contains computational data supporting the results discussed in the manuscript.

en: Enhancing chemoselectivity of on-surface reactions on reconstructed Au(110): the case of a pentalene-bridged polyacene analogue

annulation

display_name: annulation

dft

display_name: dft

hydrocarbons

display_name: hydrocarbons

on-surface-synthesis

display_name: on-surface-synthesis

scanning-tunneling-microscopy

display_name: scanning-tunneling-microscopy