AndreaBozzo/ceres-open-data-index huggingface dataset:
First-principles calculation of electron-phonon coupling in doped KTaO₃

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.

Tags: ab-initio, electron-phonon-coupling, ktao3
Ambipolar charge transfer of larger fullerenes enabled by the modulated surface potential of h-BN/Rh(111)

A detailed understanding of how molecules interact with two-dimensional materials, particularly concerning energy level alignment and charge transfer processes, is essential to incorporate functional molecular films into next-generation 2D material-organic hybrid devices. One of the major challenges in integrating molecular films in field-effect transistors is facilitating ambipolar charge transport, which is often hindered by the large electronic gap of the organic layers. In a recent work we compare the adsorption site-dependent energy level alignment of C60, C70, and C84 fullerenes induced by the spatial variation of the electrostatic surface potential of the h-BN/Rh(111) Moiré superstructure. As the size of the fullerenes increases, the HOMO-LUMO gap shrinks. In the case of C84, we find an intrinsic charge transfer from the substrate to the fullerenes adsorbed in the Moiré pore centers, rendering them negatively charged. The electric field effect-induced charging of neutral fullerenes and discharging of intrinsically negatively charged fullerenes are investigated using scanning tunneling spectroscopy, non-contact atomic force microscopy, and Kelvin probe force spectroscopy. Our findings show that on metal-supported h-BN, the LUMO level of C84 is sufficiently close to the Fermi energy that it can be neutral or 1e− negatively charged depending on slight variations of the electrostatic potential. The findings propose a path to make ambipolar charge transfer accessible and efficient by circumventing the need to overcome the fullerenes' electronic gap. This record contain data that support the results discussed in our work.

Tags: 2d-materials, ab-initio, ambipolar-charge-transfer, cscs, fullerenes, intermolecular-interactions, marvelp4, non-contact-atomic-force-microscopy, organic-hybrid-electronics, scanning-tunneling-microscopy, snsf
Range-separated hybrid functionals for accurate prediction of band gaps of extended systems

In this work, we systematically evaluate the accuracy in band gap prediction of range-separated hybrid functionals on a large set of semiconducting and insulating materials and carry out comparisons with the performance of their global counterparts. We observe that all the range-separated hybrid functionals that correctly describe the long-range dielectric screening significantly improve from standard hybrid functionals such as PBE0 and HSE06. Among this group, the choice of the short-range Fock exchange fraction and the screening length can further reduce the predicted error. We then propose a universal expression for the selection of the inverse screening parameter as a function of the short-range and long-range Fock exchange fractions, which results in a mean absolute error as small as 0.15 eV for band gap prediction.

Tags: density-functional-theory, hybrid-functional, quantum-espresso
Engineering surface dipoles on mixed conducting oxides with ultra-thin oxide decoration layers

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.

Tags: oxygen-exchange, surface-chemistry, work-function
Electronic and optical properties of the hydrogen boride sheet from the many-body perturbation theory

We study the electronic and optical properties of the hydrogen boride sheet by using the many-body perturbation theory with the perturbative GW (G₀W₀) approximation. It was found that the hydrogen boride sheet shows a semimetallic electronic structure, supporting the previous theoretical study based on the semilocal density functional theory calculations. It was also found that the optical spectrum calculated based on the quasiparticle energies agrees well with the experiments. This work suggests that G₀W₀ approximation may be useful for predicting precise electronic and optical properties of the hydrogen boride sheet and its derivatives.

Tags: density-functional-theory, dft, hydrogen-boride, many-body-perturbation-theory, nanosheet, two-dimensional
Directing CO₂ electroreduction pathways for selective C₂ product formation using single-site doped Cu catalysts

This dataset includes all the atomic coordinates of the optimized computational models from electronic structure calculations reported in the manuscript. In this work, we highlight that single-site noble metal dopants on the Cu surface can influence C–O bond dissociation in a key selectivity-determining intermediate (e.g., oxygen-bound *CH₂CHO), which in turn direct the post-C–C coupling pathways to ethylene versus ethanol. Combining theoretical and experimental analyses, we demonstrate that the favorability of C–O bond scission is controlled by the oxygen affinity of the metal dopant on the Cu catalyst. We find that the selectivity ratio of ethylene-to-ethanol displays a volcano relationship with respect to the oxygen binding strength on the doped surfaces.

