Source author record

Mark Asta

Mark Asta 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

19works
6topics
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

19 published item(s)

preprint2026arXiv

Nanocrystal Geometry Governs Phase Transformation Pathways in Palladium Hydride

Pathways and structural dynamics of phase transformations impact performance of materials in energy and information storage technologies. Palladium hydride ($\mathrm{PdH}_x$) nanocrystals are an ideal model system for studying solute-induced phase transformations, where elastic energy from lattice mismatch between $α$-$\mathrm{PdH}_x$ and $β$-$\mathrm{PdH}_x$ phases is often considered a key to determining the transformation pathways. $α/β$-$\mathrm{PdH}_x$ interfacial elastic energy is affected by the confined geometry of a nanocrystal. However, how nanocrystal geometry influences phase transformation pathways is largely unknown. Using in situ liquid phase transmission electron microscopy, we directly visualize hydrogenation in Pd nanocrystals with two geometries -- a nanocube and a hexagonal nanoplate. Both follow similar sequences of an initially curved nucleus, interface flattening, and reverse-stage nucleation; however, their evolving $α/β$-$\mathrm{PdH}_x$ interfaces exhibit geometry-dependent crystallographic alignments. In nanocubes, $\{100\}$-aligned configurations conform to static elastic energy ordering, representing a pathway that maintains a local mechanical equilibrium, whereas nanoplates display both $\{110\}$- and $\{211\}$-aligned interfaces. Theoretical simulations show that geometry determines the accessibility of alternative phase transformation pathways as the system is driven far from equilibrium during hydrogenation. These findings identify geometry as a fundamental parameter for directing phase transformation pathways, offering design principles for accessing atypical configurations and improving properties of intercalation-based devices.

preprint2022arXiv

Nonequilibrium free-energy calculation of solids using LAMMPS

This article describes nonequilibrium techniques for the calculation of free energies of solids using molecular dynamics (MD) simulations. These methods provide an alternative to standard equilibrium thermodynamic integration methods and often present superior efficiency. Here we describe the implementation in the LAMMPS (Large-scale Atomic/Molecular Massively Parallel Simulator) code of two specific nonequilibrium processes that allow the calculation of the free-energy difference between two different system Hamiltonians as well as the free-energy temperature dependence of a given Hamiltonian, respectively. The theory behind the methods is summarized, and we describe (including fragments of LAMMPS scripts) how the process parameters should be selected to obtain the best-possible efficiency in the calculations of free energies using nonequilibrium MD simulations. As an example of the application of the methods we present results related to polymorphic transitions for a classical potential model of iron.

preprint2020arXiv

Ab initio modeling of the energy landscape for screw dislocations in body-centered cubic high-entropy alloys

In traditional body-centered cubic (bcc) metals, the core properties of screw dislocations play a critical role in plastic deformation at low temperatures. Recently, much attention has been focused on refractory high-entropy alloys (RHEAs), which also possess bcc crystal structures. However, unlike face-centered cubic high-entropy alloys (HEAs), there have been far fewer investigations on bcc HEAs, specifically on the possible effects of chemical short-range order (SRO) in these multiple principal element alloys on dislocation mobility. Here, using density functional theory, we investigate the distribution of dislocation core properties in MoNbTaW RHEAs alloys, and how they are influenced by SRO. The average values of the core energies in the RHEA are found to be larger than those in the corresponding pure constituent bcc metals, and are relatively insensitive to the degree of SRO. However, the presence of SRO is shown to have a large effect on narrowing the distribution of dislocation core energies and decreasing the spatial heterogeneity of dislocation core energies in the RHEA. It is argued that the consequences for the mechanical behavior of HEAs is a change in the energy landscape of the dislocations which would likely heterogeneously inhibit their motion.

preprint2020arXiv

Kinetics of Crystallization and Orientational Ordering in Dipolar Particle Systems

