npj Computational Materials最新文献

We develop an automated high-throughput workflow for calculating lattice dynamical properties from first principles including those dictated by anharmonicity. The pipeline automatically computes interatomic force constants (IFCs) up to 4th order from perturbed training supercells, and uses the IFCs to calculate lattice thermal conductivity, coefficient of thermal expansion, and vibrational free energy and entropy. It performs phonon renormalization for dynamically unstable compounds to obtain real effective phonon spectra at finite temperatures and calculates the associated free energy corrections. The methods and parameters are chosen to balance computational efficiency and result accuracy, assessed through convergence testing and comparisons with experimental measurements. Deployment of this workflow at a large scale would facilitate materials discovery efforts toward functionalities including thermoelectrics, contact materials, ferroelectrics, aerospace components, as well as general phase diagram construction.

Superalloys are indispensable materials for the fabrication of high-temperature components in aircraft engines. The discovery of a novel class of γ/γ′ Co-Al-W alloys has ignited a surge of interest in Co-based superalloys, with the aspiration to transcend the inherent constraints of their Ni-based counterparts. However, the conventional methodologies utilized in the design and advancement of new γ/γ′ Co-based superalloys are frequently characterized by their laborious and resource-intensive nature. In this study, we employed a coupled Density Functional Theory (DFT) and machine learning (ML) approach to predict and analyze the stability of the crucial γ′ phase, which is instrumental in expediting the discovery of γ/γ′ Co-based alloys. A dataset comprised of thousands of reliable formation (H_f) and decomposition (H_d) energies was obtained through high-throughput DFT calculations. Through regression model selection and feature engineering, our trained Random Forest (RF) model achieved prediction accuracies of 98.07% for H_f and 97.05% for H_d. Utilizing the well-trained RF model, we predicted the energies of over 150,000 ternary and quaternary γ′ phases within the Co-Ni-Fe-Cr-Al-W-Ti-Ta-V-Mo-Nb system. The energy analyses revealed that the presence of Ni, Nb, Ta, Ti, and V significantly reduced the H_f and the H_d of γ′, while Mo and W deteriorate the stability by increasing both energy values. Interestingly, although Al reduces the H_f, it increases H_d, thereby adversely affecting the stability of γ′. Applying domain-specific screening based on our knowledge, we identified 1049 out of >150,000 compositions likely to form stable γ′ phases, predominantly distributed across 11 Al-containing systems and 25 Al-free systems. Combining the analysis of CALPHAD method, we experimentally synthesized two new Co-based alloys with γ/γ′ dual-phase microstructures, corroborating the reliability of our theoretical prediction model.

Lightweight refractory high-entropy alloys (LW-RHEAs) hold significant potential in the fields of aviation, aerospace, and nuclear energy due to their low density, high strength, high hardness, and corrosion resistance. However, the enormous composition space has severely hindered the development of novel LW-RHEAs with excellent comprehensive performance. In this paper, an machine learning (ML)-based alloy design strategy combined with a multi-objective optimization method was proposed and applied for a rational design of Al-Nb-Ti-V-Zr-Cr-Mo-Hf LW-RHEAs. The quantitative relation of “composition-structure-property” was first established by ML modeling. Then, feature analysis reveals that Cr content greater than 12 at.% is a key criterion for alloys with high corrosion resistance. The phase structure, density, melting point, hardness and corrosion resistance of the alloys were screened layer by layer, and finally, three LW-RHEAs with superb hard and corrosion resistance were successfully designed. Key experimental validation indicates that three target alloys have densities around 6.5 g/cm³, and all alloys are disordered bcc_A2 single-phase with the highest hardness of 593 HV and the largest pitting potential of 2.5 V_SCE, which far exceeds all the literature reports. The successful demonstration in this paper clearly demonstrates that the present design strategy driven by the ML technique should be generally applicable to other RHEA systems.

Quantum algorithms are emerging tools in the design of functional materials due to their powerful solution space search capability. How to balance the high price of quantum computing resources and the growing computing needs has become an urgent problem to be solved. We propose a novel optimization strategy based on an active learning scheme that combines the Quantum-inspired Genetic Algorithm (QGA) with machine learning surrogate model regression. Using Random Forests as the surrogate model circumvents the time-consuming physical modeling or experiments, thereby improving the optimization efficiency. QGA, a genetic algorithm embedded with quantum mechanics, combines the advantages of quantum computing and genetic algorithms, enabling faster and more robust convergence to the optimum. Using the design of planar multilayer photonic structures for transparent radiative cooling as a testbed, we show superiority of our algorithm over the classical genetic algorithm (CGA). Additionally, we show the precision advantage of the Random Forest (RF) model as a flexible surrogate model, which relaxes the constraints on the type of surrogate model that can be used in other quantum computing optimization algorithms (e.g., quantum annealing needs Ising model as a surrogate).

