Search NASA⌕ Search

SEARCH · Search NASA

Results for “free energy sampling”

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.

At least 19 records
All sourcesNASA NTRSDOE OSTI
NASA NTRS2021

Anion-Assisted Delivery of Multivalent Cations to Inert Electrodes

To understand and control key electrochemical processes - metal plating, corrosion, intercalation, etc. requires molecular-scale details of the active species at electrochemical interfaces and their mechanisms for de-solvation from the electrolyte. Using free energy sampling techniques we reveal the interfacial speciation of divalent cations in ether-based electrolytes and mechanisms for their delivery to an inert graphene electrode interface. Surprisingly, we find that anion solvophobicity drives a high population of anion-containing species to the interface that facilitate the delivery of divalent cations, even to negatively charged electrodes. Our simulations indicate that cation desolvation is greatly facilitated by cation-anion coupling. We propose anion solvophobicity as a molecular-level descriptor for rational design of electrolytes with increased efficiency for electrochemical processes limited by multivalent cation desolvation.

Electrochemical interfaces↗
NASA NTRS2005

Optimal Sampling of a Reaction Coordinate in Molecular Dynamics

Estimating how free energy changes with the state of a system is a central goal in applications of statistical mechanics to problems of chemical or biological interest. From these free energy changes it is possible, for example, to establish which states of the system are stable, what are their probabilities and how the equilibria between these states are influenced by external conditions. Free energies are also of great utility in determining kinetics of transitions between different states. A variety of methods have been developed to compute free energies of condensed phase systems. Here, I will focus on one class of methods - those that allow for calculating free energy changes along one or several generalized coordinates in the system, often called reaction coordinates or order parameters . Considering that in almost all cases of practical interest a significant computational effort is required to determine free energy changes along such coordinates it is hardly surprising that efficiencies of different methods are of great concern. In most cases, the main difficulty is associated with its shape along the reaction coordinate. If the free energy changes markedly along this coordinate Boltzmann sampling of its different values becomes highly non-uniform. This, in turn, may have considerable, detrimental effect on the performance of many methods for calculating free energies.

Pohorille, Andrew↗
NASA NTRS1990

Resonant tunneling IR detectors

Researchers propose a novel semiconductor heterojunction photodetector which would have a very low dark current and would be voltage tunable. A schematic diagram of the device and its band structure are shown. The two crucial components of the device are a cathode (InGaAs) whose condition band edge is below the conduction band edge of the quantum wells and a resonant tunneling filter (GaAs-AlGaAs). In a standard resonant tunneling device the electrodes are made of the same material as the quantum wells, and this device becomes highly conducting when the quantum levels in the wells are aligned with the Fermi level in the negatively biased electrode. In contrast, the researchers device is essentially non-conducting under the same bias conditions. This is because the Fermi Level of the cathode (InGaAs) is still well below the quantum levels so that no resonant transport occurs and the barriers (AlGaAs) effectively block current flow through the device. However, if light with the same photon energy as the conduction-band discontinuity between the cathode and the quantum wells, E sub c3-E sub c1, is shone on the sample, free carriers will be excited to an energy corresponding to the lowest quantum level in the well closest to the cathode (hv plue E sub c1 = E sub o). These electrons will resonantly tunnel through the quantum wells and be collected as a photocurrent in the anode (GaAs). To improve the quantum efficiency, the cathode (InGaAs) should be very heavily doped and capped with a highly reflective metal ohmic contact. The thickness of the device should be tailored to optimize thin film interference effects and afford the maximum absorption of light. Because the device relies on resonant tunneling, its response should be very fast, and the small voltages needed to change the responsivity should allow for very high frequency modulation of the photocurrent. In addition, the device is tuned to a specific photon energy so that it can be designed to detect a fairly narrow range of wavelengths. This selectivity is important for reducing the photocurrent due to spurious light sources.

Woodall, Jerry M.↗
NASA NTRS1985

Statistical mechanics of light elements at high pressure. VII - A perturbative free energy for arbitrary mixtures of H and He

