This research project creates a multidisciplinary intellectual partnership with the goal of designing transformative computation-based methods for the simulation of quantum dots (QDs) growth by molecular beam epitaxy. This research enables a better understanding of the kinetics and thermodynamics during the growth of quantum dots allowing for an improved control of the positioning, growth and size distribution in QD structures. The project combines expertise in materials science with state-of-the-art algorithm design for partial differential equations and moving boundary problems. This research impacts scientific and industrial communities that make use of QDs as single photon emitters or make use of the carrier spins to develop and manipulate q-bits. Outcomes help guide experiments, and the design and manufacture of new optical and electronic devices, such as QDs photonic crystal lasers, nonvolatile storage, laser scalpel and optical coherence tomography for use in medicine, and QD-based product enhancements in the energy sector. The impact on computational science is important for the simulation and design of diffusion dominated processes, moving boundary problems, and multiscale modeling where macroscopic behaviors are simulated while incorporating microscopic rates.
The research program involves the training of two students, including one student with Latino background. The research program also fosters the development of an interactive web-based computational lab for students to experiment with computational tools. In addition, advances made through this research program are disseminated through the Materials Research Laboratory at University of California Santa Barbara (UCSB), the Institute for Pure and Applied Mathematics at University of California Los Angeles (UCLA) and the California Nano Science Institute at both UCLA and UCSB.
This award is part of the Cyber-Enabled Discovery and Innovation program, and the recipients are Frederic Gibou of UCSB and Christian Ratsch of UCLA.
Normal 0 false false false false EN-US X-NONE X-NONE The goal of this research was to develop and implement a computational framework capable of simulating important processes in epitaxial growth. Epitaxial growth is a process where atoms of one material are deposited on top of another. This process is used to fabricate most modern opto-electronic devices including memory storage, lasers and optical coherence tomography, as well as to manufacture catalysts used in the energy sector, food processing, biology and environmental science. During epitaxial growth, atoms that have been deposited on a surface can diffuse on the surface, form so-called islands (or nucei) when they meet another atom, and attach to existing island.. Things get a little more complicated when atoms are deposited on top of existing islands. One of the physical phenomena that is responsible for the growth of important structures such as quantum dots and quantum posts is the so-called additional step-edge or Ehrlich-Schwoebel barrier. This barrier makes it more difficult for atoms that are on top of an island to diffuse to a lower layer (down a step), compared to diffusion on the layer. Modeling this additional barrier has been a key focus of the work funded by this NSF grant. A successful approach to the simulation of Epitaxial growth is the so-called island dynamics model [5]. Previous works have demonstrated the capabilities of this mode. However, the study of mound formation, quantum dots and in general the effect of the Ehrlich-Schwoebel barrier have not been possible because of the difficulty of imposing Robin boundary conditions in a level-set framework. As part of this NSF, we have solved that problem using a hybrid level-set/finite-volume formulation. This approach produces second-order accurate solutions and a symmetric positive definite linear system, which translates into efficient solvers. We have also extended this framework to include situations where Robin, Dirichlet and Neumann boundary conditions are imposed different parts of an irregular boundary [3]. These boundary conditions are the most general boundary conditions that can emulate virtually any epitaxial processes. We have also developed such methods in an adaptive framework designed to accurately and efficiently simulate multiscale problems. Another benefit of our approach is that the solution process directly approximates the physically correct sharp interface model. Our results demonstrate that some of the main computational components to simulate the step-edge barrier are successful in an adaptive quadtree multiscale framework [1, 4]. Our new computational approach and model has led to some very interesting new results [2]. In particular, we find that (as expected), the inclusion of a step-edge barrier leads to the formation of mounds. Figure 1 illustrates the fact that after the deposition of 20 layers, growth proceeds essentially layer-by-layer when there is no step-edge barrier present, and that mounds form when there is a step-edge barrier (middle and right panel). Moreover, the mounds become more pronounced as the step-edge barrier decreases. This new formalism has enabled us to study epitaxial growth and the formation of mounds. D'/D=.99 D'/D=.1 D'/D=.001 Figure 1: Mound formation for different values of D’/D From the point of view of computational science, communities focusing in diffusion-dominated processes and moving boundary problems will feel the impact of this research. This research will permit in a better understanding of the kinetics and thermodynamics during the growth of quantum dots and will provide a better understanding for an improved control of the positioning, growth and size distribution in Quantum Dots structures, which are used as single photon emitters or make use of the carrier spins to develop and manipulate q-bits. This grant has funded a postdoctoral scholar with Latino background (at UCLA) as well as a female graduate student (at UCSB) who obtained her PhD. This grant has also fostered the development of an interactive web-based computational lab. REFERENCES: [1] Joseph Papac, Asdis Helgadottir, Christian Ratsch, and Frederic Gibou, A level set approach for diffusion and Stefan-type problems with Robin boundary conditions on quadtree/octree adaptive Cartesian grids, Journal of Computational Physics, 233, 241-261, 2012 [2] Joseph Papac, Frederic Gibou and Christian Ratsch, Mound formations in the island dynamics model, In preparation. [3] Helgadottir, Yen-Ting Ng, Chohong Min, Christian Ratsch and Frederic Gibou, Imposing Mixed Dirichlet-Neumann-Robin Boundary Conditions in a Level-Set Framework , In preparation. [4] Frederic Gibou, Chohong Min and Ron Fedkiw, High Resolution Sharp Computational Methods for Elliptic and Parabolic Problems in Complex Geometries, Journal of Scientific Computing, 54, 369-413, 2013. [5] C. Ratsch, M.F. Gyure, R.E. Caflisch, F. Gibou, M. Petersen, M. Kang, J. Garcia, and D.D. Vvedensky, Level-Set Method for Island Dynamics in Epitaxial Growth, Phys. Rev. B 65, 195403 (2002).