Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “difference scaling”

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,189 records · Page 66
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Genome-scale phylogeny and contrasting modes of genome evolution in the fungal phylum Ascomycota

Ascomycota, the largest and most well-studied phylum of fungi, contains three subphyla: Saccharomycotina (budding yeasts), Pezizomycotina (filamentous fungi), and Taphrinomycotina (fission yeasts). Despite its importance, we lack a comprehensive genome-scale phylogeny or understanding of the similarities and differences in the mode of genome evolution within this phylum. By examining 1107 genomes from Saccharomycotina (332), Pezizomycotina (761), and Taphrinomycotina (14) species, we inferred a robust genome-wide phylogeny that resolves several contentious relationships and estimated that the Ascomycota last common ancestor likely originated in the Ediacaran period. Comparisons of genomic properties revealed that Saccharomycotina and Pezizomycotina differ greatly in their genome properties and enabled inference of the direction of evolutionary change. The Saccharomycotina typically have smaller genomes, lower guanine-cytosine contents, lower numbers of genes, and higher rates of molecular sequence evolution compared with Pezizomycotina. These results provide a robust evolutionary framework for understanding the diversity and ecological lifestyles of the largest fungal phylum.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Sparse Approximate Multifrontal Factorization with Butterfly Compression for High-Frequency Wave Equations

In this work, we present a fast and approximate multifrontal solver for large-scale sparse linear systems arising from finite-difference, finite-volume or finite-element discretization of high-frequency wave equations. The proposed solver leverages the butterfly algorithm and its hierarchical matrix extension for compressing and factorizing large frontal matrices via graph-distance guided entry evaluation or randomized matrix-vector multiplication-based schemes. Complexity analysis and numerical experiments demonstrate $\mathcal{O}(N\log^2 N)$ computation and $\mathcal{O}(N)$ memory complexity when applied to an $N\times N$ sparse system arising from 3D high-frequency Helmholtz and Maxwell problems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Scalable Predictive And Risk Technologies

The research involves developing scalable technologies for risk-informed predictive analytics to achieve condition-based monitoring and maintenance strategies to reduce overall maintenance costs. The research utilizes data (real-time data, periodic data, and institutional knowledge) related to a particular plant asset from a specific nuclear plant site to develop technologies to scale risk-informed predictive analytic algorithms across different plant assets at the plant site and across the nuclear fleet. The developed algorithms and codes are used to optimize the maintenance strategy and estimate/forecast generation costs based on the state of health of the plant asset. Developed codes specifically include 1. Parameter estimation using plant operation data 2. Federated and Transfer learning model 3. Feature group based Multi-kernel SVM 4. Three state markov model

Manjunatha, KoushikAraseethota↗
DOE OSTI2020

Cellular Response to Proton Irradiation: A Simulation Study with TOPAS-nBio

The cellular response to ionizing radiation continues to be of significant research interest in cancer radiotherapy, and DNA is recognized as the critical target for most of the biologic effects of radiation. Incident particles can cause initial DNA damages through physical and chemical interactions within a short time scale. Initial DNA damages can undergo repair via different pathways available at different stages of the cell cycle. The misrepair of DNA damage results in genomic rearrangement and causes mutations and chromosome aberrations, which are drivers of cell death. This work presents an integrated study of simulating cell response after proton irradiation with energies of 0.5–500 MeV (LET of 60–0.2 keV/µm). A model of a whole nucleus with fractal DNA geometry was implemented in TOPAS-nBio for initial DNA damage simulations. The default physics and chemistry models in TOPAS-nBio were used to describe interactions of primary particles, secondary particles, and radiolysis products within the nucleus. The initial DNA double-strand break (DSB) yield was found to increase from 6.5 DSB/Gy/Gbp at low-linear energy transfer (LET) of 0.2 keV/µm to 21.2 DSB/Gy/Gbp at high LET of 60 keV/µm. A mechanistic repair model was applied to predict the characteristics of DNA damage repair and dose response of chromosome aberrations. It was found that more than 95% of the DSBs are repaired within the first 24 h and the misrepaired DSB fraction increases rapidly with LET and reaches 15.8% at 60 keV/µm with an estimated chromosome aberration detection threshold of 3 Mbp. The dicentric and acentric fragment yields and the dose response of micronuclei formation after proton irradiation were calculated and compared with experimental results.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2020

