Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “reversible computation”

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 1,513 records · Page 84
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Multi-state Catalysts Modulated by Mechanical Force (Final Report)

The development of more efficient catalytic processes and new approaches to control catalytic activity and selectivity are central to the realization of more selective, atom economic, and energy efficient routes to value added chemicals and polymers. The reactivity and selectivity of a transition metal catalyst is intimately related to the ligand-sphere geometry and, in many cases, the ideal ligand geometry for one step of a catalytic cycle is poorly matched to the ideal ligand geometry for another, resulting in sub-optimal efficiency. Macroscopic mechanical forces are both large, potentially much larger than interatomic forces, and are directional and localized to an extent that differentiates them from other forms of energy input such as heat or light. As such, mechanical force represents a heretofore untapped approach to modulate catalyst geometry, with the potential to reversibly modulate catalyst geometry on the timescale of catalytic turnover or monomer enchainment. This project has addressed the fundamental challenges in material-to-molecule strain coupling associated with the development of a new class of mechanically responsive catalysts (mechanocatalysts) in which active organotransition metal catalysts are strategically embedded in a flexible polymer network such that application of external mechanical force (stretching or deformation) leads to modulation of catalyst geometry, and hence reactivity and selectivity. Our efforts during the tenure of this grant were directed toward the elucidation of force-reactivity relationships of elementary transformations that occur within the first coordination sphere of a transition metal complex employing stiff stilbene photoswitches tethered to a flexible bidentate phosphine ligand derived from MeOBiphep as molecular force probes which provide a range of compressive and extension forces to the coupled transition metal complex depending on the geometry of the stiff stilbene and length of the tethering chains. During the tenure of this grant, we have quantified the rate of C(sp 2 )-C(sp 2 ) reductive elimination from platinum(II) diaryl complexes containing bis(phosphine) force probe ligands as a function of mechanical force; compressive forces decreased the rate of reductive elimination whereas extension forces increased the rate relative to the strain-free MeOBiphep complex with a 3.4-fold change in rate over a ~290 pN range of restoring forces. In a similar manner, we have quantified the rate of oxidative addition of bromobenzene to low-ligated palladium(0) complexes containing force probe ligands as a function of mechanical force; compressive forces increase the rate of oxidative addition, whereas tensile forces decrease the rate with a ~6 fold change in rate across ~340 pN of force applied to the complexes. In both cases, experimental and computational analyses argue strongly against any significant force-induced perturbation of ground state geometry within the first coordination sphere of the reactant complexes. Rather, the force/rate behavior observed for these transformations across these ranges of forces is attributed to the coupling of force to the nuclear motion comprising the reaction coordinates for reductive elimination and oxidative addition. These results together inform the development of catalysts whose activity can be tuned by an external force that is adjusted within a catalytic cycle and suggest opportunities to experimentally map geometry changes associated with reactions in transition metal complexes and potential strategies for force-modulated catalysis.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2022

The reaction between the bromine atom and the water trimer: high level theoretical studies

Three different reaction pathways are found for the reaction of bromine atom (Br) with the lowest-energy structure of the water trimer [uud-(H 2 O) 3 ], initially using the MPW1K-DFT method. The three bromine pathways have closely related geometries and energetics, analogous to those found for the fluorine and chlorine reactions. Here, the lowest-energy pathway of the Br + uud-(H 2 O) 3 reaction was further investigated using the “gold standard” CCSD(T) method and the correlation-consistent basis sets up to cc-pVQZ(-PP). Based on the CCSD(T)/cc-pVQZ(-PP)//CCSD(T)/cc-pVTZ(-PP) results, the Br + (H 2 O) 3 reaction is endothermic by 33.3 kcal mol -1 . The classical barrier height is 29.0 kcal mol -1 between the reactants and the exit complex, and there is no barrier for the reverse reaction. The Br···(H 2 O) 3 entrance complex is found to lie 4.7 kcal mol -1 below the separated reactants, and the HBr···(H 2 O) 2 OH exit complex is bound by 6.4 kcal mol -1 relative to the separated products. This potential energy profile is further corrected by the zero point energies and spin–orbit coupling effects. Structurally, the Br + (H 2 O) 3 stationary points can be derived from those of the simpler Br + (H 2 O) 2 reaction by judiciously appending a H 2 O molecule. The Br + (H 2 O) 3 potential energy profile is compared with the Br + (H 2 O) 2 and Br + H 2 O reactions, as well as to the valence isoelectronic Cl + (H 2 O) 3 and F + (H 2 O) 3 systems. It is reasonable to expect that the reactions between the bromine atom and larger water clusters would be similar to the Br + (H 2 O) 3 reaction.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Electrolytic gold plating, stripping, and ion transport dynamics through a solid-state iodide perovskite