The kinetic mechanisms underlying bottom-up assembly of colloidal particles have been widely investigated in efforts to control crystallization pathways and to direct growth into targeted superstructures for applications including photonic crystals. Current work builds on recent progress in the development of kinetic theories for crystal growth of body-centered-cubic crystals in systems with short-range inter-particle interactions, accounting for a greater diversity of crystal structures and the role of the longer-ranged interactions and orientational degrees of freedom arising in polar systems. We address the importance of orientational ordering processes in influencing crystal growth in such polar systems, thus advancing the theory beyond the treatment of the translational ordering processes considered in previous investigations. The work employs comprehensive molecular-dynamics simulations that resolve key crystallization processes, and are used in the development of a quantitative theoretical framework based on ideas from time-dependent Ginzburg-Landau theory. The significant impact of orientational ordering on the crystallization kinetics could be potentially leveraged to achieve crystallization kinetics steering through external electric or magnetic fields. Our combined theory/simulation approach provides opportunities for future investigations of more complex crystallization kinetics.

preprint2019arXiv

Verification of Short-Range Order and Its Impact on the Properties of the CrCoNi Medium Entropy Alloy

Traditional metallic alloys are mixtures of elements where the atoms of minority species tend to distribute randomly if they are below their solubility limit, or lead to the formation of secondary phases if they are above it. Recently, the concept of medium/high entropy alloys (MEA/HEA) has expanded this view, as these materials are single-phase solid solutions of generally equiatomic mixtures of metallic elements that have been shown to display enhanced mechanical properties. However, the question has remained as to how random these solid solutions actually are, with the influence of chemical short-range order (SRO) suggested in computational simulations but not seen experimentally. Here we report the first direct observation of SRO in the CrCoNi MEA using high resolution and energy-filtered transmission electron microscopy. Increasing amounts of SRO give rise to both higher stacking fault energy and hardness. These discoveries suggest that the degree of chemical ordering at the nanometer scale can be tailored through thermomechanical processing, providing a new avenue for tuning the mechanical properties of MEA/HEAs.

preprint2016arXiv

Effect of strain and temperature on the threshold displacement energy in body-centered cubic iron

The threshold displacement energy (TDE) is the minimum amount of kinetic energy required to displace an atom from its lattice site. The magnitude of the TDE displays significant variance as a function of the crystallographic direction, system temperature and applied strain, among a variety of other factors. It is critically important to determine an accurate value of the TDE in order to calculate the total number of displacements due to a given irradiation condition, and thus to understand the materials response to irradiation. In this study, molecular dynamics simulations have been performed to calculate the threshold displacement energy in body-centered cubic iron as a function of strain and temperature. With applied strain, a decrease of the TDE of up to approximately 14 eV was observed. A temperature increase from 300 K to 500 K can result in an increase of the TDE of up to approximately 9 eV.

preprint2016arXiv

Predicting the lattice thermal conductivity of solids by solving the Boltzmann transport equation: AFLOW - AAPL an automated, accurate and effcient framework

One of the most accurate approaches for calculating lattice thermal conductivity, $κ_l$, is solving the Boltzmann transport equation starting from third-order anharmonic force constants. In addition to the underlying approximations of ab-initio parameterization, two main challenges are associated with this path. High computational costs and lack of automation in the frameworks using this methodology affect the discovery rate of novel materials with ad-hoc properties. Here, we present the Automatic-Anharmonic-Phonon-Library, AAPL. It efficiently computes interatomic force constants by making effective use of crystal symmetry analysis, it solves the Boltzmann transport equation to obtain $κ_l$, and allows a fully integrated operation with minimum user intervention, a rational addition to the current high-throughput accelerated materials development framework AFLOW. We show an "experiment versus theory" study of the approach, we compare accuracy and speed with respect to other available packages, and for materials characterized by strong electron localization and correlation, we demonstrate that it is possible to improve accuracy without increasing computational requirements by combining AAPL with the pseudo-hybrid functional ACBN0.

preprint2015arXiv

