Computational Design and Manufacturing of Sustainable Materials through First-Principles and Materiomics

Engineered materials are ubiquitous throughout society and are critical to the development of modern technology, yet many current material systems are inexorably tied to widespread deterioration of ecological processes. Next-generation material systems can address goals of environmental sustainability by providing alternatives to fossil fuel-based materials and by reducing destructive extraction processes, energy costs, and accumulation of solid waste. However, development of sustainable materials faces several key challenges including investigation, processing, and architecting of new feedstocks that are often relatively mechanically weak, complex, and difficult to characterize or standardize. In this review paper, we outline a framework for examining sustainability in material systems and discuss how recent developments in modeling, machine learning, and other computational tools can aid the discovery of novel sustainable materials. We consider these through the lens of materiomics, an approach that considers material systems holistically by incorporating perspectives of all relevant scales, beginning with first-principles approaches and extending through the macroscale to consider sustainable material design from the bottom-up. We follow with an examination of how computational methods are currently applied to select examples of sustainable material development, with particular emphasis on bioinspired and biobased materials, and conclude with perspectives on opportunities and open challenges.

General Graph Neural Network-Based Model To Accurately Predict Cocrystal Density and Insight from Data Quality and Feature Representation

Cocrystal engineering as an effective way to modify solid-state properties has inspired great interest from diverse material fields while cocrystal density is an important property closely correlated with the material function. In order to accurately predict the cocrystal density, we develop a graph neural network (GNN)-based deep learning framework by considering three key factors of machine learning (data quality, feature presentation, and model architecture). The result shows that different stoichiometric ratios of molecules in cocrystals can significantly influence the prediction performances, highlighting the importance of data quality. In addition, the feature complementary is not suitable for augmenting the molecular graph representation in the cocrystal density prediction, suggesting that the complementary strategy needs to consider whether extra features can sufficiently supplement the lacked information in the original representation. Based on these results, 4144 cocrystals with 1:1 stoichiometry ratio are selected as the dataset, supplemented by the data augmentation of exchanging a pair of coformers. The molecular graph is determined to learn feature representation to train the GNN-based model. Global attention is introduced to further optimize the feature space and identify important atoms to realize the interpretability of the model. Benefited from the advantages, our model significantly outperforms three competitive models and exhibits high prediction accuracy for unseen cocrystals, showcasing its robustness and generality. Overall, our work not only provides a general cocrystal density prediction tool for experimental investigations but also provides useful guidelines for the machine learning application. All source codes are freely available at https://github.com/Xiao-Gua00/CCPGraph.

Molecular Modeling in Anion Exchange Membrane Research: A Brief Review of Recent Applications

Anion Exchange Membrane (AEM) fuel cells have attracted growing interest, due to their encouraging advantages, including high power density and relatively low cost. AEM is a polymer matrix, which conducts hydroxide (OH−) ions, prevents physical contact of electrodes, and has positively charged head groups (mainly quaternary ammonium (QA) groups), covalently bound to the polymer backbone. The chemical instability of the quaternary ammonium (QA)-based head groups, at alkaline pH and elevated temperature, is a significant threshold in AEMFC technology. This review work aims to introduce recent studies on the chemical stability of various QA-based head groups and transportation of OH− ions in AEMFC, via modeling and simulation techniques, at different scales. It starts by introducing the fundamental theories behind AEM-based fuel-cell technology. In the main body of this review, we present selected computational studies that deal with the effects of various parameters on AEMs, via a variety of multi-length and multi-time-scale modeling and simulation methods. Such methods include electronic structure calculations via the quantum Density Functional Theory (DFT), ab initio, classical all-atom Molecular Dynamics (MD) simulations, and coarse-grained MD simulations. The explored processing and structural parameters include temperature, hydration levels, several QA-based head groups, various types of QA-based head groups and backbones, etc. Nowadays, many methods and software packages for molecular and materials modeling are available. Applications of such methods may help to understand the transportation mechanisms of OH− ions, the chemical stability of functional head groups, and many other relevant properties, leading to a performance-based molecular and structure design as well as, ultimately, improved AEM-based fuel cell performances. This contribution aims to introduce those molecular modeling methods and their recent applications to the AEM-based fuel cells research community.

Significance of density functional theory (DFT) calculations for electrocatalysis of N2 and CO2 reduction reactions

Density functional theory (DFT) based computational methods have shown great significance in developing high-performance electrocatalysts. In this perspective, we briefly summarized the state-of-the-art research progress of electrocatalysts for the nitrogen reduction reaction (NRR) and CO₂ reduction reaction (CO₂RR), which are important processes for the conversion of common molecules into value-added products. With the help of DFT calculations, various modulation strategies are employed to improve the catalytic activity and performance of NRR and CO₂RR electrocatalysts. DFT calculations are performed to confirm the surface catalytic sites, evaluate the catalytic activity, reveal the possible reaction mechanisms, and design novel structures with high catalytic performance. By discussing the currently applied computational methods and conditions during the calculations, we outlined our concerns on the prospects and future challenges of DFT calculations in electrocatalysis studies.

Charge-Transfer Landscape Manifesting the Structure-Rate Relationship in the Condensed Phase Via Machine Learning

In this work, we develop a machine learning (ML) strategy to map the molecular structure to condensed phase charge-transfer (CT) properties including CT rate constants, energy levels, electronic couplings, energy gaps, reorganization energies, and reaction free energies which are called CT fingerprints. The CT fingerprints of selected landmark structures covering the conformation space of an organic photovoltaic molecule dissolved in an explicit solvent are computed and used to train ML models using kernel ridge regression. The ML models show high predictive power with R² > 0.97 and both mean absolute error and root-mean-square error within chemical accuracy. The CT landscape for millions of molecular dynamics sampled structures is thus constructed, which allows for instant prediction of CT rate properties, given any conformation of the molecule. We demonstrate some immediate utilities of the CT landscape such as calculating the ensemble-averaged CT rate constant and interpreting the effects of molecular structural features on the CT rate. The unprecedented CT landscape will be useful for investigating real-time CT dynamics in nanoscale- and mesoscale-condensed phase systems and for the optimal fabrication design for homogeneous and heterogeneous optoelectronic devices.

Simultaneous Optimization of Donor/Acceptor Pairs and Device Specifications for Nonfullerene Organic Solar Cells Using a QSPR Model with Morphological Descriptors

Optimally efficient organic solar cells require not only a careful choice of new donor (D) and/or acceptor (A) molecules but also the fine-tuning of experimental fabrication conditions for organic solar cells (OSCs). Herein, a new framework for simultaneously optimizing D/A molecule pairs and device specifications of OSCs is proposed, through a quantitative structure-property relationship (QSPR) model built by machine learning. Combining the device bulk properties with structural and electronic properties, the built QSPR model achieved unprecedentedly high accuracy and consistency. Additionally, a large chemical space of 1 942 785 D/A pairs is explored to find potential synergistic ones. Favorable device bulk properties such as the root-mean-square of surfaces roughness for D/A blends and the D/A weight ratio are further screened by grid search methods. Overall, this study indicates that the simultaneous optimization of D/A molecule pairs and device specifications by theoretical calculations can accelerate the improvement of OSC efficiencies.