Tags: co2-electroreduction, selectivity, single-site-doped-catalysts
Unraveling the effects of inter-site Hubbard interactions in spinel Li-ion cathode materials

Accurate first-principles predictions of the structural, electronic, magnetic, and electrochemical properties of cathode materials can be key in the design of novel efficient Li-ion batteries. Spinel-type cathode materials LixMn2O4 and LixMn1.5Ni0.5O4 are promising candidates for Li-ion battery technologies, but they present serious challenges when it comes to their first-principles modeling. Here, we use density-functional theory with extended Hubbard functionals - DFT+U+V with on-site U and inter-site V Hubbard interactions - to study the properties of these transition-metal oxides. The Hubbard parameters are computed from first-principles using density-functional perturbation theory. We show that while U is crucial to obtain the right trends in properties of these materials, V is essential for a quantitative description of the structural and electronic properties, as well as the Li-intercalation voltages. This work paves the way for reliable first-principles studies of other families of cathode materials without relying on empirical fitting or calibration procedures.

Tags: cscs, dftu, dftuv, extended-hubbard-functionals, li-ion-batteries, lithium-intercalation-voltage, marvelosp, spinel-cathode-materials
Conversion of La₂Ti₂O₇ to LaTiO₂N via ammonolysis: An ab-initio investigation

Perovskite oxynitrides are, due to their reduced band gap compared to oxides,  promising materials for photocatalytic applications. They are most commonly synthesized from {110} layered Carpy-Galy (A₂B₂O₇) perovskites via thermal ammonolysis, i.e. the exposure to a flow of ammonia at elevated temperature. The conversion of the layered oxide to the non-layered oxynitride must involve a complex combination of nitrogen incorporation, oxygen removal and ultimately structural transition by elimination of the interlayer shear plane. Despite the process being commonly used, little is known about the microscopic mechanisms and hence factors that could ease the conversion. Here we aim to derive such insights via density functional theory calculations of the defect chemistry of the oxide and the oxynitride as well as the oxide's surface chemistry. Our results point to the crucial role of surface oxygen vacancies in forming clusters of NH₃ decomposition products and in incorporating N, most favorably substitutionally at the anion site. N then spontaneously diffuses away from the surface, more easily parallel to the surface and in interlayer regions, while diffusion perpendicular to the interlayer plane is somewhat slower. Once incorporation and diffusion lead to a local N concentration of about 70 % of the stoichiometric oxynitride composition, the nitridated oxide spontaneously transforms to a nitrogen-deficient oxynitride.

Tags: ammonolysis, defects, marvel, oxynitride, snsf
Crack path selection in 2D crystals

Even though fracture in solids can be rationalized in continuum mechanics, local processes such as the path selection and the formation of facets and kinks along the crack edges cannot be resolved in this framework. The problem is addressed by developing a high-fidelity neural network-based force field NN-F3. An active-learning approach is taken to capture the high strain and non-equilibrium nature of the crack tips in 2D crystals such as graphene and h-BN, which has not been satisfactorily addressed by existing empirical and machine-learning force fields. Atomistic simulations using NN-F3 resolve spatial complexities from lattice-scale kinks to sample-scale crack patterns, which are discussed directly with continuum mechanics predictions and experimental observation. The results show that kinking or deflection of the cracks defines the roughness of cleaved edges and is explained by the alternation of the stress intensity factors. We also find that the selection of crack paths cannot be determined by the anisotropy in the energies of relaxed edges as widely referred to in the literature. The distortion and undercoordination of the cleaved edges play critical roles in the fracture process, which have to be incorporated into the models to predict the crack paths. Measures of fracture toughness are extracted from the fracture patterns, the critical stress intensity factors, or the energies of edges in the intermediate, unrelaxed states.

Tags: 2d-materials, fracture-anisotropy, fracture-mechanics, graphene, machine-learning, neural-newrok-force-field
Non-equilibrium nature of fracture determines the crack path