Second-Nearest-Neighbor Correlations from Connection of Atomic Packing Motifs in Metallic Glasses and Liquids

Using molecular dynamics simulations, we have studied the atomic correlations characterizing the second peak in the radial distribution function (RDF) of metallic glasses and liquids. The analysis was conducted from the perspective of different connection schemes of atomic packing motifs, based on the number of shared atoms between two linked coordination polyhedra. The results demonstrate that the cluster connections by face-sharing, specifically with three common atoms, are most favored when transitioning from the liquid to glassy state, and exhibit the stiffest elastic response during shear deformation. These properties of the connections and the resultant atomic correlations are generally the same for different types of packing motifs in different alloys. Splitting of the second RDF peak was observed for the inherent structure of the equilibrium liquid, originating solely from cluster connections; this trait can then be inherited in the metallic glass formed via subsequent quenching of the parent liquid through the glass transition, in the absence of any additional type of local structural order. Increasing ordering and cluster connection during cooling, however, may tune the position and intensity of the split peaks.

preprint2015arXiv

Segregation-induced phase transformations in grain boundaries

Phase transformations in metallic grain boundaries (GBs) present significant fundamental interest in the context of thermodynamics of low-dimensional physical systems. We report on atomistic computer simulations of the Cu-Ag system that provide direct evidence that GB phase transformations in a single-component GB can continue to exist in a binary alloy. This gives rise to segregation-induced phase transformations with varying chemical composition at a fixed temperature. Furthermore, for such transformations we propose an approach to calculations of free energy differences between different GB phases by thermodynamic integration along a segregation isotherm. This opens the possibility of developing quantitative thermodynamics of GB phases, their transformations to each other, and critical phenomena in the future.

preprint2014arXiv

Effects of applied strain on radiation damage generation in body-centered cubic iron

Radiation damage in body-centered cubic (BCC) Fe has been extensively studied by computer simulations to quantify effects of temperature, impinging particle energy, and the presence of extrinsic particles. However, limited investigation has been conducted into the effects of mechanical stresses and strain. In a reactor environment, structural materials are often mechanically strained, and an expanded understanding of how this strain affects the generation of defects may be important for predicting microstructural evolution and damage accumulation under such conditions. In this study, we have performed molecular dynamics simulations in which various types of homogeneous strains are applied to BCC Fe and the effect on defect generation is examined. It is found that volume-conserving shear strains yield no statistically significant variations in the stable number of defects created via cascades in BCC Fe. However, strains that result in volume changes are found to produce significant effects on defect generation.

preprint2014arXiv

Evaluation of Constant Potential Method in Simulating Electric Double-Layer Capacitors

A major challenge in the molecular simulation of electric double layer capacitors (EDLCs) is the choice of an appropriate model for the electrode. Typically, in such simulations the electrode surface is modeled using a uniform fixed charge on each of the electrode atoms, which ignores the electrode response to local charge fluctuations induced by charge fluctuations in the electrolyte. In this work, we evaluate and compare this Fixed Charge Method (FCM) with the more realistic Constant Potential Method (CPM), [Reed, et al., J. Chem. Phys., 126, 084704 (2007)], in which the electrode charges fluctuate in order to maintain constant electric potential in each electrode. For this comparison, we utilize a simplified LiClO$_4$-acetonitrile/graphite EDLC. At low potential difference ($ΔΨ\le 2V$), the two methods yield essentially identical results for ion and solvent density profiles; however, significant differences appear at higher $ΔΨ$. At $ΔΨ\ge 4V$, the CPM ion density profiles show significant enhancement (over FCM) of "partially electrode solvated" Li$^+$ ions very close to the electrode surface. The ability of the CPM electrode to respond to local charge fluctuations in the electrolyte is seen to significantly lower the energy (and barrier) for the approach of Li$^+$ ions to the electrode surface.

preprint2014arXiv

High-Throughput Computational Screening of thermal conductivity, Debye temperature and Grüneisen parameter using a quasi-harmonic Debye Model

