Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Matched filtering”

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 145 records · Page 8

Distance-Weighted Graph Neural Networks on FPGAs for Real-Time Particle Reconstruction in High Energy Physics

Graph neural networks have been shown to achieve excellent performance for several crucial tasks in particle physics, such as charged particle tracking, jet tagging, and clustering. An important domain for the application of these networks is the FGPA-based first layer of real-time data filtering at the CERN Large Hadron Collider, which has strict latency and resource constraints. We discuss how to design distance-weighted graph networks that can be executed with a latency of less than one μs on an FPGA. To do so, we consider a representative task associated to particle reconstruction and identification in a next-generation calorimeter operating at a particle collider. We use a graph network architecture developed for such purposes, and apply additional simplifications to match the computing constraints of Level-1 trigger systems, including weight quantization. Using the hls4ml library, we convert the compressed models into firmware to be implemented on an FPGA. Performance of the synthesized models is presented both in terms of inference accuracy and resource usage.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Fox Trails

1. This software utilizes python pandas to pull data from P6 databases or XER files. The software transforms the datasets into multiple main tables by joining, filtering, iteratively flattening hierarchical structured data, and pivoting datasets to give simple flat output tables. The activity table includes all of the information related to an activity including activity codes, global, EPS, and project codes, UDFs, and WBS information as separate columns. This includes the code id, code value and sequence number for all levels in hierarchical codes. The resource table is similar to the activity table and includes all of the information related to resources on activities including UPFs and resource codes. The resource time phased table takes the resource information and time phases it for the budget, forecast, late, and actual dates/units/costs that closely matches P6's user interface's values as it implements the resource curve and calendars. The wbs table contains the WBS structure broken out by levels and includes UDFs, codes, and notebook topics. The final P6 data table is the relationships table which simply contains the relationships. 2. When a user updates the tool with data (via giving it P6 project names with database username/password information or XER files) the system creates the data in #1, then creates a networkx graph with the activity data imbedded in the node data and the relationships added as edges. Each edge also has it's float calculated (working time distance between the predecessor and successor) and attached to the edge. Activities are also tagged as a potential start of a path based on their constraints, constraint dates, remaining start date, and activity status. When a user enters an activity ID into the UI, it runs a shortest path calculation on the network graph between each node tagged as potential start to the entered activity id based on the float tagged on the edge. Each path returned by the algorithm contains all of the nodes on the path in order, as well as the total float of the edges that make the path. This data is then collected and returned to the user in the form of a gantt chart with groupings for each path that includes the total float for each group. 3. Similar to 2, if the user passes through a reference dataset each activity set in the path is checked to see if it had a path in the reference dataset, if that path was the primary path between the start and end activities, and what has changed regarding logic and durations. These changes are color coded and summarized before sent to the user to be displayed by the UI for simple discovery. 4. Utilizing the data from #1, the user can submit desired grouping code(s) and filters to the system. The system will then pull the activities, resources, and relationships and create a gantt chart based on the groupings sent and filtered based on the filters sent. 5. The system will produce a gantt chart in a similar method to #4, but allows interactivity with the data. As the user interacts with the gantt chart, the software captures the changes and stores it with the user making the change so that project controls and implement those changes in P6.

Fox, Ben↗

Optimizing structured surfaces for diffractive waveguides

We introduce universal diffractive waveguide designs that can match the performance of conventional dielectric waveguides and achieve various functionalities. Optimized using deep learning, diffractive waveguides can be cascaded to form any desired length and are comprised of transmissive diffractive surfaces that permit the propagation of desired modes with low loss and high mode purity. In addition to guiding the targeted modes through cascaded diffractive units, we also developed various waveguide components and introduced bent diffractive waveguides, rotating the direction of mode propagation, as well as spatial and spectral mode filtering and mode splitting diffractive waveguide designs, and mode-specific polarization control. This framework was experimentally validated in the terahertz spectrum to selectively pass certain spatial modes while rejecting others. Without the need for material dispersion engineering diffractive waveguides can be scaled to operate at different wavelengths, including visible and infrared spectrum, covering potential applications in, e.g., telecommunications, imaging, sensing and spectroscopy.

