Search

Tungsten (W) is a material of choice for the divertor material due to its high melting temperature, thermal conductivity, and sputtering threshold. However, W has a very high brittle-to-ductile transition temperature, and at fusion reactor temperatures (≥1000 K), it may undergo recrystallization and grain growth. Dispersion-strengthening W with zirconium carbide (ZrC) can improve ductility and limit grain growth, but much of the effects of the dispersoids on microstructural evolution and thermomechanical properties at high temperatures are still unknown. We present a machine learned Spectral Neighbor Analysis Potential for W-ZrC that can now be used to study these materials. In order to construct a potential suitable for large-scale atomistic simulations at fusion reactor temperatures, it is necessary to train on ab initio data generated for a diverse set of structures, chemical environments, and temperatures. Further accuracy and stability tests of the potential were achieved using objective functions for both material properties and high temperature stability. Validation of lattice parameters, surface energies, bulk moduli, and thermal expansion is confirmed on the optimized potential. Tensile tests of W/ZrC bicrystals show that although the W(110)-ZrC(111) C-terminated bicrystal has the highest ultimate tensile strength (UTS) at room temperature, observed strength decreases with increasing temperature. At 2500 K, the terminating C layer diffuses into the W, resulting in a weaker W-Zr interface. Meanwhile, the W(110)-ZrC(111) Zr-terminated bicrystal has the highest UTS at 2500 K.

We present a scale-bridging approach based on a multi-fidelity (MF) machine-learning (ML) framework leveraging Gaussian processes (GP) to fuse atomistic computational model predictions across multiple levels of fidelity. Through the posterior variance of the MFGP, our framework naturally enables uncertainty quantification, providing estimates of confidence in the predictions. We used density functional theory as high-fidelity prediction, while a ML interatomic potential is used as low-fidelity prediction. Practical materials' design efficiency is demonstrated by reproducing the ternary composition dependence of a quantity of interest (bulk modulus) across the full aluminum-niobium-titanium ternary random alloy composition space. The MFGP is then coupled to a Bayesian optimization procedure, and the computational efficiency of this approach is demonstrated by performing an on-the-fly search for the global optimum of bulk modulus in the ternary composition space. The framework presented in this manuscript is the first application of MFGP to atomistic materials simulations fusing predictions between density functional theory and classical interatomic potential calculations.

We consider a confined sheared active polar liquid crystal with a uniform orientation and study the effect of variations in the magnitude of polarization. Restricting our analysis to one-dimensional geometries, we demonstrate that with asymmetric boundary conditions, this system is characterized, macroscopically, by a linear shear stress vs. shear strain relationship that does not pass through the origin: At a zero strain rate, the fluid sustains a non-zero stress. Analytic solutions for the polarization, density, and velocity fields are derived for asymptotically large or small systems and are shown by comparison with precise numerical solutions to be good approximations for finite-size systems.

The Spectral Neighbor Analysis Potential (SNAP) is a classical interatomic potential that expresses the energy of each atom as a linear function of selected bispectrum components of the neighbor atoms. An extension of the SNAP form is proposed that includes quadratic terms in the bispectrum components. The extension is shown to provide a large increase in accuracy relative to the linear form, while incurring only a modest increase in computational cost. The mathematical structure of the quadratic SNAP form is similar to the embedded atom method (EAM), with the SNAP bispectrum components serving as counterparts to the two-body density functions in EAM. The effectiveness of the new form is demonstrated using an extensive set of training data for tantalum structures. Similar to artificial neural network potentials, the quadratic SNAP form requires substantially more training data in order to prevent overfitting. The quality of this new potential form is measured through a robust cross-validation analysis.

The lack of adequately predictive atomistic empirical models precludes meaningful simulations for many materials systems. We describe advances in the development of a hybrid, population based optimization strategy intended for the automated development of material specific interatomic potentials. We compare two strategies for parallel genetic programming and show that the Hierarchical Fair Competition algorithm produces better results in terms of transferability, despite a lower training set accuracy. We evaluate the use of hybrid local search and several fitness models using system energies and/or particle forces. We demonstrate a drastic reduction in the computation time with the use of a correlation-based fitness statistic. We show that the problem difficulty increases with the number of atoms present in the systems used for model development and demonstrate that vectorization can help to address this issue. Finally, we show that with the use of this method, we are able to "rediscover" the exact model for simple known two- and three-body interatomic potentials using only the system energies and particle forces from the supplied atomic configurations.

Three distinct forms are derived for the force virial contribution to the pressure and stress tensor of a collection of atoms interacting under periodic boundary conditions. All three forms are written in terms of forces acting on atoms, and so are valid for arbitrary many-body interatomic potentials. All three forms are mathematically equivalent. In the special case of atoms interacting with pair potentials, they reduce to previously published forms. (i) The atom-cell form is similar to the standard expression for the virial for a finite nonperiodic system, but with an explicit correction for interactions with periodic images. (ii) The atom form is particularly suited to implementation in modern molecular dynamics simulation codes using spatial decomposition parallel algorithms. (iii) The group form of the virial allows the contributions to the virial to be assigned to individual atoms.

The time evolution of species concentrations in biochemical reaction networks is often modeled using the stochastic simulation algorithm (SSA) [Gillespie, J. Phys. Chem. 81, 2340 (1977)]. The computational cost of the original SSA scaled linearly with the number of reactions in the network. Gibson and Bruck developed a logarithmic scaling version of the SSA which uses a priority queue or binary tree for more efficient reaction selection [Gibson and Bruck, J. Phys. Chem. A 104, 1876 (2000)]. More generally, this problem is one of dynamic discrete random variate generation which finds many uses in kinetic Monte Carlo and discrete event simulation. We present here a constant-time algorithm, whose cost is independent of the number of reactions, enabled by a slightly more complex underlying data structure. While applicable to kinetic Monte Carlo simulations in general, we describe the algorithm in the context of biochemical simulations and demonstrate its competitive performance on small- and medium-size networks, as well as its superior constant-time performance on very large networks, which are becoming necessary to represent the increasing complexity of biochemical data for pathways that mediate cell function.