Evolutionary Search and Theoretical Study of Silicene Grain Boundaries' Mechanical Properties.

Defects such as grain boundaries (GBs) are almost inevitable during the synthesis process of 2D materials. To take advantage of the fascinating properties of 2D materials, understanding the nature and impact of various GB structures on pristine 2D sheets is crucial. In this work, using an evolutionary algorithm search, we predict a wide variety of silicene GB structures with very different atomic structures compared with those found in graphene or hexagonal boron-nitride. Twenty-one GBs with the lowest energy were validated by density functional theory (DFT), a majority of which were previously unreported to our best knowledge. Based on the diversity of the GB predictions, we found that the formation energy and mechanical properties can be dramatically altered by adatom positions within a GB and certain types of atomic structures, such as four-atom rings. To study the mechanical behavior of these GBs, we apply strain to the GB structures stepwise and use DFT calculations to investigate the mechanical properties of 9 representative structures. It is observed that GB structures based on pentagon-heptagon pairs are likely to have similar or higher in-plane stiffness and strength compared to the zigzag orientation of pristine silicene. However, an adatom located at the hollow site of a heptagon ring can significantly deteriorate the mechanical strength. For all of the structures, the in-plane stiffness and strength were found to decrease with increasing formation energy. For the failure behavior of GB structures, it was found that GB structures based on pentagon-heptagon pairs have failure behavior similar to that of graphene. We also found that the GB structures with atoms positioned outside of the 2D plane tend to experience phase transitions before failure. Utilizing the evolutionary algorithm, we locate diverse silicene GBs and obtain useful information about their mechanical properties.

Accelerated Sequence Design of Star Block Copolymers: An Unbiased Exploration Strategy via Fusion of Molecular Dynamics Simulations and Machine Learning.

Star block copolymers (s-BCPs) have potential applications as novel surfactants or amphiphiles for emulsification, compatibilization, chemical transformations, and separations. s-BCPs have chain architectures where three or more linear diblock copolymer arms comprised of two chemically distinct linear polymers, e.g., solvophobic and solvophilic chains, are covalently joined at one point. The chemical composition of each of the subunit polymer chains comprising the arms, their molecular weights, and the number of arms can be varied to tailor the surface and interfacial activity of these architecturally unique molecules. This makes identification of the optimal s-BCP design nontrivial as the total number of plausible s-BCP architectures is experimentally or computationally intractable. In this work, we use molecular dynamics (MD) simulations coupled with a reinforcement learning-based Monte Carlo tree search (MCTS) to identify s-BCP designs that minimize the interfacial tension between polar and nonpolar solvents. We first validate the MCTS approach for the design of small- and medium-sized s-BCPs and then use it to efficiently identify sequences of copolymer blocks for large-sized s-BCPs. The structural origins of interfacial tension in these systems are also identified by using the configurations obtained from MD simulations. Chemical insights into the arrangement of copolymer blocks that promote lower interfacial tension were mined using machine learning (ML) techniques. Overall, this work provides an efficient approach to solve design problems via fusion of simulations and ML and provides important groundwork for future experimental investigation of s-BCPs for various applications.

Active and Transfer Learning of High-Dimensional Neural Network Potentials for Transition Metals.

