Engineering PapersSearch

SEARCH · Engineering Papers

Results for “Atomistic simulations”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records

Random Dopant Induced Threshold Voltage Lowering and Fluctuations in Sub 50 nm MOSFETs: a Statistical 3D 'Atomistic' Simulation Study

A 3D 'atomistic' simulation study of random dopant induced threshold voltage fluctuations and lowering in sub 50 nm MOSFETs is presented. The attention is focused mainly on devices with 30 nm effective channel length which represent the expected level of scaling at the end of the Silicon Roadmap. An efficient algorithm, based on a single 3D ap solution of the Poisson equation and a simplified current continuity equation, is used in the simulations. Large samples of microscopically different devices (typically 200) arc used in order to obtain statistically reliable results. The influence of different aspects of the conventional MOSFET design on the threshold voltage fluctuations and lowering are investigated. Results for fluctuation resistant device architectures based on low-doped epitaxial channel MOSFETs are also presented.

Asenov, Asen

Random Dopant Induced Threshold Voltage Lowering and Fluctuations in Sub-0.1 (micron)meter MOSFET's: A 3-D 'Atomistic' Simulation Study

A three-dimensional (3-D) "atomistic" simulation study of random dopant induced threshold voltage lowering and fluctuations in sub-0.1 microns MOSFET's is presented. For the first time a systematic analysis of random dopant effects down to an individual dopant level was carried out in 3-D on a scale sufficient to provide quantitative statistical predictions. Efficient algorithms based on a single multigrid solution of the Poisson equation followed by the solution of a simplified current continuity equation are used in the simulations. The effects of various MOSFET design parameters, including the channel length and width, oxide thickness and channel doping, on the threshold voltage lowering and fluctuations are studied using typical samples of 200 atomistically different MOSFET's. The atomistic results for the threshold voltage fluctuations were compared with two analytical models based on dopant number fluctuations. Although the analytical models predict the general trends in the threshold voltage fluctuations, they fail to describe quantitatively the magnitude of the fluctuations. The distribution of the atomistically calculated threshold voltage and its correlation with the number of dopants in the channel of the MOSFET's was analyzed based on a sample of 2500 microscopically different devices. The detailed analysis shows that the threshold voltage fluctuations are determined not only by the fluctuation in the dopant number, but also in the dopant position.

Asenov, Asen

Efficient 3D 'Atomistic' Simulation Technique for Studying of Random Dopant Induced Threshold Voltage Lowering and Fluctuations in Decanano MOSFETs

A 3D 'atomistic' simulation technique to study random dopant induced threshold voltage lowering and fluctuations in sub 0.1 micron MOSFETs is presented. It allows statistical analysis of random impurity effects down to the individual impurity level. Efficient algorithms based on a single solution of Poisson's equation, followed by the solution of a simplified current continuity equation are used in the simulations.

Asenov, Asen

Atomistic Simulations and Experimental Analysis of the Effect of Ti Additions on the Structure of NiAl

The Bozzolo-Ferrante-Smith (BFS) semiempirical method for alloy energetics is applied to the study of ternary additions to NiAl alloys. A detailed description of the method and its application to alloy design is given. Two different approaches are used in the analysis of the effect of Ti additions to NiAl. First, a thorough analytical study is performed, where the energy of formation, lattice parameter and bulk modulus are calculated for hundreds of possible atomic distributions of Ni, Al and Ti. Substitutional site preference schemes and formation of precipitates are thus predicted and analyzed. The second approach used consists of the determination of temperature effects on the final results, as obtained by performing a number of large scale numerical simulations using the Monte Carlo - Metropolis procedure and BFS for the calculation of the energy at every step in the simulation. The results indicate a sharp preference of Ti for Al sites in Ni-rich NiAl alloys and the formation of ternary Heusler precipitates beyond the predicted solubility limit of 5 at. % Ti. Experimental analysis of three NiAl+Ti alloys confirms the theoretical predictions.

Bozzolo, Guillermo

Atomistic Simulations of Ti Additions to NiAl