The pronounced electrochemical reactivity between halide perovskites and metal electrodes can introduce mobile extrinsic metal ions which can cause device instability or enable novel functionalities. Here we systematically investigate the kinetics of gold cation (Au + ) migration in indium tin oxide (ITO)/methylammonium lead triiodide (MAPbI 3 )/Au model devices under long-term potentiostatic biasing. Scanning electron microscopy (SEM), X-ray photoelectron spectroscopy (XPS), and density functional theory (DFT) analyses reveal that Au + ions, electrochemically generated at the Au anode, traverse the perovskite layer with diffusion coefficients on the order of 10 −11 to 10 −10 cm 2 s −1 and are subsequently reduced at the cathode as Au 0 clusters, resembling metal plating behavior in electrolytic cells and solid-state batteries during charging. Furthermore, reversing the applied bias strips the plated Au 0 , revealing reversibility suitable for bipolar resistive switching devices and providing direct evidence of the electrochemical and ionic nature of Au transport within the perovskite matrix. Quantitatively determining diffusion coefficients and ion concentrations provides foundational inputs for future drift-diffusion modelling opportunities and allows us to relate our findings to implications on long term operation of devices like photovoltaic modules. These results clearly demonstrate the solid-state electrochemical nature of perovskite devices, highlight methods to be more quantitative about ion transport properties, provide and emphasize the importance of disentangling electro-, photo-, photoelectrochemical processes for understanding device performance and unlocking new functionalities.

14 SOLAR ENERGY↗
DOE OSTI2021

Data Processing Package for Cyclic Integrated Reversible Bending Fatigue Testing

A data processing software package has been introduced. The package was developed using MATLAB with the aid of the Curve Fitting Toolbox. The package is made up of four modules: pre-processing, data processing for static testing, data processing for monitoring, and data processing for measurements. CIRFT data are structured with multiple levels of architecture involving group, specimen, session, and scan/block. The degree of complexity of the data structure depends on whether a test is static or cyclic.The test results are presented in figures, scatter plots, and tables. For static testing, the output in tables provides bending mechanical properties and characteristic points of moment–curvature relation: flexural rigidities in linear segments of loading and unloading stages, intersection points between characteristic segments of the curve, and equivalent stress and strain quantities. For cyclic testing, the table output lists control and fatigue life and responsive/dependent quantities including moment, curvature, flexural rigidity, flexural hysteresis, rigidity phase angle, and equivalent stresses and strains. In addition, derivatives such as half-gage length and sensor spacing correction are included. The output also provides standard deviations of the reported quantities when they are applicable or available. The data processing package can serve as a fundamental characterization tool in mechanical study of materials. The package is intended primarily for data processing for the CIRFT process and can also be used in applications for which similar testing requirements exist.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Advancing the HERO WEC Through Integrated Modeling, Testing, and Field Deployments