Classical molecular dynamics (MD) simulations represent a very popular and powerful tool for materials modeling and design. The predictive power of MD hinges on the ability of the interatomic potential to capture the underlying physics and chemistry. There have been decades of seminal work on developing interatomic potentials, albeit with a focus predominantly on capturing the properties of bulk materials. Such physics-based models, while extensively deployed for predicting the dynamics and properties of nanoscale systems over the past two decades, tend to perform poorly in predicting nanoscale potential energy surfaces (PESs) when compared to high-fidelity first-principles calculations. These limitations stem from the lack of flexibility in such models, which rely on a predefined functional form. Machine learning (ML) models and approaches have emerged as a viable alternative to capture the diverse size-dependent cluster geometries, nanoscale dynamics, and the complex nanoscale PESs, without sacrificing the bulk properties. Here, we introduce an ML workflow that combines transfer and active learning strategies to develop high-dimensional neural networks (NNs) for capturing the cluster and bulk properties for several different transition metals with applications in catalysis, microelectronics, and energy storage, to name a few. Our NN first learns the bulk PES from the high-quality physics-based models in literature and subsequently augments this learning via retraining with a higher-fidelity first-principles training data set to concurrently capture both the nanoscale and bulk PES. Our workflow departs from status-quo in its ability to learn from a sparsely sampled data set that nonetheless covers a diverse range of cluster configurations from near-equilibrium to highly nonequilibrium as well as learning strategies that iteratively improve the fingerprinting depending on model fidelity. All the developed models are rigorously tested against an extensive first-principles data set of energies and forces of cluster configurations as well as several properties of bulk configurations for 10 different transition metals. Our approach is material agnostic and provides a methodology to transfer and build upon the learnings from decades of seminal work in molecular simulations on to a new generation of ML-trained potentials to accelerate materials discovery and design.

Subnanometer Scale Mapping of Hydrogen Doping in Vanadium Dioxide.

Hydrogen donor doping of correlated electron systems such as vanadium dioxide (VO2) profoundly modifies the ground state properties. The electrical behavior of HxVO2 is strongly dependent on the hydrogen concentration; hence, atomic scale control of the doping process is necessary. It is however a nontrivial problem to quantitatively probe the hydrogen distribution in a solid matrix. As hydrogen transfers its sole electron to the material, the ionization mechanism is suppressed. In this study, a methodology mapping the doping distribution at subnanometer length scale is demonstrated across a HxVO2 thin film focusing on the oxygen-hydrogen bonds using electron energy loss spectroscopy (EELS) coupled with first-principles EELS calculations. The hydrogen distribution was revealed to be nonuniform along the growth direction and between different VO2 grains, calling for intricate hydrogenation mechanisms. Our study points to a powerful approach to quantitatively map dopant distribution in quantum materials relevant to energy and information sciences.

Infrared Nanoimaging of Hydrogenated Perovskite Nickelate Memristive Devices.

Solid-state devices made from correlated oxides, such as perovskite nickelates, are promising for neuromorphic computing by mimicking biological synaptic function. However, comprehending dopant action at the nanoscale poses a formidable challenge to understanding the elementary mechanisms involved. Here, we perform operando infrared nanoimaging of hydrogen-doped correlated perovskite, neodymium nickel oxide (H-NdNiO3, H-NNO), devices and reveal how an applied field perturbs dopant distribution at the nanoscale. This perturbation leads to stripe phases of varying conductivity perpendicular to the applied field, which define the macroscale electrical characteristics of the devices. Hyperspectral nano-FTIR imaging in conjunction with density functional theory calculations unveils a real-space map of multiple vibrational states of H-NNO associated with OH stretching modes and their dependence on the dopant concentration. Moreover, the localization of excess charges induces an out-of-plane lattice expansion in NNO which was confirmed by in situ X-ray diffraction and creates a strain that acts as a barrier against further diffusion. Our results and the techniques presented here hold great potential for the rapidly growing field of memristors and neuromorphic devices wherein nanoscale ion motion is fundamentally responsible for function.

Adhesion of impure ice on surfaces.

