Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Accounting”

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 361 records · Page 20

Relating the depth of the water table to the depth of weathering

Weathering of bedrock creates and occludes permeability, affecting subsurface water flow. Often, weathering intensifies above the water table. On the other hand, weathering can also commence below the water table. To explore relationships between weathering and the water table, a simplified weathering model for an eroding hillslope was formulated that takes into account both saturated and unsaturated subsurface water flow (but does not fully account for changes in dissolved gas chemistry). The phreatic line was calculated using solutions to mathematical treatments for both zones. In the model, the infiltration rate at the hill surface sets both the original and the eventual steady-state position of the water table with respect to the weathering reaction front. Depending upon parameters, the weathering front can locate either above or below the water table at steady state. Erosion also affects the water table position by changing porosity and permeability even when other hydrological conditions (e.g. hydraulic conductivity of parent material, infiltration rate at the surface) do not change. The total porosity in a hill (water storage capacity) was found to increase with infiltration rate (all else held constant). Furthermore, this effect was diminished by increasing the erosion rate. We also show examples of how the infiltration rate affects the position of the water table and how infiltration rate affects weathering advance.

58 GEOSCIENCES↗

Uncertainty quantification of bank vegetation impacts on the flood flow field in the American River, California, using large‐eddy simulations

Bank vegetation plays a key role in both hydrodynamics and morphodynamics of natural rivers; however, these effects are often unaccounted for in the computational flow dynamics of natural waterways. Recent studies using the large‐eddy simulation (LES), however, have attempted to gain insights into the impacts of bank vegetation on the mean flow field of the natural rivers using a vegetation model, which applies a sink term to the momentum equations of motion. This approach accounts for the effects of the vegetation and provides a practical approach to account for the complex patches of bank vegetation in large‐scale rivers. To implement the vegetation model, a drag coefficient reflecting the overall resistance of vegetal structures to the flow is needed, but due to the lack of calibrated data and range of size, density and type of vegetation, this parameter can be a significant source of uncertainty in the model results. Here, in this study, we use uncertainty quantification (UQ) to investigate the hydrodynamics and bed shear results when a bank vegetation is incorporated in an LES model. To this end, we used the polynomial chaos expansion and Monte Carlo sampling techniques to determine the uncertainties associated with the drag coefficient in the vegetation model and from uncertainties in the bed roughness and inflow discharge. The UQ analysis provided spatially varying confidence levels for the spanwise and vertical distribution of velocity magnitude and for the bed shear stress distributions. In addition, Sobol indices were computed to indicate the relative influence that each parameter had on the overall uncertainty. In general, it was found that uncertainty in flow discharge was the dominant source of uncertainty; however, the drag coefficient in the vegetation model and the bed roughness parameter also made significant contribution to the uncertainty near the banks and bed, respectively.

54 ENVIRONMENTAL SCIENCES↗

Evaluating Probability of Containment Effectiveness at a GCS Sites using integrated assessment modeling approach with Bayesian decision Networks

Improved scientific and engineering understanding of the behavior of geologic CO2 storage together with established regulatory framework and incentive structures raise the prospects for accelerated, large-scale deployment of this greenhouse gas emissions reduction approach. Incentive structures call for the establishment of appropriate verification and accounting approaches to support claims of the integrity of a geologic storage complex and to justify taking credit for long-term storage. In this study, we present a framework for assessing the probability of containment effectiveness over the lifetime of a geologic carbon storage site (e.g., after 70 years of injection and post-injection site performance) using forward stochastic model realizations based on site characterization data and using a monitoring-informed Bayesian network based on hypothetical detectability from surface seismic surveys over the site injection and post-injection phases. The National Risk Assessment Partnership’s open-source Integrated Assessment Model (NRAP-Open-IAM) was utilized to develop an ensemble of 10,000 a priori stochastic forecasts of CO2 containment. Those simulations were used to train the Bayesian network model to estimate the prior probabilities of the CO2 leakage mass into overlying, monitorable aquifers considering the uncertainties in the reservoir properties, permeability of potentially leaky wells and the overlying aquifers. The conditional probabilities in the Bayesian network were either learned from the NRAP-Open-IAM simulations or derived from the predefined detection thresholds for the monitoring method. Observations obtained from monitoring, over time during the site operation phases were then used to generate updated posterior probabilities of containment (and any loss from containment) in the Bayesian network by propagating the prior probabilities through the conditional probabilities. We demonstrate how to construct and use the Bayesian network for verifying the long-term storage complex effectiveness informed by monitoring based on the NRAP-Open-IAM simulations previously developed for the FutureGen 2.0 site. This approach may have relevance for stake holders to demonstrate secure geologic storage, provide a defensible, probabilistic approach to claim credit for geologic storage, and to estimate the likelihood that any fraction of the claimed credit may need to be refunded to the creditor based on available monitoring information.