The Hydraulic and Electric Reverse Osmosis Wave Energy Converter (HERO WEC) was developed by the National Laboratory of the Rockies as a modular platform for advancing wave-powered desalination technologies. Designed to operate in either a hydraulic or electric configuration, the system enables direct comparison of mechanical-to-water and electrical-to-water conversion pathways using a single hardware architecture. Deployments in 2022 and 2024 at the Jennette's Pier test site in North Carolina demonstrated freshwater production between 60 and 300 gallons per day, including successful operation in wave heights as low as 0.25 m. A structured evaluation approach combining numerical modeling of hydrodynamic and power take-off response, controlled laboratory testing, and in-ocean field deployments has been used to characterize and refine system performance. Specifically, insights from these efforts are now informing the second-generation HERO WEC (V2) to increase robustness and overall survivability. Changes from the original prototype include improved drivetrain components, refined geometry, enhanced control systems, and design. The HERO WEC platform continues to serve as an open-access research asset supporting broader marine energy and desalination development.

16 TIDAL AND WAVE POWER↗
DOE OSTI2017

Physiological fidelity or model parsimony? The relative performance of reverse-toxicokinetic modeling approaches

Physiologically-based toxicokinetic (PBTK) models are often developed to facilitate in vitro to in vivo extrapolation (IVIVE) using a top-down, compartmental approach, favoring architectural simplicity over physiological fidelity despite the lack of general guidelines relating model design to dynamical predictions. Here we explore the impact of design choice (high vs. low fidelity) on chemical distribution throughout an animal’s organ system. Results We contrast transient dynamics and steady states of three previously proposed PBTK models of varying complexity in response to chemical exposure. The steady states for each model were determined analytically to predict exposure conditions from tissue measurements. Steady state whole-body concentrations differ between models, despite identical environmental conditions, which originates from varying levels of physiological fidelity captured by the models. These differences affect the relative predictive accuracy of the inverted models used in exposure reconstruction to link effects-based exposure data with whole-organism response thresholds obtained from in vitro assay measurements. Conclusions Our results demonstrate how disregarding physiological fideltiy in favor of simpler models affects the internal dynamics and steady state estimates for chemical accumulation within tissues, which, in turn, poses significant challenges for the exposure reconstruction efforts that underlie many IVIVE methods. Developing standardized systems-level models for ecological organisms would not only ensure predictive consistency among future modeling studies, but also ensure pragmatic extrapolation of in vivo effects from in vitro data or modeling exposure-response relationships.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Visualizing diverse lithium growth and stripping behaviors in anode-free solid-state batteries with operando X-ray tomography

Anode-free solid-state batteries (SSBs), which eliminate the need for lithium metal use during cell assembly, have the potential to enable high energy densities and simplified manufacturing. However, the factors that control lithium growth/stripping at the anode current collector are not well understood. Here, we use operando X-ray microcomputed tomography to comprehensively image and quantify lithium deposition and stripping under various conditions in three different Li|Li 6 PS 5 Cl|current collector cells, revealing diverse behavior that depends on interface morphology, cell resistance, and solid-state electrolyte (SSE) microstructure. A cell with high resistance exhibits extensive lithium filament growth across the entire current collector interface, with filaments that grow around pre-existing pores in the SSE rather than lithium filling these pores. Lithium filament formation is partially reversible, with the cracks shrinking as lithium metal is stripped. Uniform lithium deposition is achievable at low current densities in low-resistance cells, whereas higher current densities in these cells cause an increase in interfacial roughness, which is correlated with subsequent filamentary growth at the edges of the cell. These results provide insight into filamentary vs. planar lithium growth and highlight that the evolution of lithium is sensitively dependent on SSE microstructure and electrochemical processes.

25 ENERGY STORAGE↗
DOE OSTI2016

Analysis of Reverse Flow Restriction Device to Prevent Dryout Fuel Damage during BWR Instability