The undesirable buildup of ice can compromise the operational safety of ships in the Arctic to high-flying airplanes, thereby having a detrimental impact on modern life in cold climates. The obstinately strong adhesion between ice and most functional surfaces makes ice removal an energetically expensive and dangerous affair. Hence, over the past few decades, substantial efforts have been directed toward the development of passive ice-shedding surfaces. Conventionally, such research on ice adhesion has almost always been based on ice solidified from pure water. However, in all practical situations, freezing water has dissolved contaminants; ice adhesion studies of which have remained elusive thus far. Here, we cast light on the fundamental role played by various impurities (salt, surfactant, and solvent) commonly found in natural water bodies on the adhesion of ice on common structural materials. We elucidate how varying freezing temperature & contaminant concentration can significantly alter the resultant ice adhesion strength making it either super-slippery or fiercely adherent. The entrapment of impurities in ice changes with the rate of freezing and ensuing adhesion strength increases as the cooling temperature decreases. We discuss the possible role played by the in situ generated solute enriched liquid layer and the nanometric water-like disordered ice layer sandwiched between ice and the substrate behind these observations. Our work provides useful insights into the elementary nature of impure water-to-ice transformation and contributes to the knowledge base of various natural phenomena and rational design of a broad spectrum of anti-icing technologies for transportation, infrastructure, and energy systems.

Correction: Accelerating copolymer inverse design using monte carlo tree search.

Correction for 'Accelerating copolymer inverse design using monte carlo tree search' by Tarak K. Patra et al., Nanoscale, 2020, 12, 23653-23662, https://doi.org/10.1039/D0NR06091G.

Reinforcement learning based hybrid bond-order coarse-grained interatomic potentials for exploring mesoscale aggregation in liquid-liquid mixtures.

Exploring mesoscopic physical phenomena has always been a challenge for brute-force all-atom molecular dynamics simulations. Although recent advances in computing hardware have improved the accessible length scales, reaching mesoscopic timescales is still a significant bottleneck. Coarse-graining of all-atom models allows robust investigation of mesoscale physics with a reduced spatial and temporal resolution but preserves desired structural features of molecules, unlike continuum-based methods. Here, we present a hybrid bond-order coarse-grained forcefield (HyCG) for modeling mesoscale aggregation phenomena in liquid-liquid mixtures. The intuitive hybrid functional form of the potential offers interpretability to our model, unlike many machine learning based interatomic potentials. We parameterize the potential with the continuous action Monte Carlo Tree Search (cMCTS) algorithm, a reinforcement learning (RL) based global optimizing scheme, using training data from all-atom simulations. The resulting RL-HyCG correctly describes mesoscale critical fluctuations in binary liquid-liquid extraction systems. cMCTS, the RL algorithm, accurately captures the mean behavior of various geometrical properties of the molecule of interest, which were excluded from the training set. The developed potential model along with the RL-based training workflow could be applied to explore a variety of other mesoscale physical phenomena that are typically inaccessible to all-atom molecular dynamics simulations.

Graph Neural Network Guided Evolutionary Search of Grain Boundaries in 2D Materials.

Grain boundaries (GBs) in two-dimensional (2D) materials are known to dramatically impact material properties ranging from the physical, chemical, mechanical, electronic, and optical, to name a few. Predicting a range of physically realistic GB structures for 2D materials is critical to exercising control over their properties. This, however, is nontrivial given the vast structural and configurational (defect) search space between lateral 2D sheets with varying misfits. Here, in a departure from traditional evolutionary search methods, we introduce a workflow that combines the Graph Neural Network (GNN) and an evolutionary algorithm for the discovery and design of novel 2D lateral interfaces. We use a representative 2D material, blue phosphorene (BP), and identify 2D GB structures to test the efficacy of our GNN model. The GNN was trained with a computationally inexpensive machine learning bond order potential (Tersoff formalism) and density functional theory (DFT). Systematic downsampling of the training data sets indicates that our model can predict structural energy under 0.5% mean absolute error with sparse (<2000) DFT generated energy labels for training. We further couple the GNN model with a multiobjective genetic algorithm (MOGA) and demonstrate strong accuracy in the ability of the GNN to predict GBs. Our method is generalizable, is material agnostic, and is anticipated to accelerate the discovery of 2D GB structures.

Selective area doping for Mott neuromorphic electronics.