Bayesian network, Risk assessment, Monitoring, car↗

A New Estimator to Correct for Bias from Tag Rate Expansion on Natural-Origin Fish Attributes in Mixed-Stock Analysis Using Parentage-Based Tagging

Abstract In fisheries where hatchery- and natural-origin conspecifics occur as mixed stocks, it is often important to estimate both the natural-origin proportion of the mixture and the composition of attributes within the natural-origin portion (e.g., genetic stock, sex, and age-class). These estimates are facilitated by parentage-based tagging, which allows large numbers of hatchery fish to be efficiently tagged and later identified. When tag rates are less than 100%, the untagged fish in the mixture include both untagged hatchery fish and natural-origin fish. Unbiased estimation of the abundance and attribute composition of the natural-origin portion requires an estimator that accounts for tag rates. Two estimators are described: one “accounting-style” estimator similar to previously described approaches and one maximum likelihood method. These estimators were evaluated and compared using simulations mimicking estimation of the composition of fish migrating past a dam. The two estimators performed similarly at reducing bias when tag rates were high, but the maximum likelihood method had smaller mean square error when tag rates were low. We provide an R package to allow usage of this estimator in a wide variety of fisheries applications.

Delomas, Thomas A. (ORCID:000000015154759X)↗

Improving Abundance Estimates of Spring–Summer Snake River Chinook Salmon for Fisheries Management

Abstract The Columbia River basin is home to a run of spring–summer Chinook Salmon Oncorhynchus tshawytscha that returns to the Snake River drainage of Idaho, Oregon, and Washington in the Pacific Northwest. Historically, the run was one of the more productive throughout the Columbia River basin. However, Snake River spring–summer Chinook Salmon have experienced declines in abundance due to overfishing, habitat degradation, and dams. Several stocks are listed as threatened under the U.S. Endangered Species Act and are supported by mitigation hatcheries funded by Idaho Power Company, the Lower Snake River Compensation Plan, and the Bonneville Power Administration. To maximize tribal and state harvest of returning hatchery adults, minimize impacts on wild fish, and ensure that enough hatchery fish return to meet broodstock needs, careful fisheries management is required. Since 2008, managers have used hatchery adults, PIT-tagged as juveniles and detected at Lower Granite Dam, to generate adult abundance estimates. In season, these estimates inform state and tribal harvest shares and ensure that broodstock needs are met. Postseason, they provide smolt-to-adult survival and return rates. Since 2012, parentage-based tagging (PBT) has provided an alternative method to estimate stock- and age-specific returns at Lower Granite Dam, since returning hatchery adults sampled at Lower Granite Dam can be assigned to their parents. We compared stock-specific abundance estimates between PIT- and PBT-derived methodologies for return years 2016–2019. Across all years, PIT tag estimates accounted for 65% of the PBT-based estimates at Lower Granite Dam across all age-groups and release sites combined. This underrepresentation across all groups equated to 49,833 fish that were not accounted for in PIT tag abundance estimates. It is clear that PBT-based estimates should aide in-season harvest management and postseason run reconstruction to avoid the known bias of estimates from PIT tags, especially during years of low returns when increased accuracy is critical.

Coykendall, D. Katharine (ORCID:0000000211482397)↗

A phase-field model for hydraulic fracture nucleation and propagation in porous media

