Enabling Solar Thermal Propulsion with Computational Materials Design
No abstract available
Haskins, Justin B.↗
Engineering topics
Bauschlicher, Charles W.
Publications and source records attributed to Bauschlicher, Charles W..
The Convergence of the Coupled Cluster Approach for MgO
The X(exp 1)Sigma(+) state of MgO shows significant differences between the CCSD(T) and CCSDT results and even some significant differences between the CCSDT and CCSDTQ results. The CCSDT(Q) results for this state do not agree well with the CCSDTQ results, as has been found for other systems. Unlike the X state, the state, which is well described by a single configuration, shows much less dependence on level of correlation treatment. Despite the slow convergence of the valence treatment of the X state, the effect of core correlation converges rapidly with level of correlation treatment.
Bauschlicher, Charles W.↗
Ionic Liquids at Electrified Interfaces: From Double Layers to Decomposition
Ionic liquids are versatile electrolytes whose properties at electrified interfaces have the potential to enable technologies such as supercapacitors and Li-metal battery anodes. At electrified carbon surfaces, ionic liquids form an electric double layer that stores energy and provides the foundation for supercapacitors. At electrified lithium surfaces, ionic liquids decompose to form a solid electrolyte interphase that has the potential to stabilize Li-metal anodes in rechargable batteries. The behavior of two ionic liquids of technological importance, [pyr14][TFSI] and [EMIM][BF4], are examined at these electrified interfaces through molecular dynamics and ab initio techniques.
Haskins, Justin B.↗
Phenolic Polymer Interactions with Water and Ethylene Glycol Solvents
Interactions between pre-cured phenolic polymer chains and a solvent have a significant impact on the structure and properties of the final post-cured phenolic resin. Developing an understanding of the nature of these interactions is important and will aid in the selection of the proper solvent that will lead to the desired final product. Here, we investigate the role of the phenolic chain structure and the solvent type on the overall solvation performance of the system through ab initio techniques and molecular dynamics simulations. Two types of solvents are considered: ethylene glycol (EGL) and H2O. Three phenolic chain structures are considered, including two novolac-type chains with either an ortho-ortho (OON) or an ortho-para (OPN) backbone network and a resole-type (RES) chain with an ortho-ortho network. Each system is characterized through a structural analysis of the solvation shell and the hydrogen-bonding environment as well as through a quantification of the solvation free energy along with partitioned interaction energies between specific molecular species. The combination of simulations and the analyses indicate that EGL provides a higher solvation free energy than H2O due to more energetically favorable hydrophilic interactions as well as favorable hydrophobic interactions between CH element groups. In addition, the phenolic chain structure significantly affects the solvation performance, with OON having limited intermolecular hydrogen-bond formations, while OPN and RES interact more favorably with the solvent molecules. The results suggest that a resole-type phenolic chain with an ortho-para network should have the best solvation performance in EGL, H2O, and other similar solvents.
Haskins, Justin B.↗
Molecular Dynamics Simulations of Liquid and Polymer Electrolytes for Energy Storage Devices
Advancing beyond current lithium-ion technology is necessary in order to enable energy storage devices for electric airplanes. Electrolyte stability is a key limiting factor, yet the design of improved electrolytes remains a formidable challenge. Molecular dynamics (MD) simulations are a powerful tool for studying electrolytes, since they can be used to evaluate structural, thermodynamic, and transport properties, and can provide molecular-level detail often inaccessible to experimental techniques. Our computational materials groups at the NASA Ames Research Center has developed models and methods to accurately simulate both liquid and polymer electrolytes.We report the results from atomistic MD simulations of several electrolyte materials, with lithium salts dissolved in ionic liquids, dimethoxyethane (DME), and polyethylene oxide (PEO). For improved accuracy, we employ polarizable models, where each atom is given an environment-dependent atomic dipole. The simulations accurately predict bulk transport properties, including viscosity, diffusion, and ionic conductivity, in quantitative agreement with available experimental data. Moreover, the simulations provide important insights into the solvation structure of the lithium ions.We also report the results from coarse-grained MD simulations of polyanion electrolytes. In order to more efficiently capture the longer length- and time-scales of these systems, we employ a generic bead-spring model. These simulations provide important insight into how the polymer chain architecture and ionic interaction strengths affect the ionic aggregation behavior and cation dynamics. Despite the simplicity of the model, the simulations yield qualitative agreement with experimental data for similar systems.
Abbott, Lauren J.↗
Yield and Failure Behavior Investigated for Cross-linked Phenolic Resins Using Molecular Dynamics
Molecular dynamics simulations were conducted to fundamentally evaluate the yield and failure behavior of cross-linked phenolic resins at temperatures below the glass transition. Yield stress was investigated at various temperatures, strain rates, and degrees of cross-linking. The onset of non-linear behavior in the cross-linked phenolic structures was caused by localized irreversible molecular rearrangements through the rotation of methylene linkers followed by the formation or annihilation of neighboring hydrogen bonds. The yield stress results, with respect to temperature and strain rate, could be fit by existing models used to describe yield behavior of amorphous glasses. The degree of cross-linking only indirectly influences the maximum yield stress through its influence on the glass transition temperature (Tg), however there is a strong relationship between the degree of cross-linking and the failure mechanism. Low cross-linked samples were able to separate through void formation, whereas the highly cross-linked structures exhibited bond scission.
Stress-strain curves↗
Investigation of Structure and Transport in Li-Doped Ionic Liquid Electrolytes: [pyr14][TFSI], [pyr13][FSI] and [EMIM][BF4]
Ionic liquid electrolytes have been proposed as a means of improving the safety and cycling behavior of advanced lithium batteries; however, the properties of these electrolytes under high lithium doping are poorly understood. Here, we employ both polarizable molecular dynamics simulation and experiment to investigate the structure, thermodynamics and transport of three potential electrolytes, N-methyl-Nbutylpyrrolidinium bis(trifluoromethylsufonyl)imide ([pyr14][TFSI]), N- methyl-Npropylpyrrolidinium bis(fluorosufonyl)imide ([pyr13][FSI]), and 1-ethyl-3-- methylimidazolium boron tetrafluoride ([EMIM][BF4]), as a function of Li-salt concentration and temperature. Structurally, Li(+) is shown to be solvated by three anion neighbors in [pyr14][TFSI] and four anion neighbors in both [pyr13][FSI] and [EMIM][BF4], and at all levels of x(sub Li) we find the presence of lithium aggregates. Furthermore, the computed density, diffusion, viscosity, and ionic conductivity show excellent agreement with experimental data. While the diffusion and viscosity exhibit a systematic decrease and increase, respectively, with increasing x(sub Li), the contribution of Li(+) to ionic conductivity increases until reaching a saturation doping level of x(sub Li) is approximately 0.10. Comparatively, the Li(+) conductivity of [pyr14][TFSI] is an order of magnitude lower than that of the other liquids, which range between 0.1 - 0.3 mS/cm. The differences in Li(+) transport are reflected in the residence times of Li(+) with the anions, which are revealed to be much larger for [pyr14][TFSI] (up to 100 ns at the highest doping levels) than in either [EMIM][BF4] or [pyr13][FSI]. Finally, we comment on the relative kinetics of Li(+) transport in each liquid and we present strong evidence for transport through anion exchange (hopping) as opposed to the net motion of Li(+) with its solvation shell (vehicular).
lithium batteries↗
Computational and Experimental Investigation of Li-doped Ionic Liquid Electrolytes: [pyr14][tfsi], [pyr13][fsi], and [EMIM][BF4]
We employ molecular dynamics (MD) simulation and experiment to investigate the structure, thermodynamics, and transport of N-methyl-N-butylpyrrolidinium bis(trifluoromethylsufonyl)imide ([pyr14][TFSI]), N -methyl-N-propylpyrrolidinium bis(fluorosufonyl)imide ([pyr13][FSI]), and 1-ethyl-3-methylimidazolium boron tetrafluoride ([EMIM][BF4]), as a function of Li-salt mole fraction (0.05 xLi+ 0.33) and temperature (298 K T 393 K). Structurally, Li+ is shown to be solvated by three anion neigh- bors in [pyr14][TFSI] and four anion neighbors in both [pyr13][FSI] and [EMIM][BF4], and at all levels of xLi+ we find the presence of lithium aggregates. Pulsed field gradient spin-echo NMR measurements of diffusion and electrochemical impedance spectroscopy measurements of ionic conductivity are made for the neat ionic liquids as well as 0.5 molal solutions of Li-salt in the ionic liquids. Bulk ionic liquid properties (density, diffusion, viscosity, and ionic conductivity) are obtained with MD and show excellent agreement with experiment. While the diffusion exhibits a systematic decrease with increasing xLi+, the contribution of Li+ to ionic conductivity increases until reach- ing a saturation doping level of xLi+ 0.10. Comparatively, the Li+ conductivity of [pyr14][TFSI] is an order of magnitude lower than that of the other liquids, which range between 0.1-0.3 mScm. Our transport results also demonstrate the necessity of long MD simulation runs ( 200 ns) required to converge transport properties at room T. The differences in Li+ transport are reflected in the residence times of Li+ with the anions (Li), which are revealed to be much larger for [pyr14][TFSI] (up to 100 ns at the highest doping levels) than in either [EMIM][BF4] or [pyr13][FSI]. Finally, to comment on the relative kinetics of Li+ transport in each liquid, we find that while the net motion of Li+ with its solvation shell (vehicular) significantly contributes to net diffusion in all liquids, the importance of transport through anion exchange (hopping) increases at high xLi+ and in liquids with large anions.
electrolytes↗
Multiscale Modeling of UHTC: Thermal Conductivity
We are developing a multiscale framework in computational modeling for the ultra high temperature ceramics (UHTC) ZrB2 and HfB2. These materials are characterized by high melting point, good strength, and reasonable oxidation resistance. They are candidate materials for a number of applications in extreme environments including sharp leading edges of hypersonic aircraft. In particular, we used a combination of ab initio methods, atomistic simulations and continuum computations to obtain insights into fundamental properties of these materials. Ab initio methods were used to compute basic structural, mechanical and thermal properties. From these results, a database was constructed to fit a Tersoff style interatomic potential suitable for atomistic simulations. These potentials were used to evaluate the lattice thermal conductivity of single crystals and the thermal resistance of simple grain boundaries. Finite element method (FEM) computations using atomistic results as inputs were performed with meshes constructed on SEM images thereby modeling the realistic microstructure. These continuum computations showed the reduction in thermal conductivity due to the grain boundary network.
Lawson, John W.↗
Lattice Thermal Conductivity from Atomistic Simulations: ZrB2 and HfB2
Ultra high temperature ceramics (UHTC) including ZrB2 and HfB2 have a number of properties that make them attractive for applications in extreme environments. One such property is their high thermal conductivity. Computational modeling of these materials will facilitate understanding of fundamental mechanisms, elucidate structure-property relationships, and ultimately accelerate the materials design cycle. Progress in computational modeling of UHTCs however has been limited in part due to the absence of suitable interatomic potentials. Recently, we developed Tersoff style parameterizations of such potentials for both ZrB2 and HfB2 appropriate for atomistic simulations. As an application, Green-Kubo molecular dynamics simulations were performed to evaluate the lattice thermal conductivity for single crystals of ZrB2 and HfB2. The atomic mass difference in these binary compounds leads to oscillations in the time correlation function of the heat current, in contrast to the more typical monotonic decay seen in monoatomic materials such as Silicon, for example. Results at room temperature and at elevated temperatures will be reported.
Lawson, John W.↗
Lattice Thermal Conductivity of Ultra High Temperature Ceramics (UHTC) ZrB2 and HfB2 from Atomistic Simulations
Ultra high temperature ceramics (UHTC) including ZrB2 and HfB2 have a number of properties that make them attractive for applications in extreme environments. One such property is their high thermal conductivity. Computational modeling of these materials will facilitate understanding of fundamental mechanisms, elucidate structure-property relationships, and ultimately accelerate the materials design cycle. Progress in computational modeling of UHTCs however has been limited in part due to the absence of suitable interatomic potentials. Recently, we developed Tersoff style parameterizations of such potentials for both ZrB2 and HfB2 appropriate for atomistic simulations. As an application, Green-Kubo molecular dynamics simulations were performed to evaluate the lattice thermal conductivity for single crystals of ZrB2 and HfB2. The atomic mass difference in these binary compounds leads to oscillations in the time correlation function of the heat current, in contrast to the more typical monotonic decay seen in monoatomic materials such as Silicon, for example. Results at room temperature and at elevated temperatures will be reported.
Lawson, John W.↗
Multiscale Modeling of Ultra High Temperature Ceramics (UHTC) ZrB2 and HfB2: Application to Lattice Thermal Conductivity
We are developing a multiscale framework in computational modeling for the ultra high temperature ceramics (UHTC) ZrB2 and HfB2. These materials are characterized by high melting point, good strength, and reasonable oxidation resistance. They are candidate materials for a number of applications in extreme environments including sharp leading edges of hypersonic aircraft. In particular, we used a combination of ab initio methods, atomistic simulations and continuum computations to obtain insights into fundamental properties of these materials. Ab initio methods were used to compute basic structural, mechanical and thermal properties. From these results, a database was constructed to fit a Tersoff style interatomic potential suitable for atomistic simulations. These potentials were used to evaluate the lattice thermal conductivity of single crystals and the thermal resistance of simple grain boundaries. Finite element method (FEM) computations using atomistic results as inputs were performed with meshes constructed on SEM images thereby modeling the realistic microstructure. These continuum computations showed the reduction in thermal conductivity due to the grain boundary network.
Lawson, John W.↗
Ab initio Computations of the Electronic, Mechanical, and Thermal Properties of Ultra High Temperature Ceramics (UHTC) ZrB2 and HfB2
Refractory materials such as metallic borides, often considered as ultra high temperature ceramics (UHTC), are characterized by high melting point, high hardness, and good chemical inertness. These materials have many applications which require high temperature materials that can operate with no or limited oxidation. Ab initio, first principles methods are the most accurate modeling approaches available and represent a parameter free description of the material based on the quantum mechanical equations. Using these methods, many of the intrinsic properties of these material can be obtained. We performed ab initio calculations based on density functional theory for the UHTC materials ZrB2 and HfB2. Computational results are presented for structural information (lattice constants, bond lengths, etc), electronic structure (bonding motifs, densities of states, band structure, etc), thermal quantities (phonon spectra, phonon densities of states, specific heat), as well as information about point defects such as vacancy and antisite formation energies.
Lawson, John W.↗
Lattice Thermal Conductivity of Ultra High Temperature Ceramics (UHTC) ZrB2 and HfB2 from Atomistic Simulations
Ultra high temperature ceramics (UHTC) including ZrB2 and HfB2 are candidate materials for applications in extreme environments because of their high melting point, good mechanical properties and reasonable oxidation resistance. Unlike many ceramics, these materials have high thermal conductivity which can be advantageous, for example, to reduce thermal shock. Recently, we developed Tersoff style interatomic potentials for both ZrB2 and HfB2 appropriate for atomistic simulations. As an application, Green-Kubo molecular dynamics simulations were performed to evaluate the lattice thermal conductivity for single crystals of ZrB2 and HfB2. The atomic mass difference in these binary compounds leads to oscillations in the time correlation function of the heat current. Results at room temperature and at elevated temperatures will be reported.
Lawson, JOhn W.↗
Variations in the Peak Position of the 6.2 micron Interstellar Emission Feature: A Tracer of N in the Interstellar Polycyclic Aromatic Hydrocarbon Population
This paper presents the results of an investigation of the molecular characteristics that underlie the observed peak position and profile of the nominal 6.2 micron interstellar emission band generally attributed to the CC stretching vibrations of polycyclic aromatic hydrocarbons (PAHs). It begins with a summary of recent experimental and theoretical studies ofthe spectroscopic properties of large (>30 carbon atoms) PAH cations as they relate to this aspect of the astrophysical problem. It then continues with an examination of the spectroscopic properties of a number of PAH variants within the context of the interstellar 6.2 micron emission, beginning with a class of compounds known as polycyclic aromatic nitrogen heterocycles (PANHs; PAHs with one or more nitrogen atoms substituted into their carbon skeleton). In this regard, we summarize the results of recent relevant experimental studies involving a limited set of small PANHs and their cations and then report the results of a comprehensive computational study that extends that work to larger PANH cations including many nitrogen-substituted variants of coronene(+) (C24H12(+)), ovalene(+) (C32H14(+)), circumcoronene(+) (C54H18(+)), and circum-circumcoronene(+) (C96H24(+)). Finally, we report the results of more focused computational studies of selected representatives from a number of other classes of PAH variants that share one or more of the key attributes of the PANH species studied. These alternative classes of PAH variants include (1) oxygen- and silicon-substituted PAH cations; (2) PAH-metal ion complexes (metallocenes) involving the cosmically abundant elements magnesium and iron; and (3) large, asymmetric PAH cations. Overall, the studies reported here demonstrate that increasing PAH size alone is insuEcient to account for the position of the shortest wavelength interstellar 6.2 micron emission bands, as had been suggested by earlier studies. On the other hand, this work reveals that substitution of one or more nitrogen atoms within the interior of the carbon skeleton of a PAH cation induces a significant blueshift in the position of the dominant CC stretching feature of these compounds that is suf6cient to account for the position of the interstellar bands. Subsequent studies of the effects of substitution by other heteroatoms (O and Si), metal ion complexation (Fe(+), Mg(+), and Mg(2+)), and molecular symmetry variation-all of which fail to reproduce the blueshift observed in the PANH cations-indicate that N appears to be unique in its ability to accommodate the position of the interstellar 6.2 micron bands while simultaneously satisfying the other constraints of the astrophysical problem. This result implies that the peak position of the interstellar feature near 6.2 micron traces the degree of nitrogen substitution in the population, that most of the PAHs responsible for the interstellar IR emission features incorporate nitrogen within their aromatic networks, and that a lower limit of 1%-2% of the cosmic nitrogen is sequestered within the interstellar PAH population. Finally, in view of the ubiquity and abundance of interstellar PAHs and the permanent dipoles and distinctive electronic structures of these nitrogen-substituted variants, this work impacts a wide range of observational phenomena outside of the infrared region of the spectrum including the forest of unidentified molecular rotational features and the anomalous Galactic foreground emission in the microwave, and the diffuse interstellar bands (DIBs) and other structure in the interstellar extinction curve in the ulhviolet/visible. These astrophysical ramifications are discussed, and the dipole moments and rotational constants are tabulated to facilitate further investigations of the astrophysical role of nitrogen-substituted aromatic compounds.
Hudgins, Douglas M.↗
Mid-Infrared Spectroscopy of Polycyclic Aromatic Nitrogen Heterocycles (PANHS) and their Ions
In recent years, polycyclic aromatic nitrogen heterocycles (PANHs) have attracted a good deal of attention because of their potent carcinogenic and mutagenic properties, and their prevalence in our environment. Such species also play a prominent role in the chemistry of life up to and including the very nucleobases from which our DNA is constructed. Surprisingly, these compounds may even be common outside of our terrestrial environment. To wit, it is now widely accepted that polycyclic aromatic materials are abundant in space and represent a major reservoir of organic carbon in the interstellar medium and developing planetary systems. Given that nitrogen is the fourth most abundant chemically reactive element in space (surpassed only by hydrogen, carbon, and oxygen), it is entirely reasonable to suspect that PANHs may represent an important component of that organic reservoir. Motivated by their intrinsic merit and with special attention toward evaluating their exobiological significance, we have initiated a program to study the spectroscopic and chemical properties of P A " s under conditions relevant to extraterrestrial environments. Here we present the first results of that program-infrared spectroscopic measurements on a series of PANH"s in neutral and cationic forms, isolated in inert matrices at cryogenic temperatures.temperatures. The species studied include: 1 -, and 2-azabenz[a]anthracene, 1-, 2-, and 4- azachrysene, dibenz[a,h]acridine, and dibenz[a,J)acridine. The experimental measurements are also compared with theoretical spectra calculated using density functional theory. General spectroscopic trends observed in this series of compounds are discussed and the implications of these results for Astrophysics and Exobiology are considered.
Mattioda, Andrew L.↗
Infrared Spectra and Density Functional Calculations on Binary Unsaturated Transition Metal Carbonyl Cations, Neutrals and Anions
The status of theory and experiment for the determination of the CO vibrational frequencies in binary unsaturated transition metal carbonyl cations, neutrals and anions is reviewed.
Zhou, Mingfei↗
Infrared Spectra of Polycyclic Aromatic Hydrocarbons: Methyl Substitution and Loss of H
The B3LYP approach, in conjunction with the 4-31G basis set, is used to compute the harmonic frequencies of 1- and 2-methylnaphthalene, 1-, 2-, and 9-methylanthracene, and their cations. The IR spectra of the methyl substituted species are very similar to the parent spectra, except for the addition of the methyl C-H stretch at lower frequency than the aromatic C-H stretch. The loss of a single hydrogen from naphthalene, anthracene, and their cations is shown to have a very small effect on the IR spectra. Loss of a methyl hydrogen from 1- or 2-methylnaphthalene, or their cations, is shown to shift the side group C-H frequencies from below aromatic hydrogen stretching frequencies to above them. The loss of IT from 2-methylenenaphthalene shows only a small shift in the side group C-H stretching frequency.
Bauschlicher, Charles W.↗