Interactions Between Microturbulence and Magnetohydrodynamic Instabilities in Magnetized Plasmas (Final Technical Report)

The overall goal of this project was to develop a better understanding between small-scale microturbulence and large-scale magnetohydrodynamic (MHD) instabilities, particularly magnetic islands and tearing modes, in fusion relevant magnetically confined plasmas. This project advanced understanding through numerical simulations of different model equations. Questions addressed included how the magnitudes, spatial structures, and scalings of turbulent forcings and drags impact island dynamics, as well as what the "region of influence" of a tearing mode on the turbulence dynamics. A key finding of the project is that the microturbulence (and associated transport) exhibits spatial asymmetries for sufficiently large islands. The results have improved our fundamental understanding of plasma dynamics, motivated new experimental investigations, and will help inform future fusion predictive modeling efforts.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Heating and Acceleration of Solar Wind by Parametric Decay Instability Mediated Turbulence

One of the main goals of this project to understand the origin of density fluctuations in the solar winds, as revealed by the NASA’s Flagship satellite Parker Solar Probe (PSP) measurements, which motivates us to simulate and analyze magnetohydrodynamic (MHD) turbulence in the compressible regime. To achieve this goal, we have developed a set of sophisticated turbulence driving methods to study how MHD turbulence will develop under different time-correlated fluctuations at large-scales. Another key advance is that we have developed a new 3+1 (three spatial plus time domain) 4D FFT analysis of simulation data-cubes so that we can perform the spatio-temporal spectrum and mode decomposition studies on MHD turbulence.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Pumped Storage Hydropower Augmented with Pressurized Air: The Ground-Level Integrated Diverse Energy Storage (GLIDES) System — GLIDES System Configurations and Use Cases

Energy storage is essential for cost-effective integration of variable renewable energy sources to support a low-carbon grid. It is also a key enabler of a modern grid infrastructure for demand management. However, several main challenges remain for different kind of energy storage technologies in grid scale deployment. Currently, the largest source of utility-scale storage and long-duration storage in the US is pumped storage hydropower (PSH). Prospect of growth in conventional PSH faces challenges that have limited its deployment over the last three decades, including high capital costs and long deployment timelines. Batteries have high energy densities and are the primary technology of choice for small-scale energy storage. Compressed air energy storage (CAES) is another large-scale energy storage technology, but there are few plants deployed worldwide. They suffer from their low round trip efficiency (RTE) due to the use of high-pressure air compressors. To address some of the challenges associated with these various storage technologies, the Ground-Level Integrated Diverse Energy Storage (GLIDES) is a modular PSH technology that was invented in 2015 at Oak Ridge National Laboratory. It utilizes gas compression to store electric energy. GLIDES stores energy by compressing gas using a liquid piston in high-pressure vessels. In doing so the vessels act as the upper reservoir in conventional PSH. Initially, the vessels are filled with gas to a prescribed pressure. To store energy, GLIDES uses a hydraulic piston pump to pump water into the pressurized vessels. As the water volume increases inside the vessels, water acts as a hydraulic piston compressing the gas on top of it. This process can be thought of as pumping water from the lower reservoir to the higher reservoir in PSH, increasing the water head. To dispatch the stored energy, the high-head water in the vessel is discharge through a high head Pelton hydraulic turbine that is connected to an electric generator. Employing high-pressure vessels enables GLIDES to reach water heads ~10-80 times higher than conventional PSH, achieving ~40 times higher energy densities, and overcomes the geographic limitation of conventional PSH. Although its energy density is much lower than that of batteries, GLIDES holds the potential advantages of having long service life, ease of system integration and being less hazardous over batteries. GLIDES prospective scalability could make it suitable for wide range of applications from behind the meter storage in buildings to grid-scale storage. It also makes it suitable for installations in densely populated urban areas where energy storage is most needed and real estate is limited. Over the last 5 years, work has focused on increasing GLIDES’ energy density, decreasing its initial capital cost of the system, and increasing its revenue potential. Several designs were developed and prototyped to verify and demonstrate the improvement in energy density. The latest prototype achieved energy density of 1.21 kWh/m 3 . Our analysis showed that it could achieve up to 1.7 kWh/m 3 with a mixture of air and carbon dioxide as the gas being compressed.