Many geo-engineering applications, for example, enhanced geothermal systems, rely on hydraulic fracturing to enhance the permeability of natural formations and allow for sufficient fluid circulation. Over the past few decades, the phase-field method has grown in popularity as a valid approach to modeling hydraulic fracturing because of the ease of handling complex fracture propagation geometries. However, existing phase-field methods cannot appropriately capture nucleation of hydraulic fractures because their formulations are solely energy-based and do not explicitly take into account the strength of the material. Thus, in this work, we propose a novel phase-field formulation for hydraulic fracturing with the main goal of modeling fracture nucleation in porous media, for example, rocks. Built on the variational formulation of previous phase-field methods, the proposed model incorporates the material strength envelope for hydraulic fracture nucleation through two important steps: (i) an external driving force term, included in the damage evolution equation, that accounts for the material strength; (ii) a properly designed damage function that defines the fluid pressure contribution on the crack driving force. In conclusion, the comparison of numerical results for two-dimensional test cases with existing analytical solutions demonstrates that the proposed phase-field model can accurately model both nucleation and propagation of hydraulic fractures. Additionally, we present the simulation of hydraulic fracturing in a three-dimensional domain with various stress conditions to demonstrate the applicability of the method to realistic scenarios.

15 GEOTHERMAL ENERGY↗

Shortest path network interdiction with asymmetric uncertainty

Abstract This paper considers an extension of the shortest path network interdiction problem that incorporates robustness to account for parameter uncertainty. The shortest path interdiction problem is a game of two players with conflicting agendas and capabilities: an evader, who traverses the arcs of a network from a source node to a sink node using a path of shortest length, and an interdictor, who maximizes the length of the evader's shortest path by interdicting arcs on the network. It is usually assumed that the parameters defining the network are known exactly by both players. We consider the situation where the evader assumes the nominal parameter values while the interdictor uses robust optimization techniques to account for parameter uncertainty or sensor degradation. We formulate this problem as a nonlinear mixed‐integer semi‐infinite bilevel program and show that it can be converted into a mixed‐integer linear program with a second order cone constraint. We use random geometric networks and transportation networks to perform computational studies and demonstrate the unique decision strategies that our variant produces. Solving the shortest path interdiction problem with asymmetric uncertainty protects the interdictor from investing in a strategy that hinges on key interdictions performing as promised. It also provides an alternate strategy that mitigates the risk of these worst‐case possibilities.

Punla‐Green, She'ifa Z.↗

A multiscale anisotropic polymer network model coupled with phase field fracture

The study of polymers has continued to gain substantial attention due to their expanding range of applications, spanning essential engineering fields to emerging domains like stretchable electronics, soft robotics, and implantable sensors. These materials exhibit remarkable properties, primarily stemming from their intricate polymer chain network, which, in turn, increases the complexity of precisely modeling their behavior. Especially for modeling elastomers and their fracture behavior, accurately accounting for the deformations of the polymer chains is vital for predicting the rupture in highly stretched chains. Despite the importance, many robust multiscale continuum frameworks for modeling elastomer fracture tend to simplify network deformations by assuming uniform behavior among chains in all directions. Recognizing this limitation, our study proposes a multiscale fracture model that accounts for the anisotropic nature of elastomer network responses. At the microscale, damage in the chains is assumed to be driven by both the chain's entropy and the internal energy due to molecular bond distortions. In order to bridge the stretching in the chains to the macroscale deformation, we employ the maximal advance path constraint network model, inherently accommodating anisotropic network responses. As a result, chains oriented differently can be predicted to exhibit varying stretch and, consequently, different damage levels. To drive macroscale fracture based on damages in these chains, we utilize the micromorphic regularization theory, which involves the introduction of dual local-global damage variables at the macroscale. The macroscale local damage variable is obtained through the homogenization of the chain damage values, resulting in the prediction of an isotropic material response. The macroscale global damage variable is subjected to nonlocal effects and boundary conditions in a thermodynamically consistent phase field continuum formulation. Moreover, the total dissipation in the system is considered to be mainly due to the breaking of the molecular bonds at the microscale. To validate our model, we employ the double-edge notched tensile test as a benchmark, comparing simulation predictions with existing experimental data. Additionally, to enhance our understanding of the fracturing process, we conduct uniaxial tensile experiments on a square film made up of polydimethylsiloxane (PDMS) rubber embedded with a hole and notches and then compare our simulation predictions with the experimental observations. Furthermore, we visualize the evolution of stretch and damage values in chains oriented along different directions to assess the predictive capacity of the model. In conclusion, the results are also compared with another existing model to evaluate the utility of our model in accurately simulating the fracture behavior of rubber-like materials.

42 ENGINEERING↗

Data‐driven variational method for discrepancy modeling: Dynamics with small‐strain nonlinear elasticity and viscoelasticity

