A Simple Palladium Hydride Embedded Atom Method Potential for Hydrogen Energy Applications

2019 ◽  
Vol 141 (6) ◽  
Author(s):  
Iyad Hijazi ◽  
Yang Zhang ◽  
Robert Fuller
Keyword(s):  
Elevated Temperatures ◽  
Clean Energy ◽  
Md Simulations ◽  
Miscibility Gap ◽  
Embedded Atom Method ◽  
Hydrogen Energy ◽  
Hydrogen Purification ◽  
Palladium Hydride ◽  
Embedded Atom ◽  
Eam Potential

When hydrogen is produced from a biomass or coal gasifier, it is necessary to purify it from syngas streams containing components such as CO, CO2, N2, CH4, and other products. Therefore, a challenge related to hydrogen purification is the development of hydrogen-selective membranes that can operate at elevated temperatures and pressures, provide high fluxes, long operational lifetime, and resistance to poisoning while still maintaining reasonable cost. Palladium-based membranes have been shown to be well suited for these types of high-temperature applications and have been widely utilized for hydrogen separation. Palladium's unique ability to absorb a large quantity of hydrogen can also be applied in various clean energy technologies, like hydrogen fuel cells. In this paper, a fully analytical interatomic embedded atom method (EAM) potential for the Pd-H system has been developed, that is easily extendable to ternary Palladium-based hydride systems, such as Pd-Cu-H and Pd-Ag-H. The new potential has fewer fitting parameters than previously developed EAM Pd-H potentials and is able to accurately predict the cohesive energy, lattice constant, bulk modulus, elastic constants, melting temperature, and the stable Pd-H structures in molecular dynamics (MD) simulations with various hydrogen concentrations. The EAM potential also well predicts the miscibility gap, the segregation of the palladium hydride system into dilute (α), and concentrated (β) phases.

A Simple EAM Potential for Hydrogen-Selective Palladium Based Membranes for Biomass Derived Syngas Processing

2018 ◽  
Author(s):  
Iyad Hijazi ◽  
Yang Zhang ◽  
Robert Fuller
Keyword(s):  
Elevated Temperatures ◽  
Clean Energy ◽  
Md Simulations ◽  
Miscibility Gap ◽  
Hydrogen Purification ◽  
Embedded Atom ◽  
Energy Technologies ◽  
Renewable Feedstock ◽  
Operational Lifetime ◽  
Eam Potential

Biomass offers the potential to economically produce hydrogen via gasification from an abundant and renewable feedstock. When hydrogen is produced from a biomass gasifier, it is necessary to purify it from syngas streams containing components such as CO, CO2, N2, CH4, and other products. Therefore, a challenge related to hydrogen purification is the development of hydrogen-selective membranes that can operate at elevated temperatures and pressures, provide high fluxes, long operational lifetime, and resistance to poisoning while still maintaining reasonable cost. Palladium based membranes have been shown to be well suited for these types of high-temperature applications and have been widely utilized for hydrogen separation. Palladium’s unique ability to absorb a large quantity of hydrogen can also be applied in various clean energy technologies, like hydrogen fuel cells. In this paper, a fully analytical interatomic Embedded Atom Potential (EAM) for the Pd-H system has been developed, that is easily extendable to ternary Palladium based hydride systems such as Pd-Cu-H and Pd-Ag-H. The new potential has fewer fitting parameters than previously developed EAM Pd-H potentials and is able to accurately predict the cohesive energy, lattice constant, bulk modulus, elastic constants, melting temperature, and the stable Pd-H structures in molecular dynamics (MD) simulations with various hydrogen concentrations. The EAM potential also well predicts the miscibility gap, the segregation of the palladium hydride system into dilute (α) and concentrated (β) phases.

scholarly journals Quaternary Hydrides Pd1-y-zAgyCuzHx Embedded Atom Method Potentials for Hydrogen Energy Applications

2021 ◽  
Vol 3 (1) ◽  
Author(s):  
Chaonan Zhang ◽  
Robert Fuller ◽  
Iyad Hijazi ◽  
Keyword(s):  
Mechanical Properties ◽  
Carbon Monoxide Poisoning ◽  
Hydrogen Permeability ◽  
Sulfur Poisoning ◽  
Ternary Alloys ◽  
Miscibility Gap ◽  
Embedded Atom Method ◽  
Hydrogen Energy ◽  
Energy Applications ◽  
Embedded Atom