13 HYDRO ENERGY↗
DOE OSTI2022

Nek5000/RS Performance on Advanced GPU Architectures

We demonstrate NekRS performance results on various advanced GPU architectures. NekRS is a GPU-accelerated version of Nek5000 that is targeting high performance on forthcoming exascale platforms. It is being developed in DOE’s Center of Efficient Exascale Discretizations (CEED), which is one of the co-design centers under the Exascale Computing Project (ECP). In this report, we consider Frontier, Crusher, Spock, Polaris, Perlmutter, ThetaGPU and Summit. Simulations are performed using ExaSMR’s 17x17 rod-bundle geometries with different problem sizes. The report focuses on strong scaling performance and analysis. Many of the results shown in this report are the outcome from participation in the ALCF GPU Hackathon on Polaris, which was held on 7/19/22, and 7/26–7/28/22. Members of the Nek5000/RS team included Misun Min, Yu-Hsiang Lan, Paul Fischer and Thilina Rathnayake. Mentors were Kris Rowe (ALCF) and Peng Wang (NVIDIA). Also presented in this report are results on Frontier, which were obtained in collaboration with John Holmen at OLCF.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Ducted Fuel Injection and Cooled Spray Technologies for Particulate Control in Heavy-Duty Diesel Engines

Heavy-duty diesel engine manufacturers are continually seeking simple, low-cost technologies to reduce emissions. Ducted Fuel Injection (DFI) and Cooled Spray (CS) technologies have emerged as promising solutions, offering the potential for nearly sootless operation. These innovations could significantly decrease harmful Particulate Matter (PM) emissions while enabling further optimization to reduce Nitrogen Oxides (NOx) and increase efficiency. While initial experiments across various engine types have shown promise, uncertainties remain regarding the ideal geometry, scaling properties, and effectiveness of these technologies across different operating conditions. This project aims to address these knowledge gaps through experiments in both optical and metal single-cylinder engines.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2024

Low Power, Radiation Resilient Synchronous Edge Processing for Remote Monitoring

Next-generation space remote sensing systems may be equipped with imaging arrays that sense data at a rate that outstrips the processing capability of any computing hardware that can operate within a satellite’s power budget. This project developed novel convolutional and recurrent neural networks to detect and estimate point-like events amid clutter, and investigated their efficient and accurate implementation on analog in-memory computing systems that are 10-1000× more energy-efficient than digital processors. This project leveraged two memory devices at different levels of technological maturity: a large-scale analog computing prototype using commercial SONOS charge-trap memory, and electrochemical memory (ECRAM) with intrinsic radiation hardness. We experimentally demonstrated end-to-end analog processing of our neural networks on SONOS and characterized the radiation response of both SONOS and ECRAM. We advanced the state-of-the-art in ECRAM precision and reliability, and developed co-design methods to enable accurate long-term operation of SONOS analog accelerators in space radiation environments.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2020

GCAM v4.3 SSP-RCP-GCM Output Products

Supplemental dataset to: Graham N.T., M.I. Hejazi, M. Chen, E. Davies, J.A. Edmonds, S.H. Kim, and S. Turner, et al. 2020. "Humans drive future water scarcity changes across all Shared Socioeconomic Pathways." Environmental Research Letters 15, no. 1:Article No. 014007. PNNL-SA-151297. doi:10.1088/1748-9326/ab639b This analysis makes deterministic classifications of how, when water resources are limited and constrained, the human and climate systems interact. The relative effects of both systems on water scarcity are quantified at global and basin scales across 15 global futures that include five different socioeconomic conditions (the Shared Socioeconomic Pathways, SSPs) and four different climatic conditions (the Representative Concentration Pathways, RCPs). These 15 scenarios are used to first analyze a 'Human Alone' component by isolating the human impact on future scarcity. This is accomplished by holding all climate variables to their 2005 levels while altering socioeconomic growth and technological change. General circulation model (GCM) derived climate impacts for five models are applied to establish 75 'Human and Climate' scenarios which allows for the quantification of climate impacts. Climatic impacts to water supply, agricultural productivity and change, hydropower availability, and building energy demands are applied from 5 different bias-corrected GCMs to make a suite of 75 climate runs. By subtracting the human derived impacts from the 'Human Alone' scenarios, from the 'Human and Climate' scenarios the climate impact is isolated. First, the drivers of future water scarcity are evaluated by isolating the impacts that both humans and climate have while accounting for feedbacks between humans, energy, and land. Secondly, the simultaneous impacts that human and climate systems have on water scarcity are analyzed by determining whether each system is increasing or decreasing scarcity in all global water basins. Below we describe the GCAM model, the scenario components, climate derived impacts, and the calculations of scarcity changes and attribution. This study uses the global change assessment model (GCAM), with inclusions of water constraints to both renewable and nonrenewable sources of water, to investigate the relative contributions of climate and human systems on water scarcity regionally and globally under a wide range of scenarios. GCAM links socioeconomics, the energy system, land-use change, climate, and the water sector. GCAM is a market-equilibrium model that allows for prices to be adjusted within each time step to ensure that the supply and demand of goods and services remains equilibrated at each time step allowing for simultaneous market clearing across sectors. This study accounts for a limited supply of water by employing cost resource curves across all 235 basins that follow a logit formulation to determine the share of each water source (renewable surface water, nonrenewable groundwater, and desalinated water) needed to meet the water demands within all basins. As depletion of various water sources increases the extraction price increases, which leads to compounding price increases on the goods and services that require higher-priced water sources.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2023

