Researcher profile

Dallas R. Trinkle

Dallas R. Trinkle contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
14works
0followers
3topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

14 published item(s)

preprint2025arXiv

Investigating the effects of local environment on nitrogen vacancies in high entropy metal nitrides

High entropy metal nitrides are an important material class in a variety of applications, and the role of nitrogen vacancies is of great importance for understanding their stability and mechanical properties. We study six different high entropy nitrides with eight different metal species to build a predictive model of the nitrogen vacancy formation energy. We construct sets of supercells that maximize the number of unique nitrogen environments for a given chemistry, and then use density-functional theory to calculate the energy density for all nitrogen sites, and the vacancy formation energies for the highest, lowest, and a median subset based on the energy densities. The energy density of nitrogen sites correlates with the vacancy formation energies, for binary, ternary and high entropy nitrides. A linear regression model predicts the vacancy formation energies using only the nearest-neighbor composition; across our eight metals, we find the largest vacancy formation energies next to Hf, then Zr, Ti, V, Cr, Ta, Nb, and the lowest near Mo. Additionally, we see that binary nitride data shows qualitatively similar vacancy formation energy trends for high entropy nitrides; however, the binary data alone is insufficient to predict the complex nitride behavior. Our model is both predictive and easily interpretable, and correlates with experimental data.

preprint2016arXiv

Diffusivity and derivatives for interstitial solutes: Activation energy, volume, and elastodiffusion tensors

Computational atomic-scale methods continue to provide new information about geometry, energetics, and transition states for interstitial elements in crystalline lattices. This data can be used to determine the diffusivity of interstitials by finding steady-state solutions to the master equation. In addition, atomic-scale computations can provide not just the site energy, but also the stress in the cell due to the introduction of the defect to compute the elastic dipole. We derive a general expression for the fully anistropic diffusivity tensor from site and transition state energies, and three derivatives of the diffusivity: the elastodiffusion tensor (derivative of diffusivity with respect to strain), the activation barrier tensor (logarithmic derivative of diffusivity with respect to inverse temperature) and activation volume tensor (logarithmic derivative of diffusivity with respect to pressure). Computation of these quantities takes advantage of crystalline symmetry, and we provide an open-source implementation of the algorithm. We provide analytic results for octahedral-tetrahedral networks in face-centered cubic, body-centered cubic, and hexagonal closed-packed lattices, and conclude with numerical results for C in Fe.

preprint2016arXiv

Light element diffusion in Mg using first principles calculations: Anisotropy and elastodiffusion

The light elemental solutes B, C, N, and O can penetrate the surface of Mg alloys and diffuse during heat treatment or high temperature application, forming undesirable compounds. We investigate the diffusion of these solutes by determining their stable interstitial sites and the inter-penetrating network formed by these sites. We use density functional theory (DFT) to calculate the site energies, migration barriers, and attempt frequencies for these networks to inform our analytical model for bulk diffusion. Due to the nature of the networks, O diffuses isotropically, while B, C, and N diffuse anisotropically. We compute the elastodiffusion tensor which quantifies changes in diffusivity due to small strains that perturb the diffusion network geometry and the migration barriers. The DFT-computed elastic dipole tensor which quantifies the change in site energies and migration barriers due to small strains is used as an input to determine the elastodiffusion tensor. We employ the elastodiffusion tensor to determine the effect of thermal strains on interstitial diffusion and find that B, C, and N diffusivity increases on crystal expansion, while O diffusivity decreases. From the elastodiffusion and compliance tensors we calculate the activation volume of diffusion and find that it is positive and anisotropic for B, C and N diffusion, whereas it is negative and isotropic for O diffusion.

preprint2012arXiv

Direct calculation of lattice Green function with arbitrary interactions for general crystals

Efficient computation of lattice defect geometries such as point defects, dislocations, disconnections, grain boundaries, interfaces and free surfaces requires accurate coupling of displacements near the defect to the long-range elastic strain. Flexible boundary condition methods embedded a defect in infinite harmonic bulk through the lattice Green function. We demonstrate an efficient and accurate calculation of the lattice Green function from the force-constant matrix for general crystals with an arbitrary basis by extending a method for Bravais lattices. New terms appear due to the presence of optical modes and the possible loss of inversion symmetry. By separately treating poles and discontinuities in reciprocal space, numerical accuracy is controlled at all distances. We compute the lattice Green function for a two-dimensional model with broken symmetry to elucidate the role of different coupling terms. The algorithm is generally applicable in two and three dimensions, to crystals with arbitrary number of atoms in the unit cell, symmetry, and interactions.