To advance the development of high-strength polycrystalline metallic materials towards achieving carbon neutrality, it is essential to design materials in which the atomic level control of general grain boundaries (GGBs), which govern the material properties, is achieved. However, owing to the complex and diverse structures of GGBs, there have been no reports on interatomic potentials capable of reproducing them. This accuracy is essential for conducting molecular dynamics analyses to derive material design guidelines. In this study, we constructed a machine learning interatomic potential (MLIP) with density functional theory (DFT) accuracy to model the energy, atomic structure, and dynamics of arbitrary grain boundaries (GBs), including GGBs, in Œ±-Fe. Specifically, we employed a training dataset comprising diverse atomic structures generated based on crystal space groups. The GGB accuracy was evaluated by directly comparing with DFT calculations performed on cells cut near GBs from nano-polycrystals, and extrapolation grades of the local atomic environment based on active learning methods for the entire nano-polycrystal. Furthermore, we analyzed the GB energy and atomic structure in Œ±-Fe polycrystals through large-scale molecular dynamics analysis using the constructed MLIP. The average GB energy of Œ±-Fe polycrystals calculated by the constructed MLIP is 1.57‚ÄâJ/m², exhibiting good agreement with experimental predictions. Our findings demonstrate the methodology for constructing an MLIP capable of representing GGBs with high accuracy, thereby paving the way for materials design based on computational materials science for polycrystalline materials.

Recent advances in deep learning generative models (GMs) have created high capabilities in accessing and assessing complex high-dimensional data, allowing superior efficiency in navigating vast material configuration space in search of viable structures. Coupling such capabilities with physically significant data to construct trained models for materials discovery is crucial to moving this emerging field forward. Here, we present a universal GM for crystal structure prediction (CSP) via a conditional crystal diffusion variational autoencoder (Cond-CDVAE) approach, which is tailored to allow user-defined material and physical parameters such as composition and pressure. This model is trained on an expansive dataset containing over 670,000 local minimum structures, including a rich spectrum of high-pressure structures, along with ambient-pressure structures in Materials Project database. We demonstrate that the Cond-CDVAE model can generate physically plausible structures with high fidelity under diverse pressure conditions without necessitating local optimization, accurately predicting 59.3% of the 3547 unseen ambient-pressure experimental structures within 800 structure samplings, with the accuracy rate climbing to 83.2% for structures comprising fewer than 20 atoms per unit cell. These results meet or exceed those achieved via conventional CSP methods based on global optimization. The present findings showcase substantial potential of GMs in the realm of CSP.

Molecular dynamics (MD) is an indispensable atomistic-scale computational tool widely-used in various disciplines. In the past decades, nearly all ab initio MD and machine-learning MD have been based on the general-purpose central/graphics processing units (CPU/GPU), which are well-known to suffer from their intrinsic “memory wall” and “power wall” bottlenecks. Consequently, nowadays MD calculations with ab initio accuracy are extremely time-consuming and power-consuming, imposing serious restrictions on the MD simulation size and duration. To solve this problem, here we propose a special-purpose MD processing unit (MDPU), which could reduce MD time and power consumption by about 10³ times (10⁹ times) compared to state-of-the-art machine-learning MD (ab initio MD) based on CPU/GPU, while keeping ab initio accuracy. With significantly-enhanced performance, the proposed MDPU may pave a way for the accurate atomistic-scale analysis of large-size and/or long-duration problems which were impossible/impractical to compute before.

In metal additive manufacturing (AM), processing parameters can affect the probability of macroscopic defect formation (lack-of-fusion, keyholing, balling), which can, in turn, jeopardize the final product’s integrity. A printability map classifies regions in the processing space where an alloy can be printed with or without porosity defects. However, the creation of these printability maps is resource-intensive. Previous efforts to generate printability maps have required single-track experiments on pre-alloyed powder, limiting the utilization of these printability maps for the high-throughput design of printable alloys. We address these challenges in the case of Laser Powder Bed Fusion AM (L-PBF-AM) by introducing a fully computational, predictive approach to create printability maps for arbitrary alloys. Our framework uses physics-based thermal models and a variety of defect formation criteria. We benchmark the predictive ability of the proposed framework against literature data for the following commonly printed alloys: 316 Stainless Steel, Inconel 718, Ti-6Al-4V, AF96, and Ni-5Nb. Furthermore, we deploy the framework on NiTi-based Shape Memory Alloys (SMAs) as a case study. We scrutinize the accuracy of various sets of defect criteria and use these accuracy measurements to create an uncertainty-aware probabilistic framework capable of predicting the printability maps of arbitrary alloys. This framework has the potential to guide alloy designers to potentially easy-to-print alloys, enabling the co-design of high-performing printable alloys.

Upgrades to advanced scientific user facilities such as next-generation x-ray light sources, nanoscience centers, and neutron facilities are revolutionizing our understanding of materials across the spectrum of the physical sciences, from life sciences to microelectronics. However, these facility and instrument upgrades come with a significant increase in complexity. Driven by more exacting scientific needs, instruments and experiments become more intricate each year. This increased operational complexity makes it ever more challenging for domain scientists to design experiments that effectively leverage the capabilities of and operate on these advanced instruments. Large language models (LLMs) can perform complex information retrieval, assist in knowledge-intensive tasks across applications, and provide guidance on tool usage. Using x-ray light sources, leadership computing, and nanoscience centers as representative examples, we describe preliminary experiments with a Context-Aware Language Model for Science (CALMS) to assist scientists with instrument operations and complex experimentation. With the ability to retrieve relevant information from facility documentation, CALMS can answer simple questions on scientific capabilities and other operational procedures. With the ability to interface with software tools and experimental hardware, CALMS can conversationally operate scientific instruments. By making information more accessible and acting on user needs, LLMs could expand and diversify scientific facilities’ users and accelerate scientific output.

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.