In-silico screening of novel biofuel molecules based on chemical and fuel properties is a critical first step in the biofuel evaluation process due to the significant volumes of samples required for experimental testing, the destructive nature of engine tests, and the costs associated with bench-scale synthesis of novel fuels. Predictive models are limited by training sets of few existing measurements, often containing similar classes of molecules that represent just a subset of the potential molecular fuel space. Software tools can be used to generate every possible molecular descriptor for use as input features, but most of these features are largely irrelevant and training models on datasets with higher dimensionality than size tends to yield poor predictive performance. Feature selection has been shown to improve machine learning models, but correlation-based feature selection fails to provide scientific insight into the underlying mechanisms that determine structure–property relationships. The implementation of causal discovery in feature selection could potentially inform the biofuel design process while also improving model prediction accuracy and robustness to new data. In this study, we investigate the benefits causal-based feature selection might have on both model performance and identification of key molecular substructures. We found that causal-based feature selection performed on par with alternative filtration methods, and that a structural causal model provides valuable scientific insights into the relationships between molecular substructures and fuel properties.
Image-based simulation, the use of 3D images to calculate physical quantities, relies on image segmentation for geometry creation. However, this process introduces image segmentation uncertainty because different segmentation tools (both manual and machine-learning-based) will each produce a unique and valid segmentation. First, we demonstrate that these variations propagate into the physics simulations, compromising the resulting physics quantities. Second, we propose a general framework for rapidly quantifying segmentation uncertainty. Through the creation and sampling of segmentation uncertainty probability maps, we systematically and objectively create uncertainty distributions of the physics quantities. We show that physics quantity uncertainty distributions can follow a Normal distribution, but, in more complicated physics simulations, the resulting uncertainty distribution can be surprisingly nontrivial. We establish that bounding segmentation uncertainty can fail in these nontrivial situations. While our work does not eliminate segmentation uncertainty, it improves simulation credibility by making visible the previously unrecognized segmentation uncertainty plaguing image-based simulation.
Journal of Water Resources Planning and Management
Haxton, Terranna; Klise, Katherine A.; Laky, Daniel; Murray, Regan; Laird, Carl D.; Burkhardt, Jonathan B.
Drinking water systems commonly use manual or grab sampling to monitor water quality, identify or confirm issues, and verify that corrective or emergency response actions have been effective. In this paper, the effectiveness of regulatory sampling locations for emergency response is explored. An optimization formulation based on the literature was used to identify manual sampling locations to maximize overall nodal coverage of the system. Results showed that sampling locations could be effective in confirming incidents for which they were not designed. When evaluating sampling locations optimized for emergency response against regulatory scenarios, the average performance was reduced by 3%-4%, while using optimized regulatory sampling locations for emergency response reduced performance by 7%-10%. Secondary constraints were also included in the formulation to ensure geographical and water age diversity with minimal impact on the performance. This work highlighted that regulatory sampling locations provide value in responding to an emergency for these networks.
Shu, Yi; Galles, Daniel; Tertuliano, Ottman A.; Mcwilliams, Brandon A.; Yang, Nancy; Cai, Wei; Lew, Adrian J.
The study of microstructure evolution in additive manufacturing of metals would be aided by knowing the thermal history. Since temperature measurements beneath the surface are difficult, estimates are obtained from computational thermo-mechanical models calibrated against traces left in the sample revealed after etching, such as the trace of the melt pool boundary. Here we examine the question of how reliable thermal histories computed from a model that reproduces the melt pool trace are. To this end, we perform experiments in which one of two different laser beams moves with constant velocity and power over a substrate of 17-4PH SS or Ti-6Al-4V, with low enough power to avoid generating a keyhole. We find that thermal histories appear to be reliably computed provided that (a) the power density distribution of the laser beam over the substrate is well characterized, and (b) convective heat transport effects are accounted for. Poor control of the laser beam leads to potentially multiple three-dimensional melt pool shapes compatible with the melt pool trace, and therefore to multiple potential thermal histories. Ignoring convective effects leads to results that are inconsistent with experiments, even for the mild melt pools here.
Photon detection at microwave frequency is of great interest due to its application in quantum computation information science and technology. Herein are results from studying microwave response in a topological superconducting quantum interference device (SQUID) realized in Dirac semimetal Cd3As2. The temperature dependence and microwave power dependence of the SQUID junction resistance are studied, from which we obtain an effective temperature at each microwave power level. It is observed the effective temperature increases with the microwave power. This observation of large microwave response may pave the way for single photon detection at the microwave frequency in topological quantum materials.
Point defects in SiC are an attractive platform for quantum information and sensing applications because they provide relatively long spin coherence times, optical spin initialization, and spin-dependent fluorescence readout in a fabrication-friendly semiconductor. The ability to precisely place these defects at the optimal location in a host material with nano-scale accuracy is desirable for integration of these quantum systems with traditional electronic and photonic structures. Here, we demonstrate the precise spatial patterning of arrays of silicon vacancy (VSi) emitters in an epitaxial 4H-SiC (0001) layer through mask-less focused ion beam implantation of Li+. We characterize these arrays with high-resolution scanning confocal fluorescence microscopy on the Si-face, observing sharp emission lines primarily coming from the V1 ′ zero-phonon line (ZPL). The implantation dose is varied over 3 orders of magnitude, leading to VSi densities from a few per implantation spot to thousands per spot, with a linear dependence between ZPL emission and implantation dose. Optically-detected magnetic resonance (ODMR) is also performed, confirming the presence of V2 VSi. Our investigation reveals scalable and reproducible defect generation.
We demonstrate the generation of a cold-atom ensemble within a sub-millimeter diameter hole in a transparent membrane, a so-called “membrane MOT”. With a sub-Doppler cooling process, the atoms trapped by the membrane MOT are cooled down to 10 μ K. The atom number inside the unbridged/bridged membrane hole is about 10 4 to 10 5, and the 1 / e2-diameter of the MOT cloud is about 180 μ m for a 400 μ m-diameter membrane hole. Such a membrane device can, in principle, efficiently load cold atoms into the evanescent-field optical trap generated by the suspended membrane waveguide for strong atom-light interaction and provide the capability of sufficient heat dissipation at the waveguide. This represents a key step toward the photonic atom trap integrated platform (ATIP).
Kim, Sung E.; Seo, Yongwon; Hwang, Junshik; Yoon, Hongkyu; Lee, Jonghyun
Stochastic network modeling is often limited by high computational costs to generate a large number of networks enough for meaningful statistical evaluation. In this study, Deep Convolutional Generative Adversarial Networks (DCGANs) were applied to quickly reproduce drainage networks from the already generated network samples without repetitive long modeling of the stochastic network model, Gibb’s model. In particular, we developed a novel connectivity-informed method that converts the drainage network images to the directional information of flow on each node of the drainage network, and then transforms it into multiple binary layers where the connectivity constraints between nodes in the drainage network are stored. DCGANs trained with three different types of training samples were compared; (1) original drainage network images, (2) their corresponding directional information only, and (3) the connectivity-informed directional information. A comparison of generated images demonstrated that the novel connectivity-informed method outperformed the other two methods by training DCGANs more effectively and better reproducing accurate drainage networks due to its compact representation of the network complexity and connectivity. This work highlights that DCGANs can be applicable for high contrast images common in earth and material sciences where the network, fractures, and other high contrast features are important.
The solid-state joining of oxide-dispersion-strengthened (ODS) austenitic steels was achieved using a pulsed electric current joining (PECJ) process. Microstructures of the austenitic grain structures and oxide dispersions in the joint areas were characterized using electron microscopy. Negligible grain growth was observed in austenitic grain structures, while slight coarsening of oxide dispersions occurred at a short holding time. The mechanisms of the PECJ process may involve three steps that occur simultaneously, including the sintering of mechanical alloying powders in the bonding layer, formation of oxide dispersions, and bonding of the mechanical alloying powders with the base alloy. The high hardness and irradiation resistance of ODS alloys were retained in the joint areas. This research revealed the fundamental mechanisms during the PECJ process, which is beneficial for its potential applications during the advanced manufacturing of ODS alloys.
Conventional engineering methods oftentimes have challenges in the quantification of crack nucleation processes from manufacturing defects that are relevant for engineering component lifing. We present peridynamic simulation framework for the description of the crack nucleation process from cluster of non-metallic inclusions in aluminum alloys. Our non-local simulation framework characterizes crack nucleation process as multiple micro-crack nucleation events from individual inclusions, and eventually one micro-crack dominates. We define individual stages of the crack nucleation process, i.e., nucleation, micro-crack, technical, and crack initiation, that allow a quantification and meta model development of both the individual stages and the entire crack nucleation process.
Organophosphorus hydrolase (OPH) is a metalloenzyme that can hydrolyze organophosphorus agents resulting in products that are generally of reduced toxicity. The best OPH substrate found to date is diethyl p-nitrophenyl phosphate (paraoxon). Most structural and kinetic studies assume that the binding orientation of paraoxon is identical to that of diethyl 4-methylbenzylphosphonate, which is the only substrate analog co-crystallized with OPH. In the current work, we used a combined docking and molecular dynamics (MD) approach to predict the likely binding mode of paraoxon. Then, we used the predicted binding mode to run MD simulations on the wild type (WT) OPH complexed with paraoxon, and OPH mutants complexed with paraoxon. Additionally, we identified three hot-spot residues (D253, H254, and I255) involved in the stability of the OPH active site. We then experimentally assayed single and double mutants involving these residues for paraoxon binding affinity. The binding free energy calculations and the experimental kinetics of the reactions between each OPH mutant and paraoxon show that mutated forms D253E, D253E-H254R, and D253E-I255G exhibit enhanced substrate binding affinity over WT OPH. Interestingly, our experimental results show that the substrate binding affinity of the double mutant D253E-H254R increased by 19-fold compared to WT OPH.
Eye tracking is a useful tool for studying human cognition, both in the laboratory and in real-world applications. However, there are cases in which eye tracking is not possible, such as in high-security environments where recording devices cannot be introduced. After facing this challenge in our own work, we sought to test the effectiveness of using artificial foveation as an alternative to eye tracking for studying visual search performance. Two groups of participants completed the same list comparison task, which was a computer-based task designed to mimic an inventory verification process that is commonly performed by international nuclear safeguards inspectors. We manipulated the way in which the items on the inventory list were ordered and color coded. For the eye tracking group, an eye tracker was used to assess the order in which participants viewed the items and the number of fixations per trial in each list condition. For the artificial foveation group, the items were covered with a blurry mask except when participants moused over them. We tracked the order in which participants viewed the items by moving their mouse and the number of items viewed per trial in each list condition. We observed the same overall pattern of performance for the various list display conditions, regardless of the method. However, participants were much slower to complete the task when using artificial foveation and had more variability in their accuracy. Our results indicate that the artificial foveation method can reveal the same pattern of differences across conditions as eye tracking, but it can also impact participants’ task performance.
Cerjan, Alexander W.; Jorg, Christina; Vaidya, Sachin; Augustine, Shyam; Benalcazar, Wladimir A.; Hsu, Chia W.; Von Freymann, Georg; Rechtsman, Mikael C.
In the past decade, symmetry-protected bound states in the continuum (BICs) have proven to be an important design principle for creating and enhancing devices reliant upon states with high-quality (Q) factors, such as sensors, lasers, and those for harmonic generation. However, as we show, current implementations of symmetry-protected BICs in photonic crystal slabs can only be found at the center of the Brillouin zone and below the Bragg diffraction limit, which fundamentally restricts their use to single-frequency applications. By microprinting a three-dimensional (3D) photonic crystal structure using two-photon polymerization, we demonstrate that this limitation can be overcome by altering the radiative environment surrounding the slab to be a 3D photonic crystal. This allows for the protection of a line of BICs by embedding it in a symmetry bandgap of the crystal. This concept substantially expands the design freedom available for developing next-generation devices with high-Q states.
Big Data in the area of Remote Sensing has been growing rapidly. Remote sensors are used in surveillance, security, traffic, environmental monitoring, and autonomous sensing. Real-time detection of small moving targets using a remote sensor is an ongoing, challenging problem. Since the object is located far away from the sensor, the object often appears too small. The object’s signal-to-noise-ratio (SNR) is often very low. Occurrences such as camera motion, moving backgrounds (e.g., rustling leaves), low contrast and resolution of foreground objects makes it difficult to segment out the targeted moving objects of interest. Due to the limited appearance of the target, it is tough to obtain the target’s characteristics such as its shape and texture. Without these characteristics, filtering out false detections can be a difficult task. Detecting these targets, would often require the detector to operate under a low detection threshold. However, lowering the detection threshold could lead to an increase of false alarms. In this paper, the author will introduce a new method that improves the probability to detect low SNR objects, while decreasing the number of false alarms as compared to using the traditional baseline detection technique.
Kirby, James; Geiselman, Gina M.; Yaegashi, Junko; Kim, Joonhoon; Zhuang, Xun; Tran-Gyamfi, Mary; Prahl, Jan P.; Sundstrom, Eric R.; Gao, Yuqian; Munoz, Nathalie; Burnum-Johnson, Kristin E.; Benites, Veronica T.; Baidoo, Edward E.K.; Fuhrmann, Anna; Seibel, Katharina; Webb-Robertson, Bobbie J.M.; Zucker, Jeremy; Nicora, Carrie D.; Tanjore, Deepti; Magnuson, Jon K.; Skerker, Jeffrey M.; Gladden, John M.
Background: Mitigation of climate change requires that new routes for the production of fuels and chemicals be as oil-independent as possible. The microbial conversion of lignocellulosic feedstocks into terpene-based biofuels and bioproducts represents one such route. This work builds upon previous demonstrations that the single-celled carotenogenic basidiomycete, Rhodosporidium toruloides, is a promising host for the production of terpenes from lignocellulosic hydrolysates. Results: This study focuses on the optimization of production of the monoterpene 1,8-cineole and the sesquiterpene α-bisabolene in R. toruloides. The α-bisabolene titer attained in R. toruloides was found to be proportional to the copy number of the bisabolene synthase (BIS) expression cassette, which in turn influenced the expression level of several native mevalonate pathway genes. The addition of more copies of BIS under a stronger promoter resulted in production of α-bisabolene at 2.2 g/L from lignocellulosic hydrolysate in a 2-L fermenter. Production of 1,8-cineole was found to be limited by availability of the precursor geranylgeranyl pyrophosphate (GPP) and expression of an appropriate GPP synthase increased the monoterpene titer fourfold to 143 mg/L at bench scale. Targeted mevalonate pathway metabolite analysis suggested that 3-hydroxy-3-methyl-glutaryl-coenzyme A reductase (HMGR), mevalonate kinase (MK) and phosphomevalonate kinase (PMK) may be pathway bottlenecks are were therefore selected as targets for overexpression. Expression of HMGR, MK, and PMK orthologs and growth in an optimized lignocellulosic hydrolysate medium increased the 1,8-cineole titer an additional tenfold to 1.4 g/L. Expression of the same mevalonate pathway genes did not have as large an impact on α-bisabolene production, although the final titer was higher at 2.6 g/L. Furthermore, mevalonate pathway intermediates accumulated in the mevalonate-engineered strains, suggesting room for further improvement. Conclusions: This work brings R. toruloides closer to being able to make industrially relevant quantities of terpene from lignocellulosic biomass.
Curwen, Christopher A.; Addamane, Sadhvikas J.; Reno, John L.; Shahili, Mohammad; Kawamura, Jonathan H.; Briggs, Ryan M.; Karasik, Boris S.; Williams, Benjamin S.
We compare the performance of 10 and 5 μm thick metal-metal waveguide terahertz quantum-cascade laser ridges operating around 2.7 THz and based on a 4-well phonon depopulation active region design. Thanks to reduced heat dissipation and lower thermal resistance, the 5 μm thick material shows an 18 K increase in continuous wave operating temperature compared to the 10 μm material, despite a lower maximum pulsed-mode operating temperature and a larger input power density. A maximum continuous wave operating temperature of 129 K is achieved using the 5 μm thick material and a 15 μm wide ridge waveguide, which lased up to 155 K in the pulsed mode. The use of thin active regions is likely to become increasingly important to address the increasing input power density of emerging 2- and 3-well active region designs that show the highest pulsed operating temperatures.
Sandia provided technical assistance to Energy Cloud to support the assessment of Energy Cloud’s technology for air filtration of the SARS-CoV-2 virus. The Energy Cloud system was developed to filter the SARS-CoV-2 virus from commercial HVAC systems, ducts, and enclosed rooms without generating ozone. The project included a 3-week experimental series with aerosolized SARS-CoV-2 virus at the University of Southern California’s biosafety level 3 (BSL-3) laboratory. USC provided the facility as well as experimentally derived measurements and assessment of the Energy Cloud system efficacy for SARSCoV-2 filtration from air handling systems.
This report discusses the progress on the collaboration between Sandia National Laboratories (Sandia) and Japan Atomic Energy Agency (JAEA) on the sodium fire research in fiscal year (FY) 2021 and is a continuation of the FY2020 progress report. In this report, we only report the changes made to the current sodium pool fire model in MELCOR. We modified and corrected many control functions to enhance the current model from the FY2020 report. This year’s enhancements relate to better agreement of the suspended aerosol measurement from JAEA’s F7 series tests. Staff from Sandia and JAEA conducted the validation studies of the sodium pool fire model in MELCOR. To validate this pool fire model with the latest enhancement, JAEA sodium pool fire experiments (F7-1 and F7-2) were used. The results of the calculation, including the code-to-code comparisons are discussed as well as suggestions for further model improvement. Finally, recommendations are made for new MELCOR simulations for FY2022.
A 3 foot x 3 foot x 3 foot aluminum solar collector was manufactured using computer numerical control. The interior of the device included six triangular dimpled fins for enhanced heat transfer. The interior vertical wall on the south side was also dimpled. The solar collector working fluid was based on water, and the collector consisted solely of passive heat transfer mechanisms (no moving parts), making it ideal for off-the-grid and rural applications. Two types of heat transfer experiments were conducted. One experiment had external flat heaters attached on the top and the front side, while the other four sides were insulated. Except for the bottom surface, the second experiment had all its exterior surfaces sprayed with black solar paint to collect as much solar heat as possible. Temperature data as a function of time was collected using 14 thermocouples spread strategically throughout the solar collector. In addition, computational fluid dynamics (CFD) simulations were conducted using the dynamic Smagorinsky large eddy simulation turbulence model. The first simulation considered that both the top and front surfaces were exposed to a fixed temperature of 313.7 K (105 °F), while the remaining four surfaces were insulated. For the second simulation, all conditions were the same, except that the temperature for both heated surfaces was raised to 350 K (170.3 °F). The two temperatures are expected to bound the solar collector operational temperature during the late- Spring, Summer, and early-Fall months. The solar collector design, experimental data, CFD output, and a discussion of five manufacturing approaches and costs are documented in this report.
Local-to-nonlocal (LtN) coupling refers to a class of methods aimed at combining nonlocal and local modeling descriptions of a given system into a unified coupled representation. This allows to consolidate the accuracy of nonlocal models with the computational expediency of their local counterparts, while often simultaneously removing nonlocal modeling issues such as surface effects. The number and variety of proposed LtN coupling approaches have significantly grown in recent years, yet the field of LtN coupling continues to grow and still has open challenges. This review provides an overview of the state of the art of LtN coupling in the context of nonlocal diffusion and nonlocal mechanics, specifically peridynamics. Furthermore, we present a classification of LtN coupling methods and discuss common features and challenges. The goal of this review is not to provide a preferred way to address LtN coupling but to present a broad perspective of the field, which would serve as guidance for practitioners in the selection of appropriate LtN coupling methods based on the characteristics and needs of the problem under consideration.
