Source author record

Danny Perez

Danny Perez appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

11works
7topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

11 published item(s)

preprint2026arXiv

Active Learning of A Crystal Plasticity Flow Rule From Discrete Dislocation Dynamics Simulations

Continuum-scale material deformation models, such as crystal plasticity, can significantly enhance their predictive accuracy by incorporating input from lower-scale (i.e., mesoscale) models. The procedure to generate and extract the relevant information is however typically complex and ad hoc, involving decision and intervention by domain experts, leading to long development times. In this study, we develop a principled approach for calibration of continuum-scale models using lower scale information by representing a crystal plasticity flow rule as a Gaussian process model. This representation allows for efficient parameter space exploration, guided by the uncertainty embedded in the model through a process known as Bayesian optimization. We demonstrate a semi-autonomous Bayesian optimization loop which instantiates discrete dislocation dynamics simulations whose initial conditions are automatically chosen to optimize the uncertainty of a model crystal plasticity flow rule. Our self-guided computational pipeline efficiently generated a dataset and corresponding model whose error, uncertainty, and physical feature sensitivities were validated with comparison to an independent dataset four times larger, demonstrating a valuable and efficient active learning implementation readily transferable to similar material systems.

preprint2026arXiv

LAMDA: Aiding Visual Exploration of Atomic Displacements in Molecular Dynamics Simulations

Contemporary materials science research is heavily conducted in silico, involving massive simulations of the atomic-scale evolution of materials. Cataloging basic patterns in the atomic displacements is key to understanding and predicting the evolution of physical properties. However, the combinatorial complexity of the space of possible transitions coupled with the overwhelming amount of data being produced by high-throughput simulations make such an analysis extremely challenging and time-consuming for domain experts. The development of visual analytics systems that facilitate the exploration of simulation data is an active field of research. While these systems excel in identifying temporal regions of interest, they treat each timestep of a simulation as an independent event without considering the behavior of the atomic displacements between timesteps. We address this gap by introducing LAMDA, a visual analytics system that allows domain experts to quickly and systematically explore state-to-state transitions. In LAMDA, transitions are hierarchically categorized, providing a basis for cataloging displacement behavior, as well as enabling the analysis of simulations at different resolutions, ranging from very broad qualitative classes of transitions to very narrow definitions of unit processes. LAMDA supports navigating the hierarchy of transitions, enabling scientists to visualize the commonalities between different transitions in each class in terms of invariant features characterizing local atomic environments, and LAMDA simplifies the analysis by capturing user inputs through annotations. We evaluate our system through a case study and report on findings from our domain experts.

preprint2022arXiv

Training Data Selection for Accuracy and Transferability of Interatomic Potentials

Advances in machine learning (ML) techniques have enabled the development of interatomic potentials that promise both the accuracy of first principles methods and the low-cost, linear scaling, and parallel efficiency of empirical potentials. Despite rapid progress in the last few years, ML-based potentials often struggle to achieve transferability, that is, to provide consistent accuracy across configurations that significantly differ from those used to train the model. In order to truly realize the promise of ML-based interatomic potentials, it is therefore imperative to develop systematic and scalable approaches for the generation of diverse training sets that ensure broad coverage of the space of atomic environments. This work explores a diverse-by-construction approach that leverages the optimization of the entropy of atomic descriptors to create a very large ($>2\cdot10^{5}$ configurations, $>7\cdot10^{6}$ atomic environments) training set for tungsten in an automated manner, i.e., without any human intervention. This dataset is used to train polynomial as well as multiple neural network potentials with different architectures. For comparison, a corresponding family of potentials were also trained on an expert-curated dataset for tungsten. The models trained to entropy-optimized data exhibited vastly superior transferability compared to the expert-curated models. Furthermore, while the models trained with heavy user input (i.e., domain expertise) yield the lowest errors when tested on similar configurations, out-sample predictions are dramatically more robust when the models are trained on a deliberately diverse set of training data. Herein we demonstrate the development of both accurate and transferable ML potentials using automated and data-driven approaches for generating large and diverse training sets.

preprint2021arXiv

Reaction-drift-diffusion models from master equations: application to material defects

We present a general method to produce well-conditioned continuum reaction-drift-diffusion equations directly from master equations on a discrete, periodic state space. We assume the underlying data to be kinetic Monte Carlo models (i.e., continuous-time Markov chains) produced from atomic sampling of point defects in locally periodic environments, such as perfect lattices, ordered surface structures or dislocation cores, possibly under the influence of a slowly varying external field. Our approach also applies to any discrete, periodic Markov chain. The analysis identifies a previously omitted non-equilibrium drift term, present even in the absence of external forces, which can compete in magnitude with the reaction rates, thus being essential to correctly capture the kinetics. To remove fast modes which hinder time integration, we use a generalized Bloch relation to efficiently calculate the eigenspectrum of the master equation. A well conditioned continuum equation then emerges by searching for spectral gaps in the long wavelength limit, using an established kinetic clustering algorithm (e.g., PCCA+) to define a proper reduced state space.

preprint2020arXiv

An Entropy-Maximization Approach to Automated Training Set Generation for Interatomic Potentials