preprint2011arXiv

Au/TiO2(110) interfacial reconstruction stability from ab initio

We determine the stability and properties of interfaces of low-index Au surfaces adhered to TiO2(110), using density functional theory energy density calculations. We consider Au(100) and Au(111) epitaxies on rutile TiO2(110) surface, as observed in experiments. For each epitaxy, we consider several different interfaces: Au(111)//TiO2(110) and Au(100)//TiO2(110), with and without bridging oxygen, Au(111) on 1x2 added-row TiO2(110) reconstruction, and Au(111) on a proposed 1x2 TiO reconstruction. The density functional theory energy density method computes the energy changes on each of the atoms while forming the interface, and evaluates the work of adhesion to determine the equilibrium interfacial structure.

preprint2011arXiv

Direct diffusion through interpenetrating networks: Oxygen in titanium

How impurity atoms move through a crystal is a fundamental and recurrent question in materials. The previous understanding of oxygen diffusion in titanium relied on interstitial lattice sites that were recently found to be unstable, making the diffusion pathways for oxygen unknown. Using first-principles quantum-mechanical methods, we find three oxygen interstitial sites in titanium, and quantify the multiple interpenetrating networks for oxygen diffusion. Surprisingly, no single transition dominates, but all contribute to diffusion.

preprint2011arXiv

Nanoscale-hydride formation at dislocations in palladium: Ab initio theory and incoherent inelastic neutron scattering measurements

Hydrogen arranges at dislocations in palladium to form nanoscale hydrides, changing the vibrational spectra. An ab initio hydrogen potential energy model versus Pd neighbor distances allows us to predict the vibrational excitations for H from absolute zero up to room temperature adjacent to a partial dislocation and with strain. Using the equilibrium distribution of hydrogen with temperature, we predict excitation spectra to explain new incoherent inelastic neutron-scattering measurements. At 0K, dislocation cores trap H to form nanometer-sized hydrides, while increased temperature dissolves the hydrides and disperses H throughout bulk Pd.

preprint2011arXiv

Prediction of thermal cross-slip stress in magnesium alloys from direct first principles data

We develop a first-principles model of thermally-activated cross-slip in magnesium in the presence of a random solute distribution. Electronic structure methods provide data for the interaction of solutes with prismatic dislocation cores and basal dislocation cores. Direct calculations of interaction energies are possible for solutes---K, Na, and Sc---that lower the Mg prismatic stacking fault energy to improve formability. To connect to thermally activated cross-slip, we build a statistical model for the distribution of activation energies for double kink nucleation, barriers for kink migration, and roughness of the energy landscape to be overcome by an athermal stress. These distributions are calculated numerically for a range of concentrations, as well as alternate approximate analytic expressions for the dilute limit. The analytic distributions provide a simplified model for the maximum cross-slip softening for a solute as a function of temperature. The direct interaction calculations predict lowered forming temperatures for Mg-0.7at.%Sc, Mg-0.4at.%K, and Mg-0.6at.%Na of approximately 250C.

preprint2010arXiv

Energy density in density functional theory: Application to crystalline defects and surfaces

We propose a method to decompose the total energy of a supercell containing defects into contributions of individual atoms, using the energy density formalism within density functional theory. The spatial energy density is unique up to a gauge transformation, and we show that unique atomic energies can be calculated by integrating over Bader and charge-neutral volumes for each atom. Numerically, we implement the energy density method in the framework of the Vienna ab initio simulation package (VASP) for both norm-conserving and ultrasoft pseudopotentials and the projector augmented wave method, and use a weighted integration algorithm to integrate the volumes. The surface energies and point defect energies can be calculated by integrating the energy density over the surface region and the defect region, respectively. We compute energies for several surfaces and defects: the (110) surface energy of GaAs, the mono-vacancy formation energies of Si, the (100) surface energy of Au, and the interstitial formation energy of O in the hexagonal close-packed Ti crystal. The surface and defect energies calculated using our method agree with size-converged calculations of the difference between the total energies of the system with and without the defect. Moreover, the convergence of the defect energies with size can be found from a single calculation.