Boiling Water Reactors (BWRs) are subjected to instabilities due to thermal hydraulics feedback or coupled Neutronics-thermal hydraulics feedback. Feedback is affected by many reactor parameters that could destabilize the core such as distribution of pressure drop across the core, axial power shape, subcooling level at the core inlet, radial power peaking factor, reactivity coefficients, and others. There are three common modes of instability in BWRs: single channel instability, core-wide instability, and regional instability. Each instability mode could damage the fuel if the power oscillations are left without an intervention like reactor scram. However, the instability problem becomes more dangerous if the automatic reactor scram is lost, like in Anticipated Transient without Scram (ATWS) accidents. A reverse flow restriction device (RFRD) would mitigate this problem. Core simulation using TRACE/PARCS code has been used to examine the potential effect of restricting the flow in the downward direction on mitigation of the power and flow oscillations to increase the safety of BWRs. The presented results demonstrate the ability of the device on limiting flow oscillation amplitude in the one bundle where it is used, on an adjacent bundle without it and reducing the power peak when the flow is low and power is high during BWR power oscillations. This exercise demonstrates the potential of previously described hardware device on preventing clad temperature excursion to keep the fuel safe from dryout damage and provide encouraging results ahead of the planned more detailed investigation where the device is to be used in more bundles up to full core. In the future work it is planned to examine the effect of device installed in a fraction of the bundles on the flow rate in other bundles. Also, the effect on the characteristics of dryout and rewetting and the resulting temperature excursions with and without the device will be given special attention. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2023

Re-Examination of the N 2 O + O Reaction

The reaction of N 2 O with O is a key step in consumption of nitrous oxide in thermal processes. It has two product channels, NO + NO (R2) and N 2 + O 2 (R3). The rate constant for R2 has been measured both in the forward and the reverse direction at elevated temperature and is well established. However, the rate constant for the N 2 + O 2 channel (R3) has been difficult to quantify and has significant error limits. Here, the direct reaction on the triplet surface has a barrier of around 40 kcal mol –1 , and it is too slow for the N 2 + O 2 channel to have any practical significance. Recently, Pham and Lin (2022) suggested an alternative low activation energy reaction path that involves intersystem crossing and reaction on the singlet surface. In the present work, we re-examined a wide range of experiments relevant for the N 2 O + O reaction through kinetic modeling, paying attention to the impact of artifacts such as impurities and surface reactions. Experimental results from shock tubes and batch reactors on the final NO yield in N 2 O decomposition, covering temperatures of 973–2200 K and pressures of 0.013–11.5 atm, support k 3 ~ 0, consistent with the high activation energy for reaction on the triplet surface and a low probability of ISC.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Temperature uncertainty modelling with proxy structural data as geostatistical constraints for well siting: an example applied to Granite Springs Valley, NV, USA

Utilizing existing temperature and structural geology information around Granite Springs Valley, Nevada, we build 3D stochastic temperature models with the aims of evaluating the 3D uncertainty of temperature and choosing between candidate exploration well locations. The data used to support the modelling are measured temperatures and structural proxies from 3D geologic modelling (distance to fault, distance to fault intersections and terminations, Coulomb stress change and dilation tendency), the latter considered ‘secondary’ data. Two stochastic geostatistical techniques are explored for incorporating the structural proxies: cosimulation and local varying mean. With both the cosimulation and local varying mean methods, many equally-likely temperature models (i.e. realizations) are produced, from which temperature probability profiles are calculated at candidate well locations. To aid in choosing between the candidate locations, two quantities summarize the temperature probabilities: V prior and entropy. V prior quantifies the likelihood for economic temperatures at each candidate location, whereas entropy identifies where new information has the most potential to reduce uncertainty. In general, the cosimulation realizations have smoother spatial structure, and extrapolate high temperatures at candidate locations that are located along the direction of the longest spatial correlation, which are down dip from existing temperature logs. The smooth realizations result in tight temperature probability profiles that are easier to interpret, but they have unrealistic temperature reversals in some locations because of the dipping ellipsoid shape created and that the cosimulation technique does not enforce a conductive geothermal gradient as a baseline (i.e. linearly increasing temperature with depth). The local varying mean results produce realizations with more realistic geothermal gradients, with temperatures increasing downward since a depth-temperature relationship is included. However, because they have much noisier spatial nature compared to cosimulation, it is harder to interpret the temperature probability profiles. The different local varying mean results allow the geologist to determine which proxy (e.g. dilation v. distance to fault termination) should be used given the specific geothermal system. In general, V prior from local varying mean results identify locations that are close to high values for the structural proxies: areas with higher probabilities for higher temperatures. The entropy results identify where uncertainty is greatest and therefore new drilling information could be most useful. Though these techniques provide useful information, even when applied to areas of sparse data, our comparison of these two techniques demonstrates the need for new geothermal geostatistics techniques that combine the advantages of these two methods and that are tailored to the spatial uncertainty issues inherent in geothermal exploration.