Machine learning (ML)-based interatomic potentials are currently garnering a lot of attention as they strive to achieve the accuracy of electronic structure methods at the computational cost of empirical potentials. Given their generic functional forms, the transferability of these potentials is highly dependent on the quality of the training set, the generation of which is a highly labor-intensive activity. Good training sets should at once contain a very diverse set of configurations while avoiding redundancies that incur cost without providing benefits. We formalize these requirements in a local entropy maximization framework and propose an automated sampling scheme to sample from this objective function. We show that this approach generates much more diverse training sets than unbiased sampling and is competitive with hand-crafted training sets.

preprint2020arXiv

Arbitrarily accurate representation of atomistic dynamics via Markov Renewal Processes

Atomistic simulations with methods such as molecular dynamics are extremely powerful tools to understand nanoscale dynamical behavior. The resulting trajectories, by the virtue of being embedded in a high-dimensional configuration space, can however be difficult to analyze and interpret. This makes low-dimensional representations, especially in terms of discrete jump processes, extremely valuable. This simplicity however usually comes at the cost of accuracy, as tractable representations often entail simplifying assumptions that are not guaranteed to be realized in practice. In this paper, we describe a discretization scheme for continuous trajectories that enables an arbitrarily accurate representation in terms of a Markov Renewal Process over a discrete state space. The accuracy of the model converges exponentially fast as a function of a continuous parameter that has the interpretation of a local correlation time of the dynamics.

preprint2020arXiv

Automated calculation and convergence of defect transport tensors

Defect transport is a key process in materials science and catalysis, but as migration mechanisms are often too complex to enumerate a priori, calculation of transport tensors typically have no measure of convergence and require significant end user intervention. These two bottlenecks prevent high-throughput implementations essential to propagate model-form uncertainty from interatomic interactions to predictive simulations. In order to address these issues, we extend a massively parallel accelerated sampling scheme, autonomously controlled by Bayesian estimators of statewise sampling completeness, to build atomistic kinetic Monte Carlo models on a state space irreducible under exchange and space group symmetries. Focusing on isolated defects, we derive analytic expressions for defect transport tensors and provide a convergence metric by calculating the Kullback-Leiber divergence across the ensemble of diffusion processes consistent with the sampling uncertainty. The autonomy and efficacy of the method is demonstrated on surface trimers in tungsten and hexa-interstitials in magnesium oxide, both of which exhibit complex, correlated migration mechanisms.

preprint2015arXiv

Accelerating the convergence of replica exchange simulations using Gibbs sampling and adaptive temperature sets

We recently introduced a novel replica-exchange scheme in which an individual replica can sample from states encountered by other replicas at any previous time by way of a global configuration database, enabling the fast propagation of relevant states through the whole ensemble of replicas. This mechanism depends on the knowledge of global thermodynamic functions which are measured during the simulation and not coupled to the heat bath temperatures driving the individual simulations. Therefore, this setup also allows for a continuous adaptation of the temperature set. In this paper, we will review the new scheme and demonstrate its capability. The method is particularly useful for the fast and reliable estimation of the microcanonical temperature T(U) or, equivalently, of the density of states g(U) over a wide range of energies.

preprint2014arXiv

Analysis of Transition State Theory Rates upon Spatial Coarse-Graining

Spatial multiscale methods have established themselves as useful tools for extending the length scales accessible by conventional statics (i.e., zero temperature molecular dynamics). Recently, extensions of these methods, such as the finite-temperature quasicontinuum (hot-QC) or Coarse-Grained Molecular Dynamics (CGMD) methods, have allowed for multiscale molecular dynamics simulations at finite temperature. Here, we assess the quality of the long-time dynamics these methods generate by considering canonical transition rates. Specifically, we analyze the transition state theory (TST) rates in CGMD and compare them to the corresponding TST rate of the fully atomistic system. The ability of such an approach to reliably reproduce the TST rate is verified through a relative error analysis, which is then used to highlight the major contributions to the error and guide the choice of degrees of freedom. Finally, our analytical results are compared with numerical simulations for the case of a 1-D chain.

preprint2014arXiv

Diffusion and transformation kinetics of small Helium clusters in bulk Tungsten

The production of energy through nuclear fusion poses serious challenges related to the stability and performance of materials in extreme conditions. In particular, the constant bombardment of the walls of the reactor with high doses of He ions is known to lead to deleterous changes in their microstructures. These changes follow from the aggregation of He into bubbles that can grow and blister, potentially leading to the contamination of the plasma, or to the degradation of their mechanical properties. We computationally study the behavior of small clusters of He atoms in W in conditions relevant to fusion energy production. Using a wide range of techniques, we investigate the thermodynamics of the clusters and their kinetics in terms of diffusivity, growth, and breakup, as well as mutation into nano-bubbles. Our study provides the essential ingredients to model the early stages of He exposure leading up to the nucleation of He bubbles.

preprint2014arXiv

Sampling in the multicanonical ensemble: Small He clusters in W

We carry out generalized-ensemble molecular dynamics simulations of the formation of small Helium (He) clusters in bulk Tungsten (W), a process of practical relevance for fusion energy production. We calculate formation free energies of small Helium clusters at temperatures up to the melting point of W, encompassing the whole range of interest for fusion-energy production. From this, parameters like cluster break-up or formation rates can be calculated, which help to refine models of microstructure evolution in He-irradiated Tungsten.