Applied optics↗

Frequency combs in optically injected terahertz ring quantum cascade lasers

Quantum cascade lasers (QCLs) have emerged as promising candidates for generating chip-scale frequency combs in mid-infrared and terahertz wavelengths. In this work, we demonstrate frequency comb formation in ring terahertz QCLs using the injection of light from a distributed feedback (DFB) laser. The DFB design frequency is chosen to match the modes of the ring cavity (near 3.3 THz), and light from the DFB is injected into the ring QCL via a bus waveguide. By controlling the power and frequency of the optical injection, we show that combs can be selectively formed and controlled in the ring cavity. Numerical modeling suggests that this comb is primarily frequency-modulated in character, with the injection serving to trigger comb formation. We also show that the ring can be used as a filter to control the output of the DFB QCL, potentially being of interest in terahertz photonic integrated circuits. Our work demonstrates that waveguide couplers are a compelling approach for injecting and extracting radiation from ring terahertz combs and offer exciting possibilities for the generation of new comb states in terahertz, such as frequency-modulated waves, solitons, and more.

47 OTHER INSTRUMENTATION↗

Numerical discreteness errors in multispecies cosmological N -body simulations

ABSTRACT We present a detailed analysis of numerical discreteness errors in two-species, gravity-only, cosmological simulations using the density power spectrum as a diagnostic probe. In a simple set-up where both species are initialized with the same total matter transfer function, biased growth of power forms on small scales when the solver force resolution is finer than the mean interparticle separation. The artificial bias is more severe when individual density and velocity transfer functions are applied. In particular, significant large-scale offsets in power are measured between simulations with conventional offset grid initial conditions when compared against converged high-resolution results where the force resolution scale is matched to the interparticle separation. These offsets persist even when the cosmology is chosen so that the two particle species have the same mass, indicating that the error is sourced from discreteness in the total matter field as opposed to unequal particle mass. We further investigate two mitigation strategies to address discreteness errors: the frozen potential method and softened interspecies short-range forces. The former evolves particles under the approximately ‘frozen’ total matter potential in linear theory at early times, while the latter filters cross-species gravitational interactions on small scales in low-density regions. By modelling closer to the continuum limit, both mitigation strategies demonstrate considerable reductions in large-scale power spectrum offsets.

79 ASTRONOMY AND ASTROPHYSICS↗

Influence of Fuel Properties on Gasoline Direct Injection Particulate Matter Emissions Over First 200 s of World-Harmonized Light-Duty Test Procedure Using an Engine Dynamometer and Novel “Virtual Drivetrain” Software

Abstract The influence of fuel properties on particulate matter (PM) emissions from a catalytic gasoline particulate filter (GPF) equipped gasoline direct injection (GDI) engine was investigated using novel “virtual drivetrain” software and an engine mated to an engine dynamometer. The virtual drivetrain software was developed in labview to operate the engine on an engine dynamometer as if it were in a vehicle undergoing a driving cycle. The software uses a physics-based approach to determine vehicle acceleration and speed based on engine load and a programed “shift” schedule to control engine speed. The software uses a control algorithm to modulate engine load and braking to match a calculated vehicle speed with the prescribed speed trace of the driving cycle of choice. The first 200 s of the World-harmonized Light-duty Test Procedure (WLTP) driving cycle was tested using six different fuel formulations of varying volatility, aromaticity, and ethanol concentration. The first 200 s of the WLTP was chosen as the test condition because it is the most problematic section of the driving cycle for controlling PM emissions due to the cold start and cold drive off. It was found that there was a strong correlation between aromaticity of the fuel and the engine-out PM emissions, with the highest emitting fuel producing more than double the mass emissions of the low PM production fuel. However, the post-GPF PM emissions depended greatly on the soot loading state of the GPF. The fuel with the highest engine-out PM emissions produced comparable post-GPF emissions to the lowest PM producing fuel over the driving cycle when the GPF was loaded over three cycles with the respective fuels. These results demonstrate the importance of GPF loading state when aftertreatment systems are used for PM reduction. It also shows that GPF control may be more important than fuel properties, and that regulatory compliance for PM can be achieved with proper GPF control calibration irrespective of fuel type.

