DENSITY OF AMORPHOUS CARBON BY USING DENSITY FUNCTIONAL THEORY

Amorphous carbon has been investigated for a long time. Since it has the random orientation of carbon atoms, its density depends on the position of each carbon atom. It is important to know the density of amorphous carbon to use it for modeling advance carbon materials in the future. Two methods were used to create the initial structures of amorphous carbon. One is the random placement method by randomly locating 100 carbon atoms in a cubic lattice. Another method is the liquid-quench method by using reactive force field (ReaxFF) to rapidly decrease the system of 100 carbon atoms from the melting temperature. Density functional theory (DFT) was used to refine the position of each carbon atom and the dimensions of the boundaries to minimize the ground energy of the structure. The average densities of amorphous carbon structures created by the random placement method and the liquid-quench method are 2.59 and 2.44 g/cm3, respectively. Both densities have a good agreement with previous works. In addition, the final structure of amorphous carbon generated by the liquid-quench method has lower energy.

Synthesis of novel solid materials from carbon dioxide and carbon monoxide

The increase of atmospheric CO2 has been identified as the primary cause for the observed global warming over the past century. The geological and oceanic sequestration of CO2 has issues, such as cost and leakage as well as effects on sea biota. The ideal solution should be the conversion of CO2 into useful materials. However, most processes require high energy input. Therefore, it is necessary to explore novel processes with low energy demands to convert CO2 to useful solid materials. Amorphous carbon nitride and graphone received much attention due to their unusual structures and properties as well as their potential applications. However, to date there has been no attempt to synthesize those solid materials from CO2. Lithium nitride (Li3N) and lithium imide (Li2NH) are important hydrogen storage materials. However, their optical properties and reactivity has not yet studied. This dissertation research is aimed at the synthesis of carbon nitrides and graphone from CO2 and CO via their reaction with Li3N and Li2NH. The research was focused on (1) the evaluation of Li3N and Li2NH properties, (2) thermodynamic analysis of conversion of carbon dioxide and carbon monoxide into carbon nitride and other solid materials, (3) synthesis of carbon nitride from carbon dioxide, and (4) synthesis of graphone from carbon monoxide. First, the properties of Li3N, Li2NH, and LiNH2 were investigated. The X-ray diffraction measurements revealed that heat-treatment at 500°C introduce a phase transformation of β-Li3N to α-Li3N. Furthermore, the UV-visible absorption evaluation showed that the energy gaps of α-Li3N and β-Li3N are 1.81 and 2.14 eV, respectively. The UV-visible absorption measurements also revealed that energy gaps are 3.92 eV for Li2NH and 3.93 eV for LiNH2. This thermodynamic analysis was performed to predict the reactions. It was demonstrated that the reaction between carbon dioxide and lithium nitride is thermodynamically favorable and exothermic, which can generate carbon nitride and lithium cyanamide. Furthermore, the thermodynamic calculation indicated that the reaction between carbon monoxide and lithium imide can produce graphone and lithium cyanamide along with releasing heat. Based on the above thermodynamic analysis, the experiment of CO2 and Li3N reaction and CO and Li2NH were carried out. It was found that the reaction between CO2 and Li3N is very fast and exothermic. The XRD and element analysis revealed that the products are crystal lithium cyanamide and amorphous carbon nitrides with Li2O and Li2CO3. Furthermore, TEM images showed that carbon nitrides possess layer-structure, namely, it is graphene-structured carbon nitride. It was found that the reaction between Li2NH and CO was also exothermic, which produced graphone instead of carbon nitride. The composition and structures of graphone were evaluated by XRD, element analysis, TEM observation, and Raman spectra.

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.

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.

Synthesis, structure, and stability of metal-organic frameworks

Metal-organic frameworks (MOFs) obtained much attention because of their unusual structures and properties as well as their potential applications. This dissertation research was focused on (1) the effects of synthesis conditions on the structures of MOFs, (2) the thermal stability of MOFs, (3) pressure-induced amorphization, and (4) the effect of high-valent ions on the structure of a MOF. This research demonstrated that the crystal structure of MOF-5 could be controlled by drying solvents. If the vacuum solvent is dimethylformamide (DMF), the crystal structure of MOF-5 is tetragonal. In contrast, if the DMF is displaced by CH2Cl2 before the vacuum, the obtained MOF-5 occupies a cubic structure. Furthermore, it was found that the tetragonal MOF-5 exhibited a mediate surface area (300-1000 m2/g). The surface area of tetragonal MOF-5 is also dependent on Zn(NO3)2/H2BDC (H2BDC: terephthalic acid) molar ratios used for its synthesis. The optimum ratio is 1.38, at which synthesized tetragonal MOF-5 exhibits the highest crystallinity and surface area (1297 m2/g). The thermal stability and decomposition of MOF-5 were systematically investigated. The thermal decomposition of cubic and tetragonal MOF-5s resulted in the same products: CO2, benzene, amorphous carbon, and crystal ZnO. The thermal decomposition is due to breaking carboxylic bridges between benzene rings and Zn4O clusters. Identifying structural relationships between crystalline and noncrystalline states is of fundamental interest in materials research. Currently, amorphization of solid materials at ambient temperature requires an ultra-high pressure (several GPa). However, this research demonstrated that MOF-5 and IRMOF-8 can be irreversibly amorphized at ambient temperature by employing a low compressing pressure of 3.5 MPa, which is 100 times lower than that required for amorphization of other solids. Furthermore, the pressure-induced amorphization (PIA) of MOFs is strongly dependent on the changeability of bond angles. If the geometric structure of a MOF can allow bond angles to be changed without breaking bonds, it can easily be amorphized by compression. This can explain why MOF-5 and IRMOF-8 can easily be amorphized via compression than Cu-BTC. It is generally recognized that zeolitic imidazolate frameworks (ZIFs) occupy much higher stability than other types of MOFs. The representative of ZIFs is Zn(2-methylimidazole)2 (ZIF-8) exhibiting high-decomposition temperature and high chemical resistance to various solvents. However, so far, it is still unknown whether the high stability of ZIF-8 can be challenged by ions, which is important for its modification by doping ions. In this research, we performed aqueous salt solution treatment on ZIF-8, and the results showed that anions (Cl¯ and NO3¯) in a solution exhibited no effect on the crystal structure of ZIF-8. However, the effect of cations (in a solution) on structure of ZIF-8 strongly depends on the cation valences. The univalent metal cations showed no effect on the structure of ZIF-8, whereas the bivalent or higher-valent metal cations caused the collapse of ZIF-8 crystal structure. Therefore, structure stability of ZIF-8 is considered when it is subjected to the application, in which high-valent metal cations are involved.