No Arabic abstract
Finding amorphous polymers with higher thermal conductivity is important, as they are ubiquitous in heat transfer applications. With recent progress in material informatics, machine learning approaches have been increasingly adopted for finding or designing materials with desired properties. However, relatively limited effort has been put into finding thermally conductive polymers using machine learning, mainly due to the lack of polymer thermal conductivity databases with reasonable data volume. In this work, we combine high-throughput molecular dynamics (MD) simulations and machine learning to explore polymers with relatively high thermal conductivity (> 0.300 W/m-K). We first randomly select 365 polymers from the existing PolyInfo database and calculate their thermal conductivity using MD simulations. The data are then employed to train a machine learning regression model to quantify the structure-thermal conductivity relation, which is further leveraged to screen polymer candidates in the PolyInfo database with thermal conductivity > 0.300 W/m-K. 133 polymers with MD-calculated thermal conductivity above this threshold are eventually identified. Polymers with a wide range of thermal conductivity values are selected for re-calculation under different simulation conditions, and those polymers found with thermal conductivity above 0.300 W/m-K are mostly calculated to maintain values above this threshold despite fluctuation in the exact values. A classification model is also constructed, and similar results were obtained compared to the regression model in predicting polymers with thermal conductivity above or below 0.300 W/m-K. The strategy and results from this work may contribute to automating the design of polymers with high thermal conductivity.
The bulk and surface dynamics of Cu50Zr50 metallic glass were studied using classical molecular dynamics (MD) simulations. As the alloy undergoes cooling, it passes through liquid, supercooled, and glassy states. While bulk dynamics showed a marked slowing down prior to glass formation, with increasing activation energy, the slowdown in surface dynamics was relatively subtle. The surface exhibited a lower glass transition temperature than the bulk, and the dynamics preceding the transition were accurately described by a temperature-independent activation energy. Surface dynamics were much faster than bulk at a given temperature in the supercooled state, but surface and bulk dynamics were found to be very similar when compared at their respective glass transition temperatures. The manifestation of dynamical heterogeneity, as characterized by the non-Gaussian parameter and breakdown of the Stokes-Einstein equation, was also similar between bulk and surface for temperatures scaled by their respective glass transition temperatures. Individual atom motion was dominated by a cage and jump mechanism in the glassy state for both the bulk and surface. We utilize this cage and jump mechanisms to separate the activation energy for diffusion into two parts: (i) cage-breaking barrier (Q1), associated with the rearrangement of neighboring atoms to free up space and (ii) the subsequent jump barrier (Q2). It was observed that Q1 dominates Q2 for both bulk and surface diffusion, and the difference in activation energies for bulk and surface diffusion mainly arose from the differences in cage-breaking barrier Q1.
Recently, high-entropy alloys (HEAs) have attracted wide attention due to their extraordinary materials properties. A main challenge in identifying new HEAs is the lack of efficient approaches for exploring their huge compositional space. Ab initio calculations have emerged as a powerful approach that complements experiment. However, for multicomponent alloys existing approaches suffer from the chemical complexity involved. In this work we propose a method for studying HEAs computationally. Our approach is based on the application of machine-learning potentials based on ab initio data in combination with Monte Carlo simulations. The high efficiency and performance of the approach are demonstrated on the prototype bcc NbMoTaW HEA. The approach is employed to study phase stability, phase transitions, and chemical short-range order. The importance of including local relaxation effects is revealed: they significantly stabilize single-phase formation of bcc NbMoTaW down to room temperature. Finally, a so-far unknown mechanism that drives chemical order due to atomic relaxation at ambient temperatures is discovered.
We develop a Python-based open-source package to analyze the results stemming from ab initio molecular-dynamics simulations of fluids. The package is best suited for applications on natural systems, like silicate and oxide melts, water-based fluids, various supercritical fluids. The package is a collection of Python scripts that include two major libraries dealing with file formats and with crystallography. All the scripts are run at the command line. We propose a simplified format to store the atomic trajectories and relevant thermodynamic information of the simulations, which is saved in UMD files, standing for Universal Molecular Dynamics. The UMD package allows the computation of a series of structural, transport and thermodynamic properties. Starting with the pair-distribution function it defines bond lengths, builds an interatomic connectivity matrix, and eventually determines the chemical speciation. Determining the lifetime of the chemical species allows running a full statistical analysis. Then dedicated scripts compute the mean-square displacements for the atoms as well as for the chemical species. The implemented self-correlation analysis of the atomic velocities yields the diffusion coefficients and the vibrational spectrum. The same analysis applied on the stresses yields the viscosity. The package is available via the GitHub website and via its own dedicated page of the ERC IMPACT project as open-access package.
High-temperature alloy design requires a concurrent consideration of multiple mechanisms at different length scales. We propose a workflow that couples highly relevant physics into machine learning (ML) to predict properties of complex high-temperature alloys with an example of the 9-12 wt.% Cr steels yield strength. We have incorporated synthetic alloy features that capture microstructure and phase transformations into the dataset. Identified high impact features that affect yield strength of 9Cr from correlation analysis agree well with the generally accepted strengthening mechanism. As part of the verification process, the consistency of sub-datasets has been extensively evaluated with respect to temperature and then refined for the boundary conditions of trained ML models. The predicted yield strength of 9Cr steels using the ML models is in excellent agreement with experiments. The current approach introduces physically meaningful constraints in interrogating the trained ML models to predict properties of hypothetical alloys when applied to data-driven materials.
We have investigated thermal conductivity of graphene laminate films deposited on polyethylene terephthalate substrates. Two types of graphene laminate were studied - as deposited and compressed - in order to determine the physical parameters affecting the heat conduction the most. The measurements were performed using the optothermal Raman technique and a set of suspended samples with the graphene laminate thickness from 9 to 44 micrometers. The thermal conductivity of graphene laminate was found to be in the range from 40 W/mK to 90 W/mK at room temperature. It was found unexpectedly that the average size and the alignment of graphene flakes are more important parameters defining the heat conduction than the mass density of the graphene laminate. The thermal conductivity scales up linearly with the average graphene flake size in both uncompressed and compressed laminates. The compressed laminates have higher thermal conductivity for the same average flake size owing to better flake alignment. The possibility of up to 600X enhancement of the thermal conductivity of plastic materials by coating them with the thin graphene laminate films has important practical implications.