Local processes during fracture, such as the selection of crack paths and the formation of facets and kinks along the paths, are non-equilibrium and discrete in nature and cannot be resolved in the framework of continuum thermodynamics. Atomistic models using existing empirical or machine-learning force fields also fail to offer satisfactory solutions. The problem is addressed here by developing a high-fidelity neural network-based force field (NN-F³) following an active-learning approach that covers the space spanned by strain states up to material failure and the non-equilibrium, intermediate states of atomic-level structures during cracking. Atomistic simulations using NN-F³ identify spatial complexities from lattice-scale kinks to sample-scale crack patterns in 2D crystals such as graphene and h-BN. The selection of crack paths is found to be associated with alternation in the stress intensity factors and anisotropy in the fracture toughness, which defines the roughness of cleaved edges. We find that the non-uniform lattice distortion and undercoordination of cleaved edges at the crack front play critical roles in the fracture process, which cannot be described by the energy densities of relaxed edges as widely used in the literature. Instead, the fracture patterns, the critical stress intensity factors corresponding to the lattice kinks, and the energy densities of edges in the intermediate, unrelaxed states offer reasonable measures for the fracture toughness and its anisotropy, which can be determined from experiments or simulations with the atomic-level resolution and well explain the experimental observation.

Tags: 2d-materials, fracture-anisotropy, fracture-mechanics, graphene, machine-learning, neural-newrok-force-field
Non-equilibrium nature of fracture determines the crack path

A high-fidelity neural network-based force field (NN-F³) is developed to cover the space of strain states up to material failure and the non-equilibrium, intermediate nature of fracture. Simulations of fracture in 2D crystals using NN-F³ reveal spatial complexities from lattice-scale kinks to sample-scale patterns. We find that the fracture resistance cannot be captured by the energy densities of relaxed edges as used in the literature. Instead, the fracture patterns, critical stress intensity factors at the kinks, and energy densities of edges in the intermediate, unrelaxed states offer reasonable measures for the fracture toughness and its anisotropy.

Tags: 2d-materials, fracture-anisotropy, fracture-mechanics, graphene, machine-learning, neural-newrok-force-field
Crystal structure validation of verinurad via proton-detected ultra-fast MAS NMR and machine learning

The recent development of ultra-fast MAS (>100 kHz) provides new opportunities for structural characterization in solids. Here we use NMR crystallography to validate the structure of verinurad, a microcrystalline active pharmaceutical ingredient. To do this, we take advantage of 1H resolution improvement at ultra-fast MAS and use solely 1H-detected experiments and machine learning methods to assign all the experimental proton and carbon chemical shifts. This framework provides a new tool for elucidating chemical information from crystalline samples with limited sample volume and yields remarkably faster acquisition times compared to 13C-detected experiments, without the need to employ dynamic nuclear polarization.

Tags: nmr-crystallography, solid-state-nmr, verinurad
Non-relativistic ferromagnetotriakontadipolar order and spin splitting in hematite

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.

Tags: erc, h2020, multipoles
Hubbard U through polaronic defect states

Since the preliminary work of Anisimov and co-workers, the Hubbard corrected DFT+U functional has been used for predicting properties of correlated materials by applying on-site effective Coulomb interactions to specific orbitals. However, the determination of the Hubbard U parameter has remained under intense discussion despite the multitude of approaches proposed. Here, we define a selection criterion based on the use of polaronic defect states for the enforcement of the piecewise linearity of the total energy upon electron occupation. A good agreement with results from piecewise-linear hybrid functionals is found for the electronic and structural properties of polarons, including the formation energies. The values of U determined in this way are found to give a robust description of the polaron energetics upon variation of the considered state. In particular, we also address a polaron hopping pathway, finding that the determined value of U leads to accurate energetics without requiring a configurational-dependent U. It is emphasized that the selection of U should be based on physical properties directly associated with the orbitals to which U is applied, rather than on more global properties such as band gaps and band widths.  For comparison, we also determine U through a well-established linear-response scheme finding noticeably different values of U and consequently different formation energies. Possible origins of these discrepancies are discussed. As case studies, we consider the self-trapped electron in BiVO₄, the self-trapped hole in MgO, the Li-trapped hole in MgO, and the Al-trapped hole in 𝛼-SiO₂.

Tags: density-functional-theory, dftu, hybrid-functionals, polarons, self-interaction
On the sign of the linear magnetoelectric coefficient in Cr₂O₃