15 GEOTHERMAL ENERGY↗
DOE OSTI2017

Targeted energy transfer in laminar vortex-induced vibration of a sprung cylinder with a nonlinear dissipative rotator

We computationally investigate the dynamics of a linearly-sprung circular cylinder immersed in an incompressible flow and undergoing transverse vortex-induced vibration (VIV), to which is attached a rotational nonlinear energy sink (NES) consisting of a mass that freely rotates at constant radius about the cylinder axis, and whose motion is restrained by a rotational linear viscous damper. The inertial coupling between the rotational motion of the attached mass and the rectilinear motion of the cylinder is “essentially nonlinear”, which, in conjunction with dissipation, allows for one-way, nearly irreversible targeted energy transfer (TET) from the oscillating cylinder to the nonlinear dissipative attachment. At the intermediate Reynolds number , the NES-equipped sprung cylinder undergoes repetitive cycles of slowly decaying oscillations punctuated by intervals of chaotic instabilities. During the slowly decaying portion of each cycle, the dynamics of the cylinder is regular and, for large enough values of the ratio of the NES mass to the total mass (i.e., NES mass plus cylinder mass), can lead to significant vortex street elongation with partial stabilization of the wake. As approaches zero, no such vortex elongation is observed and the wake patterns appear similar to that for a sprung cylinder with no NES. Here we apply proper orthogonal decomposition (POD) to the velocity flow field during a slowly decaying portion of the solution and show that, in situations where vortex elongation occurs, the NES, though not in direct contact with the surrounding fluid, has a drastic effect on the underlying flow structures, imparting significant and continuous passive redistribution of energy among POD modes. We construct a POD-based reduced-order model for the lift coefficient to characterize energy transactions between the fluid and the cylinder throughout the slowly decaying cycle. We introduce a quantitative signed measure of the work done by the fluid on the cylinder over one quasi-period of the slowly decaying response and find that vortex elongation is associated with a sign change of that measure, indicating that a reversal of the direction of energy transfer, with the cylinder “leaking energy back” to the flow, is responsible for partial stabilization and elongation of the wake. We interpret these findings in terms of the spatial structure and energy distribution of the POD modes, and relate them to the mechanism of transient resonance capture into a slow invariant manifold of the fluid–structure interaction dynamics.

42 ENGINEERING↗
DOE OSTI2013

From Vibrational Spectroscopy to Force Fields and Structures of Saccharides: New Computational Algorithms and Applications