The Pd-H system has attracted extensive attention. Pd can absorb considerable amount of H at room temperature, this ability is reversible, so it is suitable for multiple energy applications. Pd-Ag alloys possess higher H permeability, solubility and narrower miscibility gap with better mechanical properties than pure Pd, but sulfur poisoning remains an issue. Pd-Cu alloys have excellent resistance to sulfur and carbon monoxide poisoning and hydrogen embrittlement, good mechanical properties, and broader temperature working environments over pure Pd, but relatively lower hydrogen permeability and solubility than pure Pd and Pd-Ag alloys. This suggests that alloying Pd with Ag and Cu to create Pd-Ag-Cu ternary alloys can optimize the overall performance and substantially lowers the cost. Thus, in this paper, we provide the first embedded atom method potentials for the quaternary hydrides Pd1-y-zAgyCuzHx. The fully analytical potentials are fitted utilizing the central atom method without performing time-consuming molecular dynamics simulations.

An embedded-atom method interatomic potential for Pd–H alloys

2008 ◽  
Vol 23 (3) ◽  
pp. 704-718 ◽  
Author(s):  
X.W. Zhou ◽  
J.A. Zimmerman ◽  
B.M. Wong ◽  
J.J. Hoyt
Keyword(s):  
Computational Models ◽  
Mechanical Response ◽  
Hydrogen Diffusion ◽  
Interatomic Potential ◽  
Diffusion Mechanism ◽  
Miscibility Gap ◽  
Embedded Atom Method ◽  
Embedded Atom ◽  
Formidable Challenge ◽  
Embedded Atom Method Potential

Palladium hydrides have important applications. However, the complex Pd–H alloy system presents a formidable challenge to developing accurate computational models. In particular, the separation of a Pd–H system to dilute (α) and concentrated (β) phases is a central phenomenon, but the capability of interatomic potentials to display this phase miscibility gap has been lacking. We have extended an existing palladium embedded-atom method potential to construct a new Pd–H embedded-atom method potential by normalizing the elemental embedding energy and electron density functions. The developed Pd–H potential reasonably well predicts the lattice constants, cohesive energies, and elastic constants for palladium, hydrogen, and PdHx phases with a variety of compositions. It ensures the correct hydrogen interstitial sites within the hydrides and predicts the phase miscibility gap. Preliminary molecular dynamics simulations using this potential show the correct phase stability, hydrogen diffusion mechanism, and mechanical response of the Pd–H system.

scholarly journals Molecular Dynamics Study on Deformation Mechanism of Grain Boundaries in Magnesium Crystal: Based on Coincidence Site Lattice Theory

2018 ◽  
Vol 2018 ◽  
pp. 1-10 ◽  
Author(s):  
Ken-ichi Saitoh ◽  
Kohei Kuramitsu ◽  
Tomohiro Sato ◽  
Masanori Takuma ◽  
Yoshimasa Takahashi
Keyword(s):  
Molecular Dynamics ◽  
Deformation Mechanism ◽  
Coincidence Site Lattice ◽  
Md Simulations ◽  
Uniaxial Tensile ◽  
High Angle ◽  
Site Lattice ◽  
Embedded Atom ◽  
Significant Difference ◽  
Eam Potential

As for magnesium (Mg) alloys, it has been noted that they are inferior to plastic deformation, but improvement in the mechanical properties by further refinement of grain size has been recently suggested. It means the importance of atomistic view of polycrystalline interface of Mg crystal. In this study, to discuss the deformation mechanism of polycrystalline Mg, atomistic grain boundary (GB) models by using coincidence site lattice (CSL) theory are constructed and are simulated for their relaxed and deformatted structures. First, GB structures in which the axis of rotation is in [11¯00] direction are relaxed at 10 Kelvin, and the GB energies are evaluated. Then, the deformation mechanism of each GB model under uniaxial tensile loading is observed by using the molecular dynamics (MD) method. The present MD simulations are based on embedded atom method (EAM) potential for Mg crystal. As a result, we were able to observe atomistically a variety of GB structures and to recognize significant difference in deformation mechanism between low-angle GBs and high-angle GBs. A close scrutiny is made on phenomena of dislocation emission processes peculiar to each atomistic local structure in high-angle GBs.

scholarly journals Graphene Adhesion Mechanics on Iron Substrates: Insight from Molecular Dynamic Simulations