The cointegration of artificial neuronal and synaptic devices with homotypic materials and structures can greatly simplify the fabrication of neuromorphic hardware. We demonstrate experimental realization of vanadium dioxide (VO2) artificial neurons and synapses on the same substrate through selective area carrier doping. By locally configuring pairs of catalytic and inert electrodes that enable nanoscale control over carrier density, volatility or nonvolatility can be appropriately assigned to each two-terminal Mott memory device per lithographic design, and both neuron- and synapse-like devices are successfully integrated on a single chip. Feedforward excitation and inhibition neural motifs are demonstrated at hardware level, followed by simulation of network-level handwritten digit and fashion product recognition tasks with experimental characteristics. Spatially selective electron doping opens up previously unidentified avenues for integration of emerging correlated semiconductors in electronic device technologies.

Complex Oxides for Brain-Inspired Computing: A Review.

The fields of brain-inspired computing, robotics, and, more broadly, artificial intelligence (AI) seek to implement knowledge gleaned from the natural world into human-designed electronics and machines. In this review, the opportunities presented by complex oxides, a class of electronic ceramic materials whose properties can be elegantly tuned by doping, electron interactions, and a variety of external stimuli near room temperature, are discussed. The review begins with a discussion of natural intelligence at the elementary level in the nervous system, followed by collective intelligence and learning at the animal colony level mediated by social interactions. An important aspect highlighted is the vast spatial and temporal scales involved in learning and memory. The focus then turns to collective phenomena, such as metal-to-insulator transitions (MITs), ferroelectricity, and related examples, to highlight recent demonstrations of artificial neurons, synapses, and circuits and their learning. First-principles theoretical treatments of the electronic structure, and in situ synchrotron spectroscopy of operating devices are then discussed. The implementation of the experimental characteristics into neural networks and algorithm design is then revewed. Finally, outstanding materials challenges that require a microscopic understanding of the physical mechanisms, which will be essential for advancing the frontiers of neuromorphic computing, are highlighted.

Atomic Cross-Talk at the Interface: Enhanced Lubricity and Wear and Corrosion Resistance in Sub 2 nm Hybrid Overcoats via Strengthened Interface Chemistry.

Friction, wear, and corrosion remain the major causes of premature failure of diverse systems including hard-disk drives (HDDs). To enhance the areal density of HDDs beyond 1 Tb/in2, the necessary low friction and high wear and corrosion resistance characteristics with sub 2 nm overcoats remain unachievable. Here we demonstrate that atom cross-talk not only manipulates the interface chemistry but also strengthens the tribological and corrosion properties of sub 2 nm overcoats. High-affinity (HA) atomically thin (∼0.4 nm) interlayers (ATIs, XHA), namely Ti, Si, and SiNx, are sandwiched between the hard-disk media and 1.5 nm thick carbon (C) overlayer to develop interface-enhanced sub 2 nm hybrid overcoats that consistently outperform a thicker conventional commercial overcoat (≥2.7 nm), with the C/SiNx bilayer overcoat bettering all other <2 nm thick overcoats. These hybrid overcoats can enable the development of futuristic 2-4 Tb/in2 areal density HDDs and can transform various moving-mechanical-system based technologies.

Protein-Ligand Binding Free-Energy Calculations with ARROW─A Purely First-Principles Parameterized Polarizable Force Field.

Protein-ligand binding free-energy calculations using molecular dynamics (MD) simulations have emerged as a powerful tool for in silico drug design. Here, we present results obtained with the ARROW force field (FF)─a multipolar polarizable and physics-based model with all parameters fitted entirely to high-level ab initio quantum mechanical (QM) calculations. ARROW has already proven its ability to determine solvation free energy of arbitrary neutral compounds with unprecedented accuracy. The ARROW FF parameterization is now extended to include coverage of all amino acids including charged groups, allowing molecular simulations of a series of protein-ligand systems and prediction of their relative binding free energies. We ensure adequate sampling by applying a novel technique that is based on coupling the Hamiltonian Replica exchange (HREX) with a conformation reservoir generated via potential softening and nonequilibrium MD. ARROW provides predictions with near chemical accuracy (mean absolute error of ∼0.5 kcal/mol) for two of the three protein systems studied here (MCL1 and Thrombin). The third protein system (CDK2) reveals the difficulty in accurately describing dimer interaction energies involving polar and charged species. Overall, for all of the three protein systems studied here, ARROW FF predicts relative binding free energies of ligands with a similar accuracy level as leading nonpolarizable force fields.

