Predicting the Mechanical Properties of Carbon-Based Materials Using Molecular Dynamics

The development of innovative carbon-based materials can be greatly facilitated by molecular modeling techniques. Although molecular modeling has been used extensively to predict elastic properties of materials, modeling of more complex phenomenon such as fracture has only recently been possible with the development of new force fields such as ReaxFF, which is used in this work. It is not fully understood what molecular modeling parameters such as thermostat type, thermostat coupling, time step, system size, and strain rate are required for accurate modeling of fracture. Selection of modeling parameters to model fracture can be difficult and non-intuitive compared to modeling elastic properties using traditional force fields, and the errors generated by incorrect parameters may be non-obvious. These molecular modeling parameters are systematically investigated and their effects on the fracture of well-known carbon materials are analyzed. It is determined that for coupling coefficients of 250 fs and greater do not result in substantial differences in the stress-strain response of the materials using any thermostat type. A time step of 0.5 fs of smaller is required for accurate results. Strain rates greater than 2.2 ns-1 are sufficient to obtain repeatable results with slower strain rates for the materials studied. The results of this study indicate that further refinement of the Chenoweth parameter set is required to accurately predict the mechanical response of carbon-based systems. The ReaxFF has been used extensively to model systems in which bond breaking and formation occur. In particular ReaxFF has been used to model reactions of small molecules. Some elastic and fracture properties have been successfully modeled using ReaxFF in materials such as silicon and some metals. However, it is not clear if current parameterizations for ReaxFF are able to accurately reproduce the elastic and fracture properties of carbon materials. The stress-strain response of a new ReaxFF parameterization is compared to the previous parameterization and density functional theory results for well-known carbon materials. The new ReaxFF parameterization makes xv substantial improvements to the predicted mechanical response of carbon materials, and is found to be suitable for modeling the mechanical response of carbon materials. Finally, a new material composed of carbon nanotubes within an amorphous carbon (AC) matrix is modeled using the ReaxFF. Various parameters that may be experimentally controlled are investigated such as nanotube bundling, comparing multi-walled nanotube with single-walled nanotubes, and degree of functionalization of the nanotubes. Elastic and fracture properties are investigated for the composite systems and compared to results of pure-nanotube and pure-AC models. It is found that the arrangement of the nanotubes and degree of crosslinking may substantially affect the properties of the systems, particularly in the transverse directions.

Integration of computational models and experimental characterization to study internal frost damage in cementitious materials

The objective of this doctoral research is to investigate the internal frost damage due to crystallization pore pressure in porous cement-based materials by developing computational and experimental characterization tools. As an essential component of the U.S. infrastructure system, the durability of concrete has significant impact on maintenance costs. In cold climates, freeze-thaw damage is a major issue affecting the durability of concrete. The deleterious effects of the freeze-thaw cycle depend on the microscale characteristics of concrete such as the pore sizes and the pore distribution, as well as the environmental conditions. Recent theories attribute internal frost damage of concrete is caused by crystallization pore pressure in the cold environment. The pore structures have significant impact on freeze-thaw durability of cement/concrete samples. The scanning electron microscope (SEM) and transmission X-ray microscopy (TXM) techniques were applied to characterize freeze-thaw damage within pore structure. In the microscale pore system, the crystallization pressures at sub-cooling temperatures were calculated using interface energy balance with thermodynamic analysis. The multi-phase Extended Finite Element Modeling (XFEM) and bilinear Cohesive Zone Modeling (CZM) were developed to simulate the internal frost damage of heterogeneous cement-based material samples. The fracture simulation with these two techniques were validated by comparing the predicted fracture behavior with the captured damage from compact tension (CT) and single-edge notched beam (SEB) bending tests. The study applied the developed computational tools to simulate the internal frost damage caused by ice crystallization with the two dimensional (2-D) SEM and three dimensional (3-D) reconstructed SEM and TXM digital samples. The pore pressure calculated from thermodynamic analysis was input for model simulation. The 2-D and 3-D bilinear CZM predicted the crack initiation and propagation within cement paste microstructure. The favorably predicted crack paths in concrete/cement samples indicate the developed bilinear CZM techniques have the ability to capture crack nucleation and propagation in cement-based material samples with multiphase and associated interface. By comparing the computational prediction with the actual damaged samples, it also indicates that the ice crystallization pressure is the main mechanism for the internal frost damage in cementitious materials.