Fine-scale evaluation of two standard 16S rRNA gene amplicon primer pairs for analysis of total prokaryotes and archaeal nitrifiers in differently managed soils

The advance of high-throughput molecular biology tools allows in-depth profiling of microbial communities in soils, which possess a high diversity of prokaryotic microorganisms. Amplicon-based sequencing of 16S rRNA genes is the most common approach to studying the richness and composition of soil prokaryotes. To reliably detect different taxonomic lineages of microorganisms in a single soil sample, an adequate pipeline including DNA isolation, primer selection, PCR amplification, library preparation, DNA sequencing, and bioinformatic post-processing is required. Besides DNA sequencing quality and depth, the selection of PCR primers and PCR amplification reactions arguably have the largest influence on the results. This study tested the performance and potential bias of two primer pairs, i.e., 515F (Parada)-806R (Apprill) and 515F (Parada)-926R (Quince) in the standard pipelines of 16S rRNA gene Illumina amplicon sequencing protocol developed by the Earth Microbiome Project (EMP), against shotgun metagenome-based 16S rRNA gene reads. The evaluation was conducted using five differently managed soils. We observed a higher richness of soil total prokaryotes by using reverse primer 806R compared to 926R, contradicting to in silico evaluation results. Both primer pairs revealed various degrees of taxon-specific bias compared to metagenome-derived 16S rRNA gene reads. Nonetheless, we found consistent patterns of microbial community variation associated with different land uses, irrespective of primers used. Total microbial communities, as well as ammonia oxidizing archaea (AOA), the predominant ammonia oxidizers in these soils, shifted along with increased soil pH due to agricultural management. In the unmanaged low pH plot abundance of AOA was dominated by the acid-tolerant NS-Gamma clade, whereas limed agricultural plots were dominated by neutral-alkaliphilic NS-Delta/NS-Alpha clades. This study stresses how primer selection influences community composition and highlights the importance of primer selection for comparative and integrative studies, and that conclusions must be drawn with caution if data from different sequencing pipelines are to be compared.

16S rRNA gene amplicon Illumina sequencing↗
DOE OSTI2022

Persistent, extensive channelized drainage modeled beneath Thwaites Glacier, West Antarctica

