SEARCH · Engineering Papers
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.
A new computational screening pipeline accelerates the process of identifying potential Li-ion cathodes, with a demonstration case upon its preliminary application.
Macroscopic phase equilibria between solutions define the functionality of many biological and industrial processes, yet they are challenging to predict due to the inherent complexity of liquids containing large molecules. This work introduces an approach for the purely predictive calculation of such phase equilibria in temperature-composition space from molecular dynamics (MD) simulations at one temperature in the single-phase region. We use an approach developed previously to obtain the entropic and enthalpic contributions to the free energy of mixing from the atomic-scale information given by MD simulations via Kirkwood–Buff theory. This allows us to accurately estimate the free energy of mixing as a function of temperature, and thus obtain liquid–liquid phase equilibria, including liquid–liquid critical points, associated binodal and spinodal lines, and composition fluctuations across a region of temperature and composition. Results for binary malonamide–alkane systems are validated by comparison to a direct experimental probe of the fluctuations: the small angle X-ray scattering intensity near zero wavenumber. The MDKB → Phase method demonstrated here provides a significant improvement in predicting liquid–liquid equilibria and free energy as a function of temperature for our systems of interest compared to conventional thermodynamic models. The accurate performance of this purely predictive approach lies in its preservation of atomistic details when determining thermodynamic properties. Furthermore, its inherent extensibility to multi-component systems will likely make the MDKB → Phase approach a valuable general tool for connecting molecular interactions to macroscopic phase equilibria and for the computational screening of materials for targeted thermodynamic behavior.
We report solid-state electrolyte materials with superior lithium ionic conductivities are vital to the next-generation Li-ion batteries. Molecular dynamics could provide atomic scale information to understand the diffusion process of Li-ion in these superionic conductor materials. Here, we implement the deep potential generator to set up an efficient protocol to automatically generate interatomic potentials for Li 10 GeP 2 S 12 -type solid-state electrolyte materials (Li 10 GeP 2 S 12 , Li 10 SiP 2 S 12 , and Li 10 SnP 2 S 12 ). The reliability and accuracy of the fast interatomic potentials are validated. With the potentials, we extend the simulation of the diffusion process to a wide temperature range (300 K–1000 K) and systems with large size (~1000 atoms). Important technical aspects such as the statistical error and size effect are carefully investigated, and benchmark tests including the effect of density functional, thermal expansion, and configurational disorder are performed. The computed data that consider these factors agree well with the experimental results, and we find that the three structures show different behaviors with respect to configurational disorder. Our work paves the way for further research on computation screening of solid-state electrolyte materials.
While computational screening with density functional theory (DFT) is frequently employed for the screening of metal–organic frameworks (MOFs) for gas separation and storage, commonly applied generalized gradient approximations (GGAs) exhibit self-interaction errors, which hinder the predictions of adsorption energies. We investigate the Hubbard U parameter to augment DFT calculations for full periodic MOFs, targeting a more precise modeling of gas molecule–MOF interactions, specifically for N2, CO2, and O2. We introduce a calibration scheme for the U parameter, which is tailored for each MOF, by leveraging higher-level calculations on the secondary building unit (SBU) of the MOF. When applied to the full periodic MOF, the U parameter calibrated against hybrid HSE06 calculations of SBUs successfully reproduces hybrid-quality calculations of the adsorption energy of the periodic MOF. The mean absolute deviation of adsorption energies reduces from 0.13 eV for a standard GGA treatment to 0.06 eV with the calibrated U, demonstrating the utility of the calibration procedure when applied to the full MOF structure. Furthermore, attempting to use coupled cluster singles and doubles with perturbative triples calculations of isolated SBUs for this calibration procedure shows varying degrees of success in predicting the experimental heat of adsorption. It improves accuracy for N2 adsorption for cases of overbinding, whereas its impact on CO2 is minimal, and ambiguities in spin state assignment hinder consistent improvements of O2 adsorption. Our findings emphasize the limitations of cluster models and advocate the use of full periodic MOF systems with a calibrated U parameter, providing a more comprehensive understanding of gas adsorption in MOFs.
This paper presents a fully ab initio many-body photoemission framework that includes coherent three-body electron-photon-phonon scattering to predict the transverse momentum distributions and the mean transverse energies (MTEs) of bulk photoelectrons from single-crystal photocathodes. The need to develop such a theory stems from the lack of studies that provide complete understanding of the underlying fundamental processes governing the transverse momentum distribution of photoelectrons emitted from single crystals. For example, initial predictions based on density-functional theory calculations of effective electron masses suggested that the (111) surface of PbTe would produce very small MTEs (≤ 15 meV), whereas our experiments yielded MTEs ten to twenty times larger than these predictions and also exhibited a lower photoemission threshold than predicted. The ab initio framework presented in this paper correctly reproduces both the magnitude of the MTEs from our measurements in PbTe(111) and the observed photoemission below the predicted threshold. Our results show that both photoexcitations into states that propagate in the bulk of the material and coherent many-body electron-photon-phonon scattering processes, which initial predictions ignored, play surprisingly important roles in photoemission from PbTe(111). Lastly, from the lessons learned, we recommend a procedure for rapid computational screening of potential single-crystal photocathodes for applications in next-generation ultrafast electron diffraction and X-ray free-electron lasers, which will enable significant advances in condensed matter research.
The remaining goals of the project were the design OSDAs for zeolites for 1) ethylene/ethane separation, 2) CO2 /nitrogen separation, 3) CO2 /methane separation. These separation methods are recognized as steps in strategies to mitigate global warming. We have deployed targeted de novo design to design OSDAs directed towards zeolite structures that have been identified as hits in computational screens for zeolites that are predicted to be effective in the above mentioned gas separations. In addition, the impact on zeolite science of the methodology developed and refined in the course of the project has been recognized by the international zeolite research community. It will be the subject of an invited book chapter to be released in 2022. It is also illustrated by the co-recipient ship of M.W. Deem and F. Daeyaert of the 2019 Donald W. Breck Award in Molecular Sieve Science.
Computation screening is crucial for designing efficient electrochemical catalysts for carbon dioxide (CO2R) reduction to valuable hydrocarbons and oxygenates. Herein, leveraging density functional theory calculations of the CO adsorption energy ΔE_CO on seventeen Cu terminations, we discover a strong linear correlation between ΔE_CO and the recently experimentally measured CO2R electrochemical currents (ACS Catal. 2022, 12, 11, 6578–6588). Examining the ab initio thermodynamics of the early critical intermediates CO*, COH*, and CHO*, we find that CO* → CHO* is the thermodynamically preferred step, and notably shows a volcano trend with the experimental currents where the maximum CO2R current corresponds to the moderate CHO* formation energy. Importantly, we show that increasing the step and kink density of the Cu termination not only enhances CO adsorption strength but also modulates the CO* → CHO* pathway, as respectively exemplified in the (941) and (741) facets. We also explain why (741) is exceptional with high CO2R activity as measured experimentally due to its relatively low activity toward the hydrogen evolution reaction compared with the other Cu surfaces. Beyond the general CO adsorption energy that only shows a linear trend with CO2R activity, we show that the reaction CO* → CHO* free energy is a descriptor that displays a volcano relationship with the overall CO2R activity on Cu facets.
This report summarizes key findings from DOE Award DE‑SC0007004 – Singlet Fission, documenting more than a decade of research led by Professor Josef Michl. The project advanced foundational understanding of singlet fission (SF) through integrated synthetic, spectroscopic, crystallographic, and computational studies. Across a diverse set of chromophores—including tetracenes, cibalackrot derivatives, BODIPY dimers, TDPPs, and fluorinated isobenzofurans—the team identified structural and packing features that enable or inhibit SF, clarified the roles of charge‑transfer intermediates and excimer formation, and developed simplified theoretical models and computational screening tools for predicting optimal geometries and excited‑state-energetics. The work significantly broadened design principles for SF‑active materials and provided insights relevant to next‑generation solar energy conversion and molecular photophysics.
Researchers worldwide are seeking to repurpose existing drugs or discover new drugs to counter the disease caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). A promising source of candidates for such studies is molecules that have been reported in the scientific literature to be drug-like in the context of viral research. However, this literature is too large for human review and features unusual vocabularies for which existing named entity recognition (NER) models are ineffective. We report here on a project that leverages both human and artificial intelligence to detect references to such molecules in free text. We present 1) a iterative model-in-the-loop method that makes judicious use of scarce human expertise in generating training data for a NER model, and 2) the application and evaluation of this method to the problem of identifying drug-like molecules in the COVID-19 Open Research Dataset Challenge (CORD-19) corpus of 198,875 papers. We show that by repeatedly presenting human labelers only with samples for which an evolving NER model is uncertain, our human-machine hybrid pipeline requires only modest amounts of non-expert human labeling time (tens of hours to label 1778 samples) to generate an NER model with an F-1 score of 80.5%—on par with that of non-expert humans—and when applied to CORD’19, identifies 10,912 putative drug-like molecules. This enriched the computational screening team’s targets by 3,591 molecules, of which 18 ranked in the top 0.1% of all 6.6 million molecules screened for docking against the 3CLPro protein.
Human Interleukin-6 (hIL-6) is a pro inflammatory cytokine that binds to its receptor, IL-6Rα followed by binding to gp130 and subsequent dimerization to form a hexamer signaling complex. As a critical inflammation mediator, hIL-6 is associated with a diverse range of diseases and monoclonal antibodies in clinical use that either target IL-6Rα or hIL-6 to inhibit signaling. Here, we perform high-throughput structure-based computational screening using ensemble docking for small-molecule antagonists for which the target conformations were taken from 600 ns long molecular dynamics simulations of the apo protein. Prior knowledge of the contact sites from binary complex studies and experimental work was incorporated into the docking studies. The top 20 scoring ligands from the in silico studies after post analysis were subjected to in vitro functional assays. Among these compounds, the ligand with the second-highest calculated binding affinity experimentally showed an ~84% inhibitory effect on IL6-induced STAT3 reporter activity at 10 μM concentration. This finding may pave the way for designing small-molecule inhibitors of hIL-6 of therapeutic significance.
Practical methods of computer screening cloud-contaminated pixels from data of various satellite systems are proposed. Examples are given of the location of clouds and representative landscape features in HCMM spectral space of reflectance (VIS) vs emission (IR). Methods of screening out cloud affected HCMM are discussed. The character of subvisible absorbing-emitting atmospheric layers (subvisible cirrus or SCi) in HCMM data is considered and radiosonde soundings are examined in relation to the presence of SCi. The statistical characteristics of multispectral meteorological satellite data in clear and SCi affected areas are discussed. Examples in TIROS-N and NOAA-7 data from several states and Mexico are presented. The VIS-IR cluster screening method for removing clouds is applied to a 262, 144 pixel HCMM scene from south Texas and northeast Mexico. The SCi that remain after cluster screening are sited out by applying a statistically determined IR limit.
In rotating machinery dynamics an orbit (Lissajous curve) represents the dynamic path of the shaft centerline motion during shaft rotation and resulting precession. The orbit can be observed with an oscilloscope connected to XY promixity probes. The orbits can also be simulated by a computer. The software for HP computer simulates orbits for two cases: (1) Symmetric orbit with four frequency components with different radial amplitudes and relative phase angles; and (2) Nonsymmetric orbit with two frequency components with two different vertical/horizontal amplitudes and two different relative phase angles. Each orbit carries a Keyphasor mark (one-per-turn reference). The frequencies, amplitudes, and phase angles, as well as number of time steps for orbit computation, have to be chosen and introduced to the computer by the user. The orbit graphs can be observed on the computer screen.
A workstation was developed which uses the operator's eye movements and position to determine the placement of the cursor on a computer screen. A brainwave sensing technology overview and an introduction to the known rhythms or signals generated by the brain are given. This is followed by a descriptive explanation of the Ocular Attention Interface System (OASIS) and its intended integration into the proposed testbed.
Displays are now being used extensively throughout the society. More and more time is spent watching television, movies, computer screens, etc. Furthermore, in an increasing number of cases, the observer interacts with the display and plays the role of operator as well as observer. To a large extent, the normal behavior in the normal environment can also be thought of in these same terms. Taking liberties with Shakespeare, it might be said, all the world's a display and all the individuals in it are operators in and on the display. Within this general context of interactive display systems, a discussion is began with a conceptual overview of a particular class of such systems, namely, teleoperator systems. The notion is considered of telepresence and the factors that limit telepresence, including decorrelation between the: (1) motor output of the teleoperator as sensed directly via the kinesthetic/tactual system, and (2) the motor output of the teleoperator as sensed indirectly via feedback from the slave robot, i.e., via a visual display of the motor actions of the slave robot. Finally, the deleterious effect of time delay (a particular decorrelation) on sensory-motor adaptation (an important phenomenon related to telepresence) is examined.
Land managers increasingly are becoming dependent upon remote sensing and automated analysis techniques for information gathering and synthesis. Remote sensing and geographic information system (GIS) techniques provide quick and economical information gathering for large areas. The outputs of remote sensing classification and analysis are most effective when combined with a total natural resources data base within the capabilities of a computerized GIS. Some examples are presented of the successes, as well as the problems, in integrating remote sensing and geographic information systems. The need to exploit remotely sensed data and the potential that geographic information systems offer for managing and analyzing such data continues to grow. New microcomputers with vastly enlarged memory, multi-fold increases in operating speed and storage capacity that was previously available only on mainframe computers are a reality. Improved raster GIS software systems have been developed for these high performance microcomputers. Vector GIS systems previously reserved for mini and mainframe systems are available to operate on these enhanced microcomputers. One of the more exciting areas that is beginning to emerge is the integration of both raster and vector formats on a single computer screen. This technology will allow satellite imagery or digital aerial photography to be presented as a background to a vector display.
The effects of competition on performance of a video-formatted task were examined in a series of experiments. Two rhesus monkeys (Macaca mulatta) were trained to manipulate a joystick to shoot at moving targets on a computer screen. The task was made competitive by requiring both animals to shoot at the same target and by rewarding only the animal that hit the target first each trial. The competitive task produced a significant and robust speed-accuracy trade-off in performance. The monkeys hit the target in significantly less time on contested than on uncontested trials. However, they required significantly more shots to hit the target on contested trials in relation to uncontested trials. This effect was unchanged when various schedules of reinforcement were introduced in the uncontested trials. This supports the influence of competition qua competition on performance, a point further bolstered by other findings of behavioral contrast presented here.
The Automatic Gore Mapping System is being developed to reduce the time and labor costs associated with manufacturing the External Tank. The present chem-milling processes and procedures are discussed. The down loading of the simulation of the system has to be performed to verify that the simulation package will translate the simulation code into robot code. Also a simulation of this system has to be programmed for a gantry robot instead of the articulating robot that is presently in the system. It was discovered using the simulation package that the articulation robot cannot reach all the point on some of the panels, therefore when the system is ready for production, a gantry robot will be used. Also a hydrosensor system is being developed to replace the point-to-point contact probe. The hydrosensor will allow the robot to perform a non-contact continuous scan of the panel. It will also provide a faster scan of the panel because it will eliminate the in-and-out movement required for the present end effector. The system software is currently being modified so that the hydrosensor will work with the system. The hydrosensor consists of a Krautkramer-Branson transducer encased in a plexiglass nozzle. The water stream pumped through the nozzle is the couplant for the probe. Also, software is being written so that the robot will have the ability to draw the contour lines on the panel displaying the out-of-tolerance regions. Presently the contour lines can only be displayed on the computer screens. Research is also being performed on improving and automating the method of scribing the panels. Presently the panels are manually scribed with a sharp knife. The use of a low power laser or water jet is being studied as a method of scribing the panels. The contour drawing pen will be replaced with scribing tool and the robot will then move along the contour lines. With these developments the Automatic Gore Mapping Systems will provide a reduction in time and labor costs associated with manufacturing the External Task. The system also has the potential of inspecting other manufactured parts.
The objective is to make available a method of easily recording the microstructure of a sample by means of a computer. The method requires a minimum investment and little or no instruction on the operation of a computer. An outline of the setup involving a black and white TV camera, a digitizer control box, a metallurgical microscope and a computer screen, printer, and keyboard is shown.