Molybdenum disulfide (MoS2) is a lamellar solid lubricant often used in aerospace applications because of its extremely low friction coefficient (~0.01) in inert environments. The lubrication performance of MoS2 is significantly impaired by exposure to even small amounts of water and oxygen, and the mechanisms behind this remain poorly understood. Here we use density functional theory calculations to study the binding of water on MoS2 sheets with and without defects. In general, we find that pristine MoS2 is slightly hydrophilic but that defects greatly increase the binding affinity for water. Intercalated water disrupts the crystal structure of bulk MoS2 due to the limited space between lamellae (~3.4 Å), and this leads to generally unfavorable adsorption, except in the cases where water molecules are located on the sites of sulfur vacancies. We also find that water adsorption is more favorable directly below a surface layer of MoS2 compared to in the bulk.
The friction behavior of metals is directly linked to the mechanisms that accommodate deformation. We examine the links between mechanisms of strengthening, deformation, and the wide range of friction behaviors that are exhibited by shearing metal interfaces. Specifically, the focus is on understanding the shear strength of nanocrystalline and nanostructured metals, and conditions that lead to low friction coefficients. Grain boundary sliding and the breakdown of Hall–Petch strengthening at the shearing interface are found to generally and predictably explain the low friction of these materials. While the following is meant to serve as a general discussion of the strength of metals in the context of tribological applications, one important conclusion is that tribological research methods also provide opportunities for probing the fundamental properties and deformation mechanisms of metals.
In this work, scratch and nanoindentation testing was used to determine hardness, fracture toughness, strain rate sensitivity, and activation volumes on additively manufactured graded and uniform Ni-Nb bulk specimens. Characterization showed the presence of a two phase system consisting of Ni3Nb and Ni6Nb7 intermetallics. Intermetallics were multimodal in nature, having grain and cell sizes spanning from a few nanometers to 10s of micrometers. The unique microstructure resulted in impressively high hardness, up to 20 GPa in the case of the compositionally graded sample. AM methods with surface deformation techniques are a useful way to rapidly probe material properties and alloy composition space.
In this work, scratch and nanoindentation testing was used to determine hardness, fracture toughness, strain rate sensitivity, and activation volumes on additively manufactured graded and uniform Ni-Nb bulk specimens. Characterization showed the presence of a two phase system consisting of Ni3Nb and Ni6Nb7 intermetallics. Intermetallics were multimodal in nature, having grain and cell sizes spanning from a few nanometers to 10s of micrometers. The unique microstructure resulted in impressively high hardness, up to 20 GPa in the case of the compositionally graded sample. AM methods with surface deformation techniques are a useful way to rapidly probe material properties and alloy composition space.
Heterojunctions of semiconductors and metals are the fundamental building blocks of modern electronics. Coherent heterostructures between dissimilar materials can be achieved by composition, doping, or heteroepitaxy of chemically different elements. Here, we report the formation of coherent single-layer 1H-1T MoS2 heterostructures by mechanical exfoliation on Au(111), which are chemically homogeneous with matched lattices but show electronically distinct semiconducting (1H phase) and metallic (1T phase) character, with the formation of these heterojunctions attributed to a combination of lattice strain and charge transfer. The exfoliation approach employed is free of tape residues usually found in many exfoliation methods and yields single-layer MoS2 with millimeter (mm) size on the Au surface. Raman spectroscopy, X-ray photoelectron spectroscopy (XPS), atomic force microscopy (AFM), scanning tunneling microscopy (STM), and scanning tunneling spectroscopy (STS) have collectively been employed to elucidate the structural and electronic properties of MoS2 monolayers on Au substrates. Bubbles in the MoS2 formed by the trapping of ambient adsorbates beneath the single layer during deposition, have also been observed and characterized. Our work here provides a basis to produce two-dimensional heterostructures which represent potential candidates for future electronic devices.
We present evidence of inverse Hall-Petch behavior for a single-phase high entropy alloy (CoCrFeMnNi) in ultra-high vacuum and show that it is associated with low friction coefficients (~0.3). Grain size measurements by STEM validate a recently proposed dynamic amorphization model that accurately predicts grain size-dependent shear strength in the inverse Hall-Petch regime. Wear rates in the initially soft (coarse grained) material were shown to be remarkably low (~10–6 mm3/N-m), the lowest for any HEA tested in an inert environment where oxidation and the formation of mixed metal-oxide films is mitigated. The combined high wear resistance and low friction are linked to the formation of an ultra-nanocrystalline near-surface layer. The dynamic amorphization model was also used to predict an average high angle grain boundary energy (0.87 J/m2). This value was used to explain cavitation-induced nanoporosity found in the highly deformed surface layer, a phenomenon that has been linked to superplasticity.
Low friction is demonstrated with pure polycrystalline tantalum sliding contacts in both molecular dynamics simulations and ultrahigh vacuum experiments. This phenomenon is shown to be correlated with deformation occurring primarily through grain boundary sliding and can be explained using a recently developed predictive model for the shear strength of metals. Specifically, low friction is associated with grain sizes at the interface being smaller than a critical, material-dependent value, where a crossover from dislocation mediated plasticity to grain-boundary sliding occurs. Low friction is therefore associated with inverse Hall-Petch behavior and softening of the interface. Direct quantitative comparisons between experiments and atomistic calculations are used to illustrate the accuracy of the predictions.
We present a theoretical model that predicts the peak strength of polycrystalline metals based on the activation energy (or stress) required to cause deformation via amorphization. Building on extensive earlier work, this model is based purely on materials properties, requires no adjustable parameters, and is shown to accurately predict the strength of four exemplar metals (fcc, bcc, and hcp, and an alloy). This framework reveals new routes for design of more complex high-strength materials systems, such as compositionally complex alloys, multiphase systems, nonmetals, and composite structures.
Additive Manufacturing (AM) offers the opportunity to transform design, manufacturing, and qualification with its unique capabilities. AM is a disruptive technology, allowing the capability to simultaneously create part and material while tightly controlling and monitoring the manufacturing process at the voxel level, with the inherent flexibility and agility in printing layer-by-layer. AM enables the possibility of measuring critical material and part parameters during manufacturing, thus changing the way we collect data, assess performance, and accept or qualify parts. It provides an opportunity to shift from the current iterative design-build-test qualification paradigm using traditional manufacturing processes to design-by-predictivity where requirements are addressed concurrently and rapidly. The new qualification paradigm driven by AM provides the opportunity to predict performance probabilistically, to optimally control the manufacturing process, and to implement accelerated cycles of learning. Exploiting these capabilities to realize a new uncertainty quantification-driven qualification that is rapid, flexible, and practical is the focus of this paper.
Intermetallic alloys possess exceptional soft magnetic properties, including high permeability, low coercivity, and high saturation induction, but exhibit poor mechanical properties that make them impractical to bulk process and use at ideal compositions. We used laser-based Additive Manufacturing to process traditionally brittle Fe–Co and Fe–Si alloys in bulk form without macroscopic defects and at near-ideal compositions for electromagnetic applications. The binary Fe–50Co, as a model material, demonstrated simultaneous high strength (600–700 MPa) and high ductility (35%) in tension, corresponding to a ∼300% increase in strength and an order-of-magnitude improvement in ductility relative to conventionally processed material. Atomic-scale toughening and strengthening mechanisms, based on engineered multiscale microstructures, are proposed to explain the unusual combination of mechanical properties. This work presents an instance in which metal Additive Manufacturing processes are enabling, rather than limiting, the development of higher-performance alloys.
Stress intensity factors (SIFs) are used in continuum fracture mechanics to quantify the stress fields surrounding a crack in a homogeneous material in the linear elastic regime. Critical values of the SIFs define an intrinsic measure of the resistance of a material to propagate a crack. At atomic scales, however, fracture occurs as a series of atomic bonds breaking, differing from the continuum description. As a consequence, a formal analog of the continuum SIFs calculated from atomistic simulations can have spatially localized, microstructural contributions that originate from varying bond configurations. The ability to characterize fracture at the atomic scale in terms of the SIFs offers both an opportunity to probe the effects of chemistry, as well as how the addition of a microstructural component affects the accuracy. We present a novel numerical method to determine SIFs from molecular dynamics (MD) simulations. The accuracy of this approach is first examined for a simple model, and then applied to atomistic simulations of fracture in amorphous silica. MD simulations provide time and spatially dependent SIFs, with results that are shown to be in good agreement with experimental values for fracture toughness in silica glass.
Nanocrystalline metals offer significant improvements in structural performance over conventional alloys. However, their performance is limited by grain boundary instability and limited ductility. Solute segregation has been proposed as a stabilization mechanism, however the solute atoms can embrittle grain boundaries and further degrade the toughness. In the present study, we confirm the embrittling effect of solute segregation in Pt-Au alloys. However, more importantly, we show that inhomogeneous chemical segregation to the grain boundary can lead to a new toughening mechanism termed compositional crack arrest. Energy dissipation is facilitated by the formation of nanocrack networks formed when cracks arrested at regions of the grain boundaries that were starved in the embrittling element. This mechanism, in concert with triple junction crack arrest, provides pathways to optimize both thermal stability and energy dissipation. A combination of in situ tensile deformation experiments and molecular dynamics simulations elucidate both the embrittling and toughening processes that can occur as a function of solute content.
Diamond-like carbon (DLC) films were tribochemically formed from ambient hydrocarbons on the surface of a highly stable nanocrystalline Pt-Au alloy. A sliding contact between an alumina sphere and Pt-Au coated steel exhibited friction coefficients as low as μ = 0.01 after dry sliding in environments containing trace (ppb) organics. Ex situ analysis indicated that the change in friction coefficient was due to the formation of amorphous carbon films, and Raman spectroscopy and elastic recoil analysis showed that these films consist of sp2/sp3 amorphous carbon with as much as 20% hydrogen. Transmission electron microscopy indicated these films had thicknesses exceeding 100 nm, and were enhanced by the incorporation of worn Pt-Au nanoparticles. The result was highly wear-resistant, low-friction DLC/Pt-Au nanocomposites. Atomistic simulations of hydrocarbons under shear between rigid Pt slabs using a reactive force field showed stress-induced changes in bonding through chain scission, a likely route towards the formation of these coatings. This novel demonstration of in situ tribochemical formation of self-lubricating films has significant impact potential in a wide range of engineering applications.
This SAND report fulfills the final report requirement for the Born Qualified Grand Challenge LDRD. Born Qualified was funded from FY16-FY18 with a total budget of ~$13M over the 3 years of funding. Overall 70+ staff, Post Docs, and students supported this project over its lifetime. The driver for Born Qualified was using Additive Manufacturing (AM) to change the qualification paradigm for low volume, high value, high consequence, complex parts that are common in high-risk industries such as ND, defense, energy, aerospace, and medical. AM offers the opportunity to transform design, manufacturing, and qualification with its unique capabilities. AM is a disruptive technology, allowing the capability to simultaneously create part and material while tightly controlling and monitoring the manufacturing process at the voxel level, with the inherent flexibility and agility in printing layer-by-layer. AM enables the possibility of measuring critical material and part parameters during manufacturing, thus changing the way we collect data, assess performance, and accept or qualify parts. It provides an opportunity to shift from the current iterative design-build-test qualification paradigm using traditional manufacturing processes to design-by-predictivity where requirements are addressed concurrently and rapidly. The new qualification paradigm driven by AM provides the opportunity to predict performance probabilistically, to optimally control the manufacturing process, and to implement accelerated cycles of learning. Exploiting these capabilities to realize a new uncertainty quantification-driven qualification that is rapid, flexible, and practical is the focus of this effort.
Recent work suggests that thermally stable nanocrystallinity in metals is achievable in several binary alloys by modifying grain boundary energies via solute segregation. The remarkable thermal stability of these alloys has been demonstrated in recent reports, with many alloys exhibiting negligible grain growth during prolonged exposure to near-melting temperatures. Pt–Au, a proposed stable alloy consisting of two noble metals, is shown to exhibit extraordinary resistance to wear. Ultralow wear rates, less than a monolayer of material removed per sliding pass, are measured for Pt–Au thin films at a maximum Hertz contact stress of up to 1.1 GPa. This is the first instance of an all-metallic material exhibiting a specific wear rate on the order of 10−9 mm3 N−1 m−1, comparable to diamond-like carbon (DLC) and sapphire. Remarkably, the wear rate of sapphire and silicon nitride probes used in wear experiments are either higher or comparable to that of the Pt–Au alloy, despite the substantially higher hardness of the ceramic probe materials. High-resolution microscopy shows negligible surface microstructural evolution in the wear tracks after 100k sliding passes. Mitigation of fatigue-driven delamination enables a transition to wear by atomic attrition, a regime previously limited to highly wear-resistant materials such as DLC.
We demonstrate inverse Hall-Petch behavior (softening) in pure copper sliding contacts at cryogenic temperatures. By kinetically limiting grain growth, it is possible to generate a quasi-stable ultra-nanocrystalline surface layer with reduced strength. In situ electrical contact resistance measurements were used to determine grain size evolution at the interface, in agreement with reports of softening in highly nanotwinned copper. We also show evidence of a direct correlation between surface grain size and friction coefficient, validating a model linking friction in pure metals and the transition from dislocation mediated plasticity to grain boundary sliding.
A correlation is established between the macro-scale friction regimes of metals and a transition between two dominant atomistic mechanisms of deformation. Metals tend to exhibit bi-stable friction behavior—low and converging or high and diverging. These general trends in behavior are shown to be largely explained using a simplified model based on grain size evolution, as a function of contact stress and temperature, and are demonstrated for self-mated pure copper and gold sliding contacts. Specifically, the low-friction regime (where µ < 0.5) is linked to the formation of ultra-nanocrystalline surface films (10–20 nm), driving toward shear accommodation by grain boundary sliding. Above a critical combination of stress and temperature—demonstrated to be a material property—shear accommodation transitions to dislocation dominated plasticity and high friction, with µ > 0.5. We utilize a combination of experimental and computational methods to develop and validate the proposed structure–property relationship. This quantitative framework provides a shift from phenomenological to mechanistic and predictive fundamental understanding of friction for crystalline materials, including engineering alloys.
The friction behavior of ultra-nanocrystalline Ni-W coatings was investigated. A critical stress threshold was identified below which friction remained low, and above which a time-dependent evolution toward higher friction behavior occurred. Founded on established plasticity models we propose a correlation between surface grain size and applied stress that can be used to predict the critical stress separating the two friction regimes. This interpretation of plasticity models suggests that macro-scale low and high friction regimes are respectively associated with the nano-scale mechanisms of grain boundary and dislocation-mediated plasticity.
Aerosol deposition (AD) is a solid-state deposition technology that has been developed to fabricate ceramic coatings nominally at room temperature. Sub-micron ceramic particles accelerated by pressurized gas impact, deform, and consolidate on substrates under vacuum. Ceramic particle consolidation in AD coatings is highly dependent on particle deformation and bonding; these behaviors are not well understood. In this work, atomistic simulations and in situ micro-compressions in the scanning electron microscope, and the transmission electron microscope (TEM) were utilized to investigate fundamental mechanisms responsible for plastic deformation/fracture of particles under applied compression. Results showed that highly defective micron-sized alumina particles, initially containing numerous dislocations or a grain boundary, exhibited no observable shape change before fracture/fragmentation. Simulations and experimental results indicated that particles containing a grain boundary only accommodate low strain energy per unit volume before crack nucleation and propagation. In contrast, nearly defect-free, sub-micron, single crystal alumina particles exhibited plastic deformation and fracture without fragmentation. Dislocation nucleation/motion, significant plastic deformation, and shape change were observed. Simulation and TEM in situ micro-compression results indicated that nearly defect-free particles accommodate high strain energy per unit volume associated with dislocation plasticity before fracture. The identified deformation mechanisms provide insight into feedstock design for AD.
This work presents molecular dynamics and Monte Carlo simulations aimed at developing an understanding of the formation of core-shell Cu-Ag nanoparticles. The effects of surface and interfacial energies were considered and used to form a phenomenological model that calculates the energy gained upon the formation of a core-shell structure from two previously distinct, non-interacting nanoparticles. In the majority of cases, the core-shell structure was found to be energetically favored. Specifically, the difference in energy as a function of the radii of the individual Cu and Ag particles was examined, with the assumption that a core-shell structure forms. In general, it was found that the energetic gain from forming such a structure increased with increasing size of the initial Ag particle. This result was interpreted as a result of the reduction in surface energy. For two separate particles, both Cu and Ag contribute to the surface energy; however, for a core-shell structure, the only contribution to the surface energy is from the Ag shell and the Cu contribution is changed to a Cu-Ag interfacial energy, which is always smaller.
The ability to integrate ceramics with other materials has been limited due to high temperature (>800degC) ceramic processing. Recently, researchers demonstrated a novel process , aerosol deposition (AD), to fabricate ceramic films at room temperature (RT). In this process, sub - micro n sized ceramic particles are accelerated by pressurized gas, impacted on the substrate, plastically deformed, and form a dense film under vacuum. This AD process eliminates high temperature processing thereby enabling new coatings and device integration, in which ceramics can be deposited on metals, plastics, and glass. However, k nowledge in fundamental mechanisms for ceramic particle s to deform and form a dense ceramic film is still needed and is essential in advancing this novel RT technology. In this wo rk, a combination of experimentation and atomistic simulation was used to determine the deformation behavior of sub - micron sized ceramic particle s ; this is the first fundamental step needed to explain coating formation in the AD process . High purity, singl e crystal, alpha alumina particles with nominal size s of 0.3 um and 3.0 um were examined. Particle characterization, using transmission electron microscopy (TEM ), showed that the 0.3 u m particles were relatively defect - free single crystals whereas 3.0 u m p articles were highly defective single crystals or particles contained low angle grain boundaries. Sub - micron sized Al 2 O 3 particles exhibited ductile failure in compression. In situ compression experiments showed 0.3um particles deformed plastically, fractured, and became polycrystalline. Moreover, dislocation activit y was observed within the se particles during compression . These sub - micron sized Al 2 O 3 particles exhibited large accum ulated strain (2 - 3 times those of micron - sized particles) before first fracture. I n agreement with the findings from experimentation , a tomistic simulation s of nano - Al 2 O 3 particles showed dislocation slip and significant plastic deformation during compressi on . On the other hand, the micron sized Al 2 O 3 particles exhibited brittle f racture in compression. In situ compression experiments showed 3um Al 2 O 3 particles fractured into pieces without observable plastic deformation in compression. Particle deformation behaviors will be used to inform Al 2 O 3 coating deposition parameters and particle - particle bonding in the consolidated Al 2 O 3 coatings.
We found that the wetting and spreading of molten filler materials (pure Al, pure Ag, and AgAl alloys) on a Kovar ™ (001) substrate was studied with molecular dynamics simulations. A suite of different simulations was used to understand the effects on spreading rates due to alloying as well as reactions with the substrate. Moreover, the important conclusion is that the presence of Al in the alloy enhances the spreading of Ag, while the Ag inhibits the spreading of Al.
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.
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.
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.
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.
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.
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.