A model free energy is presented which accurately represents results from 45 high-precision Monte Carlo calculations of the thermodynamics of hydrogen-helium mixtures at pressures of astrophysical and planetophysical interest. The free energy is calculated using free-electron perturbation theory (dielectric function theory), and is an extension of the expression given in an earlier paper in this series. However, it fits the Monte Carlo results more accurately, and is valid for the full range of compositions from pure hydrogen to pure helium. Using the new free energy, the phase diagram of mixtures of liquid metallic hydrogen and helium is calculated and compared with earlier results. Sample results for mixing volumes are also presented, and the new free energy expression is used to compute a theoretical Jovian adiabat and compare the adiabat with results from three-dimensional Thomas-Fermi-Dirac theory. The present theory gives slightly higher densities at pressures of about 10 megabars.

Hubbard, W. B.↗
NASA NTRS2000

Rate Constant and RRKM Product Study for the Reaction Between CH3 and C2H3 at T = 298K

The total rate constant k1 has been determined at P = 1 Torr nominal pressure (He) and at T = 298 K for the vinyl-methyl cross-radical reaction CH3 + C2H3 yields products. The measurements were performed in a discharge flow system coupled with collision-free sampling to a mass spectrometer operated at low electron energies. Vinyl and methyl radicals were generated by the reactions of F with C2H4 and CH4, respectively. The kinetic studies were performed by monitoring the decay of C2H3 with methyl in excess, 6 < |CH3|(sub 0)/|C2H3|(sub 0) < 21. The overall rate coefficient was determined to be k1(298 K) = (1.02 +/- 0.53)x10(exp -10) cubic cm/molecule/s with the quoted uncertainty representing total errors. Numerical modeling was required to correct for secondary vinyl consumption by reactions such as C2H3 + H and C2H3 + C2H3. The present result for k1 at T = 298 K is compared to two previous studies at high pressure (100-300 Torr He) and to a very recent study at low pressure (0.9-3.7 Torr He). Comparison is also made with the rate constant for the similar reaction CH3 + C2H5 and with a value for k1 estimated by the geometric mean rule employing values for k(CH3 + CH3) and k(C2H3 + C2H3). Qualitative product studies at T = 298 K and 200 K indicated formation of C3H6, C2H2, and C2H5 as products of the combination-stabilization, disproportionation, and combination-decomposition channels, respectively, of the CH3 + C2H3 reaction. We also observed the secondary C4H8 product of the subsequent reaction of C3H5 with excess CH3; this observation provides convincing evidence for the combination-decomposition channel yielding C3H5 + H. RRKM calculations with helium as the deactivator support the present and very recent experimental observations that allylic C-H bond rupture is an important path in the combination reaction. The pressure and temperature dependencies of the branching fractions are also predicted.

Thorn, R. Peyton, Jr.↗
NASA NTRS2013

Good Practices in Free-energy Calculations

As access to computational resources continues to increase, free-energy calculations have emerged as a powerful tool that can play a predictive role in drug design. Yet, in a number of instances, the reliability of these calculations can be improved significantly if a number of precepts, or good practices are followed. For the most part, the theory upon which these good practices rely has been known for many years, but often overlooked, or simply ignored. In other cases, the theoretical developments are too recent for their potential to be fully grasped and merged into popular platforms for the computation of free-energy differences. The current best practices for carrying out free-energy calculations will be reviewed demonstrating that, at little to no additional cost, free-energy estimates could be markedly improved and bounded by meaningful error estimates. In energy perturbation and nonequilibrium work methods, monitoring the probability distributions that underlie the transformation between the states of interest, performing the calculation bidirectionally, stratifying the reaction pathway and choosing the most appropriate paradigms and algorithms for transforming between states offer significant gains in both accuracy and precision. In thermodynamic integration and probability distribution (histogramming) methods, properly designed adaptive techniques yield nearly uniform sampling of the relevant degrees of freedom and, by doing so, could markedly improve efficiency and accuracy of free energy calculations without incurring any additional computational expense.

Pohorille, Andrew↗
NASA NTRS1990

Thermodynamic properties and crystallization kinetics at high liquid undercooling

The heat capacities of liquid and crystalline Au-Pb-Sb alloys in the glass-forming composition range were measured with droplet emulsion and bulk samples. Based on the measured C(sub p) data, the entropy, enthalpy, and Gibbs free energy functions of the eutectic, solid mixture, and undercooled liquid were determined as a function of undercooling and compared with theoretical predictions. The results indicate an isentropic temperature at 313 + or - 5 K, which agrees well with experimental data for the glass transition. A kinetics analysis of the nucleation undercooling response suggests that the proper choice for the Gibbs free energy change during crystallization is most important in analyzing the nucleation kinetics. By classical nucleation theory, the prefactors obtained, based on a variety of theoretical predictions for the driving force, can differ by six orders of magnitude. If the nucleation rates are extrapolated to high undercooling, the extrapolations based on measured heat capacity data show agreement, whereas the predicted nucleation rates are inconsistent with results from drop tower experiments. The implications for microg experiments are discussed.