This work was undertaken with the main objective to investigate basic reactions that take place in relatively simple saccharides (mono-saccharides and cellobiose - the building block of cellulose) , in isolation and in cluster with few water molecules or with (gas-phase) clusters of few waters and ionic compounds (salt, isolated ions like H + or OH - ). Within the context of this work, different potentials were investigated; among them, were the PM3 semi empirical potential, DFT/BLYP and a new hybrid potential constructed from MP2 for the harmonic part and from adjusted Hartree-Fock anharmonic interactions (VSCF-PT2). These potentials were evaluated by comparison with experimental data from published sources and from several collaborating groups. The findings show excellent agreement between experiments and predictions with the hybrid VSCF-PT2 potential and very good agreement with predictions obtained from dynamics with dispersion corrected DFT/BLYP potential. Investigation of hydration of cellobiose, was another topic of interest. Guided by a hydration motif demonstrated by our experimental collaborators (team of Prof J.P. Simons), we demonstrated large energetic and structural differences between the two species of cellobiose: cis and trans. The later, which is dominant in solid and liquid phases, is higher in energy in the gas-phase and compared to pure water, it does not disturb as much the network of H bonds. In contrast, the cis species exhibits asymmetric hydration in cluster with up to 25 waters, indicating that it has surfactant properties. Another highlight of this research effort was the successful first time spectrometric and spectroscopic study of a gas-phase protonated sugar derivative (alpha-D-Galactopyranoside) and its interpretation by Ab Initio molecular dynamics (AIMD) simulations. The findings demonstrate the formation of a motif in which a proton bridges between two Oxygen atoms (belonging to OH groups) at the sugar; The vibrational bands involving the shared proton were strongly shifted to lower frequencies ( by about ~ 500 cm -1 for the symmetric mode, in this case). A similar motif was also observed recently by us in protonated cellobiose, indicating that this might be a common mechanism for interaction of a proton with sugars, perhaps similar to the proton wires observed in proteins. The simulations with protonated sugars also shed light on different mechanisms of interaction of a sugar with a proton, including formation of a carboxonium ion, mutarotation events, ring puckering and in the disaccharide cellobiose, the breaking of the glycosidic bond (in both forms of cis and trans). One final highlight to note in this summary, is the finding that in the monosaccharide beta-D-Galactose a hydroxyl ion abstracts a proton (forming water) in a barrierless process at room temperature, but the water remains bound to the sugar backbone, though it migrates around it; actual degradation occurs at ~500 K when water leaves the sugar. However, the study also shows that the water abstraction reaction can be reversed in the presence of 2 additional water molecules complexed with the sugar.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2013

Fluorescence Imaging of Single-Molecule Retention Trajectories in Reversed-Phase Chromatographic Particles

Due to its high specific surface area and chemical stability, porous silica is used as a support structure in numerous applications, including heterogeneous catalysis, biomolecule immobilization, sensors, and liquid chromatography. Reversed-phase liquid chromatography (RPLC), which uses porous silica support particles, has become an indispensable separations tool in quality control, pharmaceutics, and environmental analysis requiring identification of compounds in mixtures. For complex samples, the need for higher resolution separations requires an understanding of the time scale of processes responsible for analyte retention in the stationary phase. In the present work, single-molecule fluorescence imaging is used to observe transport of individual molecules within RPLC porous silica particles. This technique allows direct measurement of intraparticle molecular residence times, intraparticle diffusion rates, and the spatial distribution of molecules within the particle. On the basis of the localization uncertainty and characteristic measured diffusion rates, statistical criteria were developed to resolve the frame-to-frame behavior of molecules into moving and stuck events. The measured diffusion coefficient of moving molecules was used in a Monte Carlo simulation of a random-walk model within the cylindrical geometry of the particle diameter and microscope depth-of-field. Furthermore, the simulated molecular transport is in good agreement with the experimental data, indicating transport of moving molecules in the porous particle is described by a random-walk. Histograms of stuck-molecule event times, locations, and their contributions to intraparticle residence times were also characterized.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2000

Primate translational vestibuloocular reflexes. III. Effects of bilateral labyrinthine electrical stimulation

The effects of functional, reversible ablation and potential recruitment of the most irregular otolith afferents on the dynamics and sensitivity of the translational vestibuloocular reflexes (trVORs) were investigated in rhesus monkeys trained to fixate near and far targets. Translational motion stimuli consisted of either steady-state lateral and fore-aft sinusoidal oscillations or short-lasting transient lateral head displacements. Short-duration (usually <2 s) anodal (inhibitory) and cathodal (excitatory) currents (50-100 microA) were delivered bilaterally during motion. In the presence of anodal labyrinthine stimulation, trVOR sensitivity and its dependence on viewing distance were significantly decreased. In addition, anodal currents significantly increased phase lags. During transient motion, anodal stimulation resulted in significantly lower initial eye acceleration and more sluggish responses. Cathodal currents tended to have opposite effects. The main characteristics of these results were simulated by a simple model where both regularly and irregularly discharging afferents contribute to the trVORs. Anodal labyrinthine currents also were found to decrease eye velocity during long-duration, constant velocity rotations, although results were generally more variable compared with those during translational motion.