The quasi-harmonic Debye approximation has been implemented within the AFLOW and Materials Project frameworks for high-throughput computational science (Automatic Gibbs Library, AGL), in order to calculate thermal properties such as the Debye temperature and the thermal conductivity of materials. We demonstrate that the AGL method, which is significantly cheaper computationally compared to the fully ab initio approach, can reliably predict the ordinal ranking of the thermal conductivity for several different classes of semiconductor materials. We also find that for the set of 182 materials investigated in this work the Debye temperature, calculated with the AGL, is often a better predictor of the ordinal ranking of the experimental thermal conductivities than the calculated thermal conductivity. The Debye temperature is thus a potential descriptor for high-throughput screening of the thermal properties of materials.

preprint2013arXiv

Structural phase transformations in metallic grain boundaries

Structural transformations at interfaces are of profound fundamental interest as complex examples of phase transitions in low-dimensional systems. Despite decades of extensive research, no compelling evidence exists for structural transformations in high-angle grain boundaries in elemental systems. Here we show that the critical impediment to observations of such phase transformations in atomistic modeling has been rooted in inadequate simulation methodology. The proposed new methodology allows variations in atomic density inside the grain boundary and reveals multiple grain boundary phases with different atomic structures. Reversible first-order transformations between such phases are observed by varying temperature or injecting point defects into the boundary region. Due to the presence of multiple metastable phases, grain boundaries can absorb significant amounts of point defects created inside the material by processes such as irradiation. We propose a novel mechanism of radiation damage healing in metals which may guide further improvements in radiation resistance of metallic materials through grain boundary engineering.

preprint2012arXiv

Phase-field-crystal study of grain boundary premelting and shearing in bcc iron

We use the phase-field-crystal (PFC) method to investigate the equilibrium premelting and nonequilibrium shearing behaviors of $[001]$ symmetric tilt grain boundaries (GBs) at high homologous temperature over the complete range of misorientation $0<θ<90^\circ$ in classical models of bcc Fe. We characterize the dependence of the premelted layer width $W$ as a function of temperature and misorientation and compute the thermodynamic disjoining potential whose derivative with respect to $W$ represents the structural force between crystal-melt interfaces due to the spatial overlap of density waves. The disjoining potential is also computed by molecular dynamics (MD) simulations, for quantitative comparison with PFC simulations, and coarse-grained amplitude equations (AE) derived from PFC that provide additional analytical insights. We find that, for GBs over an intermediate range of misorientation ($θ_{\rm min}<θ<θ_{\rm max}$), $W$ diverges as the melting temperature is approached from below, corresponding to a purely repulsive disjoining potential, while for GBs outside this range ($θ<θ_{\rm min}$ or $θ_{\rm max}<θ<90^\circ$), $W$ remains finite at the melting point, with its value corresponding to a shallow attractive minimum of the disjoining potential. In response to a shear stress parallel to the GB plane, GBs in PFC simulations exhibit coupled motion normal to this plane, with a discontinuous change of the coupling factor as a function of $θ$ that reflects a transition between two coupling modes, and/or sliding (shearing of the two grains). The coupling factor for the two coupling modes is in excellent quantitative agreement with previous theoretical predictions [J. W. Cahn, Y. Mishin, and A. Suzuki, Acta Mater. 54, 4953 (2006)].

preprint2012arXiv

Solid-liquid interfacial premelting

We report the observation of a premelting transition at chemically sharp solid-liquid interfaces using molecular-dynamics simulations. The transition is observed in the solid-Al/liquid-Pb system and involves the formation of a liquid interfacial film of Al with a width that grows logarithmically as the bulk melting temperature is approached from below, consistent with current theories of premelting. The premelting behavior leads to a sharp change in the temperature dependence of the diffusion coefficient in the interfacial region, and could have important consequences for phenomena such as particle coalescence and shape equilibration, which are governed by interfacial kinetic processes.

preprint2012arXiv