PARAMETRIC STUDY OF REAXFF SIMULATION PARAMETERS FOR MOLECULAR DYNAMICS MODELING OF REACTIVE CARBON GASES

Abstract The development of innovative carbon-based materials can be greatly facilitated by molecular modeling techniques. Although the Reax Force Field (ReaxFF) can be used to simulate the chemical behavior of carbon-based systems, the simulation settings required for accurate predictions have not been fully explored. Using the ReaxFF, molecular dynamics (MD) simulations are used to simulate the chemical behavior of pure carbon and hydrocarbon reactive gases that are involved in the formation of carbon structures such as graphite, buckyballs, amorphous carbon, and carbon nanotubes. It is determined that the maximum simulation time step that can be used in MD simulations with the ReaxFF is dependent on the simulated temperature and selected parameter set, as are the predicted reaction rates. It is also determined that different carbon-based reactive gases react at different rates, and that the predicted equilibrium structures are generally the same for the different ReaxFF parameter sets, except in the case of the predicted formation of large graphitic structures with the Chenoweth parameter set under specific conditions.

Processing and characterization of PbSnTe-based thermoelectric materials made by mechanical alloying

The research reported in this dissertation investigates the processes required to mechanically alloy Pb1-xSnxTe and AgSbTe2 and a method of combining these two end compounds to result in (y)(AgSbTe2)–(1 - y)(Pb1-xSnxTe) thermoelectric materials for power generation applications. In general, traditional melt processing of these alloys has employed high purity materials that are subjected to time and energy intensive processes that result in highly functional material that is not easily reproducible. This research reports the development of mechanical alloying processes using commercially available 99.9% pure elemental powders in order to provide a basis for the economical production of highly functional thermoelectric materials. Though there have been reports of high and low ZT materials fabricated by both melt alloying and mechanical alloying, the processing-structure-properties-performance relationship connecting how the material is made to its resulting functionality is poorly understood. This is particularly true for mechanically alloyed material, motivating an effort to investigate bulk material within the (y)(AgSbTe2)–(1 - y)(Pb1-xSnx- Te) system using the mechanical alloying method. This research adds to the body of knowledge concerning the way in which mechanical alloying can be used to efficiently produce high ZT thermoelectric materials. The processes required to mechanically alloy elemental powders to form Pb1-xSnxTe and AgSbTe2 and to subsequently consolidate the alloyed powder is described. The composition, phases present in the alloy, volume percent, size and spacing of the phases are reported. The room temperature electronic transport properties of electrical conductivity, carrier concentration and carrier mobility are reported for each alloy and the effect of the presence of any secondary phase on the electronic transport properties is described. An mechanical mixing approach for incorporating the end compounds to result in (y)(AgSbTe2)–(1-y)(Pb1-xSnxTe) is described and when 5 vol.% AgSbTe2 was incorporated was found to form a solid solution with the Pb1-xSnxTe phase. An initial attempt to change the carrier concentration of the Pb1-xSnxTe phase was made by adding excess Te and found that the carrier density of the alloys in this work are not sensitive to excess Te. It has been demonstrated using the processing techniques reported in this research that this material system, when appropriately doped, has the potential to perform as highly functional thermoelectric material.

PROPERTIES AND STRUCTURES OF Li-N BASED HYDROGEN STORAGE MATERIALS