NASA Program Biomedical Research and Countermeasur↗
DOE OSTI2022

Quantum model learning agent: characterisation of quantum systems through machine learning

Accurate models of real quantum systems are important for investigating their behaviour, yet are difficult to distil empirically. Here, we report an algorithm—the quantum model learning agent (QMLA)—to reverse engineer Hamiltonian descriptions of a target system. We test the performance of QMLA on a number of simulated experiments, demonstrating several mechanisms for the design of candidate Hamiltonian models and simultaneously entertaining numerous hypotheses about the nature of the physical interactions governing the system under study. QMLA is shown to identify the true model in the majority of instances, when provided with limited a priori information, and control of the experimental setup. Our protocol can explore Ising, Heisenberg and Hubbard families of models in parallel, reliably identifying the family which best describes the system dynamics. We demonstrate QMLA operating on large model spaces by incorporating a genetic algorithm to formulate new hypothetical models. The selection of models whose features propagate to the next generation is based upon an objective function inspired by the Elo rating scheme, typically used to rate competitors in games such as chess and football. In all instances, our protocol finds models that exhibit F 1 score ≥ 0.88 when compared with the true model, and it precisely identifies the true model in 72% of cases, whilst exploring a space of over 250 000 potential models. By testing which interactions actually occur in the target system, QMLA is a viable tool for both the exploration of fundamental physics and the characterisation and calibration of quantum devices.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Hydrogen-Bonded Organic Frameworks: A Rising Class of Porous Molecular Materials

Hydrogen-bonded organic frameworks (HOFs) are a class of porous molecular materials that rely on the assembly of organic building blocks by means of hydrogen-bonding interactions to form two-dimensional (2D) and three-dimensional (3D) crystalline networks. The reversible nature of the hydrogen-bond formation endows HOFs with the attributes of solution processability and simple regeneration. High-quality single crystals of HOFs can be grown easily for unambiguous superstructure determination by single-crystal X-ray diffraction, which is crucial for the elucidation of superstructure–property relationships. During the past decade, considerable progress has been achieved in realizing stable HOFs with permanent porosities by focusing on the design of molecular building blocks in order to introduce rigidity, auxiliary [π···π] interactions, and interpenetration of their frameworks to sustain the extended networks. The applications of HOFs are far-reaching, spanning catalysis, energy, and biomedical products as well as the storage and separation of fine chemicals. In this paper, we, first of all, provide an overview of the chronological development of HOFs, starting from the seminal work by Marsh and Duchamp in 1969 on the crystal superstructure of the hydrogen-bonded networks of trimesic acid. We identify the development of novel hydrogen-bonding motifs such as diaminotriazine (DTA), the introduction of the concept of molecular tectonics, and the establishment of permanent porosity in HOFs as being some of the milestones, which incentivized the current burgeoning research endeavors on developing HOFs as multifunctional materials. This Account is focused primarily on surveying the strategies for constructing porous 3D HOFs based on organic building blocks with peripheral carboxyl groups. These strategies are presented in the following categories: (1) the polycatenation of 2D networks by trigonal building blocks to form global 3D frameworks, (2) the utilization of building blocks with 3D geometries—tetrahedral and trigonal prismatic—that are predisposed to form 3D networks, and (3) the docking by shape-fitting of geometrically labile building blocks. We emphasize how the molecular geometry of the building blocks plays an important role in modulating the superstructures of extended frameworks so as to address specific applications. Recognizing that the in silico design of HOFs is the ultimate goal of researchers in this field, we also discuss the recent advances in superstructure prediction that lead to the formation of porous supramolecular crystals and assess the complications in implementing computational methods for HOFs with complex superstructures. We hope this Account will inspire the development of new supramolecular designs and creative approaches to crystal engineering that aid and abet the assembly of multifunctional HOFs with customizable properties.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Explaining and predicting human behavior and social dynamics in simulated virtual worlds: reproducibility, generalizability, and robustness of causal discovery methods