Fecht, Hans J.↗
NASA NTRS2017

Testing the Accuracy of Data-Driven MHD Simulations of Active Region Evolution

Models for the evolution of the solar coronal magnetic field are vital for understanding solar activity, yet the best measurements of the magnetic field lie at the photosphere, necessitating the development of coronal models which are "data-driven" at the photosphere. We present an investigation to determine the feasibility and accuracy of such methods. Our validation framework uses a simulation of active region (AR) formation, modeling the emergence of magnetic flux from the convection zone to the corona, as a ground-truth data set, to supply both the photospheric information and to perform the validation of the data-driven method. We focus our investigation on how the accuracy of the data-driven model depends on the temporal frequency of the driving data. The Helioseismic and Magnetic Imager on NASA's Solar Dynamics Observatory produces full-disk vector magnetic field measurements at a 12-minute cadence. Using our framework we show that ARs that emerge over 25 hours can be modeled by the data-driving method with only approximately 1 percent error in the free magnetic energy, assuming the photospheric information is specified every 12 minutes. However, for rapidly evolving features, under-sampling of the dynamics at this cadence leads to a strobe effect, generating large electric currents and incorrect coronal morphology and energies. We derive a sampling condition for the driving cadence based on the evolution of these small-scale features, and show that higher-cadence driving can lead to acceptable errors. Future work will investigate the source of errors associated with deriving plasma variables from the photospheric magnetograms as well as other sources of errors, such as reduced resolution, instrument bias, and noise.

Leake, James E.↗
NASA NTRS2022

Free Energy Landscapes for Elucidating the Structural Consequences of Exon-20 mutations on the ErbB Family of Protein Kinases

The ErbB family of protein kinases plays an important role in major cellular functions and consequently mutations in the functional regions of these proteins are implicated in several types of cancer growths. To envision rational design of small molecule drugs that target the diseased proteins it is important to quantify the structural effects of the mutations, as some of these mutants render the protein resistant to tyrosine kinase inhibitors (TKIs). Herein we use advanced sampling techniques and long-timescale molecular dynamics simulations to predict the effect of major exon 20 mutations on the ErbB family, specifically EGFR and HER2 proteins. Exon 20 mutations have been clinically known to induce TKI resistance, though the mechanisms of such an effect is poorly understood. By mapping out the free energy landscape of the mutants and comparing them against the wild-type, we elucidate the structural differences in the binding pocket region that alter the nature of drug-protein interactions. We believe that these insights will play a pivotal role in developing small molecule drugs that overcome the TKI resistance.

Ashwin Ravichandran↗
NASA NTRS2016

Cyclotron Resonant Scattering Feature Simulations I. Thermally Averaged Cyclotron Scattering Cross Sections, Mean Free Photon-Path Tables, and Electron Momentum Sampling

Electron cyclotron resonant scattering features (CRSFs) are observed as absorption-like lines in the spectra of X-ray pulsars. A significant fraction of the computing time for Monte Carlo simulations of these quantum mechanical features is spent on the calculation of the mean free path for each individual photon before scattering, since it involves a complex numerical integration over the scattering cross section and the (thermal) velocity distribution of the scattering electrons. We aim to numerically calculate interpolation tables which can be used in CRSF simulations to sample the mean free path of the scattering photon and the momentum of the scattering electron. The tables also contain all the information required for sampling the scattering electron's final spin. The tables were calculated using an adaptive Simpson integration scheme. The energy and angle grids were refined until a prescribed accuracy is reached. The tables are used by our simulation code to produce artificial CRSF spectra. The electron momenta sampled during these simulations were analyzed and justified using theoretically determined boundaries. We present a complete set of tables suited for mean free path calculations of Monte Carlo simulations of the cyclotron scattering process for conditions expected in typical X-ray pulsar accretion columns (0.01 B/B(sub crit) is less than or equal to 0.12, where B(sub crit) = 4.413 x 1013 G, and 3 keV kBT 15 keV). The sampling of the tables is chosen such that the results have an estimated relative error of at most 1/15 for all points in the grid. The tables are available online (see link in footnote, page 1).

