Xiaowang Zhou

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

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.

Molecular Dynamics Studies of Dislocations in CdTe Crystals from a New Bond Order Potential

Cd1-xZnxTe (CZT) crystals are the leading semiconductors for radiation detection, but their application is limited by the high cost of detector-grade materials. High crystal costs primarily result from property nonuniformity that causes low manufacturing yield. Although tremendous efforts have been made in the past to reduce Te inclusions/precipitates in CZT, this has not resulted in an anticipated improvement in material property uniformity. Moreover, it is recognized that in addition to Te particles, dislocation cells can also cause electric field perturbations and the associated property nonuniformities. Further improvement of the material, therefore, requires that dislocations in CZT crystals be understood and controlled. Here, we use a recently developed CZT bond order potential to perform representative molecular dynamics simulations to study configurations, energies, and mobilities of 29 different types of possible dislocations in CdTe (i.e., x = 1) crystals. An efficient method to derive activation free energies and activation volumes of thermally activated dislocation motion will be explored. Our focus gives insight into understanding important dislocations in the material and gives guidance toward experimental efforts for improving dislocation network structures in CZT crystals.

Influence of crystallographic orientation and anisotropy on Kapitza conductance via classical molecular dynamics simulations

We investigate the influence of crystallographic orientation and anisotropy on local phonon density of states, phonon transmissivity, and Kapitza conductance at interfaces between Lennard-Jones solids via classical molecular dynamics simulations. In agreement with prior works, we find that the Kapitza conductance at an interface between two face-centered cubic materials is independent of crystallographic orientation. On the other hand, at an interface between a face-centered cubic material and a tetragonal material, the Kapitza conductance is strongly dependent on the relative orientation of the tetragonal material, albeit this dependence is subject to the overlap in vibrational spectra of the cubic and tetragonal materials. Furthermore, we show that interactions between acoustic phonons in the cubic material and optical phonons in the tetragonal material can lead to the interface exhibiting greater “thermal anisotropy” as compared to that of the constituent materials. Finally, it is noted that the relati...

A refined parameterization of the analytical Cd–Zn–Te bond-order potential

This paper reports an updated parameterization for a CdTe bond order potential. The original potential is a rigorously parameterized analytical bond order potential for ternary the Cd–Zn–Te systems. This potential effectively captures property trends of multiple Cd, Zn, Te, CdZn, CdTe, ZnTe, and Cd1-xZnxTe phases including clusters, lattices, defects, and surfaces. It also enables crystalline growth simulations of stoichiometric compounds/alloys from non-stoichiometric vapors. However, the potential over predicts the zinc-blende CdTe lattice constant compared to experimental data. Here, we report a refined analytical Cd–Zn–Te bond order potential parameterization that predicts a better CdTe lattice constant. Characteristics of the second potential are given based on comparisons with both literature potentials and the quantum mechanical calculations.

Atomistic Potentials for Palladium-Silver Hydrides.

New embedded-atom method potentials for the ternary palladium–silver–hydrogen system are developed by extending a previously developed palladium–hydrogen potential. The ternary potentials accurately capture the heat of mixing and structural properties associated with solid solution alloys of palladium–silver. Stable hydrides are produced with properties that smoothly transition across the compositions. Additions of silver to palladium are predicted to alter the properties of the hydrides by decreasing the miscibility gap and increasing the likelihood of hydrogen atoms occupying tetrahedral interstitial sites over octahedral interstitial sites.

An analytical bond-order potential for carbon

Carbon is the most widely studied material today because it exhibits special properties not seen in any other materials when in nano dimensions such as nanotube and graphene. Reduction of material defects created during synthesis has become critical to realize the full potential of carbon structures. Molecular dynamics (MD) simulations, in principle, allow defect formation mechanisms to be studied with high fidelity, and can, therefore, help guide experiments for defect reduction. Such MD simulations must satisfy a set of stringent requirements. First, they must employ an interatomic potential formalism that is transferable to a variety of carbon structures. Second, the potential needs to be appropriately parameterized to capture the property trends of important carbon structures, in particular, diamond, graphite, graphene, and nanotubes. Most importantly, the potential must predict the crystalline growth of the correct phases during direct MD simulations of synthesis to achieve a predictive simulation of defect formation. Because an unlimited number of structures not included in the potential parameterization are encountered, the literature carbon potentials are often not sufficient for growth simulations. We have developed an analytical bond order potential for carbon, and have made it available through the public MD simulation package LAMMPS. We demonstrate that our potential reasonably captures the property trends of important carbon phases. Stringent MD simulations convincingly show that our potential accounts not only for the crystalline growth of graphene, graphite, and carbon nanotubes but also for the transformation of graphite to diamond at high pressure.