preprint2010arXiv

First-principles calculation of H vibrational excitations at a dislocation core of Pd

Palladium is an ideal system for understanding the behavior of hydrogen in metals. In Pd, H is located both in octahedral sites and in dislocation cores, which act as nanoscale H traps and form Cottrell atmospheres. Adjacent to a dislocation core, H experiences the largest possible distortion in alpha-Pd. Ab initio density-functional theory computes the potential energy for a hydrogen in an octahedral site in alpha-Pd and in a trap site at the core of a partial of an edge dislocation. The Pd partial dislocation core changes the environment for H, distorting the H-Pd bonding which changes the local potential, vibrational spectra, and inelastic form factor for an isolated H atom. The decrease in excitation energy is consistent with experiments, and the calculations predict distortions to the H wavefunction.

preprint2010arXiv

First-principles data for solid-solution strengthening of magnesium: From geometry and chemistry to properties

Solid-solution strengthening results from solutes impeding the glide of dislocations. Existing theories of strength rely on solute-dislocation interactions, but do not consider dislocation core structures, which need an accurate treatment of chemical bonding. Here, we focus on strengthening of Mg, the lightest of all structural metals and a promising replacement for heavier steel and aluminum alloys. Elasticity theory, which is commonly used to predict the requisite solute-dislocation interaction energetics, is replaced with quantum-mechanical first-principles calculations to construct a predictive mesoscale model for solute strengthening of Mg. Results for 29 different solutes are displayed in a "strengthening design map" as a function of solute misfits that quantify volumetric strain and slip effects. Our strengthening model is validated with available experimental data for several solutes, including Al and Zn, the two most common solutes in Mg. These new results highlight the ability of quantum-mechanical first-principles calculations to predict complex material properties such as strength.

preprint2010arXiv

Misfit-Dislocation-Mediated Heteroepitaxial Island Diffusion

Scanning tunneling microscopy combined with molecular dynamics simulations reveal a dislocation-mediated island diffusion mechanism for Cu on Ag(111), a highly mismatched system. Cluster motion is tracked with atomic precision at multiple temperatures and diffusion barriers and prefactors are determined from direct measurements of hop rates. The non-monotonic size dependence of the diffusion barrier is in good agreement with simulations and can lead to enhanced mass transport upon coarsening, in surprising contrast to the traditional island diffusion models where diffusivity reduces with cluster size.

preprint2009arXiv

Cu/Ag EAM Potential Optimized for Heteroepitaxial Diffusion from ab initio Data

A binary embedded-atom method (EAM) potential is optimized for Cu on Ag(111) by fitting to ab initio data. The fitting database consists of DFT calculations of Cu monomers and dimers on Ag(111), specifically their relative energies, adatom heights, and dimer separations. We start from the Mishin Cu-Ag EAM potential and first modify the Cu-Ag pair potential to match the FCC/HCP site energy difference then include Cu-Cu pair potential optimization for the entire database. The optimized EAM potential reproduce DFT monomer and dimer relative energies and geometries correctly. In trimer calculations, the potential produces the DFT relative energy between FCC and HCP trimers, though a different ground state is predicted. We use the optimized potential to calculate diffusion barriers for Cu monomers, dimers, and trimers. The predicted monomer barrier is the same as DFT, while experimental barriers for monomers and dimers are both lower than predicted here. We attribute the difference with experiment to the overestimation of surface adsorption energies by DFT and a simple correction is presented. Our results show that the optimized Cu-Ag EAM can be applied in the study of larger Cu islands on Ag(111).

preprint2009arXiv

Island Shape Controls Magic-Size Effect for Heteroepitaxial Diffusion

Lattice mismatch of Cu on Ag(111) produces fast diffusion for special "magic sizes" of islands. A size- and shape-dependent reptation mechanism is responsible for low diffusion barriers. Initiating the reptation mechanism requires a suitable island shape, a property not considered in previous studies of 1D island chains and 2D closed-shell islands. Shape determines the dominant diffusion mechanism and leads to multiple clearly identifiable magic-size trends for diffusion depending on the number of atoms whose bonds are shortened during diffusion.