Abstract. Subglacial hydrology is a leading control on basal friction and the dynamics of glaciers and ice sheets. At low discharge, subglacial water flows through high-pressure, sheet-like systems that lead to low effective pressures. However, at high discharge, subglacial water melts the overlying ice into localized channels that efficiently remove water from the bed, thereby increasing effective pressure and basal friction. Recent observations suggest channelized subglacial flow exists beneath Thwaites Glacier, yet it remains unclear if stable channelization is feasible in West Antarctica, where surface melting is nonexistent and water at the bed is limited. Here, we use the MPAS-Albany Land Ice model to run a suite of over 130 subglacial hydrology simulations of Thwaites Glacier across a wide range of physical parameter choices to assess the likelihood of channelization. We then narrow our range of viable simulations by comparing modeled water thicknesses to previously observed radar specularity content, which indicates flat, spatially extensive water bodies at the bed. In all of our data-compatible simulations, stable channels reliably form within 100–200 km of the grounding line and reach individual discharge rates of 35–110 m3 s−1 at the ice–ocean boundary. While only one to two channels typically form across the 200 km width of the glacier in our simulations, their high efficiency drains water across the entire lateral extent of the glacier. We posit the large catchment size of Thwaites Glacier, its funnel-like geometry, and high basal melt rates together accumulate enough water to form stable channels. No simulations resembled observed specularity content when channelization was disabled. Our results suggest channelized subglacial hydrology has two consequences for Thwaites Glacier dynamics: (i) amplifying submarine melting of the terminus and ice shelf while (ii) simultaneously raising effective pressure within 100 km of the grounding line and increasing basal friction. The distribution of effective pressure implied from our modeling differs from parameterizations typically used in large-scale ice sheet models, suggesting the development of more process-based parameterizations may be necessary.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Autonomous Inversion of In Situ Deformation Measurement Data for Injection-Induced Stress Change

Geologic carbon storage (GCS) is likely to play a key part of the global effort to dramatically reduce CO2 emissions and perhaps even reduce atmospheric CO2 concentrations through carbon negative operations. A critical part of effort to commercialize and widely deploy this technology is developing the capability to rapidly assimilate real-time monitoring data into a form that will enable site operators to make decisions to manage the safe and efficient operations. Two of the risks associate with GCS are the risk of inducing fractures in the sealing formations that can create leakage pathways and the risk of inducing earthquakes of sufficient magnitude to cause public concern, property damage, or safety risks. To properly manage these risks the site operator needs to know the initial state of stress, the change in stress induced by injection, and the relationship between operational parameters such as injection rate and pressure and the change in stress. Current methods of estimating the change in stress require choosing the type of constitutive model and the model parameters based on core, log, and geophysical data during the characterization phase, with little feedback from operational observations to validate or refine these choices. These characterization methods interrogate the geologic formations using length scales, loading rates or magnitudes that are quite different from those encountered by the actual storage system. It is shown that errors in the assumed constitutive response, even when informed by laboratory tests on core samples, are likely to be common, large, and underestimate the magnitude of stress change caused by injection. Recent advances in borehole-based strain instruments and borehole and surface-based tilt and displacement instruments have now enabled monitoring of the deformation of the storage system throughout its operational lifespan. This data can enable validation and refinement of the knowledge of the geomechanical properties and state of the system, but brings with it a challenge to transform the raw data into actionable knowledge. We demonstrate a method that uses automatic differentiation and a finite-element based geomechanical model perform a gradient-based deterministic inversion of geomechanical monitoring data. This approach allows autonomous integration of the instrument data without the need for time consuming manual interpretation and selection of updated model parameters. Furthermore, only isotropic linear elasticity is considered in this paper, the approach presented is very flexible as to what type of geomechanical constitutive response can be used. The approach is easily adaptable to nonlinear physics-based constitutive models to account for common rock behaviors such as creep and plasticity. The approach also enables training of machine learning-based constitutive models by allowing back propagation of errors through the finite element calculations. This enables strongly enforcing known physics, such as conservation of momentum and continuity, while allowing data-driven models to learn the truly unknown physics such as the constitutive or petrophysical responses.

Burghardt, Jeffrey A.↗
DOE OSTI2024

Building a Prospective LCA Framework to Analyze Emerging Technologies in a Dynamic System Context

