BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Alejandro Strachan

Alejandro Strachan contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
cond-mat.mtrl-sciphysics.comp-phComputational Engineering, Finance, and Sciencecond-mat.softDigital LibrariesMachine Learningnucl-thphysics.app-phphysics.chem-ph

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
15works
0followers
9topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

15 published item(s)

preprint2026arXiv

Multi-Fidelity Predictive Model for Shock Response of Energetic Materials Using Conditional U-Net

Mapping microstructure to properties is central to materials science. Perhaps most famously, the Hall-Petch relationship relates average grain size to strength. More challenging has been deriving relationships for properties that depend on subtle microstructural features and not average properties. One such example is the initiation of energetic materials under dynamical loading, dominated by energy localization on microstructural features such as pores, cracks, and interfaces. We propose a conditional convolutional neural network to predict the shock-induced temperature field as a function of shock strength, for a wide range of microstructures, and obtained via two different simulation methods. The proposed model, denoted MISTnet2, significantly extends prior work that was limited to a single shock strength, model, and type of microstructure. MISTnet2 can contribute to bridging atomistics with coarse-grain simulations and enable first principles predictions of detonation initiation and safety of this class of materials.

0 citations0 reviewsNo code linkedOpen access
preprint2026arXiv

Physics-constrained Gaussian Processes for Predicting Shockwave Hugoniot Curves

A physics-constrained Gaussian Process regression framework is developed for predicting shocked material states along the Hugoniot curve using data from a small number of shockwave simulations. The proposed Gaussian process employs a probabilistic Taylor series expansion in conjunction with the Rankine-Hugoniot jump conditions between the various shocked material states to construct a thermodynamically consistent covariance function. This leads to the formulation of an optimization problem over a small number of interpretable hyperparameters and enables the identification of regime transitions, from a leading elastic wave to trailing plastic and phase transformation waves. This work is motivated by the need to investigate shock-driven material response for materials discovery and for offering mechanistic insights in regimes where experimental characterizations and simulations are costly. The proposed methodology relies on large-scale molecular dynamics which are an accurate but expensive computational alternative to experiments. Under these constraints, the proposed methodology establishes Hugoniot curves from a limited number of molecular dynamics simulations. We consider silicon carbide as a representative material and atomic-level simulations are performed using a reverse ballistic approach together with appropriate interatomic potentials. The framework reproduces the Hugoniot curve with satisfactory accuracy while also quantifying the uncertainty in the predictions using the Gaussian Process posterior.

0 citations0 reviewsNo code linkedOpen access
preprint2026arXiv

Predictive models for strain energy in condensed phase reactions

Molecular modeling of thermally activated chemistry in condensed phases is essential to understand polymerization, depolymerization, and other processing steps of molecular materials. Current methods typically combine molecular dynamics (MD) simulations to describe short-time relaxation with a stochastic description of predetermined chemical reactions. Possible reactions are often selected on the basis of geometric criteria, such as a capture distance between reactive atoms. Although these simulations have provided valuable insight, the approximations used to determine possible reactions often lead to significant molecular strain and unrealistic structures. We show that the local molecular environment surrounding the reactive site plays a crucial role in determining the resulting molecular strain energy and, in turn, the associated reaction rates. We develop a graph neural network capable of predicting the strain energy associated with a cyclization reaction from the pre-reaction, local, molecular environment surrounding the reactive site. The model is trained on a large dataset of condensed-phase reactions during the activation of polyacrylonitrile (PAN) obtained from MD simulations and can be used to adjust relative reaction rates in condensed systems and advance our understanding of thermally activated chemical processes in complex materials

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Deviatoric Stress Driven Transient Melting Below the Glass Transition Temperature in Shocked Polymers

The relaxation of polymers around and below their glass transition temperature is governed by a range of correlated unit processes with a wide range of timescales. The fast deformation rates of shock loading can negate a significant fraction of these processes resulting in the dynamical glass transition in rubbers. In this letter we report the inverse, a transient melting of glassy polymer under shock loading. The large deviatoric stresses near the shock front induce fast transitions in backbone dihedral angles and a stress relaxation characteristic of polymer melts. This is followed by the slower relaxation expected for glasses.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Extemporaneous Mechanochemistry: Shockwave Induced Ultrafast Chemical Reactions Due to Intramolecular Strain Energy

Regions of energy localization referred to as hotspots are known to govern shock initiation and the run-to-detonation in energetic materials. Mounting computational evidence points to accelerated chemistry in hotspots from large intramolecular strains induced via the interactions between the shockwave and microstructure. However, definite evidence mapping intramolecular strain to accelerated or altered chemical reactions has so far been elusive. From a large-scale reactive molecular dynamics simulation of the energetic material TATB, we map molecular temperature and intramolecular strain energy prior to reaction to decomposition kinetics. Both temperature and intramolecular strain are shown to accelerate chemical kinetics. A detailed analysis of the atomistic trajectory shows that intramolecular strain can induce a mechanochemical alteration of decomposition mechanisms. The results in this paper can inform continuum-level chemistry models to account for a wide range of mechanochemical effects.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Role of electronic thermal transport in amorphous metal recrystallization: a molecular dynamics study