Energy & Fuels↗

Ion Exchange Processing of AN-107 Hanford Tank Waste through Crystalline Silicotitanate in a Staged 2- then 3-Column System

The Hanford Site stores an estimated 56 million gallons of mixed radioactive and chemically hazardous waste in large underground tanks. In support of the Direct Feed Low-Activity Waste (DFLAW) Program for expediting Hanford tank waste supernate treatment, laboratory-scale ion exchange processing using prototypic unit operations was conducted on AN-107 tank waste at the Pacific Northwest National Laboratory Radiochemical Processing Laboratory. This report describes the small-scale ion exchange testing with 13.7 L of diluted and filtered supernate from Tank 241-AN-107 (hereafter referred to as AN-107) at 16 °C (62 °F). One of the waste acceptance criteria (WAC) for the Waste Treatment Plant (WTP) Low-Activity Waste Facility is that the waste must contain less than 3.18×10 -5 Ci 137 Cs per mole of Na. For the AN-107 tank waste to meet this criterion, only 0.147% of the influent 137 Cs concentration may be delivered to the WTP; this requires a Cs decontamination factor of 678. Testing with AN-107 matched current Tank Side Cesium Removal (TSCR) facility prototypic operations where a lead-lag configuration was used until the lag column reached the WAC limit, then a polish column was brought online for continued processing in a lead-lag-polish column configuration. Feed was processed at 1.9 bed volumes (BVs) per hour; the flowrate, in terms of contact time with the crystalline silicotitanate (CST) bed, matched the expected flowrate at TSCR. The Cs-decontaminated product was retained for vitrification testing (to be reported separately). The lead column reached 40% Cs breakthrough after processing ~1700 BVs of feed; the 50% Cs breakthrough was extrapolated from the breakthrough data to occur at 1873 BVs. Testing compared to previous AP-101 and AP-107 testing at 16 °C showed ~300 BV increases in volume processed to reach the WAC limit for both lead and lag columns. The increase in capacity was determined to be due to the significantly lower K concentration in the AN-107 compared to the other tank waste matrices. A comparison in breakthrough curves for the three tests indicated slightly slower kinetic behavior in the AN-107, with variations in feed matrices (high organic complexants) likely responsible for the deviation. The Cs effluent from the lag column reached the WAC limit after processing 1097 BVs. Anticipating this breakthrough point, the polish column was preemptively installed around 900 BVs. Cs breakthrough from the lag column began at 500 BVs, reaching 3.06×10 0 µCi/mL, or 2.6 % Cs breakthrough, after processing all 1700 BVs of feed. Table S.1 and Figure S.1 summarize the observed column performance and relevant Cs loading characteristics.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Ion Exchange Processing of AW-105 Hanford Tank Waste through Crystalline Silicotitanate in a Staged 2- then 3-Column System

