Publications

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Large-scale, high-throughput production of nano-structured materials (i.e. nanomanufacturing) is a strategic area in manufacturing, with markets projected to exceed $1T by 2015. Nanomanufacturing is still in its infancy; process/product developments are costly and only touch on potential opportunities enabled by growing nanoscience discoveries. The greatest promise for high-volume manufacturing lies in age-old coating and imprinting operations. For materials with tailored nm-scale structure, imprinting/embossing must be achieved at high speeds (roll-to-roll) and/or over large areas (batch operation) with feature sizes less than 100 nm. Dispersion coatings with nanoparticles can also tailor structure through self- or directed-assembly. Layering films structured with these processes have tremendous potential for efficient manufacturing of microelectronics, photovoltaics and other topical nano-structured devices. This project is designed to perform the requisite R and D to bring Sandia's technology base in computational mechanics to bear on this scale-up problem. Project focus is enforced by addressing a promising imprinting process currently being commercialized.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

The kinetic Monte Carlo method and its variants are powerful tools for modeling materials at the mesoscale, meaning at length and time scales in between the atomic and continuum. We have completed a 3 year LDRD project with the goal of developing a parallel kinetic Monte Carlo capability and applying it to materials modeling problems of interest to Sandia. In this report we give an overview of the methods and algorithms developed, and describe our new open-source code called SPPARKS, for Stochastic Parallel PARticle Kinetic Simulator. We also highlight the development of several Monte Carlo models in SPPARKS for specific materials modeling applications, including grain growth, bubble formation, diffusion in nanoporous materials, defect formation in erbium hydrides, and surface growth and evolution.

We have conducted a molecular dynamics (MD) simulation study of water confined between methyl-terminated and carboxyl-terminated alkylsilane self-assembled monolayers (SAMs) on amorphous silica substrates. In doing so, we have investigated the dynamic and structural behavior of the water molecules when compressed to loads ranging from 20 to 950 MPa for two different amounts of water (27 and 58 water molecules/nm{sup 2}). Within the studied range of loads, we observe that no water molecules penetrate the hydrophobic region of the carboxyl-terminated SAMs. However, we observe that at loads larger than 150 MPa water molecules penetrate the methyl-terminated SAMs and form hydrogen-bonded chains that connect to the bulk water. The diffusion coefficient of the water molecules decreases as the water film becomes thinner and pressure increases. When compared to bulk diffusion coefficients of water molecules at the various loads, we found that the diffusion coefficients for the systems with 27 water molecules/nm{sup 2} are reduced by a factor of 20 at low loads and by a factor of 40 at high loads, while the diffusion coefficients for the systems with 58 water molecules/nm{sup 2} are reduced by a factor of 25 at all loads.

Abstract not provided.

Abstract not provided.

We apply density functional theory (DFT) and the DFT+U technique to study the adsorption of transition metal porphine molecules on atomistically flat Au(111) surfaces. DFT calculations using the Perdew?Burke?Ernzerhof exchange correlation functional correctly predict the palladium porphine (PdP) low-spin ground state. PdP is found to adsorb preferentially on gold in a flat geometry, not in an edgewise geometry, in qualitative agreement with experiments on substituted porphyrins. It exhibits no covalent bonding to Au(111), and the binding energy is a small fraction of an electronvolt. The DFT+U technique, parametrized to B3LYP-predicted spin state ordering of the Mn d-electrons, is found to be crucial for reproducing the correct magnetic moment and geometry of the isolated manganese porphine (MnP) molecule. Adsorption of Mn(II)P on Au(111) substantially alters the Mn ion spin state. Its interaction with the gold substrate is stronger and more site-specific than that of PdP. The binding can be partially reversed by applying an electric potential, which leads to significant changes in the electronic and magnetic properties of adsorbed MnP and 0.1 {angstrom} changes in the Mn-nitrogen distances within the porphine macrocycle. We conjecture that this DFT+U approach may be a useful general method for modeling first-row transition metal ion complexes in a condensed-matter setting.

In this study, we perform molecular dynamics simulations of adhesive contact and friction between alkylsilane Si(OH){sub 3}(CX{sub 2}){sub 10}CX{sub 3} and alkoxylsilane Si(OH){sub 2}(CX{sub 2}){sub 10}CX{sub 3} (where X = H or F) self-assembled monolayers (SAMs) on an amorphous silica substrate. The alkylsilane SAMs are primarily hydrogen-bonded or physisorbed to the surface. The alkoxylsilane SAMs are covalently bonded or chemisorbed to the surface. Previously, we studied the chemisorbed systems. In this work, we study the physisorbed systems and compare the tribological properties with the chemisorbed systems. Furthermore, we examine how water at the interface of the SAMs and substrate affects the tribological properties of the physisorbed systems. When less than a third of a monolayer is present, very little difference in the microscopic friction coefficient {mu} or shear stresses is observed. For increasing amounts of water, the values of {mu} and the shear stresses decrease; this effect is somewhat more pronounced for fluorocarbon alkylsilane SAMs than for the hydrocarbon SAMs. The observed decrease in friction is a consequence of a slip plane that occurs in the water as the amount of water is increased. We studied the frictional behavior using relative shear velocities ranging from v = 2 cm/s to 2 m/s. Similar to previously reported results for alkoxylsilane SAMs, the values of the measured stress and {mu} for the alkylsilane SAM systems decrease monotonically with v.