Traditional transportation fuel, petroleum, is limited and nonrenewable, and it also causes pollutions. Hydrogen is considered one of the best alternative fuels for transportation. The key issue for using hydrogen as fuel for transportation is hydrogen storage. Lithium nitride (Li3N) is an important material which can be used for hydrogen storage. The decompositions of lithium amide (LiNH2) and lithium imide (Li2NH) are important steps for hydrogen storage in Li3N. The effect of anions (e.g. Cl-) on the decomposition of LiNH2 has never been studied. Li3N can react with LiBr to form lithium nitride bromide Li13N4Br which has been proposed as solid electrolyte for batteries. The decompositions of LiNH2 and Li2NH with and without promoter were investigated by using temperature programmed decomposition (TPD) and X-ray diffraction (XRD) techniques. It was found that the decomposition of LiNH2 produced Li2NH and NH3 via two steps: LiNH2 into a stable intermediate species (Li1.5NH1.5) and then into Li2NH. The decomposition of Li2NH produced Li, N2 and H2 via two steps: Li2NH into an intermediate species --- Li4NH and then into Li. The kinetic analysis of Li2NH decomposition showed that the activation energies are 533.6 kJ/mol for the first step and 754.2 kJ/mol for the second step. Furthermore, XRD demonstrated that the Li4NH, which was generated in the decomposition of Li2NH, formed a solid solution with Li2NH. In the solid solution, Li4NH possesses a similar cubic structure as Li2NH. The lattice parameter of the cubic Li4NH is 0.5033nm. The decompositions of LiNH2 and Li2NH can be promoted by chloride ion (Cl-). The introduction of Cl- into LiNH2 resulted in the generation of a new NH3 peak at low temperature of 250 °C besides the original NH3 peak at 330 °C in TPD profiles. Furthermore, Cl- can decrease the decomposition temperature of Li2NH by about 110 °C. The degradation of Li3N was systematically investigated with techniques of XRD, Fourier transform infrared (FT-IR) spectroscopy, and UV-visible spectroscopy. It was found that O2 could not affect Li3N at room temperature. However, H2O in air can cause the degradation of Li3N due to the reaction between H2O and Li3N to LiOH. The produced LiOH can further react with CO2 in air to Li2CO3 at room temperature. Furthermore, it was revealed that Alfa-Li3N is more stable in air than Beta-Li3N. The chemical stability of Li13N4Br in air has been investigated by XRD, TPD-MS, and UV-vis absorption as a function of time. The aging process finally leads to the degradation of the Li13N4Br into Li2CO3, lithium bromite (LiBrO2) and the release of gaseous NH3. The reaction order n = 2.43 is the best fitting for the Li13N4Br degradation in air reaction. Li13N4Br energy gap was calculated to be 2.61 eV.

An Investigation of Waste Glass-Based Geopolymers Supplemented with Alumina

An increased consideration of sustainability throughout society has resulted in a surge of research investigating sustainable alternatives to existing construction materials. A new binder system, called a geopolymer, is being investigated to supplement ordinary portland cement (OPC) concrete, which has come under scrutiny because of the CO2 emissions inherent in its production. Geopolymers are produced from the alkali activation of a powdered aluminosilicate source by an alkaline solution, which results in a dense three-dimensional matrix of tetrahedrally linked aluminosilicates. Geopolymers have shown great potential as a building construction material, offering similar mechanical and durability properties to OPC. Additionally, geopolymers have the added value of a considerably smaller carbon footprint than OPC. This research considered the compressive strength, microstructure and composition of geopolymers made from two types of waste glass with varying aluminum contents. Waste glass shows great potential for mainstream use in geopolymers due to its chemical and physical homogeneity as well as its high content of amorphous silica, which could eliminate the need for sodium silicate. However, the lack of aluminum is thought to negatively affect the mechanical performance and alkali stability of the geopolymer system. Mortars were designed using various combinations of glass and metakaolin or fly ash to supplement the aluminum in the system. Mortar made from the high-Al glass (12% Al2O3) reached over 10,000 psi at six months. Mortar made from the low-Al glass (<1% Al2O3) did not perform as well and remained sticky even after several weeks of curing, most likely due to the lack of Al which is believed to cause hardening in geopolymers. A moderate metakaolin replacement (25-38% by mass) was found to positively affect the compressive strength of mortars made with either type of glass. Though the microstructure of the mortar was quite indicative of mechanical performance, composition was also found to be important. The initial stoichiometry of the bulk mixture was maintained fairly closely, especially in mixtures made with fine glass. This research has shown that glass has great potential for use in geopolymers, when care is given to consider the compositional and physical properties of the glass in mixture design.