The Hanford Site stores an estimated 56 million gallons of mixed radioactive and chemically hazardous waste in large underground tanks. In support of the Direct Feed Low-Activity Waste (DFLAW) Program for expediting Hanford tank waste supernate treatment, laboratory-scale ion exchange processing using prototypic unit operations was conducted on AW-105 tank waste at the Pacific Northwest National Laboratory Radiochemical Processing Laboratory. This report describes the small-scale ion exchange testing with 9.2 L of diluted and filtered supernate from Tank 241-AW-105 (hereafter referred to as AW-105) at 16 °C (62 °F). One of the waste acceptance criteria (WAC) for the Waste Treatment Plant (WTP) Low-Activity Waste Facility is that the waste must contain less than 3.18×10 -5 Ci 137 Cs per mole of Na. For the AW-105 tank waste to meet this criterion, only 0.225% of the influent 137 Cs concentration may be delivered to the WTP; this requires a Cs decontamination factor of 445. Testing with AW-105 matched current Tank Side Cesium Removal (TSCR) facility prototypic operations where a lead-lag configuration was used until the lag column reached the WAC limit, then a polish column was brought online for continued processing in a lead-lag-polish column configuration. Feed was processed at 1.9 bed volumes (BVs) per hour; the flowrate, in terms of contact time with the crystalline silicotitanate (CST) bed, matched the expected flowrate at TSCR. The Cs-decontaminated product was retained for vitrification testing (to be reported separately). The lead column reached 83% Cs breakthrough after processing ~1500 BVs of feed; the 50% Cs breakthrough was interpolated from the breakthrough data and occurred at 1041 BVs. Despite the AW-105 having a significantly higher K concentration (0.55 M compared to 0.10 M), testing compared to previous AP-107 ion exchange column testing at 16 °C showed no difference in BVs processed to reach the WAC on the lead column and only an approximate ~20 BV decrease in volume processed to reach the WAC limit on the lag column. The negligible differences in capacity despite the 5x concentration differences in K was determined to be due to the significantly lower NO3 concentration in the AW-105 supernate compared to the AP-107 tank waste matrix. A comparison in breakthrough curves for the two tests also indicated slightly faster kinetic behavior in the AW-105, with the variations in feed matrices (lower NO3 concentration) likely responsible for the deviation. The Cs effluent from the lag column reached the WAC limit after processing 772 BVs. Anticipating this breakthrough point, the polish column was preemptively installed around 675 BVs. Cs breakthrough from the lag column began at 300 BVs, reaching 1.10×10 1 µCi/mL, or 14.13 % Cs breakthrough, after processing all 1500 BVs of feed. The polish column processed nominally 830 BVs and reached 2.10×10 -1 µCi/mL, or 0.27 % Cs breakthrough at the conclusion of the test. Table S.1 and Figure S.1 summarize the observed column performance and relevant Cs loading characteristics.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

High-dimensional Data-driven Energy optimization for Multi-Modal Transit Agencies (HD-EMMA) (Final Technical Report)

Public bus transit services in the U.S. are responsible for at least 19.7 million metric tons of CO 2 emission annually. Electric vehicles (EVs) can have a much lower environmental impact than comparable internal combustion engine vehicles (ICEVs), especially in urban areas. Unfortunately, EVs are also much more expensive than ICEVs. As a result, many public transit agencies can afford only mixed fleets of transit vehicles, consisting of EVs, hybrids (HEVs), and ICEVs. Transit agencies that operate such mixed fleets of vehicles face a challenging optimization problem: these agencies need to decide which vehicles are assigned to serving which transit trips. Since the advantage of EVs over ICEVs varies depending on the route and time of day (e.g., the benefit of EVs is higher in slower traffic with frequent stops and lower on highways), the assignment can have a significant effect on energy use and, hence, environmental impact. Through this project, we have developed reference data about energy collections and constructed a set of machine learning models that can accurately predict the energy consumption for the whole fleet at the level of each trip. We have used these models to develop a scheduling and assignment strategy that can rotate the different vehicle types across the transit agencies’ routes. The optimization algorithm ensures that the vehicles are matched to trips considering weather patterns, expected congestion, and road gradients to minimize the overall energy usage. We list the key observations from our project for other practitioners below. Details are available in the report, and the list of source code and our publications are included in the appendix. 1. We have demonstrated the feasibility of collecting, merging and analyzing large volumes of high-resolution real-world telemetry data from a mixed vehicle fleet. To mitigate the inherent noise of the recorded GPS points, the team developed an algorithm that filters data and maps the points onto a street. The algorithm considers previous and subsequent location measurements and different characteristics of nearby streets to determine how likely the vehicle travels on them. Then, the team segmented the time series into disjoint contiguous samples based on adjacent road segments and repeated the outlier detection and removal. For each data point, the team added features corresponding to elevation changes within the samples, weather features, such as temperature, and traffic data, such as speed ratio between actual speed and free-flow speed. 2. We have developed two forms of machine learning models that be used to understand and analyze the energy operations of a mixed vehicle transit fleet. The micro prediction model provides estimates of instantaneous energy prediction for all types of buses (diesel, hybrid, and electric). Such a model is important in evaluating the energy impacts of real-time bus operation strategies, but it is challenging due to diversified driving cycles of transit buses. The model can help the drivers understand the impact of their driving behaviors and short-term congestions. The macro prediction models estimate average energy consumption across the whole trip considering the features: distance traveled, various road-type features, elevation change, day of the week, time of day, various weather features (temperature, humidity, etc.), and traffic features (speed ratio and jam factor). 3. We have demonstrated that it is possible to transfer the machine learning models we have developed in this project to other teams and cities by using inductive transfer learning. We also showed that the performance of the macro energy prediction models can be improved using a multi-task learning approach where the learning parameters are shared between the models being developed for different vehicle types. The advantage of this approach is improved learning performance as the models can exploit common spatio-temporal and environmental characteristics. 4. Finally, we have developed trip and vehicle assignment and scheduling algorithms that use the energy prediction models and develop a trip to vehicle type (diesel, electric, hybrid) assignment for the whole operation to reduce overall emissions and cost. We have shown through simulations that the proposed algorithms can save $\$$ 48,910 in energy costs and 175 metric tons of CO 2 emission annually for CARTA.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Accurate Identification of Deamidation and Citrullination from Global Shotgun Proteomics Data Using a Dual-Search Delta Score Strategy

