npj Computational Materials最新文献

The quantification of interphase properties between metals and their corresponding hydrides is crucial for modeling the thermodynamics and kinetics of the hydrogenation processes in solid-state hydrogen storage materials. In particular, interphase boundary energies assume a pivotal role in determining the kinetics of nucleation, growth, and coarsening of hydrides, alongside accompanying morphological evolution during hydrogenation. The total interphase energy arises from both chemical bonding and mechanical strains in these solid-state systems. Since these contributions are usually coupled, it is challenging to distinguish via conventional computational approaches. Here, a comprehensive atomistic modeling methodology is developed to decouple chemical and mechanical energy contributions using first-principles calculations, of which feasibility is demonstrated by quantifying chemical and elastic strain energies of key interfaces within the FeTi metal-hydride system. Derived materials parameters are then employed for mesoscopic micromechanical analysis, predicting crystallographic orientations in line with experimental observations. The multiscale approach outlined verifies the importance of the chemo-mechanical interplay in the morphological evolution of growing hydride phases, and can be generalized to investigate other systems. In addition, it can streamline the design of atomistic models for the quantitative evaluation of interphase properties between dissimilar phases and allow for efficient predictions of their preferred phase boundary orientations.

The self-consistent field (SCF) generation of the three-dimensional (3D) electron density distribution (ρ) represents a fundamental aspect of density functional theory (DFT) and related first-principles calculations, and how one can shorten or bypass the SCF loop represents a critical question in electronic structure theory from both practical and fundamental standpoints. Herein, a machine learning strategy, DeepSCF, is presented in which the map between the SCF ρ and the initial guess density (ρ₀) constructed by the summation of neutral atomic densities is learned using 3D convolutional neural networks (CNNs). High accuracy and transferability of DeepSCF are achieved by first encoding ρ₀ on a 3D grid and then expanding the input features to include atomic fingerprints beyond ρ₀. The prediction of the residual density (δρ) rather than ρ itself is targeted, and given that δρ is indicative of chemical bonding information, a dataset of small-sized organic molecules featuring diverse bonding characters is adopted. The fidelity of DeepSCF is finally enhanced by subjecting the atomic geometries of the dataset to random rotations and strains. The effectiveness of DeepSCF is demonstrated using a complex carbon nanotube-based DNA sequencer model. This work evidences that the nearsightedness in electronic structure can be optimally represented via the spatial locality in CNNs, offering insight into the success of various machine learning-based atomistic materials simulations.

The conversion efficiency from charge current to spin current via the spin Hall effect is evaluated by the spin Hall ratio (SHR). Through state-of-the-art ab initio calculations involving both charge conductivity and spin Hall conductivity, we report the SHRs of the III-V monolayer family, revealing an ultrahigh ratio of 0.58 in the hole-doped GaAs monolayer. In order to find more promising 2D materials, a descriptor for high SHR is proposed and applied to a high-throughput database, which provides the fully relativistic band structures and Wannier Hamiltonians of 216 exfoliable monolayer semiconductors and has been released to the community. Among potential candidates for high SHR, the MXene monolayer Sc₂CCl₂ is identified with the proposed descriptor and confirmed by computation, demonstrating the descriptor validity for high SHR materials discovery.

Single-atom catalysts (SACs) with multiple active sites exhibit high activity for a wide range of sluggish reactions, but identifying optimal multimetallic SAC is challenging due to the vast design space. Here, we present a self-driving computational strategy that combines first-principles calculations and equivariant graph neural network (GNN) to explore over 30,000 binary metallic sites with varying combinations of 3d transition metals and different ligand environments for oxygen reduction and evolution reactions (ORR/OER). Active learning facilitates the investigation of the search space by balancing the exploration of unseen atomic structures with the exploitation of the active ones. The GNN learns the chemical environments to capture composition-structure-property relationships for ORR/OER activity and selectivity. The computational predictions of promising Co-Fe, Co-Co, and Co-Zn metal pairs are consistent with the state-of-the-art results of experimental measurements reported in the literature. This approach can be extended to a broader class of multi-element high entropic materials systems.

Materials datasets usually contain many redundant (highly similar) materials due to the tinkering approach historically used in material design. This redundancy skews the performance evaluation of machine learning (ML) models when using random splitting, leading to overestimated predictive performance and poor performance on out-of-distribution samples. This issue is well-known in bioinformatics for protein function prediction, where tools like CD-HIT are used to reduce redundancy by ensuring sequence similarity among samples greater than a given threshold. In this paper, we survey the overestimated ML performance in materials science for material property prediction and propose MD-HIT, a redundancy reduction algorithm for material datasets. Applying MD-HIT to composition- and structure-based formation energy and band gap prediction problems, we demonstrate that with redundancy control, the prediction performances of the ML models on test sets tend to have relatively lower performance compared to the model with high redundancy, but better reflect models’ true prediction capability.