Crystals ◽  
2019 ◽  
Vol 9 (11) ◽  
pp. 579 ◽  
Author(s):  
Wang ◽  
Jin ◽  
Yang ◽  
Zong ◽  
Peng
Keyword(s):  
Molecular Dynamic ◽  
Surface Engineering ◽  
Graphene Nanosheets ◽  
Pair Potential ◽  
Md Simulations ◽  
Adhesion Energy ◽  
Dynamic Simulations ◽  
Embedded Atom ◽  
Adhesion Mechanics ◽  
Eam Potential

The adhesion feature of graphene on metal substrates is important in graphene synthesis, transfer and applications, as well as for graphene-reinforced metal matrix composites. We investigate the adhesion energy of graphene nanosheets (GNs) on iron substrate using molecular dynamic (MD) simulations. Two Fe–C potentials are examined as Lennard–Jones (LJ) pair potential and embedded-atom method (EAM) potential. For LJ potential, the adhesion energies of monolayer GN are 0.47, 0.62, 0.70 and 0.74 J/m2 on the iron {110}, {111}, {112} and {100} surfaces, respectively, compared to the values of 26.83, 24.87, 25.13 and 25.01 J/m2 from EAM potential. When the number of GN layers increases from one to three, the adhesion energy from EAM potential increases. Such a trend is not captured by LJ potential. The iron {110} surface is the most adhesive surface for monolayer, bilayer and trilayer GNs from EAM potential. The results suggest that the LJ potential describes a weak bond of Fe–C, opposed to a hybrid chemical and strong bond from EAM potential. The average vertical distances between monolayer GN and four iron surfaces are 2.0–2.2 Å from LJ potential and 1.3–1.4 Å from EAM potential. These separations are nearly unchanged with an increasing number of layers. The ABA-stacked GN is likely to form on lower-index {110} and {100} surfaces, while the ABC-stacked GN is preferred on higher-index {111} surface. Our insights of the graphene adhesion mechanics might be beneficial in graphene growing, surface engineering and enhancement of iron using graphene sheets.

The Infidence of Defect Concentrations on Migration Energies in AgZn Alloys

1988 ◽  
Vol 128 ◽  
Author(s):  
T. D. Andreadis ◽  
M. Rosen ◽  
J. M. Eridon ◽  
D. J. Rosen
Keyword(s):  
Linear Approximation ◽  
Embedded Atom Method ◽  
Local Defect ◽  
Migration Energy ◽  
New Model ◽  
Impurity Atoms ◽  
Embedded Atom ◽  
Defect Reactions ◽  
Eam Potential ◽  
Formation Energies

ABSTRACTMigration energies in Ag of vacancies, interstitials, Zn impurity atoms, interstitial-iipurity cumplexes, and vacancy-impurity complexes were calculated using Embedded Atom. Method (EAM) potentials in Molecular Statics calculations. A new Zn EAM potential was determined and used in these calculations. The dependence of the migration energies on local defect concentrations was determined in a linear approximation. Binding and formation energies of defects are also presented. A new model for the migration energy appropriate for defect reactions is introduced.

scholarly journals A numerical fitting procedure for the embedded atom method interatomic potential and a bridged finite element-molecular dynamics method for large atomic systems

2021 ◽  
Author(s):  
Karthik Narayan
Keyword(s):  
Molecular Dynamics ◽  
Md Simulation ◽  
Embedded Atom Method ◽  
Discretization Scheme ◽  
Element Discretization ◽  
Atomic Systems ◽  
Embedded Atom ◽  
Fitting Procedure ◽  
Numerical Fitting ◽  
Eam Potential

This thesis presents a powerful numerical fitting procedure for generating Embedded Atom Method (EAM) inter-atomic potentials for pure Face Centred Cubic (FCC) and Body Centred Cubic (BCC) metals. The numerical fitting procedure involves assuming a reasonable parameterized form for a portion of the EAM potential, and then fitting the remaining portion to select thermal and elastic properties of the metal. Molecular Dynamics (MD) simulation is used to effect the fitting procedure. The procedure is used to generate an EAM potential for copper, an FCC metal. This resulting EAM potential is used to conduct MD simulations of perfect copper crystals containing voids of different geometries. Following this, a bridged Finite Element-Molecular Dynamics (FE-MD) method is presented, which can be used to simulate large atomic systems much more efficiently than MD simulation alone. The method implements a novel element discretization scheme proposed by the author that is so general that it can be applied to any system of objects interacting with each other via any potential (simple or complex, EAM or otherwise). This bridged FE-MD method is used to reanalyze the voids in the copper crystal lattice. The resulting virial stress increment patterns are found to agree remarkably with the earlier MD simulation results. Furthermore, the bridged FE-MD method is much quicker than the pure MD simulation. These two facts prove the power and usefulness of the bridged FE-MD method, and validate the proposed element discretization scheme

