Matching Items (9)
Filtering by
- All Subjects: Molecular Dynamics
- Creators: Dai, Lenore
- Creators: Muhich, Christopher
- Resource Type: Text
Description
Metal-organic frameworks (MOFs) are a new set of porous materials comprised of metals or metal clusters bonded together in a coordination system by organic linkers. They are becoming popular for gas separations due to their abilities to be tailored toward specific applications. Zirconium MOFs in particular are known for their high stability under standard temperature and pressure due to the strength of the Zirconium-Oxygen coordination bond. However, the acid modulator needed to ensure long range order of the product also prevents complete linker deprotonation. This leads to a powder product that cannot easily be incorporated into continuous MOF membranes. This study therefore implemented a new bi-phase synthesis technique with a deprotonating agent to achieve intergrowth in UiO-66 membranes. Crystal intergrowth will allow for effective gas separations and future permeation testing. During experimentation, successful intergrown UiO-66 membranes were synthesized and characterized. The degree of intergrowth and crystal orientations varied with changing deprotonating agent concentration, modulator concentration, and ligand:modulator ratios. Further studies will focus on achieving the same results on porous substrates.
ContributorsClose, Emily Charlotte (Author) / Mu, Bin (Thesis director) / Shan, Bohan (Committee member) / Chemical Engineering Program (Contributor) / Barrett, The Honors College (Contributor)
Created2016-12
Description
A hybrid molecular dynamics (MD) simulation framework is developed to emulate mechanochemical reaction of mechanophores in epoxy-based nanocomposites. Two different force fields, a classical force field and a bond order based force field are hybridized to mimic the experimental processes from specimen preparation to mechanical loading test. Ultra-violet photodimerization for mechanophore synthesis and epoxy curing for thermoset polymer generation are successfully simulated by developing a numerical covalent bond generation method using the classical force field within the framework. Mechanical loading tests to activate mechanophores are also virtually conducted by deforming the volume of a simulation unit cell. The unit cell deformation leads to covalent bond elongation and subsequent bond breakage, which is captured using the bond order based force field. The outcome of the virtual loading test is used for local work analysis, which enables a quantitative study of mechanophore activation. Through the local work analysis, the onset and evolution of mechanophore activation indicating damage initiation and propagation are estimated; ultimately, the mechanophore sensitivity to external stress is evaluated. The virtual loading tests also provide accurate estimations of mechanical properties such as elastic, shear, bulk modulus, yield strain/strength, and Poisson’s ratio of the system. Experimental studies are performed in conjunction with the simulation work to validate the hybrid MD simulation framework. Less than 2% error in estimations of glass transition temperature (Tg) is observed with experimentally measured Tgs by use of differential scanning calorimetry. Virtual loading tests successfully reproduce the stress-strain curve capturing the effect of mechanophore inclusion on mechanical properties of epoxy polymer; comparable changes in Young’s modulus and yield strength are observed in experiments and simulations. Early damage signal detection, which is identified in experiments by observing increased intensity before the yield strain, is captured in simulations by showing that the critical strain representing the onset of the mechanophore activation occurs before the estimated yield strain. It is anticipated that the experimentally validated hybrid MD framework presented in this dissertation will provide a low-cost alternative to additional experiments that are required for optimizing material design parameters to improve damage sensing capability and mechanical properties.
In addition to the study of mechanochemical reaction analysis, an atomistic model of interphase in carbon fiber reinforced composites is developed. Physical entanglement between semi-crystalline carbon fiber surface and polymer matrix is captured by introducing voids in multiple graphene layers, which allow polymer matrix to intertwine with graphene layers. The hybrid MD framework is used with some modifications to estimate interphase properties that include the effect of the physical entanglement. The results are compared with existing carbon fiber surface models that assume that carbon fiber has a crystalline structure and hence are unable to capture the physical entanglement. Results indicate that the current model shows larger stress gradients across the material interphase. These large stress gradients increase the viscoplasticity and damage effects at the interphase. The results are important for improved prediction of the nonlinear response and damage evolution in composite materials.
In addition to the study of mechanochemical reaction analysis, an atomistic model of interphase in carbon fiber reinforced composites is developed. Physical entanglement between semi-crystalline carbon fiber surface and polymer matrix is captured by introducing voids in multiple graphene layers, which allow polymer matrix to intertwine with graphene layers. The hybrid MD framework is used with some modifications to estimate interphase properties that include the effect of the physical entanglement. The results are compared with existing carbon fiber surface models that assume that carbon fiber has a crystalline structure and hence are unable to capture the physical entanglement. Results indicate that the current model shows larger stress gradients across the material interphase. These large stress gradients increase the viscoplasticity and damage effects at the interphase. The results are important for improved prediction of the nonlinear response and damage evolution in composite materials.
ContributorsKoo, Bonsung (Author) / Chattopadhyay, Aditi (Thesis advisor) / Dai, Lenore (Committee member) / Jiang, Hanqing (Committee member) / Jiao, Yang (Committee member) / Rajadas, John (Committee member) / Arizona State University (Publisher)
Created2017
Description
Amine-modified solid sorbents and membrane separation are promising technologies for separation and capture of carbon dioxide (CO2) from combustion flue gas. Amine absorption processes are mature, but still have room for improvement. This work focused on the synthesis of amine-modified aerogels and metal-organic framework-5 (MOF-5) membranes for CO2 separation. A series of solid sorbents were synthesized by functionalizing amines on the surface of silica aerogels. This was done by three coating methods: physical adsorption, magnetically assisted impact coating (MAIC) and atomic layer deposition (ALD). CO2 adsorption capacity of the sorbents was measured at room temperature in a Cahn microbalance. The sorbents synthesized by physical adsorption show the largest CO2 adsorption capacity (1.43-1.63 mmol CO2/g). An additional sorbent synthesized by ALD on hydrophilic aerogels at atmospheric pressures shows an adsorption capacity of 1.23 mmol CO2/g. Studies on one amine-modified sorbent show that the powder is of agglomerate bubbling fluidization (ABF) type. The powder is difficult to fluidize and has limited bed expansion. The ultimate goal is to configure the amine-modified sorbents in a micro-jet assisted gas fluidized bed to conduct adsorption studies. MOF-5 membranes were synthesized on α-alumina supports by two methods: in situ synthesis and secondary growth synthesis. Characterization by scanning electron microscope (SEM) imaging and X-ray diffraction (XRD) show that the membranes prepared by both methods have a thickness of 14-16 μm, and a MOF-5 crystal size of 15-25 μm with no apparent orientation. Single gas permeation results indicate that the gas transport through both membranes is determined by a combination of Knudsen diffusion and viscous flow. The contribution of viscous flow indicates that the membranes have defects.
ContributorsRosa, Teresa M (Author) / Lin, Jerry (Thesis advisor) / Pfeffer, Robert (Thesis advisor) / Dai, Lenore (Committee member) / Nielsen, David (Committee member) / Arizona State University (Publisher)
Created2010
Description
In order to better understand the physical properties of polyethylene, an extremely common plastic used mostly in packaging, many scientists and engineers use olecular dynamics. To reduce the computational expense associated with traditional atomistic molecular dynamics, coarse-grained molecular dynamics is often used. Coarse-grained molecular dynamics groups multiple atoms into single beads, reducing the number of degrees of freedom in a system and eliminating smaller atoms with faster kinematics. However, even coarse-grained methods have their limitations, one of which is timestep duration, which is limited by the maximum vibrational frequency in the coarse-grained system. To study this limitation, a coarse-grained model of polyethylene was created such that every C 2 H 4 unit was replaced with a bead. Coarse-grained potentials for bond-stretching, bond-bending, and non-bonded interaction were generated using the iterative Boltzmann inversion method, which matches coarse-grained distribution functions to atomistic distribution functions. After the creation of the model, the coarse-grained potentials were rescaled by a constant so that they were less stiff, decreasing the maximum vibrational frequency of the system. It is found that by diminishing the bond-stretching potential to 6.25% of its original value, the maximum stable timestep can be increased 85% over that of the unmodified potential functions. The results of this work suggest that it may be possible to simulate lengthy processes, such as the crystallization of polyethylene, in less time with adjusted coarse-grained potentials. Additionally, the large discrepancies in the speed of bond-stretching, bond-bending, and non- bonded interaction dynamics suggest that a multi-timestep method may be worth investigating in future work.
ContributorsWiles, Christian Scott (Author) / Oswald, Jay (Thesis director) / Dai, Lenore (Committee member) / Mechanical and Aerospace Engineering Program (Contributor) / Barrett, The Honors College (Contributor)
Created2015-12
Description
Intelligent engineering designs require an accurate understanding of material behavior, since any uncertainties or gaps in knowledge must be counterbalanced with heightened factors of safety, leading to overdesign. Therefore, building better structures and pushing the performance of new components requires an improved understanding of the thermomechanical response of advanced materials under service conditions. This dissertation provides fundamental investigations of several advanced materials: thermoset polymers, a common matrix material for fiber-reinforced composites and nanocomposites; aluminum alloy 7075-T6 (AA7075-T6), a high-performance aerospace material; and ceramic matrix composites (CMCs), an advanced composite for extreme-temperature applications. To understand matrix interactions with various interfaces and nanoinclusions at their fundamental scale, the properties of thermoset polymers are studied at the atomistic scale. An improved proximity-based molecular dynamics (MD) technique for modeling the crosslinking of thermoset polymers is carefully established, enabling realistic curing simulations through its ability to dynamically and probabilistically perform complex topology transformations. The proximity-based MD curing methodology is then used to explore damage initiation and the local anisotropic evolution of mechanical properties in thermoset polymers under uniaxial tension with an emphasis on changes in stiffness through a series of tensile loading, unloading, and reloading experiments. Aluminum alloys in aerospace applications often require a fatigue life of over 109 cycles, which is well over the number of cycles that can be practically tested using conventional fatigue testing equipment. In order to study these high-life regimes, a detailed ultrasonic cycle fatigue study is presented for AA7075-T6 under fully reversed tension-compression loading. The geometric sensitivity, frequency effects, size effects, surface roughness effects, and the corresponding failure mechanisms for ultrasonic fatigue across different fatigue regimes are investigated. Finally, because CMCs are utilized in extreme environments, oxidation plays an important role in their degradation. A multiphysics modeling methodology is thus developed to address the complex coupling between oxidation, mechanical stress, and oxygen diffusion in heterogeneous carbon fiber-reinforced CMC microstructures.
ContributorsSchichtel, Jacob (Author) / Chattopadhyay, Aditi (Thesis advisor) / Dai, Lenore (Committee member) / Ghoshal, Anindya (Committee member) / Huang, Huei-Ping (Committee member) / Jiao, Yang (Committee member) / Oswald, Jay (Committee member) / Arizona State University (Publisher)
Created2022
Description
Polyurea is a highly versatile material used in coatings and armor systems to protect against extreme conditions such as ballistic impact, cavitation erosion, and blast loading. However, the relationships between microstructurally-dependent deformation mechanisms and the mechanical properties of polyurea are not yet fully understood, especially under extreme conditions. In this work, multi-scale coarse-grained models are developed to probe molecular dynamics across the wide range of time and length scales that these fundamental deformation mechanisms operate. In the first of these models, a high-resolution coarse-grained model of polyurea is developed, where similar to united-atom models, hydrogen atoms are modeled implicitly. This model was trained using a modified iterative Boltzmann inversion method that dramatically reduces the number of iterations required. Coarse-grained simulations using this model demonstrate that multiblock systems evolve to form a more interconnected hard phase, compared to the more interrupted hard phase composed of distinct ribbon-shaped domains found in diblock systems. Next, a reactive coarse-grained model is developed to simulate the influence of the difference in time scales for step-growth polymerization and phase segregation in polyurea. Analysis of the simulated cured polyurea systems reveals that more rapid reaction rates produce a smaller diameter ligaments in the gyroidal hard phase as well as increased covalent bonding connecting the hard domain ligaments as evidenced by a larger fraction of bridging segments and larger mean radius of gyration of the copolymer chains. The effect that these processing-induced structural variations have on the mechanical properties of the polymer was tested by simulating uniaxial compression, which revealed that the higher degree of hard domain connectivity leads to a 20% increase in the flow stress. A hierarchical multiresolution framework is proposed to fully link coarse-grained molecular simulations across a broader range of time scales, in which a family of coarse-grained models are developed. The models are connected using an incremental reverse–mapping scheme allowing for long time scale dynamics simulated at a highly coarsened resolution to be passed all the way to an atomistic representation.
ContributorsLiu, Minghao (Author) / Oswald, Jay (Thesis advisor) / Muhich, Christopher (Committee member) / Jiang, Hanqing (Committee member) / Peralta, Pedro (Committee member) / Jiao, Yang (Committee member) / Arizona State University (Publisher)
Created2020
Description
Functional materials can be characterized as materials that have tunable properties and are attractive solutions to the improvement and optimization of processes that require specific physiochemical characteristics. Through tailoring and altering these materials, their characteristics can be fine-tuned for specific applications. Computational modeling proves to be a crucial methodology in the design and optimization of such materials. This dissertation encompasses the utilization of molecular dynamics simulations and quantum calculations in two fields of functional materials: electrolytes and semiconductors. Molecular dynamics (MD) simulations were performed on ionic liquid-based electrolyte systems to identify molecular interactions, structural changes, and transport properties that are often reflected in experimental results. The simulations aid in the development process of the electrolyte systems in terms of concentrations of the constituents and can be invoked as a complementary or predictive tool to laboratory experiments. The theme of this study stretches further to include computational studies of the reactivity of atomic layer deposition (ALD) precursors. Selected aminosilane-based precursors were chosen to undergo density functional theory (DFT) calculations to determine surface reactivity and viability in an industrial setting. The calculations were expanded to include the testing of a semi-empirical tight binding program to predict growth per cycle and precursor reactivity with a high surface coverage model. Overall, the implementation of computational methodologies and techniques within these applications improves materials design and process efficiency while streamlining the development of new functional materials.
ContributorsGliege, Marisa Elise (Author) / Dai, Lenore (Thesis advisor) / Derecskei-Kovacs, Agnes (Thesis advisor) / Muhich, Christopher (Committee member) / Emady, Heather (Committee member) / Zhuang, Houlong (Committee member) / Arizona State University (Publisher)
Created2021
Description
Ionic liquids are salts with low melting temperatures that maintain their liquid form below 100 °C, or even at ambient temperature. Ionic liquids are conductive, electrochemically stable, non-volatile, and have a low vapor pressure, making them a class of excellent candidate materials for electrolytes in energy storage, electrodeposition, batteries, fuel cells, and supercapacitors. Due to their multiple advantages, the use of ionic liquids on Earth has been widely studied; however, further research must be done before their implementation in space. The extreme temperatures encountered during space travel and extra-terrestrial deployment have the potential to greatly affect the liquid electrolyte system. Examples of low temperature planetary bodies are the permanently shadowed sections of the moon or icy surfaces of Jupiter’s moons. Recent studies have explored the limits of glass transition temperatures for ionic liquid systems. The project is centered around the development of an ionic liquid system for a molecular electronic transducer seismometer that would be deployed on the low temperature system of Europa. For this project, molecular dynamics simulations used input intermolecular and intramolecular parameters that then simulated molecular interactions. Molecular dynamics simulations are based around the statistical mechanics of chemistry and help calculate equilibrium properties that are not easily calculated by hand. These simulations will give insight into what interactions are significant inside a ionic liquid solution. The simulations aim to create an understanding how ionic liquid electrolyte systems function at a molecular level. With this knowledge one can tune their system and its contents to adapt the systems properties to fit all environments the seismometers will experience.
ContributorsDavis, Vincent Champneys (Author) / Dai, Lenore (Thesis director) / Gliege, Marisa (Committee member) / Chemical Engineering Program (Contributor) / Barrett, The Honors College (Contributor)
Created2020-05
Description
This research seeks to answer the question if there is a singular relationship between stishovite nucleation and the atomistic structure of the preshocked amorphous SiO$_2$. To do this a stishovite manufacturing method is developed in which 1,152 samples were produced. The majority of these samples did crystallize. The method was produced through two rounds of experiments and fine-tuning with the pressure damp, temperature damp, shock pressure using an NPHug fix, and sample origin. A new random atomic insertion method was used to generate a new and different SiO$_2$ amorphous structure not before seen within the research literature. The optimal values for shock were found to be 60~GPa for randomly atom insertion samples and 55~GPa for quartz origin samples. Temperature damp appeared to have a slight effect optimizing at 0.05~ps and the pressure damp had no visible effect, testing was done with temperature damp from 0.05 to 0.5~ps and pressure damp from 0.1 to 10.0~ps. There appeared to be significant randomness in crystallization behavior. The preshocked and postnucleated samples were transformed into Gaussian fields of crystal, mass, and charge. These fields were divided and classified using a cut-off method taking the number of crystals produced in portions of each simulation and classifying each potion as nucleated or non-nucleated. Data in which some nucleation but not a critical amount was present was removed constituting 2.6\% to 20.3\% of data in all tests. A max method was also used which takes only the maximum portions of each simulation to classify as nucleating. There are three other variables tested within this work, a sample size of 18,000 or 72,728~atoms, Gaussian variance of 1 or 4~\AA, and Convolutional neural network (CNN) architecture of a garden verity or all convolution along with the portioning classification method, sample origination, and Gaussian field type. In total 64 tests were performed to try every combination of variable. No significant classifications were made by the CNNs to nucleation or non-nucleation portions. The results clearly confirmed that the data was not abstracting to atomistic structure and was random by all classifications of the CNNs. The all convolution CNN testing did show smoother outcomes in training with less fluctuations. 59\% of all validation accuracy was held at 0.5 for a random state and 84\% was within $\pm0.02$ of 0.5. It is conclusive that prenucleation structure is not the sole predictor of nucleation behavior. It is not conclusive if prenucleation structure is a partial or non-factor within nucleation of stishovite from amorphous SiO$_2$.
ContributorsChristen, Jonathan Scorr (Author) / Oswald, Jay (Thesis advisor) / Muhich, Christopher (Committee member) / Zhuang, Houlong (Committee member) / Arizona State University (Publisher)
Created2021