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.

Charge induced enhancement of adsorption for hydrogen storage materials

The rising concerns about environmental pollution and global warming have facilitated research interest in hydrogen energy as an alternative energy source. To apply hydrogen for transportations, several issues have to be solved, within which hydrogen storage is the most critical problem. Lots of materials and devices have been developed; however, none is able to meet the DOE storage target. The primary issue for hydrogen physisorption is a weak interaction between hydrogen and the surface of solid materials, resulting negligible adsorption at room temperature. To solve this issue, there is a need to increase the interaction between the hydrogen molecules and adsorbent surface. In this study, intrinsic electric dipole is investigated to enhance the adsorption energy. The results from the computer simulation of single ionic compounds with hydrogen molecules to form hydrogen clusters showed that electrical charge of substances plays an important role in generation of attractive interaction with hydrogen molecules. In order to further examine the effects of static interaction on hydrogen adsorption, activated carbon with a large surface area was impregnated with various ionic salts including LiCl, NaCl, KCl, KBr, and NiCl and their performance for hydrogen storage was evaluated by using a volumetric method. Corresponding computer simulations have been carried out by using DFT (Density Functional Theory) method combined with point charge arrays. Both experimental and computational results prove that the adsorption capacity of hydrogen and its interaction with the solid materials increased with electrical dipole moment. Besides the intrinsic dipole, an externally applied electric field could be another means to enhance hydrogen adsorption. Hydrogen adsorption under an applied electric field was examined by using porous nickel foil as electrodes. Electrical signals showed that adsorption capacity increased with the increasing of gas pressure and external electric voltage. Direct measurement of the amount of hydrogen adsorption was also carried out with porous nickel oxides and magnesium oxides using the piezoelectric material PMN-PT as the charge supplier due to the pressure. The adsorption enhancement from the PMN-PT generated charges is obvious at hydrogen pressure between 0 and 60 bars, where the hydrogen uptake is increased at about 35% for nickel oxide and 25% for magnesium oxide. Computer simulation reveals that under the external electric field, the electron cloud of hydrogen molecules is pulled over to the adsorbent site and can overlap with the adsorbent electrons, which in turn enhances the adsorption energy Experiments were also carried out to examine the effects of hydrogen spillover with charge induced enhancement. The results show that the overall storage capacity in nickel oxide increased remarkably by a factor of 4.

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.

Functional ZnO nanostructures for electronic and energy applications