NREL's Lifecycle Analysis Integration into Opensource Numerical models (LiAISON) framework computes temporally explicit life cycle impacts and resource uses for specific technologies (foreground) in a dynamic system context (background). LiAISON computes results for a series of environmental mid-points enabling an analysis of prospective tradeoffs of emerging technologies toward 2100. This prospective feature is of critical importance when analyzing present-day emerging technologies whose large-scale impacts during deployment phases will occur in different, future system contexts. LiAISON systematically accounts for dynamic system changes by applying an integrated background of the future energy-economy-land-climate system, generated by exogenous integrated assessment models (IAMs). Using IAM scenarios, LiAISON generates a time-series of life cycle inventory (LCI) databases, which are then used to calculate the impacts per functional unit per time step. This expands current practice of using static, future system assumptions, e.g., a specific grid-mix each year. Further, IAM scenarios are provided in a standardized format of shared-socioeconomic pathways (SSP) and representative concentration pathways (RCP) combinations. These are coherent, regularly, published, and peer-reviewed scenario combinations that establish a reproducible and standardized societal and climate mitigation futures context. They are comparable across IAMs and expand the system boundary of the traditional LCA by including dimensions such as societal and behavioral changes. We apply the framework to assess two emerging Power-to-Hydrogen processes, high temperature electrolysis using solid oxide fuel cell (HT-SOE) and polymer electrolyte membrane electrolysis (PEME). We compare the technologies to a baseline Hydrogen production process via steam methane reforming. Despite the decarbonized electricity systems' beneficial effects on the PtH2 processes' carbon intensities, we find environmental tradeoffs, which require technology improvements via learning-by-doing to be alleviated. Future work via ongoing collaborations will focus on linking the framework to other energy-economy-land-climate models and open-source life cycle inventory databases.

ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2022

The role of breakup and coalescence in fine-scale bubble-induced turbulence. I. Dynamics

We study the effect of bubble breakup and coalescence on fine-scale dynamics of bubbly turbulent flows using direct numerical simulations. We perform two different simulations of dilute bubbly flows of void fraction 0.5%: one with bubbles breaking up and coalescing and the other without these physical processes. The volume of the fluid method is used for simulating bubbles undergoing breakup and coalescence while the bubbles are treated as rigid spheres in the immersed boundary method simulation. The energy spectrum in both types of simulation, consistent with previous studies, exhibits a -3 slope. We follow a single infinitesimal fluid element as it evolves to understand velocity gradient dynamics using conditional mean trajectories. We note finite-time divergence when the fluid element evolves under the action of inertial and pressure dynamics. The inertial, pressure, and viscous velocity gradient dynamics, when considered individually, produce the same results for bubble-induced turbulence (BIT) as with the classical homogeneous isotropic turbulence (HIT). Yet when the overall velocity gradient dynamics is considered, BIT results in non-cyclic trajectories moving toward stable node and unstable saddle while classical HIT shows cyclic behavior in their trajectories that move toward the origin. Interestingly, both the volume of fluid and immersed boundary simulations produce similar results. Therefore, there are two main takeaways from this research. First, new velocity gradient models are needed for BIT as their velocity gradients behavior is entirely different from the HIT. Second, we can neglect the bubble topology, breakup, and coalescence while studying or modeling the fine-scale dynamics of BIT.

42 ENGINEERING↗
DOE OSTI2022

Large-scale application of free energy perturbation calculations for antibody design

Abstract Alchemical free energy perturbation (FEP) is a rigorous and powerful technique to calculate the free energy difference between distinct chemical systems. Here we report our implementation of automated large-scale FEP calculations, using the Amber software package, to facilitate antibody design and evaluation. In combination with Hamiltonian replica exchange, our FEP simulations aim to predict the effect of mutations on both the binding affinity and the structural stability. Importantly, we incorporate multiple strategies to faithfully estimate the statistical uncertainties in the FEP results. As a case study, we apply our protocols to systematically evaluate variants of the m396 antibody for their conformational stability and their binding affinity to the spike proteins of SARS-CoV-1 and SARS-CoV-2. By properly adjusting relevant parameters, the particle collapse problems in the FEP simulations are avoided. Furthermore, large statistical errors in a small fraction of the FEP calculations are effectively reduced by extending the sampling, such that acceptable statistical uncertainties are achieved for the vast majority of the cases with a modest total computational cost. Finally, our predicted conformational stability for the m396 variants is qualitatively consistent with the experimentally measured melting temperatures. Our work thus demonstrates the applicability of FEP in computational antibody design.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Nuclear effective field theory: Status and perspectives

The nuclear physics landscape has been redesigned as a sequence of effective field theories (EFTs) connected to the Standard Model through symmetries and lattice simulations of Quantum Chromodynamics (QCD). EFTs in this sequence are expansions around different low-energy limits of QCD, each with its own characteristics, scales, and ranges of applicability regarding energy and number of nucleons. We review each of the three main nuclear EFTs---Chiral, Pionless, Halo/Cluster---highlighting their similarities, differences, and connections. Here, in doing so, we survey the structural properties and reactions of nuclei that have been derived from the \abinitio solution of the few- and many-body problem built upon EFT input.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