Recrystallization of glasses is important in a wide range of applications including electronics and reactive materials. Molecular dynamics (MD) has been used to provide an atomic picture of this process, but prior work has neglected the thermal transport role of electrons, the dominant thermal carrier in metallic systems. We characterize the role of electronic thermal conductivity on the velocity of recrystallization in Ni using MD coupled to a continuum description of electronic thermal transport via a two-temperature model. Our simulations show that for strong enough coupling between electrons and ions, the increased thermal conductivity removes the heat from the exothermic recrystallization process more efficiently, leading to a lower effective temperature at the recrystallization front and, consequently, lower propagation velocity. We characterize how electron-phonon coupling strength and system size affects front propagation velocity. Interestingly, we find that initial recrystallization velocity increases with decreasing in system size due to higher overall temperatures. Overall, we show that a more accurate description of thermal transport due to the incorporation of electrons results in better agreement with experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Systematic Builder for All-Atom Simulations of Plastically Bonded Explosives

The shock to detonation transition in heterogeneous plastically bonded explosives is dominated by energy localization into hotspots that arise from the interaction of the shockwave with microstructural features and defects. The complex polycrystalline structure of these materials leads to a network of hotspot that can coalesce into deflagration and detonation waves. Significant progress has been made on the formation and potency of hotspots using atomistic simulations, but most of the work has focused on ideal and isolated defects. Hence, developed a method, denoted PBXGen, to build realistic PBX microstructures for all-atom simulations. PBXGen is generally applicable, and we demonstrate it with two systems: an RDX-polystyrene PBX with a 3D microstructure and a TATB-polystyrene with columnar grains. The resulting structure exhibit key features of PBXs, albeit at smaller scales, and are validated against experimental mechanical and shock properties.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Sim2Ls: FAIR simulation workflows and data

Just like the scientific data they generate, simulation workflows for research should be findable, accessible, interoperable, and reusable (FAIR). However, while significant progress has been made towards FAIR data, the majority of science and engineering workflows used in research remain poorly documented and often unavailable, involving ad hoc scripts and manual steps, hindering reproducibility and stifling progress. We introduce Sim2Ls (pronounced simtools) and the Sim2L Python library that allow developers to create and share end-to-end computational workflows with well-defined and verified inputs and outputs. The Sim2L library makes Sim2Ls, their requirements, and their services discoverable, verifies inputs and outputs, and automatically stores results in a globally-accessible simulation cache and results database. This simulation ecosystem is available in nanoHUB, an open platform that also provides publication services for Sim2Ls, a computational environment for developers and users, and the hardware to execute runs and store results at no cost. We exemplify the use of Sim2Ls using two applications and discuss best practices towards FAIR simulation workflows and associated data.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Error correction in multi-fidelity molecular dynamics simulations using functional uncertainty quantification

We use functional, Fréchet, derivatives to quantify how thermodynamic outputs of a molecular dynamics (MD) simulation depend on the potential used to compute atomic interactions. Our approach quantifies the sensitivity of the quantities of interest with respect to the input functions as opposed to its parameters as is done in typical uncertainty quantification methods. We show that the functional sensitivity of the average potential energy and pressure in isothermal, isochoric MD simulations using Lennard-Jones two-body interactions can be used to accurately predict those properties for other interatomic potentials (with different functional forms) without re-running the simulations. This is demonstrated under three different thermodynamic conditions, namely a crystal at room temperature, a liquid at ambient pressure, and a high pressure liquid. The method provides accurate predictions as long as the change in potential can be reasonably described to first order and does not significantly affect the region in phase space explored by the simulation. The functional uncertainty quantification approach can be used to estimate the uncertainties associated with constitutive models used in the simulation and to correct predictions if a more accurate representation becomes available.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Functional uncertainty quantification for isobaric molecular dynamics simulations and defect formation energies

Functional uncertainty quantification (FunUQ) was recently proposed to quantify uncertainties in models and simulations that originate from input functions, as opposed to parameters. This paper extends FunUQ to quantify uncertainties originating from interatomic potentials in isothermal-isobaric molecular dynamics (MD) simulations and to the calculation of defect formation energies. We derive and verify a computationally inexpensive expression to compute functional derivatives in MD based on perturbation theory. We show that this functional derivative of the quantities of interest (average internal energy, volume, and defect energies in our case) with respect to the interatomic potential can be used to predict those quantities for a different interatomic potential, without re-running the simulation. The codes and scripts to perform FunUQ in MD are freely available for download. In addition, to facilitate reproducibility and to enable use of best practices for the approach, we created Jupyter notebooks to perform FunUQ analysis on MD simulations and made them available for online simulation in nanoHUB. The tool uses cloud computing resources and users can view, edit, and run end-to-end workflows from a standard web-browser without the need to need to download or install any software.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Online simulation powered learning modules for materials science

