-
Unveiling the effect of Ni on the formation and structure of Earth's inner core
Authors:
Yang Sun,
Mikhail I. Mendelev,
Feng Zhang,
Xun Liu,
Bo Da,
Cai-Zhuang Wang,
Renata M. Wentzcovitch,
Kai-Ming Ho
Abstract:
Ni is the second most abundant element in the Earth's core. Yet, its effects on the inner core's structure and formation process are usually disregarded because of its electronic and size similarity with Fe. Using ab initio molecular dynamics simulations, we find that the bcc phase can spontaneously crystallize in liquid Ni at temperatures above Fe's melting point at inner core pressures. The melt…
▽ More
Ni is the second most abundant element in the Earth's core. Yet, its effects on the inner core's structure and formation process are usually disregarded because of its electronic and size similarity with Fe. Using ab initio molecular dynamics simulations, we find that the bcc phase can spontaneously crystallize in liquid Ni at temperatures above Fe's melting point at inner core pressures. The melting temperature of Ni is shown to be 700-800 K higher than that of Fe at 323-360 GPa. hcp, bcc, and liquid phase relation differ for Fe and Ni. Ni can be a bcc stabilizer for Fe at high temperatures and inner core pressures. A small amount of Ni can accelerate Fe's crystallization at core pressures. These results suggest Ni may substantially impact the structure and formation process of the solid inner core.
△ Less
Submitted 22 September, 2023; v1 submitted 8 August, 2023;
originally announced August 2023.
-
Ab initio melting temperatures of bcc and hcp iron under the Earth's inner core condition
Authors:
Yang Sun,
Mikhail I. Mendelev,
Feng Zhang,
Xun Liu,
Bo Da,
Cai-Zhuang Wang,
Renata M. Wentzcovitch,
Kai-Ming Ho
Abstract:
There has been a long debate on the stable phase of iron under the Earth's inner core conditions. Because of the solid-liquid coexistence at the inner core boundary, the thermodynamic stability of solid phases directly relates to their melting temperatures, which remains considerable uncertainty. In the present study, we utilized a semi-empirical potential fitted to high-temperature ab initio data…
▽ More
There has been a long debate on the stable phase of iron under the Earth's inner core conditions. Because of the solid-liquid coexistence at the inner core boundary, the thermodynamic stability of solid phases directly relates to their melting temperatures, which remains considerable uncertainty. In the present study, we utilized a semi-empirical potential fitted to high-temperature ab initio data to perform a thermodynamic integration from classical systems described by this potential to ab initio systems. This method provides a smooth path for thermodynamic integration and significantly reduces the uncertainty caused by the finite-size effect. Our results suggest the hcp phase is the stable phase of pure iron under the inner core conditions, while the free energy difference between the hcp and bcc phases is tiny, on the order of 10s meV/atom near the melting temperature.
△ Less
Submitted 9 March, 2023; v1 submitted 4 May, 2022;
originally announced May 2022.
-
Two-step nucleation of the Earth's inner core
Authors:
Yang Sun,
Feng Zhang,
Mikhail I. Mendelev,
Renata M. Wentzcovitch,
Kai-Ming Ho
Abstract:
It has long been assumed the Earth's solid inner core started to grow when molten iron cooled to its melting point. However, the nucleation mechanism, which is a necessary step of crystallization, has not been well understood. Recent studies found it requires an unrealistic degree of undercooling to nucleate the stable hexagonal close-packed (hcp) phase of iron, which can never be reached under th…
▽ More
It has long been assumed the Earth's solid inner core started to grow when molten iron cooled to its melting point. However, the nucleation mechanism, which is a necessary step of crystallization, has not been well understood. Recent studies found it requires an unrealistic degree of undercooling to nucleate the stable hexagonal close-packed (hcp) phase of iron, which can never be reached under the actual Earth's core conditions. This contradiction leads to the inner core nucleation paradox [1]. Here, using a persistent-embryo method and molecular dynamics simulations, we demonstrate that the metastable body-centered cubic (bcc) phase of iron has a much higher nucleation rate than the hcp phase under inner-core conditions. Thus, the bcc nucleation is likely to be the first step of inner core formation instead of direct nucleation of the hcp phase. This mechanism reduces the required undercooling of iron nucleation, which provides a key factor to solve the inner-core nucleation paradox. The two-step nucleation scenario of the inner core also opens a new avenue for understanding the structure and anisotropy of the present inner core.
△ Less
Submitted 7 January, 2022; v1 submitted 15 May, 2021;
originally announced May 2021.
-
Development of interatomic potential appropriate for simulation of dislocation migration in fcc Fe
Authors:
Mikhail I. Mendelev,
Valery Borovikov
Abstract:
Molecular dynamics (MD) simulation of dislocation migration requires semi-empirical potentials of the interatomic interaction. While there are many reliable semi-empirical potentials for the bcc Fe, the number of the available potentials for the fcc is very limited. In the present study we tested three EAM potentials for the fcc Fe (ABCH97 [Phil. Mag. A, 75, 713-732 (1997)], BCT13 [MSMSE 21, 08500…
▽ More
Molecular dynamics (MD) simulation of dislocation migration requires semi-empirical potentials of the interatomic interaction. While there are many reliable semi-empirical potentials for the bcc Fe, the number of the available potentials for the fcc is very limited. In the present study we tested three EAM potentials for the fcc Fe (ABCH97 [Phil. Mag. A, 75, 713-732 (1997)], BCT13 [MSMSE 21, 085004 (2013)] and ZFS18 [J. Comp. Chem. 39, 2420-2431 (2018)]) from literature. It was found that the ABCH97 potential does not provide that the fcc phase is the most stable at any temperature. On the other hand, the fcc phase is always more stable than the bcc phase for the BCT13, ZFS18 potentials. The hcp phase is the most stable phase for the BCT13 potential at any temperature. In order to fix these problems we developed two new EAM potentials (MB1 and MB2). The fcc phase is still more stable than the bcc phase for the MB1 potential but the MB2 potential provides that the bcc phase is the most stable phase from the upper fcc-bcc transformation temperature, T_gamma-delta, to the melting temperature, Tm, and the fcc phase is the most stable phase below T_gamma-delta. This potential also leads to an excellent agreement with the experimental data on the fcc elastic constants and reasonable stacking fault energy which makes it the best potential for the simulation of the dislocation migration in the fcc Fe among all semi-empirical potentials considered in the present study. The MD simulation demonstrated that only the ZFS18, MB1 and MB2 potentials are actually suitable for the simulation of the dislocation migration in the fcc Fe. They lead to the same orders of magnitude for the dislocation velocities and all of them show that the edge dislocation is faster than the screw dislocation. However, the actual values of the dislocation velocities do depend on the employed semi-empirical potential.
△ Less
Submitted 25 September, 2020;
originally announced September 2020.
-
Competitive B2 and B33 Nucleation during Solidification of Ni50Zr50 Alloy: Molecular Dynamics Simulation and Classical Nucleation Theory
Authors:
Yang Sun,
Feng Zhang,
Huajing Song,
Mikhail I. Mendelev,
Cai-Zhuang Wang,
Kai-Ming Ho
Abstract:
We investigated the homogenous nucleation of the stoichiometric B2 and B33 phases in the Ni50Zr50 alloy using the persistent embryo method and the classical nucleation theory. The two phases become very close competitors at large supercoolings, which is consistent with the experimental observations. In the case of the B2 phase, the linear temperature dependence of the solid-liquid interface (SLI)…
▽ More
We investigated the homogenous nucleation of the stoichiometric B2 and B33 phases in the Ni50Zr50 alloy using the persistent embryo method and the classical nucleation theory. The two phases become very close competitors at large supercoolings, which is consistent with the experimental observations. In the case of the B2 phase, the linear temperature dependence of the solid-liquid interface (SLI) free energy extrapolated to the melting temperature leads to the same value as the one obtained from the capillarity fluctuation method (CFM). In the case of the B33 phases, the SLI free energy is also a linear function of temperature at large supercoolings but the extrapolation to the melting temperature leads to a value which is considerably different from the CFM value. This is consistent with the large anisotropy of the SLI properties of the B33 phase nearby the melting temperature observed in the simulation of the nominally flat interface migration.
△ Less
Submitted 23 February, 2019;
originally announced February 2019.
-
Effects of dopants on the glass forming ability in Al-based metallic alloy
Authors:
Yang Sun,
Feng Zhang,
Lin Yang,
Huajing Song,
Mikhail I. Mendelev,
Cai-Zhuang Wang,
Kai-Ming Ho
Abstract:
The effect of dopants on the metallic glass forming ability is usually considered based on analysis of changes in the liquid structure or thermodynamics. What is missing in such considerations is an analysis of how a dopant changes the properties of the crystal phases which can form instead of the glass. In order to illuminate this aspect we performed molecular dynamics simulations to study the ef…
▽ More
The effect of dopants on the metallic glass forming ability is usually considered based on analysis of changes in the liquid structure or thermodynamics. What is missing in such considerations is an analysis of how a dopant changes the properties of the crystal phases which can form instead of the glass. In order to illuminate this aspect we performed molecular dynamics simulations to study the effects of Mg and Sm dopants on the crystal nucleation in Al. The simulation data were found to be consistent with the experimental observations that addition of Mg to Al does not lead to vitrification but addition of only 8% Sm does. The significant effect of Sm doping was related to the intolerance of Al to this dopant. This leads to increase in the solid-liquid interfacial free energy, and therefore, to increase in the nucleation barrier and to dramatic decrease in the nucleation rate. The intolerance mechanism also significantly affects the growth kinetics.
△ Less
Submitted 15 February, 2019; v1 submitted 7 November, 2018;
originally announced November 2018.
-
Temperature dependence of the solid-liquid interface free energy of Ni and Al from molecular dynamics simulation of nucleation
Authors:
Yang Sun,
Feng Zhang,
Huajing Song,
Mikhail I. Mendelev,
Cai-Zhuang Wang,
Kai-Ming Ho
Abstract:
The temperature dependence of the solid-liquid interfacial free energy, γ, is investigated for Al and Ni at the undercooled temperature regime based on a recently developed persistent-embryo method. The atomistic description of the nucleus shape is obtained from molecular dynamics simulations. The computed γ shows a linear dependence on the temperature. The values of γ extrapolated to the melting…
▽ More
The temperature dependence of the solid-liquid interfacial free energy, γ, is investigated for Al and Ni at the undercooled temperature regime based on a recently developed persistent-embryo method. The atomistic description of the nucleus shape is obtained from molecular dynamics simulations. The computed γ shows a linear dependence on the temperature. The values of γ extrapolated to the melting temperature agree well with previous data obtained by the capillary fluctuation method. Using the temperature dependence of γ, we estimate the nucleation free energy barrier in a wide temperature range from the classical nucleation theory. The obtained data agree very well with the results from the brute-force molecular dynamics simulations.
△ Less
Submitted 27 October, 2018;
originally announced October 2018.
-
Effect of Samarium doping on the nucleation of fcc-Aluminum in undercooled liquids
Authors:
Yang Sun,
Feng Zhang,
Lin Yang,
Zhuo Ye,
Huajing Song,
Mikhail I. Mendelev,
Cai-Zhuang Wang,
Kai-Ming Ho
Abstract:
The effect of Sm doping on the fcc-Al nucleation was investigated in Al-Sm liquids with low Sm concentrations (xSm) with molecular dynamics simulations. The nucleation in the moderately undercooled liquid is achieved by the recently developed persistent-embryo method. Systematically computing the nucleation rate with different xSm (xSm=0%, 1%, 2%, 3%, 5%) at 700 K, we found Sm dopant reduces the n…
▽ More
The effect of Sm doping on the fcc-Al nucleation was investigated in Al-Sm liquids with low Sm concentrations (xSm) with molecular dynamics simulations. The nucleation in the moderately undercooled liquid is achieved by the recently developed persistent-embryo method. Systematically computing the nucleation rate with different xSm (xSm=0%, 1%, 2%, 3%, 5%) at 700 K, we found Sm dopant reduces the nucleation rate by up to 25 orders of magnitudes with only 5% doping concentration. This effect is mostly associated with the increase in the free energy barrier with a minor contribution from suppression of the attachment to the nucleus caused by Sm doping.
△ Less
Submitted 2 May, 2018;
originally announced May 2018.
-
Effects of Grain Boundary Disorder on Yield Strength
Authors:
Valery Borovikov,
Mikhail I. Mendelev,
Alexander H. King
Abstract:
It was recently reported that segregation of Zr to grain boundaries (GB) in nanocrystalline Cu can lead to the formation of disordered intergranular films [1,2]. In this study we employ atomistic computer simulations to study how the formation of these films affects the dislocation nucleation from the GBs. We found that full disorder of the grain boundary structure leads to the suppression of disl…
▽ More
It was recently reported that segregation of Zr to grain boundaries (GB) in nanocrystalline Cu can lead to the formation of disordered intergranular films [1,2]. In this study we employ atomistic computer simulations to study how the formation of these films affects the dislocation nucleation from the GBs. We found that full disorder of the grain boundary structure leads to the suppression of dislocation emission and significant increase of the yield stress. Depending on the solute concentration and heat-treatment, however, a partial disorder may also occur and this aids dislocation nucleation rather than suppressing it, resulting in elimination of the strengthening effect.
△ Less
Submitted 22 February, 2018;
originally announced February 2018.
-
Observation of Al41Sm5: new evidence for structural hierarchy in Al-Sm alloys
Authors:
Z. Ye,
F. Meng,
F. Zhang,
Y. Sun,
L. Yang,
S. H. Zhou,
R. E. Napolitano,
M. I. Mendelev,
R. T. Ott,
M. J. Kramer,
C. Z. Wang,
K. M. Ho
Abstract:
Using an effective genetic algorithm, we uncover the structure of a metastable Al41Sm5 phase that supplements its family sharing similar short-range orders. The phase evolves upon heating an amorphous Al-9.7at.%Sm ribbon, produced by melt-spinning. The dynamical phase selection is discussed with respect to the structural connections between the short-range packing motifs in the amorphous precursor…
▽ More
Using an effective genetic algorithm, we uncover the structure of a metastable Al41Sm5 phase that supplements its family sharing similar short-range orders. The phase evolves upon heating an amorphous Al-9.7at.%Sm ribbon, produced by melt-spinning. The dynamical phase selection is discussed with respect to the structural connections between the short-range packing motifs in the amorphous precursor and those observed in the selected phases. The phase elucidated here is one of several newly discovered large-unit-cell phases found to form during devitrification from the glass in this binary system, further illustrating the power and efficiency of our approach, the important role of structural hierarchy in phase selection, and the richness of the metastable phase landscape accessible from the glassy structure.
△ Less
Submitted 27 December, 2017;
originally announced December 2017.
-
Effect of Ag and Zr solutes on dislocation emission from the Σ11(332)[110] symmetric tilt grain boundary in fcc Cu
Authors:
Valery Borovikov,
Mikhail I. Mendelev,
Alexander H. King
Abstract:
Addition of solutes is commonly used to stabilize nanocrystalline materials against grain growth. However, segregating at grain boundaries, these solutes also affect the process of dislocation nucleation from grain boundaries under applied stress. Using atomistic simulations we demonstrate that the effect of solutes on the dislocation nucleation strongly depends on the distribution of solutes at t…
▽ More
Addition of solutes is commonly used to stabilize nanocrystalline materials against grain growth. However, segregating at grain boundaries, these solutes also affect the process of dislocation nucleation from grain boundaries under applied stress. Using atomistic simulations we demonstrate that the effect of solutes on the dislocation nucleation strongly depends on the distribution of solutes at the grain boundary, which can vary dramatically depending on the solute type. In particular, our results indicate that the solutes with a smaller size mismatch can be more effective in suppressing dislocation emission from grain boundaries. Bearing in mind that dislocation slip originating from grain boundaries or their triple junctions is the dominant mechanism of plastic deformation when grain sizes are reduced to the nanoscale, we emphasize the importance of the search for the optimal solute additions, which would stabilize the nanocrystalline material against grain growth and, at the same time, effectively suppress the dislocation nucleation from the grain boundaries.
△ Less
Submitted 20 December, 2017;
originally announced December 2017.
-
A self-contained algorithm for determination of solid-liquid equilibria in an alloy system
Authors:
Lin Yang,
Yang Sun,
Zhuo Ye,
Feng Zhang,
Mikhail I. Mendelev,
Cai-Zhuang Wang,
Kai-Ming Ho
Abstract:
We describe a self-contained procedure to evaluate the free energy of liquid and solid phases of an alloy system. The free energy of a single-element solid phase is calculated with thermodynamic integration using the Einstein crystal as the reference system. Then, free energy difference between the solid and liquid phases is calculated by Gibbs-Duhem integration. The central part of our method is…
▽ More
We describe a self-contained procedure to evaluate the free energy of liquid and solid phases of an alloy system. The free energy of a single-element solid phase is calculated with thermodynamic integration using the Einstein crystal as the reference system. Then, free energy difference between the solid and liquid phases is calculated by Gibbs-Duhem integration. The central part of our method is the construction of a reversible alchemical path connecting a pure liquid and a liquid alloy to calculate the mixing enthalpy and entropy. We have applied the method to calculate the free energy of solid and liquid phases in the Al-Sm system. The driving force for fcc-Al nucleation in Al-Sm liquid and the melting curve for fcc-Al and Al3Sm are also calculated.
△ Less
Submitted 4 December, 2017;
originally announced December 2017.
-
Nucleation of stoichiometric compounds from liquid: the role of the kinetic factor
Authors:
H. Song,
Y. Sun,
F. Zhang,
C. Z. Wang,
K. M. Ho,
M. I. Mendelev
Abstract:
While the role of the free energy barrier during nucleation is a text-book subject the importance of the kinetic factor is frequently underestimated. We obtained both quantities from molecular dynamics (MD) simulations for the pure Ni and B2 phases in the Ni50Al50 and Cu50Zr50 alloys. The free-energy barrier was found to be higher in Ni but the nucleation rate is much lower in the Ni50Al50 alloy w…
▽ More
While the role of the free energy barrier during nucleation is a text-book subject the importance of the kinetic factor is frequently underestimated. We obtained both quantities from molecular dynamics (MD) simulations for the pure Ni and B2 phases in the Ni50Al50 and Cu50Zr50 alloys. The free-energy barrier was found to be higher in Ni but the nucleation rate is much lower in the Ni50Al50 alloy which was attributed to the ordered nature of the B2 phase. Since the Cu50Zr50 B2 phase can has even smaller fraction of the anti-site defects its nucleation is never observed in the MD simulation.
△ Less
Submitted 27 November, 2017;
originally announced November 2017.
-
Solute effects on interfacial dislocation emission in nanomaterials: nucleation site competition and neutralization
Authors:
Valery Borovikov,
Mikhail I. Mendelev,
Alexander H. King
Abstract:
Interfacial nucleation is the dominant process of dislocation generation during the plastic deformation of nano-crystalline materials. Solute additions intended to stabilize nano-crystalline metals against grain growth, may segregate to the grain boundaries and triple junctions where they can affect the process of the dislocation emission. In this Letter we demonstrate that the effect of solute ad…
▽ More
Interfacial nucleation is the dominant process of dislocation generation during the plastic deformation of nano-crystalline materials. Solute additions intended to stabilize nano-crystalline metals against grain growth, may segregate to the grain boundaries and triple junctions where they can affect the process of the dislocation emission. In this Letter we demonstrate that the effect of solute addition in a nano-crystalline material containing competing solute segregation sites and dislocation sources can be very complex due to different rates of segregation at different interfaces. Moreover, at large concentrations, when the solutes form clusters near the grain boundaries or triple junctions, the interfaces between these clusters and the matrix can introduce new dislocation emission sources, which can be activated under lower applied stress. Thus, the strength maximum can occur at a certain solute concentration: adding solutes beyond this optimal solute concentration can reduce the strength of the material.
△ Less
Submitted 25 October, 2017;
originally announced October 2017.
-
Overcoming the time limitation in Molecular Dynamics simulation of crystal nucleation: a persistent-embryo approach
Authors:
Yang Sun,
Huajing Song,
Feng Zhang,
Lin Yang,
Zhuo Ye,
Mikhail I. Mendelev,
Cai-Zhuang Wang,
Kai-Ming Ho
Abstract:
The crystal nucleation from liquid in most cases is too rare to be accessed within the limited timescales of the conventional molecular dynamics (MD) simulation. Here, we developed a "persistent embryo" method to facilitate crystal nucleation in MD simulations by preventing small crystal embryos from melting using external spring forces. We applied this method to the pure Ni case for a moderate un…
▽ More
The crystal nucleation from liquid in most cases is too rare to be accessed within the limited timescales of the conventional molecular dynamics (MD) simulation. Here, we developed a "persistent embryo" method to facilitate crystal nucleation in MD simulations by preventing small crystal embryos from melting using external spring forces. We applied this method to the pure Ni case for a moderate undercooling where no nucleation can be observed in the conventional MD simulation, and obtained nucleation rate in good agreement with the experimental data. Moreover, the method is applied to simulate an even more sluggish event: the nucleation of the B2 phase in a strong glass-forming Cu-Zr alloy. The nucleation rate was found to be 8 orders of magnitude smaller than Ni at the same undercooling, which well explains the good glass formability of the alloy. Thus, our work opens a new avenue to study solidification under realistic experimental conditions via atomistic computer simulation.
△ Less
Submitted 24 February, 2018; v1 submitted 31 August, 2017;
originally announced September 2017.
-
Development of an Interatomic Potential for the Simulation of Defects, Plasticity and Phase Transformations in Titanium
Authors:
M. I. Mendelev,
T. L. Underwood,
G. J. Ackland
Abstract:
New interatomic potentials describing defects, plasticity and high temperature phase transitions for Ti are presented. Fitting the martensitic hcp-bcc phase transformation temperature requires an efficient and accurate method to determine it. We apply a molecular dynamics (MD) method based on determination of the melting temperature of competing solid phases, and Gibbs-Helmholtz integration, and a…
▽ More
New interatomic potentials describing defects, plasticity and high temperature phase transitions for Ti are presented. Fitting the martensitic hcp-bcc phase transformation temperature requires an efficient and accurate method to determine it. We apply a molecular dynamics (MD) method based on determination of the melting temperature of competing solid phases, and Gibbs-Helmholtz integration, and a lattice-switch Monte Carlo method (LSMC): these agree on the hcp-bcc transformation temperatures to within 2 K. We were able to develop embedded atom potentials which give a good fit to either low or high temperature data, but not both. The first developed potential (Ti1) reproduces the hcp-bcc transformation and melting temperatures and is suitable for the simulation of phase transitions and bcc Ti. Two other potentials (Ti2 and Ti3) correctly describe defect properties, and can be used to simulate plasticity or radiation damage in hcp Ti. The fact that a single EAM potential cannot describe both low and high temperature phases may be attributed to neglect of electronic degrees of freedom, notably bcc has a much higher electronic entropy. A temperature-dependent potential obtained from the combination of potentials Ti1 and Ti2 may be used to simulate Ti properties at any temperature.
△ Less
Submitted 9 December, 2016;
originally announced December 2016.
-
Unconventional phase selection in high-driven systems: A complex metastable structure prevails over simple stable phases
Authors:
Zhuo Ye,
Feng Zhang,
Yang Sun,
Manh Cuong Nguyen,
Mikhail I. Mendelev,
Ryan T. Ott,
Eun Soo Park,
Matthew Besser,
Matthew J. Kramer,
Ze-Jun Ding,
Cai-Zhuang Wang,
Kai-Ming Ho
Abstract:
Phase selection in deeply undercooled liquids and devitrified glasses during heating involves complex interplay between the barriers to nucleation and the ability for these nuclei to grow. During the devitrification of glassy alloys, complicated metastable structures often precipitate instead of simpler, more stable compounds. Here, we access this unconventional type of phase selections by investi…
▽ More
Phase selection in deeply undercooled liquids and devitrified glasses during heating involves complex interplay between the barriers to nucleation and the ability for these nuclei to grow. During the devitrification of glassy alloys, complicated metastable structures often precipitate instead of simpler, more stable compounds. Here, we access this unconventional type of phase selections by investigating an Al-10%Sm system, where a complicated cubic structure first precipitates with a large lattice parameter of 1.4 nm. We not only solve the structure of this "big cubic" phase containing ~140 atoms but establish an explicit interconnection between the structural orderings of the amorphous alloy and the cubic phase, which provides a low-barrier nucleation pathway at low temperatures. The surprising rapid growth of the crystal is attributed to its high tolerance to point defects, which minimize the short-scale atomic rearrangements to form the crystal. Our study suggests a new scenario of devitrification, where phase transformation proceeds initially without partitioning through a complex intermediate crystal phase.
△ Less
Submitted 30 January, 2015;
originally announced February 2015.
-
Effect of Stacking Fault Energy on Mechanism of Plastic Deformation in Nanotwinned FCC Metals
Authors:
Valery Borovikov,
Mikhail I. Mendelev,
Alexander H. King,
Richard LeSar
Abstract:
Starting from a semi-empirical potential designed for Cu, we developed a series of potentials that provide essentially constant values of all significant (calculated) materials properties except for the intrinsic stacking fault energy, which varies over a range that encompasses the lowest and highest values observed in nature. These potentials were employed in molecular dynamics (MD) simulations t…
▽ More
Starting from a semi-empirical potential designed for Cu, we developed a series of potentials that provide essentially constant values of all significant (calculated) materials properties except for the intrinsic stacking fault energy, which varies over a range that encompasses the lowest and highest values observed in nature. These potentials were employed in molecular dynamics (MD) simulations to investigate how stacking fault energy affects the mechanical behavior of nanotwinned face-centered cubic (fcc) materials. The results indicate that properties such as yield strength and microstructural stability do not vary systematically with stacking fault energy, but rather fall into two distinct regimes corresponding to low and high stacking fault energies.
△ Less
Submitted 25 January, 2015;
originally announced January 2015.
-
Crystalline 'Genes' in Metallic Liquids
Authors:
Yang Sun,
Feng Zhang,
Zhuo Ye,
Xiaowei Fang,
Zejun Ding,
Cai-Zhuang Wang,
Mikhail I. Mendelev,
Ryan T. Ott,
M. J. Kramer,
Kai-Ming Ho
Abstract:
The underlying structural order that transcends the liquid, glass and crystalline states is identified using an efficient genetic algorithm (GA). GA identifies the most common energetically favorable packing motif in crystalline structures close to the alloy's Al-10 at.% Sm composition. These motifs are in turn compared to the observed packing motifs in the actual liquid structures using a cluster…
▽ More
The underlying structural order that transcends the liquid, glass and crystalline states is identified using an efficient genetic algorithm (GA). GA identifies the most common energetically favorable packing motif in crystalline structures close to the alloy's Al-10 at.% Sm composition. These motifs are in turn compared to the observed packing motifs in the actual liquid structures using a cluster-alignment method which reveals the average topology. Conventional descriptions of the short-range order, such as Voronoi tessellation, are too rigid in their analysis of the configurational poly-types when describing the chemical and topological ordering during transition from undercooled metallic liquids to crystalline phases or glass. Our approach here brings new insight into describing mesoscopic order-disorder transitions in condensed matter physics.
△ Less
Submitted 7 August, 2014;
originally announced August 2014.
-
Molecular Dynamics Study of Self-Diffusion in Zr
Authors:
M. I. Mendelev,
B. S. Bokstein
Abstract:
We employed a recently developed semi-empirical Zr potential to determine the diffusivities in the hcp and bcc Zr via molecular dynamics simulation. The point defect concentration was determined directly from MD simulation rather than from theoretical methods using T=0 calculations. We found that the diffusion proceeds via the interstitial mechanism in the hcp Zr and both the vacancy and interst…
▽ More
We employed a recently developed semi-empirical Zr potential to determine the diffusivities in the hcp and bcc Zr via molecular dynamics simulation. The point defect concentration was determined directly from MD simulation rather than from theoretical methods using T=0 calculations. We found that the diffusion proceeds via the interstitial mechanism in the hcp Zr and both the vacancy and interstitial mechanisms give contribution in diffusivity in the bcc Zr. The agreement with the experimental data is excellent for the hcp Zr and for the bcc Zr it is rather good at high temperatures but there is a considerable disagreement at low temperatures.
△ Less
Submitted 8 February, 2009;
originally announced February 2009.
-
Asperity contacts at the nanoscale: comparison of Ru and Au
Authors:
Andrea Fortini,
Mikhail I. Mendelev,
Sergey Buldyrev,
David Srolovitz
Abstract:
We develop and validate an interatomic potential for ruthenium based on the embedded atom method framework with the Finnis/Sinclair representation. We confirm that the new potential yields a stable hcp lattice with reasonable lattice and elastic constants and surface and stacking fault energies. We employ molecular dynamics simulations to bring two surfaces together; one flat and the other with…
▽ More
We develop and validate an interatomic potential for ruthenium based on the embedded atom method framework with the Finnis/Sinclair representation. We confirm that the new potential yields a stable hcp lattice with reasonable lattice and elastic constants and surface and stacking fault energies. We employ molecular dynamics simulations to bring two surfaces together; one flat and the other with a single asperity. We compare the process of asperity contact formation and breaking in Au and Ru, two materials currently in use in micro electro mechanical system switches. While Au is very ductile at 150 and 300 K, Ru shows considerably less plasticity at 300 and 600 K (approximately the same homologous temperature). In Au, the asperity necks down to a single atom thick bridge at separation. While similar necking occurs in Ru at 600 K, it is much more limited than in Au. On the other hand, at 300 K, Ru breaks by a much more brittle process of fracture/decohesion with limited plastic deformation.
△ Less
Submitted 3 July, 2008;
originally announced July 2008.
-
Solid-Liquid Phase Diagrams for Binary Metallic Alloys: Adjustable Interatomic Potentials
Authors:
H. -S. Nam,
M. I. Mendelev,
D. J. Srolovitz
Abstract:
We develop a new approach to determining LJ-EAM potentials for alloys and use these to determine the solid-liquid phase diagrams for binary metallic alloys using Kofke's Gibbs-Duhem integration technique combined with semigrand canonical Monte Carlo simulations. We demonstrate that it is possible to produce a wide-range of experimentally observed binary phase diagrams (with no intermetallic phas…
▽ More
We develop a new approach to determining LJ-EAM potentials for alloys and use these to determine the solid-liquid phase diagrams for binary metallic alloys using Kofke's Gibbs-Duhem integration technique combined with semigrand canonical Monte Carlo simulations. We demonstrate that it is possible to produce a wide-range of experimentally observed binary phase diagrams (with no intermetallic phases) by reference to the atomic sizes and cohesive energies of the two elemental materials. In some cases, it is useful to employ a single adjustable parameter to adjust the phase diagram (we provided a good choice for this free parameter). Next, we perform a systematic investigation of the effect of relative atomic sizes and cohesive energies of the elements on the binary phase diagrams. We then show that this approach leads to good agreement with several experimental binary phase diagrams. The main benefit of this approach is not the accurately reproduction of experimental phase diagrams, but rather to provide a method by which material properties can be continuously changed in simulations studies. This is one of the keys to the use of atomistic simulations to understand mechanisms and properties in a manner not available to experiment.
△ Less
Submitted 20 November, 2006; v1 submitted 2 August, 2006;
originally announced August 2006.
-
Development of an interatomic potential for phosphorus impurities in iron
Authors:
G. J. Ackland,
M. I. Mendelev,
D. Srolovitz,
S. Han,
A. V. Barashev
Abstract:
We present the derivation of an interatomic potential for the iron phosphorus system based primarily on {\it ab initio} data. Transferrability in this system is extremely problematic, and the potential is intended specifically to address the problem of radiation damage and point defects in iron containing low concentrations of phosphorus atoms. Some preliminary molecular dynamics calculations sh…
▽ More
We present the derivation of an interatomic potential for the iron phosphorus system based primarily on {\it ab initio} data. Transferrability in this system is extremely problematic, and the potential is intended specifically to address the problem of radiation damage and point defects in iron containing low concentrations of phosphorus atoms. Some preliminary molecular dynamics calculations show that P strongly affects point defect migration.
△ Less
Submitted 16 June, 2004;
originally announced June 2004.