Abstract The effective inclusion of a priori knowledge when embedding known data in physics‐based models of dynamical systems can ensure that the reconstructed model respects physical principles, while simultaneously improving the accuracy of the solution in the previously unseen regions of state space. This paper presents a physics‐constrained data‐driven discrepancy modeling method that variationally embeds known data in the modeling framework. The hierarchical structure of the method yields fine scale variational equations that facilitate the derivation of residuals which are comprised of the first‐principles theory and sensor‐based data from the dynamical system. The embedding of the sensor data via residual terms leads to discrepancy‐informed closure models that yield a method which is driven not only by boundary and initial conditions, but also by measurements that are taken at only a few observation points in the target system. Specifically, the data‐embedding term serves as residual‐based least‐squares loss function, thus retaining variational consistency. Another important relation arises from the interpretation of the stabilization tensor as a kernel function, thereby incorporating a priori knowledge of the problem and adding computational intelligence to the modeling framework. Numerical test cases show that when known data is taken into account, the data driven variational (DDV) method can correctly predict the system response in the presence of several types of discrepancies. Specifically, the damped solution and correct energy time histories are recovered by including known data in the undamped situation. Morlet wavelet analyses reveal that the surrogate problem with embedded data recovers the fundamental frequency band of the target system. The enhanced stability and accuracy of the DDV method is manifested via reconstructed displacement and velocity fields that yield time histories of strain and kinetic energies which match the target systems. The proposed DDV method also serves as a procedure for restoring eigenvalues and eigenvectors of a deficient dynamical system when known data is taken into account, as shown in the numerical test cases presented here.

Masud, Arif↗

Mechanistic Models as Framework for Understanding Biomarker Disposition: Prediction of Creatinine‐Drug Interactions

Creatinine is widely used as a biomarker of glomerular filtration, and, hence, renal function. However, transporter‐mediated secretion also contributes to its renal clearance, albeit to a lesser degree. Inhibition of these transporters causes transient serum creatinine elevation, which can be mistaken as impaired renal function. The current study developed mechanistic models of creatinine kinetics within physiologically based framework accounting for multiple transporters involved in creatinine renal elimination, assuming either unidirectional or bidirectional‐OCT2 transport (driven by electrochemical gradient). Robustness of creatinine models was assessed by predicting creatinine‐drug interactions with 10 perpetrators; performance evaluation accounted for 5% intra‐individual variability in serum creatinine. Models showed comparable predictive performances of the maximum steady‐state effect regardless of OCT2 directionality assumptions. However, only the bidirectional‐OCT2 model successfully predicted the minimal effect of ranitidine. The dynamic nature of models provides clear advantage to static approaches and most advanced framework for evaluating interplay between multiple processes in creatinine renal disposition.

Scotcher, Daniel↗

A Poisson–Nernst–Planck Model of Ion Transport and Interface Segregation in Metal–Insulator–Semiconductor Structures and Solar Cells

A numerical model that describes the transport of mobile ionic species in metal–insulator–semiconductor (MIS) and photovoltaic (PV) devices subject to temperature and voltage stress is presented. The finite element method (FEM) is used to solve the Nernst–Planck equation while imposing Poisson's equation self-consistently as a restriction for the electrostatic potential. This allows the contribution of the ionic species to the potential to be taken into account. Using a variational formulation eases the implementation of diverse boundary conditions, including the incorporation of segregation kinetics at the device interfaces. Segregation across the dielectric–semiconductor interface is relevant to modeling the electronic device degradation in systems where contamination reaches the semiconductor. The model in closed systems with no-flux boundary conditions is validated first. In the limiting case of low contamination levels with respect to the gate bias, the FEM solution matches analytically derived approximations. Then, the implementation is broadened to include an open boundary at the dielectric–semiconductor interface to account for leakage of ions. The predicted time dependence of the flatband voltage in Na-contaminated MIS test structures agrees well with measurements. This model successfully captures the role of long-range ion transport at concentrations of relevance to electronic and PV device instability and neuromorphic computing.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Bayesian Calibration of Stochastic Agent Based Model via Random Forest