Ground Truth program was designed to evaluate social science modeling approaches using simulation test beds with ground truth intentionally and systematically embedded to understand and model complex Human Domain systems and their dynamics Lazer et al. (Science 369:1060–1062, 2020). Our multidisciplinary team of data scientists, statisticians, experts in Artificial Intelligence (AI) and visual analytics had a unique role on the program to investigate accuracy, reproducibility, generalizability, and robustness of the state-of-the-art (SOTA) causal structure learning approaches applied to fully observed and sampled simulated data across virtual worlds. In addition, we analyzed the feasibility of using machine learning models to predict future social behavior with and without causal knowledge explicitly embedded. In this paper, we first present our causal modeling approach to discover the causal structure of four virtual worlds produced by the simulation teams—Urban Life, Financial Governance, Disaster and Geopolitical Conflict. Our approach adapts the state-of-the-art causal discovery (including ensemble models), machine learning, data analytics, and visualization techniques to allow a human-machine team to reverse-engineer the true causal relations from sampled and fully observed data. We next present our reproducibility analysis of two research methods team’s performance using a range of causal discovery models applied to both sampled and fully observed data, and analyze their effectiveness and limitations. We further investigate the generalizability and robustness to sampling of the SOTA causal discovery approaches on additional simulated datasets with known ground truth. Our results reveal the limitations of existing causal modeling approaches when applied to large-scale, noisy, high-dimensional data with unobserved variables and unknown relationships between them. We show that the SOTA causal models explored in our experiments are not designed to take advantage from vasts amounts of data and have difficulty recovering ground truth when latent confounders are present; they do not generalize well across simulation scenarios and are not robust to sampling; they are vulnerable to data and modeling assumptions, and therefore, the results are hard to reproduce. Finally, when we outline lessons learned and provide recommendations to improve models for causal discovery and prediction of human social behavior from observational data, we highlight the importance of learning data to knowledge representations or transformations to improve causal discovery and describe the benefit of causal feature selection for predictive and prescriptive modeling.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

SU-D-207A-02: Possible Characterization of the Brain Tumor Vascular Environment by a Novel Strategy of Quantitative Analysis in Dynamic Contrast Enhanced MR Imaging: A Combination of Both Patlak and Logan Analyses

Purpose: The majority of quantitative analyses involving dynamic contrast enhanced (DCE) MRI have been performed to obtain kinetic parameters such as Ktrans and ve. Such analyses are generally performed assuming a “reversible” tissue compartment, where the tracer is assumed to be rapidly equilibrated between the plasma and tissue compartments. However, some tumor vascular environments may be more suited for a “non-reversible” tissue compartment, where, as with FDG PET imaging, the tracer is continuously deposited into the tissue compartment (or the return back to the plasma compartment is very slow in the imaging time scale). Therefore, Patlak and Logan analyses, which represent tools for the “non-reversible” and “reversible” modeling, respectively, were performed to better characterize the brain tumor vascular environment. Methods: A voxel-by-voxel analysis was performed to generate both Patlak and Logan plots in two brain tumor patients, one with grade III astrocytoma and the other with grade IV astrocytoma or glioblastoma. The slopes of plots and the r-square were then obtained by linear fitting and compared for each voxel. Results: The 2-dimensional scatter plots of Logan (Y-axis) vs. Patlak slopes (X-axis) clearly showed increased Logan slopes for glioblastoma (Figure 3A). The scatter plots of goodness-of-fit (Figure 3B) also suggested glioblastoma, relative to grade III astrocytoma, might consist of more voxels that are kinetically Logan-like (i.e. rapidly equilibrated extravascular space and active vascular environment). Therefore, the enhanced Logan-like behavior (and the Logan slope) in glioblastoma may imply an increased fraction of active vascular environment, while the enhanced Patlak-like behavior implies the vascular environment permitting a relatively slower washout of the tracer. Conclusion: Although further verification is required, the combination of Patlak and Logan analyses in DCE MRI may be useful in characterizing the tumor vascular environment, and thus, may have implications in tumor grading and monitoring response to anti-vascular therapy.

60 APPLIED LIFE SCIENCES↗