Stability and Metastability of Liquid Water in a Machine-Learned Coarse-Grained Model with Short-Range Interactions.

Coarse-grained water models are ∼100 times more efficient than all-atom models, enabling simulations of supercooled water and crystallization. The machine-learned monatomic model ML-BOP reproduces the experimental equation of state (EOS) and ice-liquid thermodynamics at 0.1 MPa on par with the all-atom TIP4P/2005 and TIP4P/Ice models. These all-atom models were parametrized using high-pressure experimental data and are either accurate for water's EOS (TIP4P/2005) or ice-liquid equilibrium (TIP4P/Ice). ML-BOP was parametrized from temperature-dependent ice and liquid experimental densities and melting data at 0.1 MPa; its only pressure training is from compression of TIP4P/2005 ice at 0 K. Here we investigate whether ML-BOP replicates the experimental EOS and ice-water thermodynamics along all pressures of ice I. We find that ML-BOP reproduces the temperature, enthalpy, entropy, and volume of melting of hexagonal ice up to 400 MPa and the EOS of water along the melting line with an accuracy that rivals that of both TIP4P/2005 and TIP4P/Ice. We interpret that the accuracy of ML-BOP originates from its ability to capture the shift between compact and open local structures to changes in pressure and temperature. ML-BOP reproduces the sharpening of the tetrahedral peak of the pair distribution function of water upon supercooling, and its pressure dependence. We characterize the region of metastability of liquid ML-BOP with respect to crystallization and cavitation. The accessibility of ice crystallization to simulations of ML-BOP, together with its accurate representation of the thermodynamics of water, makes it promising for investigating the interplay between anomalies, glass transition, and crystallization under conditions challenging to access through experiments.

Machine learning overcomes human bias in the discovery of self-assembling peptides.

Peptide materials have a wide array of functions, from tissue engineering and surface coatings to catalysis and sensing. Tuning the sequence of amino acids that comprise the peptide modulates peptide functionality, but a small increase in sequence length leads to a dramatic increase in the number of peptide candidates. Traditionally, peptide design is guided by human expertise and intuition and typically yields fewer than ten peptides per study, but these approaches are not easily scalable and are susceptible to human bias. Here we introduce a machine learning workflow-AI-expert-that combines Monte Carlo tree search and random forest with molecular dynamics simulations to develop a fully autonomous computational search engine to discover peptide sequences with high potential for self-assembly. We demonstrate the efficacy of the AI-expert to efficiently search large spaces of tripeptides and pentapeptides. The predictability of AI-expert performs on par or better than our human experts and suggests several non-intuitive sequences with high self-assembly propensity, outlining its potential to overcome human bias and accelerate peptide discovery.

Efficient Probabilistic Computing with Stochastic Perovskite Nickelates.

Probabilistic computing has emerged as a viable approach to solve hard optimization problems. Devices with inherent stochasticity can greatly simplify their implementation in electronic hardware. Here, we demonstrate intrinsic stochastic resistance switching controlled via electric fields in perovskite nickelates doped with hydrogen. The ability of hydrogen ions to reside in various metastable configurations in the lattice leads to a distribution of transport gaps. With experimentally characterized p-bits, a shared-synapse p-bit architecture demonstrates highly parallelized and energy-efficient solutions to optimization problems such as integer factorization and Boolean satisfiability. The results introduce perovskite nickelates as scalable potential candidates for probabilistic computing and showcase the potential of light-element dopants in next-generation correlated semiconductors.

Controlled Formation of Conduction Channels in Memristive Devices Observed by X-ray Multimodal Imaging.