Agent-based models (ABM) provide an excellent framework for modeling outbreaks and interventions in epidemiology by explicitly accounting for diverse individual interactions and environments. However, these models are usually stochastic and highly parametrized, requiring precise calibration for predictive performance. When considering realistic numbers of agents and properly accounting for stochasticity, this high-dimensional calibration can be computationally prohibitive. This paper presents a random forest-based surrogate modeling technique to accelerate the evaluation of ABMs and demonstrates its use to calibrate an epidemiological ABM named CityCOVID via Markov chain Monte Carlo (MCMC). The technique is first outlined in the context of CityCOVID's quantities of interest, namely hospitalizations and deaths, by exploring dimensionality reduction via temporal decomposition with principal component analysis (PCA) and via sensitivity analysis. The calibration problem is then presented, and samples are generated to best match COVID-19 hospitalization and death numbers in Chicago from March to June in 2020. Further, these results are compared with previous approximate Bayesian calibration (IMABC) results, and their predictive performance is analyzed, showing improved performance with a reduction in computation.

60 APPLIED LIFE SCIENCES↗

Assessing the blockage effect of wind turbines and wind farms using an analytical vortex model

Wind farm blockage effects are currently neglected in the prediction of wind farm energy yield, typically leading to an overestimation of the production. This work presents a novel method to assess wind farm production, while accounting for blockage effects. We apply a vortex model, based on a cylindrical wake, to assess induction effects. We present variations of the model to account for finite wake length, finite tip-speed ratios, and the proximity to the ground. The results are applied to single rotors in aligned and yawed conditions and to different wind farm layouts. We provide far-field approximations for faster estimates of the velocity field. Further, this article includes a new methodology to couple the induction model to engineering wake models, such as the ones present in the FLOw Redirection and Induction in Steady State (FLORIS). We compare the results to actuator disk simulations for various operating conditions of a single turbine and different wind farm layouts. We found that the mean relative error of the model in the induction zone is typically around 0.2% compared with actuator disk simulations. The computational time of the velocity field using the analytical vortex model is three orders of magnitude less than the one obtained with the actuator disk simulation.

17 WIND ENERGY↗

Staying cool: The development of India's pioneering energy efficiency policy for chillers

Abstract Chillers are the key energy‐consuming piece of equipment in most large space cooling systems. They can account for more than 40% of total energy consumption in any typical commercial building and in specific industrial building installations depending on the product and process. With climate change driving rapidly rising demand for commercial and industrial cooling, especially in emerging economies, energy efficiency policies for chillers are crucial to moderate future electricity demand and thereby minimize greenhouse gas emissions. However, most developing and emerging economies lack energy efficiency policies for this significant and growing cooling equipment sector. India recently announced new energy efficiency policies for chillers, projected to avoid 28 TWh of electricity consumption and 23 million tons of CO 2 emissions cumulatively through 2030. This article presents a case study of the unique challenges and solutions required to design a new energy efficiency policy for chillers in India. First, the energy consumed in chiller systems is sensitive to building occupancy, load changes, seasonal variations, operation and maintenance, and climatic conditions. In order to address these sensitivities, the chiller testing standard in India factors local weather and temperature conditions into energy performance ratings to account for these sensitivities. Further, every large chiller is a customized product, therefore developing a viable policy and implementation process required collaboration with international and national experts in standardization, manufacturers, building consultants, and testing laboratories. This article is categorized under: Energy Efficiency > Economics and Policy Solar Energy: Thermal ‐‐> Climate and Environment

Mukherjee, P. K.↗

Coupled Boltzmann transport equations of heavy quarks and quarkonia in quark-gluon plasma

We develop a framework of coupled transport equations for open heavy flavor and quarkonium states, in order to describe their transport inside the quark-gluon plasma. Our framework is capable of studying simultaneously both open and hidden heavy flavor observables in heavy-ion collision experiments and can account for both, uncorrelated and correlated recombination. Our recombination implementation depends on real-time open heavy quark and antiquark distributions. We carry out consistency tests to show how the interplay among open heavy flavor transport, quarkonium dissociation and recombination drives the system to equilibrium. We then apply our framework to study bottomonium production in heavy-ion collisions. We include Υ(1S), Υ(2S), Υ(3S), χb (1P) and χb (2P) in the framework and take feed-down contributions during the hadronic gas stage into account. Cold nuclear matter effects are included by using nuclear parton distribution functions for the initial primordial heavy flavor production. A calibrated 2 + 1 dimensional viscous hydrodynamics is used to describe the bulk QCD medium. We calculate both the nuclear modification factor R AA of all bottomonia states and the azimuthal angular anisotropy coefficient v 2 of the Υ(1S) state and find that our results agree reasonably with experimental measurements. Our calculations indicate that correlated cross-talk recombination is an important production mechanism of bottomonium in current heavy-ion experiments. The importance of correlated recombination can be tested experimentally by measuring the ratio of R AA ( χb (1P)) and R AA (Υ(2S)).

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Low-energy radiative backgrounds in CCD-based dark-matter detectors