Schwarm, F.-W.↗
NASA NTRS1982

A two-dimensional phase separation on the spherical surface of the metallic glass Au55Pb22.5Sb22.5

Recent experiments indicate that a phase separation in a spherical sample of the metallic glass Au55Pb22.5Sb22.5 occurs near the surface of the sphere. This strongly suggests either a contribution of surface-free energy to the decomposition process or a possible influence of near surface impurities absorbed during synthesis of the sphere. The surface phase separation has been studied as a function of cooling rate of the sphere. At high cooling rates (small sphere sizes), the surface separation disappears altogether suggesting that the surface of the parent liquid droplet is initially homogeneous.

Lee, M. C.↗
NASA NTRS2018

Hydrothermal Decomposition of Amino Acids and Origins of Prebiotic Meteoritic Organic Compounds

The organic compounds found in carbonaceous chondrite meteorites provide insight into primordial solar system chemistry. Evaluating the formation and decomposition mechanisms of meteoritic amino acids may aid our understanding of the origins of life and homochirality on Earth. The amino acid glycine is widespread in meteorites and other extraterrestrial environments; other amino acids, such as isovaline, are found with enantiomeric excesses in some meteorites. The relationship between meteoritic amino acids and other compounds with similar molecular structures, such as aliphatic monoamines and monocarboxylic acids is unclear; experimental results evaluating the decomposition of amino acids have produced inconclusive results about the preferred pathways, reaction intermediates, and if the conditions applied may be compatible with those occurring inside meteoritic parent bodies. In this work, we performed extensive tandem metadynamics, umbrella sampling, and committor analysis to simulate the neutral mild hydrothermal decomposition mechanisms of glycine and isovaline and put them into context for the origins of meteoritic organic compounds. Our ab initio simulations aimed to determine free energy profiles and decomposition pathways for glycine and isovaline. We found that under our modeled conditions, methylammonium, glycolic acid, and sec-butylamine are the most likely decomposition products. These results suggest that meteoritic aliphatic monocarboxylic acids are not produced from decomposition of meteoritic amino acids. Our results also indicate that the decomposition of L-isovaline prefers an enantioselective pathway resulting in the production of (S)-sec-butylamine.

Enhanced sampling↗
NASA NTRS1988

A rapid method for the computation of equilibrium chemical composition of air to 15000 K

A rapid computational method has been developed to determine the chemical composition of equilibrium air to 15000 K. Eleven chemically reacting species, i.e., O2, N2, O, NO, N, NO+, e-, N+, O+, Ar, and Ar+ are included. The method involves combining algebraically seven nonlinear equilibrium equations and four linear elemental mass balance and charge neutrality equations. Computational speeds for determining the equilibrium chemical composition are significantly faster than the often used free energy minimization procedure. Data are also included from which the thermodynamic properties of air can be computed. A listing of the computer program together with a set of sample results are included.

Prabhu, Ramadas K.↗
NASA NTRS1972

Dielectric properties of Apollo 14 lunar samples.

Laboratory characterization of dielectric properties of lunar samples 14301,41, 14310,75, 14318,30 and 14321,163 is made. Our measurements of dielectric constants and losses were made over a range of frequency from 100 Hz to 10 MHz and temperature from 77 to 473 K by two-terminal capacitance substitution methods. The dielectric behavior of these samples is generally similar to other lunar samples from Apollo 11 and Apollo 12 sites when these samples are free from absorbed moisture. As did sample 12002,58, sample 14310,75 showed a distinctive dispersion that may be associated with the presence of water, and different values of the activation energy for conduction as the temperature was varied. The activation energies range from about 0.03 to 0.5 eV.

Chung, D. H.↗
NASA NTRS1975

Analytical and experimental evaluation of techniques for the fabrication of thermoplastic hologram storage devices

The results of an experimental investigation on recording information on thermoplastic are given. A description was given of a typical fabrication configuration, the recording sequence, and the samples which were examined. There are basically three configurations which can be used for the recording of information on thermoplastic. The most popular technique uses corona which furnishes free charge. The necessary energy for deformation is derived from a charge layer atop the thermoplastic. The other two techniques simply use a dc potential in place of the corona for deformation energy.