The development of more efficient engines and power plants for future supersonic transports depends on the advancement of new high-temperature materials with temperature capabilities exceeding those of Ni-based superalloys. Having theoretical modelling techniques to aid in the design of these alloys would greatly facilitate this development. The present paper discusses a successful attempt to correlate theoretical predictions of alloy properties with experimental confirmation for ternary NiAl-Ti alloys. The B.F.S. (Bozzolo-Ferrante-Smith) method for alloys is used to predict the solubility limit and site preference energies for Ti additions of 1 to 25 at.% to NiAl. The results show the solubility limit to be around 5% Ti, above which the formation of Heusler precipitates is favored. These results were confirmed by transmission electron microscopy performed on a series of NiAl-Ti alloys.

Bozzolo, Guillermo

Suppression of Random Dopant-Induced Threshold Voltage Fluctuations in Sub-0.1-(micron)meter MOSFET's with Epitaxial and (delta)-Doped Channels

A detailed three-dimensional (3-D) statistical 'atomistic' simulation study of fluctuation-resistant sub-0.1-(micron)meter MOSFET architectures with epitaxial channels and delta doping is presented. The need for enhancing the fluctuation resistance of the sub-0.1-(micron)meter generation transistors is highlighted by presenting summarized results from atomistic simulations of a wide range of conventional devices with uniformly doped channel. According to our atomistic results, the doping concentration dependence of the random dopant-induced threshold voltage fluctuations in conventional devices is stronger than the analytically predicted fourth-root dependence. As a result of this, the scaling of such devices will be restricted by the "intrinsic" random dopant-induced fluctuations earlier than anticipated. Our atomistic simulations confirm that the introduction of a thin epitaxial layer in the MOSFET's channel can efficiently suppress the random dopant-induced threshold voltage fluctuations in sub-0.1-(micron)meter devices. For the first time, we observe an "anomalous" reduction in the threshold voltage fluctuations with an increase in the doping concentration behind the epitaxial channel, which we attribute to screening effects. Also, for the first time we study the effect of a delta-doping, positioned behind the epitaxial layer, on the intrinsic threshold voltage fluctuations. Above a certain thickness of epitaxial layer, we observe a pronounced anomalous decrease in the threshold voltage fluctuation with the increase of the delta doping. This phenomenon, which is also associated with screening, enhances the importance of the delta doping in the design of properly scaled fluctuation-resistant sub-0.1-(micron)meter MOSFET's. Index Terms-Doping, fluctuations, MOSFET, semiconductor device simulation, silicon devices, threshold.

Asenov, Asen

Simulation studies for surfaces and materials strength

During this reporting period three investigations were carried out. The first area of research concerned the analysis of the structure-energy relationship in small clusters. This study is very closely related to the improvement of the potential energy functions which are suitable and simple enough to be used in atomistic simulation studies. Parameters obtained from ab initio calculations for dimers and trimers of Al were used to estimate energetics and global minimum energy structures of clusters continuing up to 15 Al atoms. The second research topic addressed modeling of the collision process for atoms impinging on surfaces. In this simulation study qualitative aspects of the O atom collision with a graphite surface were analyzed. Four different O/graphite systems were considered and the aftermath of the impact was analyzed. The final area of investigation was related to the simulation of thin amorphous Si films on crystalline Si substrates. Parameters obtained in an earlier study were used to model an exposed amorphous Si surface and an a-Si/c-Si interface. Structural details for various film thicknesses were investigated at an atomistic level.

Halicioglu, T.

Electronic materials high-T(sub c) superconductivity polymers and composites structural materials surface science and catalysts industry participation

The fifth year of the Center for Advanced Materials was marked primarily by the significant scientific accomplishments of the research programs. The Electronics Materials program continued its work on the growth and characterization of gallium arsenide crystals, and the development of theories to understand the nature and distribution of defects in the crystals. The High Tc Superconductivity Program continued to make significant contributions to the field in theoretical and experimental work on both bulk materials and thin films and devices. The Ceramic Processing group developed a new technique for cladding YBCO superconductors for high current applications in work with the Electric Power Research Institute. The Polymers and Composites program published a number of important studies involving atomistic simulations of polymer surfaces with excellent correlations to experimental results. The new Enzymatic Synthesis of Materials project produced its first fluorinated polymers and successfully began engineering enzymes designed for materials synthesis. The structural Materials Program continued work on novel alloys, development of processing methods for advanced ceramics, and characterization of mechanical properties of these materials, including the newly documented characterization of cyclic fatigue crack propagation behavior in toughened ceramics. Finally, the Surface Science and Catalysis program made significant contributions to the understanding of microporous catalysts and the nature of surface structures and interface compounds.