While proteins with deamidated/citrullinated amino acids play critical roles in the pathogenesis of many human diseases, identifying these modifications in complex biological samples has been an ongoing challenge. Herein we present a method to accurately identify these modifications from shotgun proteomics data from a deep proteome profiling study of human pancreatic islets obtained by laser capture microdissection. All MS/MS spectra were searched against database by MSGF+ twice with or without a +0.9840 Da mass shift on amino acids asparagine, glutamine, and arginine (NQR) as a dynamic modification. Consequently, for each spectrum the resulting two peptide-to-spectrum matches (PSM) with their respective MSGF+ scores were used for Delta Score calculation. It was observed that all PSMs with positive Delta Score values were clustered with mass errors around 0 ppm, while PSMs with negative Delta Score values were distributed nearly equally within the defined mass error range (20 ppm) for database searching. To estimate false discovery rate (FDR), a “pseudo-decoy” approach was applied whether datasets were searched against a database with a “real modification” mass shift (+0.9840 Da) and a “mock modification” mass shift (+1.0227 Da). FDR was controlled to ~2% with a Delta Score filter greater than zero. Manual inspection of spectra showed that PSMs with positive Delta Score value contained deamidated/citrullinated fragments in their MS/MS spectra. Finally, the results demonstrated that in-situ deamidated/citrullinated peptides can be accurately identified from shotgun tissue proteomics data by the Dual-search Delta Score Strategy.

59 BASIC BIOLOGICAL SCIENCES↗

From models to reality: a systematic review on simulated and measured residential heat pump energy savings

High-performance HVAC solutions are central to residential energy management. A substantial share of these are electric, reversible-cycle systems, with heat pumps representing the largest portion of current and near-term adoption. This review synthesizes peer-reviewed and grey literature on residential space heating and cooling heat pumps. The academic literature is dominated by modeling (73.8%), with limited field measurement (13.1%). Grey literature from United States serve as a supplemental resource providing measured savings. Conversions from electric-resistance heating consistently show the largest site energy reductions, while oil/propane baselines yield moderate savings, and gas baseline scenario often deliver small and region-dependent savings. This study cross-checks the grey literature measured data with simulation data filtered from the ResStock dataset. The comparison indicates a discrepancy between simulations and measured data: simulated site EUIs are typically lower than measured EUIs, but percentage energy savings fall in similar ranges, implying simulations capture directional effects while underestimating energy use. Factors associated with variability and model–measurement differences include system characterization and control representation (e.g., backup heat engagement, thermostat/setpoint strategies, commissioning/installation quality), occupant behavior, weather normalization, metering scope, and envelope characterization. This paper also outlines the proposed methodology for comparing simulation and measured data for heat pumps. It emphasizes the metrics used for comparison and units harmonization, building characteristics matching, and compact metadata are needed for simulations to match measured data. The proposed methodology is expected to improve the credibility of simulated savings as measured evidence grows.

Yu, Lili↗