Abstract not provided.

This report summarizes the work performed as part of a one-year LDRD project, 'Evolutionary Complexity for Protection of Critical Assets.' A brief introduction is given to the topics of genetic algorithms and genetic programming, followed by a discussion of relevant results obtained during the project's research, and finally the conclusions drawn from those results. The focus is on using genetic programming to evolve solutions for relatively simple algebraic equations as a prototype application for evolving complexity in computer codes. The results were obtained using the lil-gp genetic program, a C code for evolving solutions to user-defined problems and functions. These results suggest that genetic programs are not well-suited to evolving complexity for critical asset protection because they cannot efficiently evolve solutions to complex problems, and introduce unacceptable performance penalties into solutions for simple ones.

Confinement within the nanoscale pores of a zeolite strongly modifies the behavior of small molecules. Typical of many such interesting and important problems, realistic modeling of this phenomena requires simultaneously capturing the detailed behavior of chemical bonds and the possibility of collective dynamics occurring in a complex unit cell (672 atoms in the case of Zeolite-4A). Classical simulations alone cannot reliably model the breaking and formation of chemical bonds, while quantum methods alone are incapable of treating the extended length and time scales characteristic of complex dynamics. We have developed a robust and efficient model in which a small region treated with the Kohn-Sham density functional theory is embedded within a larger system represented with classical potentials. This model has been applied in concert with first-principles electronic structure calculations and classical molecular dynamics and Monte Carlo simulations to study the behavior of water, ammonia, the hydroxide ion, and the ammonium ion in Zeolite-4a. Understanding this behavior is important to the predictive modeling of the aging of Zeolite-based desiccants. In particular, we have studied the absorption of these molecules, interactions between water and the ammonium ion, and reactions between the hydroxide ion and the zeolite cage. We have shown that interactions with the extended Zeolite cage strongly modifies these local chemical phenomena, and thereby we have proven out hypothesis that capturing both local chemistry and collective phenomena is essential to realistic modeling of this system. Based on our results, we have been able to identify two possible mechanisms for the aging of Zeolite-based desiccants.

Abstract not provided.

We have developed force fields for the calculation of adsorption of NH{sub 3}, CO{sub 2}, and H{sub 2}O on zeolite 4A by performing Gibbs ensemble Monte Carlo simulations to fit experimental isotherms at 298 K. The calculated NH{sub 3} and CO{sub 2} isotherms are in excellent agreement with experimental data over a wide range of temperatures and several orders of magnitude in pressure. We have calculated isotherms for H{sub 2}O in 4A using two different models and have found that H{sub 2}O saturates zeolite 4A even at pressures as low as 0.01 kPa for the range of temperatures studied. We have studied the geometry of the adsorption sites and their dependence on loading. At low pressures, CO{sub 2} molecules adsorb with their longitudinal axis pointing toward the center of the supercage, whereas at higher pressures, the two oxygen atoms are equidistant from the Na atom in the binding site.

There is currently a great need for solid state lasers that emit in the infrared. Whether or not conjugated polymers that emit in the IR can be synthesized is an interesting theoretical challenge. We show that emission in the IR can be achieved in designer polymers in which the effective Coulomb correlation is smaller than that in existing systems. We also show that the structural requirement for having small effective Coulomb correlations is that there exist transverse {pi}--conjugation over a few bonds in addition to longitudinal conjugation with large conjugation lengths.

Abstract not provided.

Abstract not provided.

We present the results of Molecular Dynamics and Monte Carlo simulations of n-butane and isobutane in silicalite. We begin with a comparison of the bulk adsorption and diffusion properties for two different parameterizations of the interaction potential between the hydrocarbon species, both of which have been shown to reproduce experimental gas-liquid coexistence curves. We examine diffusion as a function of the loading of the zeolite, as well as the temperature dependence of the diffusion constant at loading and for infinite dilution. Both force fields give accurate descriptions of bulk properties. We continue with simulations in which interfaces are formed between single component gases and the zeolite. After reaching equilibrium, we examine the dynamics of exchange between the bulk gas and the zeolite. In particular, we examine the average time spent in the adsorption layer by molecules as they enter the zeolite from the gas in an attempt to probe the microscopic origins of the surface barrier. The microscopic barrier is found to be insignificant for experimental systems. Finally, we calculate the permeability of the zeolite for n-butane and isobutane as a function of pressure. Our results underestimate the experimental results by an order of magnitude, indicating a strong effect from the surface barrier in these simulations. Our simulations are performed for a number of different gas temperatures and pressures, covering a wide range of state points.