Grain Boundary Migration in Nanocrystalline Iron

2007 ◽  
Vol 129 ◽  
pp. 145-150 ◽  
Author(s):  
Tomasz Wejrzanowski ◽  
M. Spychalski ◽  
Roman Pielaszek ◽  
Krzysztof Jan Kurzydlowski
Keyword(s):  
Boundary Migration ◽  
3D Models ◽  
Embedded Atom Method ◽  
Interface Motion ◽  
Embedded Atom ◽  
The Matrix ◽  
Nanocrystalline Iron ◽  
Eam Potential ◽  
The Relationship ◽  
Misorientation Angles

In this study a series of 3D models for curved [100] grain boundaries (GBs) in pure α-iron have been constructed. Each model consisted of a spherical grain, with an initial size of about 9 nm, surrounded by a large single-crystal. Different orientations have been assigned to the grain and the matrix in order to obtain interfaces with misorientation angles in the range of 5-45 degrees in steps of 5 degrees. The molecular dynamics with Embedded Atom Method (EAM) potential have been performed for investigation of the temporal changes in GB migration and grain rotations at temperature of 1000 K. The relationship between GB misorientation and its mobility has been found. It was also discovered that the density of the material decreases with a reduction of GB area. The effect of a triple junction on the interface motion has been also studied by introducing a bi-crystal matrix instead of a singlecrystal one. The results are discussed in terms of grain growth investigations in nanometals.

Strengthening Behavior and Tension–Compression Strength–Asymmetry in Nanocrystalline Metal–Ceramic Composites

2012 ◽  
Vol 134 (4) ◽  
Author(s):  
A. M. Dongare ◽  
B. LaMattina ◽  
A. M. Rajendran
Keyword(s):  
Molecular Dynamics ◽  
Volume Fraction ◽  
Md Simulations ◽  
Ceramic Composites ◽  
Embedded Atom Method ◽  
Nanocrystalline Metal ◽  
Rich Region ◽  
Embedded Atom ◽  
Metal Ceramic ◽  
Strength Asymmetry

Metal–ceramic composites are an emerging class of materials for use in the next-generation high technology applications due to their ability to sustain plastic deformation and resist failure in extreme mechanical environments. Large scale molecular dynamics simulations are used to investigate the performance of nanocrystalline metal–matrix composites (MMCs) formed by the reinforcement of the nanocrystalline Al matrix with a random distribution of nanoscale ceramic particles. The interatomic interactions are defined by the newly developed angular-dependent embedded atom method (A-EAM) by combining the embedded atom method (EAM) potential for Al with the Stillinger–Weber (SW) potential for Si in one functional form. The molecular dynamics (MD) simulations are aimed to investigate the strengthening behavior and the tension–compression strength asymmetry of these composites as a function of volume fraction of the reinforcing Si phase. MD simulations suggest that the strength of the nanocomposite increases linearly with an increase in the volume fraction of Si in the Al-rich region, whereas the increase is very sharp in the Si-rich region. The higher strength of the nanocomposite is attributed to the reduced sliding/rotation between the Al/Si and the Si/Si grains as compared to the pure nanocrystalline metal.

Cascade Overlap in hcp Zirconium: Defect Accumulation and Microstructure Evolution with Radiation using Molecular Dynamics Simulations

2013 ◽  
Vol 1514 ◽  
pp. 37-42 ◽  
Author(s):  
Prithwish K. Nandi ◽  
Jacob Eapen
Keyword(s):  
Molecular Dynamics ◽  
Molecular Dynamics Simulations ◽  
Microstructure Evolution ◽  
Stacking Faults ◽  
Embedded Atom Method ◽  
Shockley Partial Dislocation ◽  
Embedded Atom ◽  
Defect Accumulation ◽  
Dynamics Simulations ◽  
Eam Potential

ABSTRACTMolecular dynamics simulations are performed to investigate the defect accumulation and microstructure evolution in hcp zirconium (Zr) – a material which is widely used as clad for nuclear fuel. Cascades are generated with a 3 keV primary knock-on atom (PKA) using an embedded atom method (EAM) potential with interactions modified for distances shorter than 0.1 Å. With sequential cascade simulations we show the emergence of stacking faults both in the basal and prism planes, and a Shockley partial dislocation on the basal plane.

Sign in / Sign up

Export Citation Format

Share Document