SEARCH · Search NASA
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Preparation of reactive cyclic compounds using cryochemical techniques, and low temperature mass spectrometric stability and molecular energetics study
Sensitive millimeter-wave C-12O line observations and near-infrared hydrogen recombination line observations of AFGL 961 are reported and analyzed. The observations show that the high-velocity flow has a maximum extent of at least six arcmin or 1.9 pc at the distance of AFGL 961. The flow is anisotropic and the emission profiles are asymmetric in velocity. The mass of the gas involved in the flow is determined to be approximately 19 solar masses, and the kinetic energy of this gas is estimated to be about 8 x 10 to the 46th ergs. These observations are interpreted as evidence that an energetic bipolar outflow of molecular gas is occurring near AFGL 961. Comparison of infrared recombination line data, radio continuum, and CO observations shows that a strong, ionized, stellar wind would not be a sufficient force to drive the outflow. The effects of the outflow on the surrounding molecular cloud was also investigated.
Complete tensors of the second- and third-order elastic constants of the organic molecular crystals acetaminophen, pentaerythritol tetranitrate (PETN), cyclotrimethylene trinitramine (RDX), cyclotetramethylene tetranitramine (HMX), 1,1-diamino-2,2-dinitroethylene (FOX-7), hexanitrohexaazaisowurtzitane (CL-20), and erythritol tetranitrate (ETN) have been calculated using dispersion-corrected density functional theory. The sets of second- and third-order elastic constants are expected to provide a more accurate and reliable description of the behavior of these materials under nonhydrostatic loads than pressure- and volume-dependent second-order elastic constants. The tensors of second-order constants have been compared with experimental data and/or other calculations when possible, and with the exception of results for CL-20 from Brillouin scattering experiments, we find good agreement. The calculated third-order elastic constants of PETN are in very good agreement with the subset of third-order constants derived from experimental wave speed measurements. The elastic anisotropies of the crystals have been estimated using the universal elastic anisotropy index, which shows that the crystals fall into three groups with low anisotropy (PETN, RDX, and CL-20), moderate anisotropy (acetaminophen, HMX, and ETN), and high elastic anisotropy (FOX-7).
When examining dynamics occurring at nonzero temperatures, both energy and entropy must be taken into account to describe activated barrier crossing events. Furthermore, good reaction coordinates need to be constructed to describe different metastable states and the transition mechanisms between them. Here we use a physics-based machine learning method called state predictive information bottleneck (SPIB) to find nonlinear reaction coordinates for three systems of varying complexity. SPIB is able to correctly predict an entropic bottleneck for an analytical flat-energy double-well system and identify the entropy- and energy-dominated pathways for an analytical four-well system. Lastly, for a simulation of benzoic acid permeation through a lipid bilayer, SPIB is able to discover the entropic and energetic barriers to the permeation process. Given these results, we thus establish that SPIB is a reasonable and robust method for finding the important entropy, energy, and enthalpy barriers in physical systems, which can then be used to enhance the understanding and sampling of different activated mechanisms.
Observations of the S255 molecular cloud in the radio lines of CO, (C-13)O, and H2O as well as in the IR region from 12 to 20 microns are presented. These results show that an extended molecular cloud is associated with the optically visible H II regions S255 - S257, both the kinetic temperature and molecular density are enhanced in an area centered on a compact near-IR source, this source is also coincident with an OH and an H2O maser (within the errors), and the IR source exhibits the silicon absorption feature. First-order techniques which allow physical properties of molecular clouds to be estimated and the energetics of the gas and dust in the clouds to be analyzed are developed and applied to S255. It is concluded that: (1) the gas cooling rate is much less than the dust cooling rate; (2) the primary energy flow is in the IR through dust emission over a range of temperatures; (3) the near-IR source is apparently responsible for the local peak in gas kinetic temperature; and (4) the exciting stars of nearby H II regions appear to be the primary heat sources for the overall molecular cloud.
The S88 molecular cloud has been observed in several molecular lines and at infrared wavelengths from 1 to 100 microns. The CO emission has a single, sharp peak which is near the H-alpha emission region S88 B and centered on a compact H II region observed in the radio continuum. The infrared observations indicate that the principal luminosity source is located near the radio continuum peak and is hidden behind substantial extinction, presumably from the molecular cloud. The molecular cloud has a size of about 6.5 pc and a mass of about 5,000 solar masses at an assumed distance of 2 kpc. Regions of very high density (more than 10,000/cu cm) do not seem to be present. Analysis of the energetics leads to a predicted dust cooling rate in good agreement with the infrared observations, which indicate L = 180,000 solar luminosities. A more detailed examination of the energetics indicates that densities may be insufficient for collisions of molecules with warm dust grains to heat the gas to the observed kinetic temperature.
The molecular cloud associated with the S235 H II region has been studied by means of molecular lines and near-infrared observations. The cloud consists of two components, one of which partly surrounds the S235 H II region. The other component contains a dense, hot region of active star formation, marked by self-reversed CO profiles, compact H II regions, masers, and infrared sources. From the molecular line data, the size and mass of the two components are estimated to be 6-8 pc and 3000-4000 solar masses. More detailed studies near the region of active star formation yield estimates of density (approximately 200,000-500,000/cu cm) and the abundances of H2CO, HCO(+), HCN, and (C-13)O. Analysis of the energetics suggests that the cloud is heated by the exciting star of S235 and by the exciting stars of the compact H II regions. Assuming that the gas is heated by collisions with warm dust grains, the far-infrared luminosity has been predicted. The observations of far-infrared emission are in reasonable agreement with predictions.
The molecular clouds of S37, which coincides with a region of visual obscuration, is observed at several molecular transition frequencies and at IR wavelengths in the 2-125 micron range. It is found that the temperature and the column density peak near the reflection nebulae vdB 118 and 119, with line splitting and possible self-absorption being observed near the latter. Although many 2-micron sources are observed in the direction of the cloud, no bright, 10-micron sources are seen, and far-IR emission is found over an extensive region in the cloud. An energetics analysis leads to the prediction of a dust cooling rate consistent with the far-IR observations, which indicate not less than 15,000-35,000 solar luminosities, and more detailed energetics examinations show that the conditions necessary for collision heating of the gas by warm dust exist.
Techniques of analysis which allow estimates of kinetic temperature, gas density, gas cooling rate, and dust cooling rate are applied to the S140 molecular cloud. Maps of this cloud in various molecular lines are examined which reveal the existence of a strong temperature and density peak near the H-alpha rim and the presence of a strong near-IR source very near that peak. This source is found to have an energy distribution which is very similar to that of the Becklin-Neugebauer object, including the silicate feature at 9.5 microns. The heating and cooling rates for the gas and dust in S140 are estimated, and the predicted cooling rate for dust is shown to be much greater than that for gas. The molecular-line observations are compared with a spherical large-velocity-gradient trapping calculation. An optimum cloud model is developed which indicates that: (1) (C-13)O, HCO(+), and H2CO are about 10 times less abundant in S140 than in dark clouds; (2) the HCN abundance is approximately 4 times the H2CO abundance in S140, compared with roughly equal values in dark clouds; and (3) the (C-13)O/C(O-18) ratio is consistent with terrestrial isotope ratios.
The present conference on astronomical studies of the Orion Nebula covers molecular clouds in Orion, the use of the nebula's cloak as a model for gas super-shells around OB associations, optical and UV data concerning the nebula's physical conditions, the presence of atomic carbon in Orion, large scale distribution of far-IR and sub-mm line emission, star formation studies in the IR, gas dynamics in the circumstellar nebula of the Becklin-Neugebauer source, maser sources and far-IR CO line emission in Orion-KL, and synthesis maps of mm molecular lines. Also discussed are Orion's star distributions, core region nebular condensations, energetic molecular flows in star-forming cloud cores, IR observations of HH objects, compact continuum radio sources, the SiO maser, shock waves, and the chemical evolution of OB associations.
Nanoindentation can be leveraged to aid in the high fidelity modeling of dislocation mediated plasticity in pentaerythritol tetranitrate (PETN), an anisotropic energetic molecular crystal. Moreover, nanoindentation tip parameters such as tip geometry, size, and degree of acuity can be utilized to target anisotropic behavior. In this work, nanoindentation was conducted across a range of orientations on the (110) face of PETN to characterize resultant yield behavior, mechanical property measurements, and resultant slip behavior and fracture initiation. Three different indentation tips were utilized: a 3-sided pyramidal Berkovich tip, a 4-sided high aspect ratio Knoop tip, and a 90° conical tip. Ultimately, indenter tip radius was documented to impact yield behavior, whereas tip geometry affected larger scale processes such as slip, and tip acuity was the dominating factor that led to fracture. The axisymmetric conical tip, serving as a baseline, showed the least amount of variation in mechanical property measurements but also the largest distribution of maximum shear stress at which initial yielding occurred. Its high degree of acuity, however, was more prone to induce fracture at higher loads. The Knoop tip was shown to be suitable for average measurements, but also for elucidation of certain anisotropic features. A distinctly higher perceived hardness at 45° was measured with the Knoop tip, indicating less dislocation motion in that direction also observed in this work via scanning probe microscopy. Lastly, the commonly used Berkovich tip was a good compromise whereby it provided a representative volume element describing the average behavior of the material. These results can be utilized to target desired anisotropic behavior in a wider range of molecular crystals, as well as to inform theoretical considerations for dislocation mediated plasticity in PETN.
The processes by which energetic electrons lose energy in a weakly ionized gas of molecular hydrogen are analyzed, and calculations are carried out taking into account the discrete nature of the excitation processes. The excitation, ionization, and heating efficiencies are computed for electrons with energies up to 100 eV absorbed in a gas with fractional ionizations up to 0.01, and the mean energy per pair of neutral hydrogen atoms is calculated.
This program explored a broad range of charge transfer (CT) states in organic solar cells with the goal of understanding how molecular energetics, structure, and morphology influence CT state energy and dynamics, and how this in turn influences solar cell performance. We identified new types of delocalized CT states at ordered molecular heterojunctions, determined the factors that dictate CT state energetic disorder in bulk heterojunctions, determined that hot CT state dissociation is negligible for most systems, that the occupation of CT states is generally non-thermal in disordered heterojunctions, and that quasi-equilibrium generally does not hold in organic solar cells. These findings allowed us to demonstrate the lowest loss in potential of any organic solar cell to date and provide valuable guidance for designing higher efficiency organic solar cells in the future.
In this study, high-performance organic photodetectors are presented which utilize a pristine chlorinated subphthalocyanine photoactive layer. Optical and optoelectronic analyses indicate that the device photocurrent is primarily generated through direct charge generation within the chlorinated subphthalocyanine layer, rather than exciton separation at layer interfaces. Molecular modelling suggests that this direct charge generation is facilitated by chlorinated subphthalocyanine high octupole moment (-80 DÅ 2 ), which generates a 200 meV shift in molecular energetics. Increasing the thickness of chlorinated subphthalocyanine leads to faster response time, correlated with a decrease in trap density. Notably, photodetectors with a 50 nm thick chlorinated subphthalocyanine photoactive layer exhibit detectivities approaching 10 13 Jones, with a dark current below 10 -7 A cm -2 up to -5 V. Based on these findings, we conclude that high octupole moment molecular semiconductors are promising materials for high-performance organic photodetectors employing single-component photoactive layer.
Using reverse non-equilibrium molecular dynamics simulations, we have determined the dependences on temperature and pressure of the thermal conductivity tensors for the monoclinic γ and ε polymorphs of hexanitrohexaazaisowurtzitane (HNIW or CL20). In this work, a recently developed non-reactive force field [X. Bidault and S. Chaudhuri, RSC Adv. 9, 39649–39661 (2019)], designed to study polymorphism and phase transitions in CL20, is employed. The effects of temperature and pressure are investigated between 200 and 500 K and up to 0.5 GPa for γ-CL20 and 2 GPa for ε-CL20. In order to obtain the full thermal conductivity tensor, $κ_{ij}$, for the monoclinic crystals, four distinct heat propagation directions are used. We find that $κ_{ij}$ for both polymorphs is more isotropic than for other energetic molecular crystals, including α- and γ-RDX, β-HMX, and PETN, with a maximum difference of 9.8% between orientations observed at 300 K and 0 GPa for γ-CL20 and a maximum difference of 4.8% for ε-CL20. The average thermal conductivity, $\bar{κ}$, of ε-CL20 is 6.4% larger than that of γ-CL20 at 300 K and 0 GPa. Analytic linear functions of the inverse temperature and the pressure are provided, which fit the data well and can be used to predict the thermal conductivity of both polymorphs for any orientation, pressure, and temperature in and around the fitting range. Our predictions agree reasonably well with the limited available experimental data, for which the polymorph type is unknown.
Many important aspects of the dynamic thermo-mechanical response of materials occur at the mesoscale, i.e. a physical scale of interactions smaller than what can be adequately described by homogenous behaviors, yet larger than the scale of the atomic lattice. Concurrent advancements in computational power, continuum theory, and experimental diagnostics are enabling unprecedented understanding of such interactions. However, we cannot develop a sufficient level of confidence in such mesoscale capability until the constitutive description of the underlying constituents is reliably representative of their actual physical behavior. Therefore, there is a strong need to combine experimental, modeling, and data-science techniques to validate models of the thermomechanical response of individual single crystals. One experimental diagnostic with high potential impact to shock physics and materials science is in-situ x-ray diffraction. This paper is primarily focused on simulation of x-ray diffraction in shock physics, but with an aim toward quantifying parametric uncertainty of simulation models. Here, we develop and demonstrate a data-science and model-driven approach to constrain the parameterization of continuum models of crystal lattice deformation associated with the shock response of crystalline materials. The framework is built around the connection between continuum hydrodynamic simulations of lattice deformation and a new Bragg diffraction simulation code, BarberShop. The dynamic deformation of a crystal lattice is modeled using the DiscoFlux model within an arbitrary Lagrangian-Eulerian hydrodynamic code, FLAG. These detailed continuum simulations of lattice deformation can be computationally slow, thus a statistical model is used to emulate the evolution of lattice deformation fields in time and across the considered model parameter space. Emulated lattice deformation fields can then be generated rapidly for any combination of physics model parameters. In turn, these fields can be fed into BarberShop to realize a rapid prediction of Bragg diffraction patterns associated with particular values of physics model parameters. The framework enables parameterization of the single crystal model to obtain Bragg diffraction patterns that most closely resemble a corresponding measurement. Furthermore, the framework naturally provides sensitivities of the lattice deformation to the physics parameters. We highlight the utility of this framework through the application to a synthetic closed-loop inverse problem leading to the parameterization of a single crystal material model. As a model problem, we consider the dynamic response of the energetic molecular crystal, cyclotrimethylenetrinitramine (or RDX), under dynamic compression induced by simulated flyer plate impact experiments.
The MSCG/FM (multiscale coarse-graining via force-matching) approach is an efficient supervised machine learning method to develop microscopically informed coarse-grained (CG) models. Here we present a theory based on the principle of maximum entropy (PME) enveloping the existing MSCG/FM approaches. This theory views the MSCG/FM method as a special case of matching the thermodynamic forces from the extended ensemble described by the set of thermodynamic (relevant) system coordinates. This set may include CG coordinates, the stress tensor, applied external fields, and so forth, and may be characterized by nonequilibrium conditions. Following the presentation of the theory, we discuss the consistent matching of both bonded and nonbonded interactions. The proposed PME formulation is used as a starting point to extend the MSCG/FM method to the constant strain ensemble, which together with the explicit matching of the bonded forces is better suited for coarse-graining anisotropic media at a submolecular resolution. The theory is demonstrated by performing the fine coarse-graining of crystalline 1,3,5-triamino-2,4,6-trinitrobenzene (TATB), a well-known insensitive molecular energetic material, which exhibits highly anisotropic mechanical properties.