BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Alejandro Strachan

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

ResearcherUnclaimed source record
cond-mat.mtrl-sciphysics.comp-phphysics.chem-phComputational Engineering, Finance, and Sciencecond-mat.mes-hallcond-mat.softDigital LibrariesMachine Learningnucl-thphysics.app-ph

Catalog footprint

What is connected

26works
10topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

26 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
preprint2016arXiv

The dynamics of copper intercalated molybdenum ditelluride

Layered transition metal dichalcogenides are emerging as key materials in nanoelectronics and energy applications. Predictive models to understand their growth, thermomechanical properties and interactions with metals are needed in order to accelerate their incorporation into commercial products. Interatomic potentials enable large-scale atomistic simulations at the device level, beyond the range of applications of first principle methods. We present a ReaxFF reactive force field to describe molybdenum ditelluride and its interactions with copper. We optimized the force field parameters to describe the properties of layered MoTe2 in various phases, the intercalation of Cu atoms and clusters within its van der Waals gap, including a proper description of energetics, charges and mechanical properties. The training set consists of an extensive set of first principle calculations computed from density functional theory. We use the force field to study the adhesion of a single layer MoTe2 on a Cu(111) surface and the results are in good agreement with density functional theory, even though such structures were not part of the training set. We characterized the mobility of the Cu ions intercalated into MoTe2 under the presence of an external electric fields via molecular dynamics simulations. The results show a significant increase in drift velocity for electric fields of approximately 0.4 V/A and that mobility increases with Cu ion concentration.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Atomistic simulations of amorphous polymers in the cloud with PolymerModeler

Molecular dynamics (MD) simulations enable the description of ma- terial properties and processes with atomistic detail by numerically solv- ing the time evolution of every atom in the system. We introduce Poly- merModeler, a general-purpose, online simulation tool to build atomistic structures of amorphous polymers and perform MD simulations on the re- sulting configurations to predict their thermo-mechanical properties. Poly- merModeler is available, free of charge, at nanoHUB.org via an interactive web interface, and the actual simulations are performed in the cloud us- ing nanoHUB.org resources. Starting from the specification of one or more monomers PolymerModeler builds the polymer chains into a simulation cell with periodic boundary conditions at the desired density. Monomers are added sequentially using the continuous configuration bias direct Monte Carlo method, and copolymers can be created describing the desired se- quence of monomers. PolymerModeler also enables users to perform MD simulations on the structures created by the builder using the Dreiding force field and Gasteiger partial atomic charges. We describe the force field im- plementation and the various options for the MD simulations that use the LAMMPS simulator. PolymerModeler MD simulations for a PMMA sam- ple show good structural agreement with experiments and are in good agree- ment with simulations obtained with commercial software.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Mesoscale simulations of shockwave energy dissipation via chemical reactions

We use a particle-based mesoscale model that incorporates chemical reactions at a coarse-grained level to study the response of materials that undergo volume-reducing chemical reactions under shockwave-loading conditions. We find that such chemical reactions can attenuate the shockwave and characterize how the parameters of the chemical model affect this behavior. The simulations show that the magnitude of the volume collapse and velocity at which the chemistry propagates are critical to weaken the shock, whereas the energetics in the reactions play only a minor role. Shock loading results in transient states where the material is away from local equilibrium and, interestingly, chemical reactions can nucleate under such non-equilibrium states. Thus, the timescales for equilibration between the various degrees of freedom in the material affect the shock-induced chemistry and its ability to attenuate the propagating shock.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Stability and Migration of Small Copper Clusters in Amorphous Dielectrics

We use density functional theory (DFT) to study the thermodynamic stability and migration of copper ions and small clusters embedded in amorphous silicon dioxide. We perform the calculations over an ensemble of statistically independent structures to quantify the role of the intrinsic atomic-level variability in the amorphous matrix affect the properties. The predicted formation energy of a Cu ion in the silica matrix is 2.7+/-2.4 eV, significantly lower the value for crystalline SiO2. Interestingly, we find that Cu clusters of any size are energetically favorable as compared to isolated ions; showing that the formation of metallic clusters does not require overcoming a nucleation barrier as is often assumed. We also find a broad distribution of activation energies for Cu migration, from 0.4 to 1.1 eV. This study provides insights into the stability of nanoscale metallic clusters in silica of interest in electrochemical metallization cell memories and optoelectronics.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Voltage equilibration for reactive atomistic simulations of electrochemical processes

We introduce EChemDID, a model to describe electrochemical driving force in reactive molecular dynamics simulations. The method describes the equilibration of external electrochemical potentials (voltage) within metallic structures and their effect on the self consistent partial atomic charges used in reactive molecular dynamics. An additional variable assigned to each atom denotes the local potential in its vicinity and we use fictitious, but computationally convenient, dynamics to describe its equilibration within not-simply connected metallic structures on-the-fly during the molecular dynamics simulation. This local electrostatic potential is used to dynamically modify the atomic electronegativities used to compute partial atomic changes via charge equilibration. Validation tests show that the method provides an accurate description of the electric fields generated by the applied voltage and the driving force for electrochemical reactions. We demonstrate EChemDID via simulations of the operation of electrochemical metallization cells. The simulations predict the switching of the device between a high-resistance to a low-resistance state as a conductive metallic bridge is formed and resistive currents that can be compared with experimental measurements. In addition to applications in nanoelectronics, EChemDID could be useful to model electrochemical energy conversion devices.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Mesodynamics with implicit degrees for freedom