Understanding misfit strain releasing mechanisms via molecular dynamics simulations of CdTe growth on {112}zinc-blende CdS

Molecular dynamics simulations have been used to analyse microstructures of CdTe films grown on {112} surfaces of zinc-blende CdS. Interestingly, CdTe films grow in ⟨331⟩ orientations as opposed to ⟨112⟩ epitaxial orientations. At the CdTe-{331}/CdS-{112} interface, however, there exists an axis that is parallel to the ⟨110⟩ orientation of both CdS and CdTe. It is the direction orthogonal to this ⟨110⟩ that becomes different, being ⟨116⟩ for CdTe and ⟨111⟩ for CdS, respectively. Missing CdTe-{110} planes are found along the ⟨110⟩ axis, suggesting that the misfit strain is released by the conventional misfit dislocation mechanism along this axis. In the orthogonal axis, the misfit strain is found to be more effectively released by the new grain orientation mechanism. Our finding is supported by literature experimental observations of the change of growth direction when Cd0.96Zn0.04Te films are deposited on GaAs. Analyses of energetics clearly demonstrate the cause for the formation of the new orientation, and the insights gained from our studies can help understand the grain structures experimentally observed in lattice mismatched systems.

Nanopatterning and bandgap grading to reduce defects in CdTe solar cells

We present simulation and experimental results proving the feasibility of a novel concept to increase efficiency of CdTe based solar cells. In order to achieve

Atomistic models for scintillator discovery

0.50/W price in CdTe based modules, higher efficiencies need to be attained. The high defect density due to lattice-mismatch between CdS and CdTe reduces lifetime, voltage, and efficiency of the cells. We propose the use of a graded composition structure and a patterned substrate to reduce defects, increase lifetime, and efficiency of the cells. Innovative simulations using high-fidelity molecular dynamics predict that defect-free films are possible if the CdTe film is graded with Zn and is constructed as nano-islands with sizes below 90 nm. Both graded structure and nano-islands reduce the lattice-mismatch stresses. Also, the graded composition creates a back surface field and an enhanced ohmic contact. We have attempted to grow ZnTe and CdTe films on CdS substrates using a template of micro and nano-islands. Selective growths on patterned substrates have shown fewer grain boundaries when the island size decreases below 300 nm. Also, larger grain sizes were obtained using a CdTe/ZnTe stack when compared to a single layer CdTe. The simulation and experimental results demonstrate for the first time the ability to use nanopatterned substrates to enhance uniformity in thin film solar cells.

Predicting growth of graphene nanostructures using high-fidelity atomistic simulations

A2BLnX6 elpasolites (A, B: alkali; Ln: lanthanide; X: halogen), LaBr3 lanthanum bromide, and AX alkali halides are three classes of the ionic compound crystals being explored for γ-ray detection applications. Elpasolites are attractive because they can be optimized from combinations of four different elements. One design goal is to create cubic crystals that have isotropic optical properties and can be grown into large crystals at lower costs. Unfortunately, many elpasolites do not have cubic crystals and the experimental trial-and-error approach to find the cubic elpasolites has been prolonged and inefficient. LaBr3 is attractive due to its established good scintillation properties. The problem is that this brittle material is not only prone to fracture during services, but also difficult to grow into large crystals resulting in high production cost. Unfortunately, it is not always clear how to strengthen LaBr3 due to the lack of understanding of its fracture mechanisms. The problem with alkali halides is that their properties decay rapidly over time especially under harsh environment. Here we describe our recent progress on the development of atomistic models that may begin to enable the prediction of crystal structures and the study of fracture mechanisms of multi-element compounds.