Rogers, J. W.↗
NASA NTRS2015

Magnetogram Forecast: An All-Clear Space Weather Forecasting System

Solar flares and coronal mass ejections (CMEs) are the drivers of severe space weather. Forecasting the probability of their occurrence is critical in improving space weather forecasts. The National Oceanic and Atmospheric Administration (NOAA) currently uses the McIntosh active region category system, in which each active region on the disk is assigned to one of 60 categories, and uses the historical flare rates of that category to make an initial forecast that can then be adjusted by the NOAA forecaster. Flares and CMEs are caused by the sudden release of energy from the coronal magnetic field by magnetic reconnection. It is believed that the rate of flare and CME occurrence in an active region is correlated with the free energy of an active region. While the free energy cannot be measured directly with present observations, proxies of the free energy can instead be used to characterize the relative free energy of an active region. The Magnetogram Forecast (MAG4) (output is available at the Community Coordinated Modeling Center) was conceived and designed to be a databased, all-clear forecasting system to support the operational goals of NASA's Space Radiation Analysis Group. The MAG4 system automatically downloads nearreal- time line-of-sight Helioseismic and Magnetic Imager (HMI) magnetograms on the Solar Dynamics Observatory (SDO) satellite, identifies active regions on the solar disk, measures a free-energy proxy, and then applies forecasting curves to convert the free-energy proxy into predicted event rates for X-class flares, M- and X-class flares, CMEs, fast CMEs, and solar energetic particle events (SPEs). The forecast curves themselves are derived from a sample of 40,000 magnetograms from 1,300 active region samples, observed by the Solar and Heliospheric Observatory Michelson Doppler Imager. Figure 1 is an example of MAG4 visual output

Barghouty, Nasser↗
NASA NTRS1994

Low encounter speed comet COMA sample return missions

Comets, being considered the most primitive bodies in the solar system, command the highest priority among solar-system objects for studying solar nebula evolution and the evolution of life through biogenic elements and compounds. The study of comets, and more especially, of material from them, provides an understanding of the physical, chemical, and mineralogical processes operative in the formation and earliest development of the solar systems. These return samples will provide valuable information on comets and serve as a rosetta stone for the analytical studies conducted on interplanetary dust particles over the past two decades, and will provide much needed extraterrestrial samples for the planetary materials community since the Apollo program. Lander sample return missions require rather complex spacecraft, intricate operations, and costly propulsion systems. By contrast, it is possible to take a highly simplified approach for sample capture and return in the case of a comet. In the past, we have considered Earth free-return trajectory to the comet, in which passive collectors intercept dust and volatiles from the cometary coma. However, standard short period cometary free-return trajectories results in the comet to the spacecraft encounter speeds in the range of 10 km/s. At these speeds the kinetic energy of the capture process can render significant modification of dust structure, change of solid phase as well as the lost of volatiles components. This paper presents a class of new missions with trajectories with significant reduction of encounter speeds by incorporating gravity assists and deep space maneuvering. Low encounter speed cometary flyby sample return will enable a marked increase in the value of the return science. Acquiring thousands of samples from a known comet and thousands of images of a comet nucleus would be space firsts. Applying new approach in flight mechanics to generate a new class of low encounter speed cometary sample return trajectories opens new possibilities in science. A systematic search of trajectories for the first decade of the twenty-first century will be made. The target encounter speed is for less than 7 km/s to short period comets.

Tsou, P.↗
NASA NTRS1995

Microstructural Characterization of Reaction-Formed Silicon Carbide Ceramics

Microstructural characterization of two reaction-formed silicon carbide ceramics has been carried out by interference layering, plasma etching, and microscopy. These specimens contained free silicon and niobium disilicide as minor phases with silicon carbide as the major phase. In conventionally prepared samples, the niobium disilicide cannot be distinguished from silicon in optical micrographs. After interference layering, all phases are clearly distinguishable. Back scattered electron (BSE) imaging and energy dispersive spectrometry (EDS) confirmed the results obtained by interference layering. Plasma etching with CF4 plus 4% O2 selectively attacks silicon in these specimens. It is demonstrated that interference layering and plasma etching are very useful techniques in the phase identification and microstructural characterization of multiphase ceramic materials.

Singh, M.↗