Mesoscale phenomena -- involving a level of description between the finest atomistic scale and the macroscopic continuum -- can be studied by a variation on the usual atomistic-level molecular dynamics (MD) simulation technique. In mesodynamics, the mass points, rather than being atoms, are mesoscopic in size, for instance representing the centers of mass of polycrystalline grains or molecules. In order to reproduce many of the overall features of fully atomistic MD, which is inherently more expensive, the equations of motion in mesodynamics must be derivable from an interaction potential that is faithful to the compressive equation of state, as well as to tensile de-cohesion that occurs along the boundaries of the mesoscale units. Moreover, mesodynamics diffrers from Newton's equations of motion in that dissipation -- the exchange of energy between mesoparticles and their internal degrees of freedom (DoFs) -- must be described, and so should the transfer of energy between the internal modes of neighboring mesoparticles. We present a formulation where energy transfer between the internal modes of a mesoparticle and its external center-of-mass DoFs occurs in the phase space of mesoparticle coordinates, rather than momenta, resulting in a Galilean invariant formulation that conserves total linear momentum and energy (including the energy internal to the mesoparticles). We show that this approach can be used to describe, in addition to mesoscale problems, conduction electrons in atomic-level simulations of metals, and we demonstrate applications of mesodynamics to shockwave propagation and thermal transport.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Role of Strain on Electronic and Mechanical Response of Semiconducting Transition-Metal Dichalcogenide Monolayers: an ab-initio study

We characterize the electronic structure and elasticity of monolayer transition-metal dichalcogenides MX2 (M=Mo, W, Sn, Hf and X=S, Se, Te) with 2H and 1T structures using fully relativistic first principles calculations based on density functional theory. We focus on the role of strain on the band structure and band alignment across the series 2D materials. We find that strain has a significant effect on the band gap; a biaxial strain of 1% decreases the band gap in the 2H structures, by as a much 0.2 eV in MoS2 and WS2, while increasing it for the 1T materials. These results indicate that strain is a powerful avenue to modulate their properties; for example, strain enables the formation of, otherwise impossible, broken gap heterostructures within the 2H class. These calculations provide insight and quantitative information for the rational development of heterostructures based on these class of materials accounting for the effect of strain.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Strain enhancement of high-k dielectric response in (La/Sc)2 O3 and LaScO3 : an ab-initio study

We use density functional theory within the generalized gradient approximation to characterize the dielectric response of rare earth oxides: (La,Sc)2 O3 bixbyite, and LaScO3 perovskite. We focus on the role of strain on the phonon contribution of the dielectric constant and find that, contrary to the classical expectation based on the Clausius-Mossotti relation, tensile volumetric strain and volume-conserving bi- axial strain on the order of +/-1% can lead to an increase in dielectric constant of up to 20%. The insight into the atomic mechanisms responsible for these effects and the quantitative results in this paper can contribute to the development and understanding of high-κ materials.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Phonon thermal transport outside of local equilibrium in nanowires via molecular dynamics

We study thermal transport through Pt nanowires that bridge planar contacts as a function of wire length and vibrational frequency of the contacts. When phonons in the contacts have lower average frequencies than those in the wires thermal transport occurs under conditions away from local equilibrium with low-frequency phonons experiencing a higher thermal gradient than high-frequency ones. This results in a size-dependent increase in the effective thermal conductivity of the wire with decreasing vibrational frequencies of the contacts. The interfacial resistivity when heat flows from the wire to the contact is also size-dependent and has the same physical origin in the lack of full equilibration in short nanowires. We develop a model based on a 1D atomic chain that captures the salient physics of the MD results.

0 citations0 reviewsNo code linkedOpen access
preprint2002arXiv

Ab initio and finite-temperature molecular dynamics studies of lattice resistance in tantalum

This manuscript explores the apparent discrepancy between experimental data and theoretical calculations of the lattice resistance of bcc tantalum. We present the first results for the temperature dependence of the Peierls stress in this system and the first ab initio calculation of the zero-temperature Peierls stress to employ periodic boundary conditions, which are those best suited to the study of metallic systems at the electron-structure level. Our ab initio value for the Peierls stress is over five times larger than current extrapolations of experimental lattice resistance to zero-temperature. Although we do find that the common techniques for such extrapolation indeed tend to underestimate the zero-temperature limit, the amount of the underestimation which we observe is only 10-20%, leaving open the possibility that mechanisms other than the simple Peierls stress are important in controlling the process of low temperature slip.

0 citations0 reviewsNo code linkedOpen access
preprint2001arXiv

Accurate Calculations of the Peierls Stress in Small Periodic Cells

The Peierls stress for a [111]-screw dislocation in bcc Tantalum is calculated using an embedded atom potential. More importantly, a method is presented which allows accurate calculations of the Peierls stress in the smallest periodic cells. This method can be easily applied to ab initio calculations, where only the smallest unit cells capable of containing a dislocation can be conviently used. The calculation specifically focuses on the case where the maximum resolved shear stress is along a {110}-plane.

0 citations0 reviewsNo code linkedOpen access