Cesium Removal from SY-101 Tank Waste Using Crystalline Silicotitanate

The newly operational Tank Side Cesium Removal (TSCR) system removes radioactive cesium-137 ( 137 Cs) and solids from tank waste supernate in the 200 east area of the Hanford site. Efforts to expand the removal capabilities to the 200 west area are underway by a system anticipated to be called the West Area Risk Management (WARM) system. Laboratory-scale ion exchange processing using expected WARM unit operations were conducted to contribute toward Washington River Protection Solutions (WRPS) establishing accurate process flowsheets for the individual feed campaigns planned for the west area supernate pretreatment. This report describes the small-scale ion exchange testing with 8.0 L of filtered supernate from tank 241-SY-101 (referred to as SY-101) at 16 °C (62 °F) to demonstrate processing conditions that would be prototypic of what the WARM system may experience. One of the waste acceptance criteria (WAC) for the Waste Treatment Plant (WTP) Low-Activity Waste Facility is that the waste must contain less than 3.18×10 -5 Ci 137 Cs per mole of Na. For the SY-101 tank waste to meet this criterion, only 0.44% of the influent 137 Cs concentration may be delivered to the WTP; this requires a Cs decontamination factor of 227. Testing with SY-101 matched current TSCR prototypic operations where a lead-lag configuration is used until the lag column reached the WAC limit. Neither the lead nor lag columns reached the WAC, so a polish column was never utilized during this test. Feed was processed at 1.9 bed volumes (BVs) per hour; the flowrate, in terms of contact time with the crystalline silicotitate (CST) bed, matched the current flowrate at TSCR (an expected flowrate to be used at WARM). The lead column only reached 0.09% Cs breakthrough after processing ~1363 BVs of feed. Cesium breakthrough from the lag column was not observed during the entire processing. Table S.1 and Figure S.1 summarize the observed column performance and relevant Cs loading characteristics. Table S.1 and Figure S.1 also summarize the measured SY-101 Cs load performance.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

CR-39 nuclear track detector response to inertial confinement fusion relevant ions

The detection properties of CR-39 were investigated for protons, deuterons, and tritons of various energies. Two models for the relationship between the track diameter and particle energy are presented and demonstrated to match experimental data for all three species. Data demonstrate that CR-39 has 100% efficiency for protons between 1 MeV and 4 MeV, deuterons between 1 MeV and 12.2 MeV, and tritons between 1 MeV and 10 MeV. The true upper bounds for deuterons and tritons exceed what could be measured in data. Simulations were developed to further explore the properties of CR-39 and suggest that the diameter–energy relationship of alpha particles cannot be captured by the conventional c-parameter model. These findings provide confidence in CR-39 track diameter based spectroscopy of all three species and provide invaluable insight for designing filtering for all CR-39 based diagnostics.

47 OTHER INSTRUMENTATION↗

CAHS: Context-Aware Homology Search

Protein homology search is foundational to bioinformatics: it supports annotation transfer, structure/function inference, and evolutionary analysis over rapidly expanding sequence repositories (e.g., UniProtKB). Profile hidden Markov models (pHMMs), as implemented in HMMER, remain the most widely trusted approach because they provide statistically calibrated E-values; however, their gap behavior is fixed once a profile is trained, despite biological evidence that insertion/deletion tolerance varies across flexible loops and intrinsically disordered regions. We present CAHS (Context-Aware Homology Search), a lightweight query-time adapter for pHMM search that incorporates learned and biologically motivated signals without changing HMMER's downstream search pipeline or its calibrated E-value reporting. Given a query sequence, CAHS computes per-residue representations from a protein language model and a disorder predictor, maps these to profile coordinates, and modulates only match-state transition rows (gap-open and gap-extension probabilities) while preserving Plan7 constraints. We comprehensively evaluate CAHS across six structurally diverse protein families and multi-domain architectures against a 570k-sequence target corpus. CAHS expands detection capability, retrieving thousands of additional remote homologs at relaxed thresholds by maintaining alignment quality through flexible regions. For multi-domain proteins, context-aware modulation resolves 94% of fragmented alignments. Crucially, CAHS preserves hit-set invariance at stringent operating points (E<10-10), demonstrating increased statistical confidence without inflating false positives. Furthermore, sharper statistical distinction between homologs and background noise during early filter stages yields up to a 3.87× acceleration in end-to-end wall-clock time on high-performance computing clusters. Overall, CAHS illustrates a practical AI-for-science design pattern: augmenting a trusted probabilistic model with query-specific learned signals to improve interpretable, reproducible inference in data-rich biology.

