Search NASA⌕ Search

Engineering topics

Samanta, Amit

Publications and source records attributed to Samanta, Amit.

33 records · Page 2

Analysis of defects in In 2 O 3 :H synthesized in presence of water vapor and hydrogen gas mixture

Using hybrid functional-based density functional theory calculations, we analyze the structure and kinetics of defects formed in two competing synthesis routes to prepare hydrogen-doped In 2 O 3 films, using a hydrogen and oxygen gas mixture vs synthesis in the presence of water vapor. For both of these synthesis routes, we find that H + is the dominant defect species: when the Fermi level is close to the conduction band, H + has a lower formation energy than other intrinsic or extrinsic defects. Our results also suggest that water molecules spontaneously split into H + (which occupies octahedral voids) and OH – interstitials (which occupies vacant oxygen lattice sites or oxygen vacancies). From the analysis of the binding energies between these different defects, we conclude that these defects do not cluster and are most likely to stay spatially distributed throughout the films. In addition, the sum of formation energies of an oxygen (i.e., O$^{2–}_i$) and a H + interstitial is close to the formation energy of a OH – interstitial, meaning that water molecules are completely split into 2H + and O$^{2–}_i$ at the synthesis conditions. Further, in the presence of H 2 + O 2 gas mixture, oxygen interstitials occupy oxygen vacancies while hydrogen interstitials occupy vacant oxygen lattice sites and form bonds with lattice oxygens. Our analysis of the defect equilibria suggests that the hydrogen content in films synthesized in the presence of water vapor is higher than films synthesized in the presence of a hydrogen gas mixture. At high dopant concentrations, a hydrogen bond network is formed in the system and this leads to large distortions in the lattice.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Intercalation of Lithium into Graphite: Insights from First-Principles Simulations

Understanding ion intercalation at electrode–electrolyte interfaces is key to the development of energy storage and water desalination. Here, we investigate Li + kinetics at a prototypical interface between graphite anodes and an organic electrolyte, and we elucidate key factors that determine ion transport, using first-principles methodology coupling ab initio molecular dynamics simulations with a solvation model. We show that surface chemical composition significantly influences the kinetics of ion intercalation from the liquid into graphite. We find that this is partly related to the ion desolvation process, which varies notably for different graphite surface chemical terminations. In addition, interfacial polarization is found to play an important role in determining energy barriers for ion transfer. We also discuss the impact of electrode potentials, which is often neglected in conventional first-principles calculations despite being a key factor in device configurations. Our study provides insights into the coupling of electronic and ionic effects of interfacial chemistry on ion transport at complex electrode–electrolyte interfaces.

25 ENERGY STORAGE↗

The role of water vapor during the synthesis of hydrogen doped In 2 O 3

Hydrogen-doped In 2 O 3 synthesized using magnetron sputtering or atomic layer deposition in the presence of water vapor has high transparency in the near infra-red region and mobility values that are almost three to four times of those in commercially used In 2 O 3 :Sn at a much lower carrier concentration. However, simple questions like what happens when water molecules enter In 2 O 3 or why the dominant charge carrier changes from an oxygen vacancy in In 2 O 3 :Sn to H + in In 2 O 3 are not clear. Using hybrid functional based density functional theory calculations, in this study we show that water molecules spontaneously split into H + and OH – . The H + interstitial bonds with an oxygen in In 2 O 3 , and the OH – interstitial occupies a vacant lattice site or is trapped by an oxygen vacancy. Thus, in agreement with experimental findings, our results suggest that the concentration of oxygen vacancies decreases as the partial pressure of water increases. The charge transition level, ϵ ( – / + ) , of a hydrogen interstitial and the Fermi level at which the formation energies of OH – and H + are equal lie above the conduction band minimum, suggesting that the H + interstitial has a lower formation energy. In addition, water molecules do not exhibit any propensity to form clusters, meaning that H + ions can diffuse and remain spatially separated from each other.

36 MATERIALS SCIENCE↗

Quantized Grain Boundary States Promote Nanoparticle Alignment During Imperfect Oriented Attachment