The reach of sub-GeV dark-matter detectors is at present severely affected by low-energy events from various origins. We present the theoretical methods to compute the single- and few-electron events that arise from secondary radiation emitted by high-energy particles as they pass through detector materials and perform a detailed simulation to quantify them at (Skipper) CCD-based experiments, focusing on the SENSEI data collected at Fermilab near the MINOS cavern. The simulations account for the generation of secondaries from Cherenkov and luminescent recombination radiation; photo-absorption in the bulk, backside layer, pitch adapter, and epoxy; the photon reflection and refraction at interfaces; thin-film interference; the roughness of the interfaces; the dynamics of charges produced in the highly doped CCD-backside-layers; and the partial charge collection on the CCD backside. We consider several systematic uncertainties, notably those stemming from the backside modeling, which we estimate with a “fiducial” and an “extreme” charge-diffusion model, with the former model being preferred due to better agreement with partial-charge collection data. We find that Cherenkov photons constitute about 30% of the observed single-electron events for both diffusion models; radiative recombination contributes negligibly to the event rate for the fiducial model, although it can dominate over Cherenkov for the extreme model. We also estimate the fraction of 2-electron events that arise from 1-electron event coincidences in the same pixel, finding that the entire 2-electron rate can be explained by coincidences of radiative events and spurious charge. Accounting for both radiative and non-radiative backgrounds, we project the sensitivity of future Skipper-CCD-based experiments to different dark-matter models. For light-mediator models with dark-matter masses of 1, 5, and 10 MeV, we find that future experiments with 10-kg-year exposures and successful background mitigation could have a sensitivity that is larger by 9, 3, and 2 orders of magnitude, respectively, when compared to an experiment without background improvements.

79 ASTRONOMY AND ASTROPHYSICS↗

Bragg-Primakoff axion photoconversion in crystal detectors

Abstract Axions and axion-like pseudoscalar particles with dimension-5 couplings to photons exhibit coherent Primakoff scattering with ordered crystals at keV energy scales, making for a natural detection technique in searches for solar axions. We find that there are large suppressive corrections, potentially greater than a factor of$$ \mathcal{O} $$ O (10 3 ), to the coherent enhancement when taking into account absorption of the final state photon. This effect has already been accounted for in light-shining-through-wall experiments through the language of Darwin classical diffraction, but is missing from the literature in the context of solar axion searches that use a matrix element approach. We extend the treatment of the event rate with a heuristic description of absorption effects to bridge the gap between these two languages. Furthermore, we explore the Borrmann effect of anomalous absorption in lifting some of the event rate suppression by increasing the coherence length of the conversion. We study this phenomenon in Ge, NaI, and CsI crystal experiments and its impact on the projected sensitivities of SuperCDMS, LEGEND, and SABRE to the solar axion parameter space. Lastly, we comment on the reach of multi-tonne scale crystal detectors and strategies to maximize the discovery potential of experimental efforts in this vein.

Physics↗

Quantum calculations of the cavity shift in electron magnetic moment measurements

The measurement of the anomalous electron magnetic moment g – 2 through quantum transitions of a single trapped electron is the most stringent test of quantum field theory. These experiments are now so precise that they must account for the effects of the cavity containing the electron. Classical calculations of this “cavity shift” must subtract the electron’s divergent self-field, and thus require knowledge of the exact Green’s function for the cavity’s electromagnetic field. We perform the first fully quantum calculation of the cavity shift in a closed cavity, which instead involves subtracting linearly divergent cavity mode sums and integrals. Using contour integration methods, we find perfect agreement with existing classical results for both spherical and cylindrical cavities, justifying their current use. Moreover, our mode-based results can be naturally generalized to account for systematic effects, necessary to push future measurements to the next order of magnitude in precision.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