Bhattaram, Swethasree [Georgia Institute of Techno↗

tapir: A tool for topologies, amplitudes, partial fraction decomposition and input for reductions

The demand for precision predictions in the field of high energy physics has dramatically increased over recent years. Experiments conducted at the LHC, as well as precision measurements at the intensity frontier such as Belle II require equally precise theoretical predictions to make full use of the acquired data. To match the experimental precision, second-, third- and, for certain quantities, even higher-order calculations in perturbative quantum field theory are required. To facilitate such calculations, computer software automating as many steps as possible is required. Yet, each calculation poses different challenges and thus, a high level of configurability is required. In this context we present tapir: a tool for identification, manipulation and minimization of Feynman integral families. It is designed to integrate in toolchains based on the computer algebra system FORM, the use of which is common practice in the field. tapir can be used to reduce the complexity of multi-loop problems with cut-filters, topology mapping, partial fraction decomposition and alike. Program Title:tapir CPC Library link to program files:https://doi.org/10.17632/ptc9t46xyn.1 Developer's repository link:https://gitlab.com/tapir-devs/tapir Licensing provisions: GPLv3 Programming language:python 3, C++ Nature of problem: Multi-loop computations require the automatization of a large number of different tasks related to Feynman integral topologies. Among them are the identification and minimization of integral topologies, partial fraction decomposition of topologies in the case of linearly dependent propagators as well as mapping scalar products of loop momenta to scalar functions. Solution method: The minimization of topologies is performed by comparison of their respective Nickel indices [1], even further minimization utilizes Pak's algorithm [2]. To efficiently map scalar products of loop momenta to scalar functions FORM [3] code is generated. Additional comments including restrictions and unusual features: Minimization based on Pak's algorithm slows down for many lines and scales. A coarser minimization using the Nickel indices, however, is still possible. [1]B. Nickel, D. Meiron, G.A.J. Baker, Compilation of 2-pt and 4-pt graphs for continuous spin model, Report, University of Guelph, 1977.[2]A. Pak, J. Phys. Conf. Ser. 368 (2012) 012049, https://doi.org/10.1088/1742-6596/368/1/012049, arXiv:1111.0868.[3]B. Ruijl, T. Ueda, J. Vermaseren, FORM version 4.2, arXiv:1707.06453, 7 2017.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Combining Remotely Sensed Evapotranspiration and an Agroecosystem Model to Estimate Center‐Pivot Irrigation Water Use at High Spatio‐Temporal Resolution

Abstract Estimating irrigation water use accurately is critical for sustainable irrigation and studying terrestrial water cycle in irrigated croplands. However, irrigation is not monitored in most places, and current estimations of irrigation water use has coarse spatial and/or temporal resolutions. This study aims to estimate irrigation water use at the daily and field scale through the proposed model‐data fusion framework, which is achieved by particle filtering with two configurations (concurrent, CON, and sequential, SEQ) by assimilating satellite‐based evapotranspiration (ET) observations into an advanced agroecosystem model, ecosys . Two types of experiments using synthetic and real ET observations were conducted to study the efficacy of the proposed framework for estimating irrigation water use at the irrigated fields in eastern and western Nebraska, United States. The experiments using synthetic ET observations indicated that, for two major sources of uncertainties of ET difference between observations and model simulations, which are bias and noise, noise had larger impacts on degrading the estimation performance of irrigation water use than bias. For the experiments using real ET observations, monthly and annual estimations of irrigation water use matched well with farmer irrigation records, with Pearson correlation coefficient ( r ) around 0.80 and 0.50, respectively. Although detecting daily irrigation records was very challenging, our method still gave a good performance with RMSE, BIAS, and r around 2.90, 0.03, and 0.4 mm/d, respectively. Our proposed model‐data fusion framework for estimating irrigation water use at high spatio‐temporal resolution could contribute to regional water management, sustainable irrigation, and better tracking terrestrial water cycle.

Environmental Sciences & Ecology↗

Tunneling magnetoresistance in MgO tunnel junctions with Fe-based leads in empirically corrected density functional theory

The minority-spin Fe/MgO interface states are at the Fermi level in density functional theory (DFT), but experimental evidence and GW calculations place them slightly higher in energy. This small shift can strongly influence tunneling magnetoresistance (TMR) in junctions with a thin MgO barrier and its dependence on the concentration of Co in the electrodes. Here, in this study, an empirical potential correction to DFT is introduced to shift the interface states up to match the tunnel spectroscopy data. With this shift, TMR in Fe/MgO/Fe junctions exceeds 800% and 3000% at 3 and 4 monolayers (ML) of MgO, respectively. We further consider the effect of alloying of the Fe electrodes with up to 30% Co or 10% V, treating them in the coherent potential approximation (CPA). Alloying with Co broadens the interface states and brings a large incoherent minority-spin spectral weight to the Fermi level. Alloying with V brings the minority-spin resonant states close to the Fermi level. However, in both cases the minority-spin spectral weight at the Fermi level resides primarily at the periphery of the Brillouin zone, which is favorable for spin filtering. Using convolutions of k II -resolved barrier densities of states calculated in CPA, it is found that TMR is strongly reduced by alloying with Co or V but still remains above 500% at 4 ML of MgO up to 30% of Co or 5% V. At 5 ML, the TMR increases above 1000% in all systems considered. However, while TMR declines sharply with increasing bias up to 0.2 eV in the tunnel junctions with pure Fe leads, it remains almost constant up to 0.5 eV if leads are alloyed with Co.

74 ATOMIC AND MOLECULAR PHYSICS↗

Energy Exascale Earth System Model v2.0

First release of version 2 of the Energy Exascale Earth System Model. The atmosphere component remains EAM. Major changes since version 1 include: all column-physics parameterizations are computed on a separate grid that has approximately half the number of points of the dynamics grid, a new nonhydrostatic dynamical core (running in hydrostatic mode) with semi-Lagrangian tracer transport, CLUBB updated from v1 to v2, a new convective trigger (dCAPE/ULL) based on the dynamic Convective Available Potential Energy (CAPE) (dCAPE) and the Unrestricted Launch Level (ULL) concepts is used in ZM. minimum cloud droplet number changed, gravity wave drag energy conservation fixed and new tunings used, dust emission size distribution changed to emit more coarse dust particles The land component is still ELM. Major changes since version 1 include: using SNICAR-AD for radiation in snow to match the sea-ice model and fixing bugs in snow compaction and water state calculation. The ocean component remains MPAS-ocean. Major change since version 1 include: Redi isopycnal mixing has been updated, tested, and tuned in combination with the Gent-McWilliams parameterization, a sign error was fixed in the 3rd-order flux routines, the mesh used in low-resolution coupled cases was modified and the time steps adjusted, new regionally refined meshes were created, one focused on North America and another on the Southern Ocean, including ice shelf cavities. The sea-ice component remains MPAS-seaice. Major changes since version 1 include: A new heat- and freshwater-conserving coupling of frazil ice, turning off of SSH filtering, addition of SNICAR-AD and snow grain aging. The land-ice component remains MPAS-Albany-landIce (MALI) and is a static ice sheet. There is more out-of-the-box support for cryosphere configurations including the new regionally refined configuration around Antarctica with ice shelf cavities. The river model is MOSART. The half-degree river mesh was redone so it no longer treats Black and Caspian seas as ocean. The coupler remains cpl7/MCT. Major changes since version 1 include handling of ice shelf melt fluxes (heat / freshwater exchange with the ocean), and data icebergs. All components allow regional refinement of their meshes and two separate example of refinement, one in and around North America and one around Antarctica and the Southern Ocean, are provided.

E3SM Project, DOE↗