Oriented attachment (OA) has become a well-recognized mechanism for the growth of metal, ceramic, and biomineral crystals. While many computational and experimental studies of OA have shown that particles can attach with some misorientation then rotate to remove adjoining grain boundaries, the underlying atomistic pathways for this “Imperfect OA” process remain the subject of debate. In this paper, molecular dynamics and in situ TEM were used to probe the crystallographic evolution of up to 30 gold and copper nanoparticles during aggregation. It was found that Imperfect OA occurs because (1) grain boundaries become quantized when their size is comparable to the separation between constituent dislocations and (2) kinetic barriers associated with the glide of grain boundary dislocations are small. In support of these findings, TEM experiments show the formation of a single crystal aggregate after annealing 9 initially misoriented, agglomerated particles with evidence of dislocation slip and twin formation during particle/grain alignment. These observations motivate future work on assembled nanocrystals with tailored defects and call for a revision of Read-Shockley models for grain boundary energies in nanocrystalline materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A novel liquid-mediated nucleation mechanism for explosive crystallization in amorphous germanium

We report a novel mechanism for explosive crystallization in amorphous germanium (a-Ge), which operates through liquid-mediated nucleation occurring under extreme thermal gradient conditions. The crystallization kinetics of sputter-deposited films with thicknesses ranging from 30 to 150 nm were characterized using in situ movie-mode dynamic transmission electron microscopy (MM-DTEM). After localized heating from a short laser pulse, explosive liquid phase nucleation (LPN) was observed to occur during the early stage (<2 μs) of crystallization in the thicker (>50 nm) films deposited on silicon nitride substrates. The crystallization front propagated at ~12–15 m/s and produced nanocrystalline microstructure with ~50 nm grains. A mechanism involving the existence of a relatively thick (>100 nm) transient liquid layer and a high nucleation rate is proposed to explain the behavior. In conclusion, the key thermodynamic and kinetic features as well as the feasibility of the mechanism are further explored by employing parametric and systematic phase-field modeling and simulations.

36 MATERIALS SCIENCE↗

Modeling Mesoscale Processes of Scalable Synthesis (Final Report)

The main objective of the center is to develop tools that can help for the understanding of equilibrium and non-equilibrium properties of nanoparticles in aqueous electrolytes. To achieve this goal, we need to develop new computational methodologies that can integrate molecular level description such as the electronic structure analysis to mesoscale description such as generalized hydrodynamics. The challenge is to develop systematic coarse graining techniques that can adaptively choose the level of description needed for the particular task. For example, if bonding breaking is involved, then one should use appropriate models from quantum mechanics to accurately model the time scale and molecular mechanism for the breaking process. However, to model transport process or the process of self-assembly, much less detailed models can be used. There are two main components to this task. One is to develop accurate models that bridge different levels of description. The classical example is the QM/MM model that bridges quantum mechanics and molecular mechanics for studying chemical reactions of large molecules. This work was awarded the Nobel prize this year. For our purpose, we also need to develop methods that bridge molecular and continuum levels of description. Another main task is to develop systematic coarse-graining methods that can give rise to accurate mesoscale and macroscale models, starting from more detailed models at the atomic or molecular scale. None of these ideas are new. By now they have been fashionable topics for quite some time. But neither problem is satisfactorily resolved. The main problem is that uncontrolled approximations are used. As a result, we have very little idea about the validity and accuracy of these models. The main purpose of the center is to exploit the Mori-Zwanzig formalism to achieve the goals stated above. The Mori-Zwanzig formalism is a general tool for eliminating degrees of freedom from existing models. What is attractive about this formalism is its generality. The flip side is that it is also somewhat of a tautology. In order to make it a practical tool, one has to make approximations and this is where the trouble comes. Instead of focusing on Mori-Zwanzig formalism, the PI proposes to study the closely related renormalization group formalism (RG). Like Mori-Zwanzig, RG is a general strategy for eliminating small scale degrees of freedom in a model. Unlike Mori-Zwanzig, RG has the advantage that at least for some problems, one can write a PDE (partial differential equation) type of model for the exact RG flow, without any approximation. This gives us a solid starting point for discussing and analyzing approximations. The PI studied approximate RG flows for turbulent diffusion and related problems. However, extending it to more practical problems proved to be hard. The second line that the PI has pursued is to develop special numerical algorithms for studying particularly important events for dynamics processes in meso-scale systems. The PI succeeded in doing so for the problem of understanding phase transition in metallic systems.

