Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userKinetic Monte Carlo simulation of shape transition in strained quantum dots
140
0
0.0
(
0
)
Authors
Chi-Hang Lam
Ask ChatGPT about the research
No Arabic abstract
The pyramid-to-dome transition in Ge$_{x}$Si$_{1-x}$ on Si(100) initiated by step bunching on pyramidal quantum dots is atomistically simulated using a novel multi-state lattice model incorporating effective surface reconstructions. Results are explained by a simple theory based on a shallow island approximation. Under given deposition conditions in $d$ dimensions, the shape transition is shown to occur at island size $n_c$ following $n_c^{1/d} propto x^{-zeta}$ independent of temperature and deposition rate, where $zetaalt 2$ and $x$ is the actual Ge concentration in the island. The transition has an energy barrier dominated by the facet interface energy. Fast deposition however can out-run and delay the transition to larger island sizes.
rate research
Read More
Accelerated algorithms for simulating the morphological evolution of strained heteroeptiaxy based on a ball and spring lattice model in three dimensions are explained. We derive exact Greens function formalisms for boundary values in the associated lattice elasticity problems. The computational efficiency is further enhanced by using a superparticle surface coarsening approximation. Atomic hoppings simulating surface diffusion are sampled using a multi-step acceptance-rejection algorithm. It utilizes quick estimates of the atomic elastic energies from extensively tabulated values modulated by the local strain. A parameter controls the compromise between accuracy and efficiency of the acceptance-rejection algorithm.
We present a new efficient method for Monte Carlo simulations of diffusion-reaction processes. First introduced by us in [Phys. Rev. Lett., 97:230602, 2006], the new algorithm skips the traditional small diffusion hops and propagates the diffusing particles over long distances through a sequence of super-hops, one particle at a time. By partitioning the simulation space into non-overlapping protecting domains each containing only one or two particles, the algorithm factorizes the N-body problem of collisions among multiple Brownian particles into a set of much simpler single-body and two-body problems. Efficient propagation of particles inside their protective domains is enabled through the use of time-dependent Greens functions (propagators) obtained as solutions for the first-passage statistics of random walks. The resulting Monte Carlo algorithm is event-driven and asynchronous; each Brownian particle propagates inside its own protective domain and on its own time clock. The algorithm reproduces the statistics of the underlying Monte-Carlo model exactly. Extensive numerical examples demonstrate that for an important class of diffusion-reaction models the new algorithm is efficient at low particle densities, where other existing algorithms slow down severely.
By carrying out Monte Carlo simulations based on the two-species atomic-scale kinetic growth model of GaAs(001) homoepitaxy and comparing the results with scanning tunneling microscope images, we show that initial growing islands undergo the structural transformation before adopting the proper beta2(2x4) reconstruction.
We describe the development of a new object kinetic Monte Carlo code where the elementary defect objects are off-lattice atomistic configurations. Atomic-level transitions are used to transform and translate objects, to split objects and to merge them together. This gradually constructs a database of atomic configurations -- a set of relevant defect objects and their possible events generated on-the-fly. Elastic interactions are handled within objects with empirical potentials at short distances, and between spatially distinct objects using the dipole tensor formalism. The model is shown to evolve mobile interstitial clusters in tungsten faster than an equivalent molecular dynamics simulation, even at elevated temperatures. We apply the model to the evolution of complex defects generated using molecular dynamics simulations of primary radiation damage in tungsten. We show that we can evolve defect structures formed in cascade simulations to experimentally observable timescales of seconds while retaining atomistic detail. We conclude that the first few nanoseconds of simulation following cascade initiation would be better performed using molecular dynamics, as this will capture some of the near-temperature-independent evolution of small highly-mobile interstitial clusters. We also conclude that, for the 20keV PKA cascades annealing simulations considered here, internal relaxations of sessile objects difficult to capture using conventional object KMC with idealised object geometries establish the conditions for long timescale evolution.
High entropy alloys (HEAs) are a series of novel materials that demonstrate many exceptional mechanical properties. To understand the origin of these attractive properties, it is important to investigate the thermodynamics and elucidate the evolution of various chemical phases. In this work, we introduce a data-driven approach to construct the effective Hamiltonian and study the thermodynamics of HEAs through canonical Monte Carlo simulation. The main characteristic of our method is to use pairwise interactions between atoms as features and systematically improve the representativeness of the dataset using samples from Monte Carlo simulation. We find this method produces highly robust and accurate effective Hamiltonians that give less than 0.1 mRy test error for all the three refractory HEAs: MoNbTaW, MoNbTaVW, and MoNbTaTiW. Using replica exchange to speed up the MC simulation, we calculated the specific heats and short-range order parameters in a wide range of temperatures. For all the studied materials, we find there are two major order-disorder transitions occurring respectively at $T_1$ and $T_2$, where $T_1$ is near room temperature but $T_2$ is much higher. We further demonstrate that the transition at $T_1$ is caused by W and Nb while the one at $T_2$ is caused by the other elements. By comparing with experiments, {color{black} the results provide insight into the role of chemical ordering in the strength and ductility of HEAs.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