Source record

Energetics and crystal chemistry of Ruddlesden-Popper type structures in high T(sub c) ceramic superconductors

The formation of Ruddlesden-Popper type layers (alternating slabs of rocksalt and perovskite structures) is seen in these oxides which is similar in many respects to what is seen in the system Sr-Ti-O. However, it was observed that there are some significant differences, for example the rocksalt and perovskite blocks in new superconducting compounds are not necessarily electrically, unlike in Sr-Ti-O systems. This will certainly render an additional coulombic bonding energy between two different types of blocks and may well lead to significant differences in their structural chemistry. In the higher order members of the various homologous series, additional Cu-O planes are inserted in the perovskite blocks. In order for the unit cell to electrically neutral the net positive charge on rocksalt block (which remains constant throughout the homologous series) should be balanced by an equal negative charge on perovskite block. It, thus becomes necessary to create oxygen vacancies in the basic perovskite structure, when width of the perovskite slab changes on addition of extra Cu-O planes. Results of atomistic simulations suggest that these missing oxygen ions allow the Cu-O planes to buckle in these compounds. This is also supported by the absence of buckling in the first member of Bi-containing compounds in which there are no missing oxygen ions and the Sr-Ti-O series of compounds. Additional results are presented on the phase stability of polytypoid structures in these crystal chemically complex systems. The studies will focus on the determination of the location of Cu(3+) in the structures of higher order members of the La-Cu-O system and whether Cu(3+) ions or oxygen vacancies are energetically more favorable charge compensating mechanism.

Dwivedi, Anurag

Diffusion on Cu surfaces

Understanding surface diffusion is essential in understanding surface phenomena, such as crystal growth, thin film growth, corrosion, physisorption, and chemisorption. Because of its importance, various experimental and theoretical efforts have been directed to understand this phenomena. The Field Ion Microscope (FIM) has been the major experimental tool for studying surface diffusion. FIM have been employed by various research groups to study surface diffusion of adatoms. Because of limitations of the FIM, such studies are only limited to a few surfaces: nickel, platinum, aluminum, iridium, tungsten, and rhodium. From the theoretical standpoint, various atomistic simulations are performed to study surface diffusion. In most of these calculations the Embedded Atom Method (EAM) along with the molecular static (MS) simulation are utilized. The EAM is a semi-empirical approach for modeling the interatomic interactions. The MS simulation is a technique for minimizing the total energy of a system of particles with respect to the positions of its particles. One of the objectives of this work is to develop the EAM functions for Cu and use them in conjunction with the molecular static (MS) simulation to study diffusion of a Cu atom on a perfect as well as stepped Cu(100) surfaces. This will provide a test of the validity of the EAM functions on Cu(100) surface and near the stepped environments. In particular, we construct a terrace-ledge-kink (TLK) model and calculate the migration energies of an atom on a terrace, near a ledge site, near a kink site, and going over a descending step. We have also calculated formation energies of an atom on the bare surface, a vacancy in the surface, a stepped surface, and a stepped-kink surface. Our results are compared with the available experimental and theoretical results.

Karimi, Majid

Analysis of Surface and Bulk Behavior in Ni-Pd Alloys

The most salient features of the surface structure and bulk behavior of Ni-Pd alloys have been studied using the BFS method for alloys. Large-scale atomistic simulations were performed to investigate surface segregation profiles as a function of temperature, crystal face, and composition. Pd enrichment of the first layer was observed in (111) and (100) surfaces, and enrichment of the top two layers occurred for (110) surfaces. In all cases, the segregation profile shows alternate planes enriched and depleted in Pd. In addition, the phase structure of bulk Ni-Pd alloys as a function of temperature and composition was studied. A weak ordering tendency was observed at low temperatures, which helps explain the compositional oscillations in the segregation profiles. Finally, based on atom-by-atom static energy calculations, a comprehensive explanation for the observed surface and bulk features will be presented in terms of competing chemical and strain energy effects.

