Source author record

Peter Smereka

Peter Smereka 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-sci math-ph math.MP

Catalog footprint

What is connected

5works

3topics

4close collaborators

Actions

Connect this record

Open graph Browse works

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.

preprint2016arXiv

Crystallization for a Brenner-like potential

Graphene is a carbon molecule with the structure of a honeycomb lattice. We show how this structure can arise in two dimensions as the minimizer of an interaction energy with two-body and three-body terms. In the engineering literature, the Brenner potential is commonly used to describe the interactions between carbon atoms. We consider a potential of Stillinger-Weber type that incorporates certain characteristics of the Brenner potential: the preferred bond angles are 180 degrees and all interactions have finite range. We show that the thermodynamic limit of the ground state energy per particle is the same as that of a honeycomb lattice. We also prove that, subject to periodic boundary conditions, the minimizers are translated versions of the honeycomb lattice.

preprint2015arXiv

Modeling the Elastic Energy of Alloys: Potential Pitfalls of Continuum Treatments

Some issues that arise when modeling elastic energy for binary alloys are discussed within the context of a Keating model and density functional calculations. The Keating model is based on atomistic modeling of elastic interactions in binary alloy using harmonic springs with species dependent equilibrium lengths. It is demonstrated that the continuum limit for the strain field are the usual equations of linear elasticity for alloys and that they correctly capture the coarse-grained displacement field. In addition, it is established that Euler-Lagrange equation of the continuum limit of the elastic energy will yield the same strain field equation. However, a direct calculation of the elastic energy of the atomistic model reveals that the continuum expression for the elastic energy is both qualitatively and quantitatively incorrect. This is because it does not take atomistic scale compositional non-uniformity into account. Importantly, we also shows that finely mixed alloys tend to have more elastic energy than segregated systems, which is the opposite of predictions by some continuum theories. It is also shown that for strained thin films the traditionally used effective misfit for alloys systematically underestimate the strain energy. In some models, this drawback is handled by including an elastic contribution to the enthalpy of mixing which is characterized in terms of the continuum concentration. The direct calculation of the atomistic model reveals that this approach suffers serious difficulties. It is demonstrated that elastic contribution to the enthalpy of mixing is non-isotropic and scale dependent. It also shown that such effects are present in density-functional theory calculations for the Si/Ge and Ag/Pt systems. This work demonstrates that it is critical to include the microscopic arrangements in any elastic model to achieve even qualitatively correct behavior.

preprint2013arXiv

A unified model of droplet epitaxy for compound semiconductor nanostructures: experiments and theory

We present a unified model of compound semiconductor growth based on kinetic Monte Carlo simulations in tandem with new experimental results that can describe and predict the mechanisms for the formation of various types of nanostructures observed during droplet epitaxy. The crucial features of the model include the explicit and independent representation of atoms with different species and the ability to treat solid and liquid phases independently. Using this model, we examine nanostructural evolution in droplet epitaxy. The model faithfully captures several of the experimentally observed structures, including compact islands and nanorings. Moreover, simulations show the presence of Ga/GaAs core-shell structures that we validate experimentally. A fully analytical model of droplet epitaxy that explains the relationship between growth conditions and the resulting nanostructures is presented, yielding key insight into the mechanisms of droplet epitaxy.

preprint2012arXiv

Mechanisms of Stranski-Krastanov Growth

Stranski-Krastanov (SK) growth is reported experimentally as the growth mode that is responsible for the transition to three dimensional islands in heteroepitaxial growth. A kinetic Monte Carlo (KMC) model is proposed that can replicate many of the experimentally observed features of this growth mode. Simulations reveal that this model effectively captures the SK transition and subsequent growth. Annealing simulations demonstrate that the wetting layer formed during SK growth is stable, with entropy playing a key role in its stability. It is shown that this model also captures the apparent critical thickness that tends to occur at higher deposition rates and for alloy films (where intermixing is significant). This work shows that the wetting layer thickness increases with increasing temperature, whereas the apparent critical thickness decreases with increasing temperature. Both of which are in agreement with experiments.

preprint2009arXiv

Kinetic Monte Carlo Simulation of Strained Heteroepitaxial Growth with Intermixing

An efficient method for the simulation of strained heteroepitaxial growth with intermixing using kinetic Monte Carlo is presented. The model used is based on a solid-on-solid bond counting formulation in which elastic effects are incorporated using a ball and spring model. While idealized, this model nevertheless captures many aspects of heteroepitaxial growth, including nucleation, surface diffusion, and long range effects due elastic interaction. The algorithm combines a fast evaluation of the elastic displacement field with an efficient implementation of a rejection-reduced kinetic Monte Carlo based on using upper bounds for the rates. The former is achieved by using a multigrid method for global updates of the displacement field and an expanding box method for local updates. The simulations show the importance of intermixing on the growth of a strained film. Further the method is used to simulate the growth of self-assembled stacked quantum dots.