Neuromorphic computing provides a means for achieving faster and more energy efficient computations than conventional digital computers for artificial intelligence (AI). However, its current accuracy is generally less than the dominant software-based AI. The key to improving accuracy is to reduce the intrinsic randomness of memristive devices, emulating synapses in the brain for neuromorphic computing. Here using a planar device as a model system, the controlled formation of conduction channels is achieved with high oxygen vacancy concentrations through the design of sharp protrusions in the electrode gap, as observed by X-ray multimodal imaging of both oxygen stoichiometry and crystallinity. Classical molecular dynamics simulations confirm that the controlled formation of conduction channels arises from confinement of the electric field, yielding a reproducible spatial distribution of oxygen vacancies across switching cycles. This work demonstrates an effective route to control the otherwise random electroforming process by electrode design, facilitating the development of more accurate memristive devices for neuromorphic computing.

Removal of Non-Specifically Bound Proteins Using Rayleigh Waves Generated on ST-Quartz Substrates.

Label-free biosensors are plagued by the issue of non-specific protein binding which negatively affects sensing parameters such as sensitivity, selectivity, and limit-of-detection. In the current work, we explore the possibility of using the Rayleigh waves in ST-Quartz devices to efficiently remove non-specifically bound proteins via acoustic streaming. A coupled-field finite element (FE) fluid structure interaction (FSI) model of a surface acoustic wave (SAW) device based on ST-Quartz substrate in contact with a liquid loading was first used to predict trends in forces related to SAW-induced acoustic streaming. Based on model predictions, it is found that the computed SAW body force is sufficient to overcome adhesive forces between particles and a surface while lift and drag forces prevent reattachment for a range of SAW frequencies. We further performed experiments to validate the model predictions and observe that the excitation of Rayleigh SAWs removed non-specifically bound (NSB) antigens and antibodies from sensing and non-sensing regions, while rinsing and blocking agents were ineffective. An amplified RF signal applied to the device input disrupted the specific interactions between antigens and their capture antibody as well. ST-quartz allows propagation of Rayleigh and leaky SH-SAW waves in orthogonal directions. Thus, the results reported here could allow integration of three important biosensor functions on a single chip, i.e., removal of non-specific binding, mixing, and sensing in the liquid phase.

Unusual High Hardness and Load-Dependent Mechanical Characteristics of Hydrogenated Carbon-Nitrogen Hybrid Films.

Mechanical components are exposed to a rigorous environment in a number of applications including engineering, aerospace, and automobiles. Thus, their service lifetime and reliability are always on the verge of risk. Protective coatings with high hardness are required to enhance their service lifetime and minimize the replacement cost and waste burden. Hydrogenated amorphous carbon including nitrogen-incorporated films, that are commonly deposited by plasma-enhanced chemical vapor deposition, are widely used for commercial protective coating applications. However, their mechanical hardness still falls into the moderate hard regime. This needs to be substantially enhanced for advanced applications. Here, we report the synthesis of very hard nanostructured hydrogenated carbon-nitrogen hybrid (n-C:H:N) films. The optimized n-C:H:N film displays a hardness of about 36 GPa, elastic modulus of 360 GPa, and reasonably good elastic recovery (ER) of 62.7%. The mechanical properties of n-C:H:N films are further tailored when nitrogen pressure is tuned during the growth. The realized remarkably improved mechanical properties are correlated with the films' structural properties and experimental growth conditions. We also conducted density functional theory calculations that show the trend for the elastic modulus of the amorphous carbon films with varying nitrogen concentrations matches well with experimentally measured values. Finally, we probed load-dependent mechanical properties of n-C:H:N films and found an anomalous behavior; some of the mechanical parameters, for instance, ER, reveal an irregular trend with indentation load, which we explain in the framework of the film-substrate composite concept. Overall, this work uncovers many unknown and exciting mechanical phenomena that could pave the way for new technological developments.