We establish the sign of the linear magnetoelectric (ME) coefficient, α, in chromia, Cr₂O₃. Cr₂O₃ is the prototypical linear ME material, in which an electric (magnetic) field induces a linearly proportional magnetization (polarization), and a single magnetic domain can be selected by annealing in combined magnetic (H) and electric (E) fields. Opposite antiferromagnetic domains have opposite ME responses, and which antiferromagnetic domain corresponds to which sign of response has previously been unclear. We use density functional theory (DFT) to calculate the magnetic response of a single antiferromagnetic domain of Cr₂O₃ to an applied in-plane electric field at 0 K. We find that the domain with nearest neighbor magnetic moments oriented away from (towards) each other has a negative (positive) in-plane ME coefficient, α⊥, at 0 K. We show that this sign is consistent with all other DFT calculations in the literature that specified the domain orientation, independent of the choice of DFT code or functional, the method used to apply the field, and whether the direct (magnetic field) or inverse (electric field) ME response was calculated. Next, we reanalyze our previously published spherical neutron polarimetry data to determine the antiferromagnetic domain produced by annealing in combined E and H fields oriented along the crystallographic symmetry axis at room temperature. We find that the antiferromagnetic domain with nearest-neighbor magnetic moments oriented away from (towards) each other is produced by annealing in (anti-)parallel E and H fields, corresponding to a positive (negative) axial ME coefficient, α∥, at room temperature. Since α⊥ at 0 K and α∥ at room temperature are known to be of opposite sign, our computational and experimental results are consistent. This dataset contains the input data to reproduce the calculation of the magnetoelectric effect as plotted in Fig. 3 of the manuscript, for Elk, Vasp, and Quantum Espresso.

Tags: cr2o3, dft, horizon-europe, magnetoelectric
Sparse tensor based nuclear gradients for periodic Hartree–Fock and low-scaling correlated wave function methods in the CP2K software package: A massively parallel and GPU accelerated implementation.

The development of novel double-hybrid density functionals offers new levels of accuracy, and is leading to fresh insights into the fundamental properties of matter. Hartree–Fock exact exchange and correlated wave function methods such as MP2 and direct RPA are usually required to build such functionals. Their high computational cost is a concern, and their application to large and periodic systems is therefore limited. In this work, low-scaling methods for HFX, SOS-MP2 and direct RPA energy gradients are developed and implemented in the CP2K software package. The use of the resolution-of-the-identity approximation with a short range metric and atom-centered basis functions lead to sparsity, allowing for sparse tensor contractions to take place. These operations are efficiently performed with the newly developed DBT and DBM libraries, which scale to hundreds of GPU nodes. The resulting methods, RI-HFX, SOS-MP2 and dRPA, were benchmarked on large supercomputers. They exhibit favorable sub-cubic scaling with system size, good strong scaling performance, and GPU acceleration up to a factor of 3. These developments will allow for double-hybrid level calculations of large and periodic condensed phase systems to take place on a more regular basis.

This record contains all CP2K input and output files used for the paper.

Tags: correlated-wave-function, gpu, hartree-fock, hpc, low-scaling, nuclear-gradients, rpa, sos-mp2, sparse-tensors
Mechanism of C-N bonds formation in electrocatalytic urea production revealed by ab initio molecular dynamics simulation

Electrosynthesis of urea from CO2 and NOX provides an exceptional opportunity for human society, given the increasingly available renewable energy. To raise the overall electrosynthesis efficiency, the most critical reaction step for such electrosynthesis, C-N coupling, needs to be significantly improved. The C-N coupling can only happen at a narrow potential window, generally in the low overpotential region, and a fundamental understanding of the C-N coupling is needed for further development of this strategy. In this regard, we performed ab initio Molecular Dynamics (AIMD) simulations to reveal the origin of C-N coupling under a small electrode potential window with both the dynamic nature of water as a solvent, and the electrode potentials considered. We explored the key reaction networks for urea formation on Cu(100) surface in neutral electrolytes. Our work shows excellent agreement with experimentally observed selectivity under different potentials on the Cu electrode. We discovered that the *NH and *CO are the key precursors for C-N bonds formation at low overpotential, while at high overpotential the C-N coupling occurs between adsorbed *NH and solvated CO. These insights provide vital information for future spectroscopic measurements and enable us to design new electrochemical systems for more value-added chemicals.

Tags: aimd, electrocatalysis, nox-reduction