Simulation tools are playing an increasingly important role in materials science and engineering and beyond their well established importance in research and development, these tools have a significant pedagogical potential. We describe a set of online simulation tools and learning modules designed to help students explore important concepts in materials science where hands-on activities with high-fidelity simulations can provide insight not easily acquired otherwise. The online tools, which involve density functional theory and molecular dynamics simulations, have been designed with non-expert end-users in mind and only a few clicks are required to perform most simulations, yet they are powered by research-grade codes and expert users can access advanced options. All tools and modules are available for online simulation in nanoHUB.org and access is open and free of charge. Importantly, instructors and students do not need to download or install any software. The learning modules cover a range of topics from electronic structure of crystals and doping, plastic deformation in metals, and physical properties of polymers. These modules have been used in several core undergraduate courses at Purdue's School of Materials Engineering, they are self contained, and are easy to incorporate into existing classes.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The thermal conductivity of "pastas" in neutron star matter

This investigation explores the phononic thermal conductivity of nuclear star matter as it undergoes the "topological" transition to the "pasta" regime, and further down to the solid-liquid phase transition. The study was carried out using molecular dynamics simulations with nuclear potentials embedded in an effective (i.e. Thomas-Fermi) Coulomb potential. The thermal conductivity experiences a dramatic change within a narrow temperature interval around T=1 MeV. This change accomplishes the "pasta" breakdown during a heating process. The thermal conductivity by flipping protons' or neutrons' velocity further shows a decoupling for asymmetric nuclear star matter.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Tunability of martensitic behavior through coherent nanoprecipitates and other nanostructures

Molecular dynamics simulations show that coherent precipitates can significantly affect the properties of martensitic transformations in Ni$_{63}$Al$_{37}$ alloys. The precipitates, consisting of non-martensitic Ni$_{50}$Al$_{50}$, modify the free energy landscape that governs the phase transformation and result in a significant reduction of the thermal hysteresis, at comparably minor expense of transformation strain, and modification of transformation temperatures. Importantly, this paper shows that free energy landscape engineering is possible with nanostructures potentially accessible through standard metallurgical processing routes. The atomistic-level nucleation and transformation mechanisms within the nanoprecipitate systems are explored and compared with epitaxial nanolaminates and nanowires. The simulations reveal three distinct regimes of transformation mechanisms and martensitic nanostructure as a function of volume fraction of the non-martensitic phase. Free energy landscape engineering is generally applicable and could contribute to the design of new shape memory alloys with novel properties, such as light weight alloys that operate at room temperature.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Tuning martensitic transformations via coherent second phases in nanolaminates using free energy landscape engineering

We explore the possibilities and limitations of using a coherent second phase to engineer the thermo-mechanical properties of a martensitic alloy by modifying the underlying free energy landscape that controls the transformation. We use molecular dynamics simulations of a model atomistic system where the properties of a coherent, nanoscale second phase can be varied systematically. With a base martensitic material that undergoes a temperature-induced transformation from a cubic austenite to a monoclinic martensite, the simulations show a significant ability to engineer the transformation temperatures, from a ~50% reduction to a ~200% increase, with 50 at. % of the cubic second phase. We establish correlations between the properties of the second phase and the transformation characteristics and microstructure, via the free energy landscape of the two-phase systems. Coherency stresses have a strong influence on the martensitic variants observed and can even cause the non-martensitic second phase to undergo a transformation. Reducing the stiffness of second phase increases the transformation strain and modifies the martensitic microstructure, increasing the volume fraction of the transformed material. This increase in transformation strain is accompanied by a significant increase in the Af and thermal hysteresis, while the Ms remains unaltered. Our findings on the tunability of martensitic transformations can be used for informed searches of second phases to achieve desired material properties, such as achieving room temperature, lightweight shape memory alloys.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Uncharacteristic second order martensitic transformation in metals via epitaxial stress fields

While most phase transformations, e.g. ferroelectric or ferromagnetic, can be first or second order depending on external applied fields, martensitic transformations in metallic alloys are nearly universally first order. We demonstrate that epitaxial stress originating from the incorporation of a tailored second phase can modify the free energy landscape that governs the phase transition and change its order from first to second. High-fidelity molecular dynamics simulations show a remarkable change in the character of the martensitic transformation in Ni-Al alloys near the critical point. We observe the continuous evolution of the transformation order parameter and scaling with power-law exponents comparable to those in other ferroic transitions exhibiting critical behavior. Our theoretical work provides a foundation to recent experimental and computational results on martensites near critical points.

0 citations0 reviewsNo code linkedOpen access