97 MATHEMATICS AND COMPUTING↗

Representing local atomic environment using descriptors based on local correlations

Statistical learning of material properties is an emerging topic of research and has been tremendously successful in areas such as representing complex energy landscapes as well as in technologically relevant areas, like identification of better catalysts and electronic materials. However, analysis of large data sets to efficiently learn characteristic features of a complex energy landscape, for example, depends on the ability of descriptors to effectively screen different local atomic environments. Thus, discovering appropriate descriptors of bulk or defect properties and the functional dependence of such properties on these descriptors remains a difficult and tedious process. To this end, we develop here a framework to generate descriptors based on many-body correlations that can effectively capture intrinsic geometric features of the local environment of an atom. These descriptors are based on the spectrum of two-body, three-body, four-body, and higher order correlations between an atom and its neighbors and are evaluated by calculating the corresponding two-body, three-body, and four-body overlap integrals. They are invariant to global translation, global rotation, reflection, and permutations of atomic indices. By systematically testing the ability to capture the local atomic environment, it is shown that the local correlation descriptors are able to successfully reconstruct structures containing 10-25 atoms which was previously not possible.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The thermodynamics of a liquid-solid interface at extreme conditions: A model close-packed system up to 100 GPa

The first experimental insight into the nature of the liquid-solid interface occurred with the pioneering experiments of Turnbull, which simultaneously demonstrated both that metals could be deeply undercooled (and therefore had relatively large barriers to nucleation) and that the inferred interfacial free energy γ was linearly proportional to the enthalpy of fusion [D. Turnbull, J. Appl. Phys. 21, 1022 (1950)]. By an atomistic simulation of a model face-centered cubic system via adiabatic free energy dynamics, we extend Turnbull’s result to the realm of high pressure and demonstrate that the interfacial free energy, evaluated along the melting curve, remains linear with the bulk enthalpy of fusion, even up to 100 GPa. This linear dependence of γ on pressure is shown to be a consequence of the entropy dominating the free energy of the interface in conjunction with the fact that the entropy of fusion does not vary greatly along the melting curve for simple monoatomic metals. Here, based on this observation, it appears that large undercoolings in liquid metals can be achieved even at very high pressure. Therefore, nucleation rates at high pressure are expected to be non-negligible, resulting in observable solidification kinetics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Predicting phase behavior of grain boundaries with evolutionary search and machine learning

The study of grain boundary phase transitions is an emerging field until recently dominated by experiments. The major bottleneck in the exploration of this phenomenon with atomistic modeling has been the lack of a robust computational tool that can predict interface structure. Here we develop a computational tool based on evolutionary algorithms that performs efficient grand-canonical grain boundary structure search and we design a clustering analysis that automatically identifies different grain boundary phases. Its application to a model system of symmetric tilt boundaries in Cu uncovers an unexpected rich polymorphism in the grain boundary structures. We find new ground and metastable states by exploring structures with different atomic densities. Our results demonstrate that the grain boundaries within the entire misorientation range have multiple phases and exhibit structural transitions, suggesting that phase behavior of interfaces is likely a general phenomenon.

36 MATERIALS SCIENCE↗

Solvation and Dynamics of Sodium and Potassium in Ethylene Carbonate from ab Initio Molecular Dynamics Simulations

The development of sodium and potassium batteries offers a promising way to meet the scaling and cost challenges of energy storage. However, compared to Li + , several intrinsic properties of Na + and K + , including their solvation and dynamics in typical organic electrolytes utilized in battery applications, are less well-understood. Here in this paper, we report a systematic investigation of Na + and K + in ethylene carbonate (EC) using first-principles molecular dynamics simulations. Our simulations reveal significant differences in the solvation structure and dynamical properties of Na + and K + compared to Li + . We find that, in contrast to Li + which exhibits a well-defined first solvation shell, the larger Na+ and K+ ions show more disordered and flexible solvation structures. These differences in solvation were found to significantly influence the ion dynamics, leading to larger diffusion coefficients of Na + and K + compared to Li + . Our simulations also reveal a clear and interesting analog in the behavior of the ions in EC and aqueous environments, particularly in the specific ion effects on the solvent dynamics. Lastly, this work provides fundamental understanding of the intrinsic properties of Na + and K + in organic electrolytes, which may ultimately influence the intercalation mechanism at the electrode–electrolyte interface and therefore battery performance, lifetime, and safety.