Ordering

New Theoretical Technique Developed for Predicting the Stability of Alloys

When alloys are being designed for aeronautical and other applications, a substantial experimental effort is necessary to make incremental changes in the desired alloy properties. A scheme to narrow the field to the most promising candidates would substantially reduce the high cost of this experimental screening. Such a method for determining alloy properties, called the BFS (Bozzolo, Ferrante, and Smith) method, has been developed at the NASA Lewis Research Center. This method was used to calculate the thermal stability and mechanical strength of 200 alloys of Ni3Al, with Cu and Au impurities forming ternary and quaternary compounds. With recent advances in the method, almost any metallic impurity and crystal structure can be addressed. In addition, thermal effects can be addressed with Monte Carlo techniques. At present, an experimental program is in progress to verify these results. The method identified a small number of the most promising candidates from the 200 alloys with the largest negative heat of formation and the highest bulk modulus. This calculation required only 5 min of CPU time on a VAX computer. It is clear that semi-empirical methods have achieved the level of development and reliability to warrant examining this new approach to the problem of alloy design. The present work was meant to demonstrate, perhaps in a rather simple way, this power. This type of application of atomistic simulation methods can narrow the gap and improve the feedback between theoretical predictions and laboratory experimentation.

Source record

Computer simulation of surface and film processes

All the investigations which were performed employed in one way or another a computer simulation technique based on atomistic level considerations. In general, three types of simulation methods were used for modeling systems with discrete particles that interact via well defined potential functions: molecular dynamics (a general method for solving the classical equations of motion of a model system); Monte Carlo (the use of Markov chain ensemble averaging technique to model equilibrium properties of a system); and molecular statics (provides properties of a system at T = 0 K). The effects of three-body forces on the vibrational frequencies of triatomic cluster were investigated. The multilayer relaxation phenomena for low index planes of an fcc crystal was analyzed also as a function of the three-body interactions. Various surface properties for Si and SiC system were calculated. Results obtained from static simulation calculations for slip formation were presented. The more elaborate molecular dynamics calculations on the propagation of cracks in two-dimensional systems were outlined.

Tiller, W. A.

Computational Modeling Method for Superalloys

Computer modeling based on theoretical quantum techniques has been largely inefficient due to limitations on the methods or the computer needs associated with such calculations, thus perpetuating the notion that little help can be expected from computer simulations for the atomistic design of new materials. In a major effort to overcome these limitations and to provide a tool for efficiently assisting in the development of new alloys, we developed the BFS method for alloys, which together with the experimental results from previous and current research that validate its use for large-scale simulations, provide the ideal grounds for developing a computationally economical and physically sound procedure for supplementing the experimental work at great cost and time savings.

Bozzolo, Guillermo

An atomistic model of slip formation

The results of an atomistic model for the simulation of the early stages of crack initiation in a two-dimensional triangular lattice are presented. In the current model, each particle in the system is treated discretely and assumed to be interacting with the surrounding particles via Lennard-Jones potentials. A uniaxial load (in incremental elongations) is applied to the rectangular two-dimensional slab in either the x or the y direction. After each incremental elongation the system is equilibrated using a static method. Initially, elastic behavior in the x and y directions is observed. Continued elongation results in plastic deformation. In lattices with point defects, the defects first move to the surface, creating vacancies which trigger plastic deformation.

Halicioglu, T.

Energetics of Single Substitutional Impurities in NiTi

Shape-memory alloys are of considerable current interest, with applications ranging from stents to Mars rover components. In this work, we present results on the energetics of single substitutional impurities in B2 NiTi. Specifically, energies of Pd, Pt, Zr and Hf impurities at both Ni and Ti sites are computed. All energies are computed using the CASTEP ab initio code, and, for comparison, using the quantum approximate energy method of Bozzolo, Ferrante and Smith. Atomistic relaxation in the vicinity of the impurities is investigated via quantum approximate Monte Carlo simulation, and in cases where the relaxation is found to be important, the resulting relaxations are applied to the ab initio calculations. We compare our results with available experimental work.

Binary Alloys