Step Free Energies at Faceted Solid-Liquid Interfaces from Equilibrium Molecular Dynamics Simulations

In this work a method is proposed for computing step free energies for faceted solid-liquid interfaces based on atomistic simulations. The method is demonstrated in an application to (111) interfaces in elemental Si, modeled with the classical Stillinger-Weber potential. The approach makes use of an adiabatic trapping procedure, and involves simulations of systems with coexisting solid and liquid phases separated by faceted interfaces containing islands with different sizes, for which the corresponding equilibrium temperatures are computed. We demonstrate that the calculated coexistence temperature is strongly affected by the geometry of the interface. We find that island radius is inversely proportional to superheating, allowing us to compute the step free energy by fitting simulation data within the formalism of classical nucleation theory. The approach outlined in this work paves the way to the calculation of step free energies relevant to the solidification of faceted crystals from liquid mixtures, as encountered in nanowire growth by the vapor-liquid-solid mechanism and in alloy casting. The present work also shows that at low undercoolings the Stillinger-Weber interatomic potential for Si tends to crystallize in the wurtzite, rather than the diamond-cubic structure.

preprint2010arXiv

First-principles computational study of defect clustering in solid solutions of ThO$_{2}$ with trivalent oxides

The energetics of mixing and defect ordering in solid solutions of fluorite-structured ThO$_{2}$ with oxides of trivalent cations (Sc, In, Y, Nd, La) are investigated by electronic density-functional-theory (DFT). Through DFT calculations of structures enumerated by lattice-algebra techniques, we identify the lowest-energy patterns for defect clustering for four separate dopant concentrations. The most stable structures are characterized by a repulsive interaction between nearest-neighbor vacancies on the oxygen sublattice. The enthalpies of formation with respect to constituent oxides are positive for all dopants considered, and show a tendency to decrease in magnitude as the size and electronegativity of the trivalent dopant decrease. Due to the small positive formation enthalpies and low oxygen-vacancy binding energy with La dopants, La$_{2}$O$_{3}$-ThO$_{2}$ solid solutions are predicted to have relatively high ionic conductivities relative to those for the other aliovalent dopants considered. Our results are compared with those for the more widely studied ZrO$_{2}$ and CeO$_{2}$ fluorite-structured solid solutions with trivalent cations.

preprint2010arXiv

Simulation of ion track ranges in uranium oxide

Direct comparisons between statistically sound simulations of ion-tracks and published experimental measurements of range densities of iodine implants in uranium dioxide have been made with implant energies in the range of 100-800 keV. Our simulations are conducted with REED-MD (Rare Event Enhanced Domain-following Molecular Dynamics) in order to account for the materials structure in both single crystalline and polycrystalline experimental samples. We find near-perfect agreement between REED-MD results and experiments for polycrystalline target materials.

preprint2010arXiv

Structural disjoining potential for grain boundary premelting and grain coalescence from molecular-dynamics simulations

We describe a molecular dynamics framework for the direct calculation of the short-ranged structural forces underlying grain-boundary premelting and grain-coalescence in solidification. The method is applied in a comparative study of (i) a Sigma 9 <115> 120 degress twist and (ii) a Sigma 9 <110> {411} symmetric tilt boundary in a classical embedded-atom model of elemental Ni. Although both boundaries feature highly disordered structures near the melting point, the nature of the temperature dependence of the width of the disordered regions in these boundaries is qualitatively different. The former boundary displays behavior consistent with a logarithmically diverging premelted layer thickness as the melting temperature is approached from below, while the latter displays behavior featuring a finite grain-boundary width at the melting point. It is demonstrated that both types of behavior can be quantitatively described within a sharp-interface thermodynamic formalism involving a width-dependent interfacial free energy, referred to as the disjoining potential. The disjoining potential for boundary (i) is calculated to display a monotonic exponential dependence on width, while that of boundary (ii) features a weak attractive minimum. The results of this work are discussed in relation to recent simulation and theoretical studies of the thermodynamic forces underlying grain-boundary premelting.