25 ENERGY STORAGE↗

In situ TEM observation of FCC Ti formation at elevated temperatures

Pure Ti traditionally exhibits the hexagonal closed packed (HCP) crystallographic structure under ambient conditions and the body centered cubic (BCC) structure at elevated temperatures. In addition to these typical structures for Ti alloys, the presence of a face centered cubic (FCC) phase associated with thin films, interfaces, or high levels of plastic deformation has occasionally been reported. Here in this paper we show that small FCC precipitates form in freestanding thin foils during in situ transmission electron microscope (TEM) heating and we discuss the potential origins of the FCC phase in light of the in situ observations. This FCC phase was found to be stable upon cooling and under ambient conditions, which allowed us to explore its mechanical properties and stability via nanomechanical in situ TEM testing. It was found that FCC platelets within the HCP matrix phase were stable under mechanical deformation and exhibited similar mechanical deformation behavior as the parent HCP phase.

36 MATERIALS SCIENCE↗

Interfacial diffusion aided deformation during nanoindentation

Nanoindentation is commonly used to quantify the mechanical response of material surfaces. Despite its widespread use, a detailed understanding of the deformation mechanisms responsible for plasticity during these experiments has remained elusive. Nanoindentation measurements often show stress values close to a material’s ideal strength which suggests that dislocation nucleation and subsequent dislocation activity dominates the deformation. However, low strain-rate exponents and small activation volumes have also been reported which indicates high temperature sensitivity of the deformation processes. Using an order parameter aided temperature accelerated sampling technique called adiabatic free energy dynamics [J. B. Abrams and M. E. Tuckerman, J. Phys. Chem. B, 112, 15742 (2008)], and molecular dynamics we have probed the diffusive mode of deformation during nanoindentation. Localized processes such as surface vacancy and ad-atom pair formation, vacancy diffusion are found to play an important role during indentation. Our analysis suggests a change in the dominant deformation mode from dislocation mediated plasticity to diffusional flow at high temperatures, slow indentation rates and small indenter tip radii.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Exploring the free energy surface using ab initio molecular dynamics

Efficient exploration of the configuration space and identification of metastable structures are challenging from both computational as well as algorithmic perspectives. Here, we extend the recently proposed orderparameter aided temperature accelerated sampling schemes to efficiently and systematically explore free energy surfaces, and search for metastable states and reaction pathways within the framework of density functional theory based molecular dynamics. The sampling method is applied to explore the relevant parts of the configuration space in prototypical materials SiO 2 and Ti to identify the different metastable structures corresponding to different phases in these materials. In addition, we use the string method in collective variables to study the melting pathways in the high pressure cotunnite phase of SiO 2 and the hcp to fcc phase transition in Ti.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

A theoretical study of the stability of anionic defects in cubic ZrO 2 at extreme conditions

Using first principles density functional theory calculations, we present a study of the structure, mobility, and the thermodynamic stability of anionic defects in the high-temperature cubic phase of ZrO 2 . Our results suggest that the local structure of an oxygen interstitial depends on the charge state and the cubic symmetry of the anionic sublattice is unstable at 0 K. In addition, the oxygen interstitials and the vacancies exhibit symmetry breaking transitions to low-energy structures with tetragonal distortion of the oxygen sublattice at 0 K. However, the vibrational entropy stabilizes the defect structures with cubic symmetry at 2600–2980 K. The formation free energies of the anionic defects and Gibbs free energy changes associated with different defect reactions are calculated by including the vibrational free energy contributions and the effect of pressure on these defect structures. By analyzing the defect chemistry, we obtain the defect concentrations at finite temperature and pressure conditions using the zero temperature ab initio results as input and find that at low oxygen partial pressures, neutral oxygen vacancies are most dominant and at high oxygen partial pressures, doubly charged anionic defects are dominant. As a result, the relevance of the results to the thermal protective coating capabilities of zirconium-based ceramic composites is elucidated.

36 MATERIALS SCIENCE↗