Crystalline solids play a fundamental role in a host of materials and technologies, ranging from pharmaceuticals to renewable energy. The thermodynamic properties of these solids are crucial determinants of their stability and therefore their behavior. The advent of large density functional theory databases with properties of solids has stimulated research on predictive methods for their thermodynamic properties, especially for the enthalpy of formation Δ_fH. Increasingly sophisticated artificial intelligence and machine learning (ML) models have primarily driven development in this field in recent years. However, these models can suffer from lack of generalizability and poor interpretability. In this work, we explore a different route and develop and evaluate a framework for the application of reaction network (RN) theory to the prediction of Δ_fH of crystalline solids. For an experimental dataset of 1550 compounds we are able to obtain a mean absolute error w.r.t Δ_fH of 29.6 meV atom⁻¹ using the RN approach. This performance is better than existing ML-based predictive methods and close to the experimental uncertainty. Moreover, we show that the RN framework allows for straightforward estimation of the uncertainty of the predictions.

By simultaneously achieving high spatial and angular sampling resolution, four dimensional scanning transmission electron microscopy (4D STEM) is enabling analysis techniques that provide great insight into the atomic structure of materials. Applying these techniques to scientifically and technologically significant beam-sensitive materials remains challenging because the low doses needed to minimise beam damage lead to noisy data. We demonstrate an unsupervised deep learning model that leverages the continuity and coupling between the probe position and the electron scattering distribution to denoise 4D STEM data. By restricting the network complexity it can learn the geometric flow present but not the noise. Through experimental and simulated case studies, we demonstrate that denoising as a preprocessing step enables 4D STEM analysis techniques to succeed at lower doses, broadening the range of materials that can be studied using these powerful structure characterization techniques.

Strain engineering is a widely used technique for enhancing the mobility of charge carriers in semiconductors, but its effect is not fully understood. In this work, we perform first-principles calculations to explore the variations of the mobility of electrons and holes in silicon upon deformation by uniaxial strain up to 2% in the [100] crystal direction. We compute the π₁₁ and π₁₂ electron piezoresistances based on the low-strain change of resistivity with temperature in the range 200 K to 400 K, in excellent agreement with experiment. We also predict them for holes which were only measured at room temperature. Remarkably, for electrons in the transverse direction, we predict a minimum room-temperature mobility about 1200 cm²V⁻¹s⁻¹ at 0.3% uniaxial tensile strain while we observe a monotonous increase of the longitudinal transport, reaching a value of 2200 cm²V⁻¹s⁻¹ at high strain. We confirm these findings experimentally using four-point bending measurements, establishing the reliability of our first-principles calculations. For holes, we find that the transport is almost unaffected by strain up to 0.3% uniaxial tensile strain and then rises significantly, more than doubling at 2% strain. Our findings open new perspectives to boost the mobility by applying a stress in the [100] direction. This is particularly interesting for holes for which shear strain was thought for a long time to be the only way to enhance the mobility.

Interpreting Raman and IR vibrational spectra in complex organic molecules lacking symmetries poses a formidable challenge. In this study, we propose an innovative approach for simulating vibrational spectra and attributing observed peaks to molecular motions, even when highly anharmonic, without the need for computationally expensive ab initio calculations. Our approach stems from the time-dependent stochastic self-consistent harmonic approximation to capture quantum nuclear fluctuations in atom dynamics while describing interatomic interaction through state-of-the-art reactive machine-learning force fields. Finally, we employ an isotropic charge model and a bond capacitor model trained on ab initio data to predict the intensity of IR and Raman signals.

Dielectric properties of the hydrogen-bonded ferroelectric crystal KH₂PO₄ (KDP) differ significantly from those of KD₂PO₄ (DKDP). It is well established that deuteration affects the interplay of hydrogen-bond switches and heavy ion displacements that underlie the emergence of macroscopic polarization, but a detailed microscopic model is missing. We show that all-atom path integral molecular dynamics simulations can predict the isotope effects, revealing the microscopic mechanism that differentiates KDP and DKDP. Proton tunneling generates phosphate configurations that do not contribute to the polarization. At low temperatures, these quantum dipolar defects are substantial in KDP but negligible in DKDP. These intrinsic defects explain why KDP has lower spontaneous polarization and transition entropy than DKDP. The prominent role of quantum fluctuations in KDP is related to the unusual strength of the hydrogen bonds and should be equally important in other crystals of the KDP family, which exhibit similar isotope effects.