ZnO has proven to be a multifunctional material with important nanotechnological applications. ZnO nanostructures can be grown in various forms such as nanowires, nanorods, nanobelts, nanocombs etc. In this work, ZnO nanostructures are grown in a double quartz tube configuration thermal Chemical Vapor Deposition (CVD) system. We focus on functionalized ZnO Nanostructures by controlling their structures and tuning their properties for various applications. The following topics have been investigated: 1. We have fabricated various ZnO nanostructures using a thermal CVD technique. The growth parameters were optimized and studied for different nanostructures. 2. We have studied the application of ZnO nanowires (ZnONWs) for field effect transistors (FETs). Unintentional n-type conductivity was observed in our FETs based on as-grown ZnO NWs. We have then shown for the first time that controlled incorporation of hydrogen into ZnO NWs can introduce p-type characters to the nanowires. We further found that the n-type behaviors remained, leading to the ambipolar behaviors of hydrogen incorporated ZnO NWs. Importantly, the detected p- and n- type behaviors are stable for longer than two years when devices were kept in ambient conditions. All these can be explained by an ab initio model of Zn vacancy-Hydrogen complexes, which can serve as the donor, acceptors, or green photoluminescence quencher, depend on the number of hydrogen atoms involved. 3. Next ZnONWs were tested for electron field emission. We focus on reducing the threshold field (Eth) of field emission from non-aligned ZnO NWs. As encouraged by our results on enhancing the conductivity of ZnO NWs by hydrogen annealing described in Chapter 3, we have studied the effect of hydrogen annealing for improving field emission behavior of our ZnO NWs. We found that optimally annealed ZnO NWs offered much lower threshold electric field and improved emission stability. We also studied field emission from ZnO NWs at moderate vacuum levels. We found that there exists a minimum Eth as we scale the threshold field with pressure. This behavior is explained by referring to Paschen’s law. 4. We have studied the application of ZnO nanostructures for solar energy harvesting. First, as-grown and (CdSe) ZnS QDs decorated ZnO NBs and ZnONWs were tested for photocurrent generation. All these nanostructures offered fast response time to solar radiation. The decoration of QDs decreases the stable current level produced by ZnONWs but increases that generated by NBs. It is possible that NBs offer more stable surfaces for the attachment of QDs. In addition, our results suggests that performance degradation of solar cells made by growing ZnO NWs on ITO is due to the increase in resistance of ITO after the high temperature growth process. Hydrogen annealing also improve the efficiency of the solar cells by decreasing the resistance of ITO. Due to the issues on ITO, we use Ni foil as the growth substrates. Performance of solar cells made by growing ZnO NWs on Ni foils degraded after Hydrogen annealing at both low (300 °C) and high (600 °C) temperatures since annealing passivates native defects in ZnONWs and thus reduce the absorption of visible spectra from our solar simulator. Decoration of QDs improves the efficiency of such solar cells by increasing absorption of light in the visible region. Using a better electrolyte than phosphate buffer solution (PBS) such as KI also improves the solar cell efficiency. 5. Finally, we have attempted p-type doping of ZnO NWs using various growth precursors including phosphorus pentoxide, sodium fluoride, and zinc fluoride. We have also attempted to create p-type carriers via introducing interstitial fluorine by annealing ZnO nanostructures in diluted fluorine gas. In brief, we are unable to reproduce the growth of reported p-type ZnO nanostructures. However; we have identified the window of temperature and duration of post-growth annealing of ZnO NWs in dilute fluorine gas which leads to suppression of native defects. This is the first experimental effort on post-growth annealing of ZnO NWs in dilute fluorine gas although this has been suggested by a recent theory for creating p-type semiconductors. In our experiments the defect band peak due to native defects is found to decrease by annealing at 300 °C for 10 – 30 minutes. One of the major future works will be to determine the type of charge carriers in our annealed ZnONWs.

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.

A THEORETICAL STUDY OF INTERACTION OF NANOPARTICLES WITH BIOMOLECULE

Many types of materials at nanoscale are currently being used in everyday life. The production and use of such products based on engineered nanomaterials have raised concerns of the possible risks and hazards associated with these nanomaterials. In order to evaluate and gain a better understanding of their effects on living organisms, we have performed first-principles quantum mechanical calculations and molecular dynamics simulations. Specifically, we will investigate the interaction of nanomaterials including semiconducting quantum dots and metallic nanoparticles with various biological molecules, such as dopamine, DNA nucleobases and lipid membranes. Firstly, interactions of semiconducting CdSe/CdS quantum dots (QDs) with the dopamine and the DNA nucleobase molecules are investigated using similar quantum mechanical approach to the one used for the metallic nanoparticles. A variety of interaction sites are explored. Our results show that small-sized Cd4Se4 and Cd4S4 QDs interact strongly with the DNA nucleobase if a DNA nucleobase has the amide or hydroxyl chemical group. These results indicate that these QDs are suitable for detecting subcellular structures, as also reported by experiments. The next two chapters describe a preparation required for the simulation of nanoparticles interacting with membranes leading to accurate structure models for the membranes. We develop a method for the molecular crystalline structure prediction of 1,2-Dimyristoyl-sn-glycero-3-phosphorylcholine (DMPC), 1,2-Dimyristoyl-sn-glycero-3-phosphorylethanolamine (DMPE) and cyclic di-amino acid peptide using first-principles methods. Since an accurate determination of the structure of an organic crystal is usually an extremely difficult task due to availability of the large number of its conformers, we propose a new computational scheme by applying knowledge of symmetry, structural chemistry and chemical bonding to reduce the sampling size of the conformation space. The interaction of metal nanoparticles with cell membranes is finally carried out by molecular dynamics simulations, and the results are reported in the last chapter. A new force field is developed which accurately describes the interaction forces between the clusters representing small-sized metal nanoparticles and the lipid bilayer molecules. The permeation of nanoparticles into the cell membrane is analyzed together with the RMSD values of the membrane modeled by a lipid bilayer. The simulation results suggest that the AgNPs could cause the same amount of deformation as the AuNPs for the dysfunction of the membrane.