Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “applied statistics”

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,621 records · Page 90

Data Processing for the Near Earth Asteroid Rendezvous (NEAR), X-Ray and Gamma-Ray Spectrometer (XRS) Ground System

An X-ray and Gamma-ray spectrometer (XGRS) is onboard the Near Earth Asteroid Rendezvous (NEAR) spacecraft to determine the elemental composition of the surface of the asteroid 433Eros. The Eros asteroid is highly non-spherical in physical shape and the development of data management and analysis methodologies are in several areas a divergence from traditional remotely sensed geographical information systems techniques. Field of view and asteroid surface geometry must be derived virtually and then combined with real measurements of solar, spectral and instrument calibration information to derive meaningful scientific results. Spatial resolution of planned geochemical maps will be improved from the initial conditions of low statistical significance per integration by repeated surface flyovers and regional spectral accumulation. This paper describes the results of a collaborative effort of design and development of the NEAR XGRS instrument ground system undertaken by participants at the Goddard Space Flight Center, University of Arizona, Cornell University, Applied Physics Laboratory, and Max Plank institute.

McClanahan, Timothy P.↗

Noise Radiation From a Leading-Edge Slat

This paper extends our previous computations of unsteady flow within the slat cove region of a multi-element high-lift airfoil configuration, which showed that both statistical and structural aspects of the experimentally observed unsteady flow behavior can be captured via 3D simulations over a computational domain of narrow spanwise extent. Although such narrow domain simulation can account for the spanwise decorrelation of the slat cove fluctuations, the resulting database cannot be applied towards acoustic predictions of the slat without invoking additional approximations to synthesize the fluctuation field over the rest of the span. This deficiency is partially alleviated in the present work by increasing the spanwise extent of the computational domain from 37.3% of the slat chord to nearly 226% (i.e., 15% of the model span). The simulation database is used to verify consistency with previous computational results and, then, to develop predictions of the far-field noise radiation in conjunction with a frequency-domain Ffowcs-Williams Hawkings solver.

Lockard, David P.↗

Missing Data and Multiple Imputation: An Unbiased Approach

The default method of dealing with missing data in statistical analyses is to only use the complete observations (complete case analysis), which can lead to unexpected bias when data do not meet the assumption of missing completely at random (MCAR). For the assumption of MCAR to be met, missingness cannot be related to either the observed or unobserved variables. A less stringent assumption, missing at random (MAR), requires that missingness not be associated with the value of the missing variable itself, but can be associated with the other observed variables. When data are truly MAR as opposed to MCAR, the default complete case analysis method can lead to biased results. There are statistical options available to adjust for data that are MAR, including multiple imputation (MI) which is consistent and efficient at estimating effects. Multiple imputation uses informing variables to determine statistical distributions for each piece of missing data. Then multiple datasets are created by randomly drawing on the distributions for each piece of missing data. Since MI is efficient, only a limited number, usually less than 20, of imputed datasets are required to get stable estimates. Each imputed dataset is analyzed using standard statistical techniques, and then results are combined to get overall estimates of effect. A simulation study will be demonstrated to show the results of using the default complete case analysis, and MI in a linear regression of MCAR and MAR simulated data. Further, MI was successfully applied to the association study of CO2 levels and headaches when initial analysis showed there may be an underlying association between missing CO2 levels and reported headaches. Through MI, we were able to show that there is a strong association between average CO2 levels and the risk of headaches. Each unit increase in CO2 (mmHg) resulted in a doubling in the odds of reported headaches.

Foy, M.↗

Statistical Study on Multispacecraft Widespread Solar Energetic Particle Events During Solar Cycle 24

We conduct a statistical study on the large three-spacecraft widespread solar energetic particle (SEP) events. Longitudinal distributions of the peak intensities, onset delays, and relation between the SEP intensity, coronal mass ejection (CME) shock speed, width, and the kinetic energy of the CME have been investigated. We apply a Gaussian fit to obtain the SEP intensity Io and distribution width and a forward-modeling fit to determine the true shock speed and true CME width. We found a good correction between and connection angel to the flare site and Io and the kinetic energy of the CME. By including the true chock speed and true CME widths, we reduce root-mean-square errors on the predicted SEP intensity by ~41% for protons compared to Richards et al.'s (2014, https://doi.org/10.1007/s11207-014-0524-8) prediction. The improved correction between the CME kinetic energy and SEP intensity provides strong evidence for the CME-shock acceleration theory of SEPs. In addition, we found that electron and proton release time delays (DTs) relative to Type II radio bursts increase with connection angles. The average electron (proton) DT is ~14 (32) min for strongly anisotropic events and ~2.5 (4.4) hr for weakly anisotropic events. Poor magnetic connectivity and large scattering effects are two main reasons to cause large delays.

Xie, H.↗

Failure Modes in Capacitors When Tested Under a Time-Varying Stress

Power-on failure has been the prevalent failure mechanism for solid tantalum capacitors in decoupling applications. A surge step stress test (SSST) has been previously applied to identify the critical stress level of a capacitor batch to give some predictability to the power-on failure mechanism [1]. But SSST can also be viewed as an electrically destructive test under a time-varying stress (voltage). It consists of rapidly charging the capacitor with incremental voltage increases, through a low resistance in series, until the capacitor under test is electrically shorted. When the reliability of capacitors is evaluated, a highly accelerated life test (HALT) is usually adopted since it is a time-efficient method of determining the failure mechanism; however, a destructive test under a time-varying stress such as SSST is even more time efficient. It usually takes days or weeks to complete a HALT test, but it only takes minutes for a time-varying stress test to produce failures. The advantage of incorporating a specific time-varying stress profile into a statistical model is significant in providing an alternative life test method for quickly revealing the failure mechanism in capacitors. In this paper, a time-varying stress that mimics a typical SSST has been incorporated into the Weibull model to characterize the failure mechanism in different types of capacitors. The SSST circuit and transient conditions for correctly surge testing capacitors are discussed. Finally, the SSST was applied for testing Ta capacitors, polymer aluminum capacitors (PA capacitors), and multi-layer ceramic (MLC) capacitors with both precious metal electrodes (PME) and base metal electrodes (BME). The test results are found to be directly associated with the dielectric layer breakdown in Ta and PA capacitors and are independent of the capacitor values, the way the capacitors were built, and the capacitors manufacturers. The test results also show that MLC capacitors exhibit surge breakdown voltages much higher than the rated voltage and that the breakdown field is inversely proportional to the dielectric layer thickness. The SSST data can also be used to comparatively evaluate the voltage robustness of capacitors for decoupling applications.

Liu, David (Donhang)↗

IRA Energy Community Data Layers

Data, geospatial data resources, and the linked mapping tool and web services reflect data for two types of potentially qualifying energy communities: 1) Census tracts and directly adjoining tracts that have had coal mine closures since 1999 or coal-fired electric generating unit retirements since 2009. These census tracts qualify as energy communities. 2) Metropolitan statistical areas (MSAs) and non-metropolitan statistical areas (non-MSAs) that are energy communities for 2023 and 2024, along with their fossil fuel employment (FFE) status. Additional information on energy communities and related tax credits can be accessed on the Interagency Working Group on Coal & Power Plant Communities & Economic Revitalization Energy Communities website (https://energycommunities.gov/energy-community-tax-credit-bonus/). Use limitations: these spatial data and mapping tool may not be relied upon by taxpayers to substantiate a tax return position or for determining whether certain penalties apply and will not be used by the IRS for examination purposes. The mapping tool does not reflect the application of the law to a specific taxpayer’s situation, and the applicable Internal Revenue Code provisions ultimately control.

Census Tract↗

Unlocking hidden information in sparse small-angle neutron scattering measurements

Hypothesis Small-Angle Neutron Scattering (SANS) is a powerful technique for studying soft matter systems such as colloids, polymers, and lyotropic phases, providing nanoscale structural insights. However, its effectiveness is limited by low neutron flux, leading to long acquisition times and noisy data. Here, we hypothesize that Bayesian statistical inference using Gaussian Process Regression (GPR) can reconstruct high-fidelity scattering data from sparse measurements by leveraging intensity smoothness and continuity. Experiments and Simulations The method was benchmarked computationally and validated through SANS experiments on various soft matter systems, including wormlike micelles, colloidal suspensions, polymeric structures, and lyotropic phases. GPR-based inference was applied to both experimental and synthetic data to evaluate its effectiveness in noise reduction and intensity reconstruction. Findings GPR significantly enhances SANS data quality and therefore reducing measurement times by up to two orders of magnitude. This cost-effective approach maximizes experimental efficiency, enabling high-throughput studies and real-time monitoring of dynamic systems. It is particularly beneficial for weakly scattering and time-sensitive studies. Beyond SANS, this framework applies to other low-SNR techniques, including laboratory-based small-angle X-ray scattering and various dynamical scattering methods. Furthermore, it offers transformative potential for compact neutron sources, enhancing their viability for structural analysis in resource-limited settings.

Small angle neutron scattering↗

The statistical significance of error probability as determined from decoding simulations for long codes

The very low error probability obtained with long error-correcting codes results in a very small number of observed errors in simulation studies of practical size and renders the usual confidence interval techniques inapplicable to the observed error probability. A natural extension of the notion of a 'confidence interval' is made and applied to such determinations of error probability by simulation. An example is included to show the surprisingly great significance of as few as two decoding errors in a very large number of decoding trials.

Massey, J. L.↗

Rapid Aero Modeling for Urban Air Mobility Aircraft in Wind-Tunnel Tests

Rapid Aero Modeling (RAM) applied to wind tunnel testing, RAM-T, is an approach to efficiently and automatically obtain aerodynamic models during testing. The approach saves time and resources by responding to the demand for experimental efficiency and model fidelity. Motivation for this demand is more acute when investigating a class of vehicles categorized as Urban Air Mobility (UAM) aircraft where many features from both aircraft and rotorcraft are present. RAM-T provides a feedback loop around the test facility to guide the test toward high-fidelity, statistically rigorous aircraft models. The general RAM approach is applicable to computational or physical experiments. It combines concepts from design of experiment theory and aircraft system identification theory that allow the user the freedom to choose, in advance of the test, a specific level of fidelity in terms of prediction error. RAM only collects data required to meet the user-specified fidelity and fidelity is only limited by the facility and test article capabilities. This paper presents results from tests conducted for development of an automated RAM-T technology. The results highlight some of the unique features of RAM applied to eVTOL configurations.

Aerodynamics↗

Quantifying the impact of climate change on erosion - 20397

The New York State Energy Research and Development Authority (NYSERDA) is the owner of the Western New York Nuclear Service Center (WNYNSC), a 1,351-ha (3,338-ac) site located approximately 48 km (30 mi) south of Buffalo, New York. In 1962, Nuclear Fuel Services, Inc. (NFS) entered into Agreements with the Atomic Energy Commission and New York State to construct the first commercial reprocessing plant of nuclear fuel in the United States. NFS, a private company, built and operated the spent fuel reprocessing plant and waste disposal facilities, processing 640 Mg (metric tons, or 705 short tons) of spent nuclear fuel from 1966 to 1972 under an Atomic Energy Commission license. Nuclear fuel reprocessing operations ended in 1972 and never reopened, leaving behind radioactive and chemical wastes. Erosion can play an important role in the fate and transport of waste at sites where disposal of long-lived waste is anticipated. The statistical characterization of key processes related to erosion is essential to the understanding of site stability through time. One of the key processes governing erosion is extreme precipitation and it is critical that trends in the distribution of extreme precipitation events through time be represented. The evidence of climate change is increasingly well documented and projected impacts on extreme precipitation events should be incorporated in performance assessment studies when relevant. Not evaluating future climate states in a performance assessment is contradictory to good modeling practice. Specifically, excluding climate change limits development of modeling information that could aid in effective decision making. The current climate literature provides both observational evidence and climate model projections of climate trends and/or climate change in the late 20. and early 21. centuries for North America and the northeast United States. In this work, this information was used to assess the impacts of potential changes in climate on erosion processes. The goal was to understand how projections of future climate relate to the performance of the WNYNSC through time. In the first stage of this work multi-temporal historical aerial images were analyzed in conjunction with orthophotography and Lidar data to develop probability distributions for variables representing important erosion processes. The information from these analyses was then used in conjunction with simulated data from the West Valley Erosion Working Group (EWG). Analysis of EWG simulations provides estimates for the change in erosion rates through time that is driven by changes in climate. The time-varying rates of erosion change were applied to the historical aerial imagery data in order to inform time-varying rates of erosion that are driven by changes in climate. Ultimately, this process resulted in the identification of locations for features like gully heads using both the Lidar dataset and projection of the estimated location from the historical aerial photo under consideration back to the Lidar dataset. The distance between the estimate of the location from the historical image and that of the Lidar dataset was the estimated distance the feature has moved. This distance was then divided by the number of years between the Lidar dataset and the year of the historical aerial image of interest to get a rate of movement through time. This was done for several dozen points on each historical aerial image. The analyses from the EWG were used to characterize the relative impact of climate change on the erosion rates. This relative impact was quantified by comparing the LEM based estimates of erosion that were derived using historical climate data with LEM based estimates of erosion that were derived using projections of future climate. The relative increase in the erosion rates was then applied to the historical estimates of erosion derived from the analysis of the historical aerial imagery. This approach was used since the LEM-based estimated of the historical erosion rates from the past century have a bias towards underestimating erosion. This underestimation is hypothesized to be a consequence of an inability of the LEMs to account for the impacts of land-use change (i.e. deforestation) which had been shown to significantly increase erosion in similar Northern hardwood forest ecosystems. In summary, gully head retreat rate and gully widening rate were characterized using statistical probability distributions from the historical aerial image analysis. Simulated data from the EWG was used to estimate the climatically-driven changes in erosion rates through time. The changes through time from the EWG were applied to the gully head retreat rate and gully widening rate characterized using the historical aerial image analysis. This information can be used to inform PPA models to project future risks from a given site. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

The Applied Meteorology Unit: Nineteen Years Successfully Transitioning Research into Operations for America's Space Program

The Applied Meteorology Unit (AMU) provides technology transition and technique development to improve operational weather support to the Space Shuttle and the entire American space program. The AMU is funded and managed by NASA and operated by a contractor that provides five meteorologists with a diverse mix of advanced degrees, operational experience, and associated skills including data processing, statistics, and the development of graphical user interfaces. The AMU's primary customers are the U.S. Air Force 45th Weather Squadron at Patrick Air Force Base, the National Weather Service Spaceflight Meteorology Group at NASA Johnson Space Center, and the National Weather Service Melbourne FL Forecast Office. The AMU has transitioned research into operations for nineteen years and worked on a wide range of topics, including new forecasting techniques for lightning probability, synoptic peak winds,.convective winds, and summer severe weather; satellite tools to predict anvil cloud trajectories and evaluate camera line of sight for Space Shuttle launch; optimized radar scan strategies; evaluated and implemented local numerical models; evaluated weather sensors; and many more. The AMU has completed 113 projects with 5 more scheduled to be completed by the end of 2010. During this rich history, the AMU and its customers have learned many lessons on how to effectively transition research into operations. Some of these lessons learned include collocating with the operational customer and periodically visiting geographically separated customers, operator submitted projects, consensus tasking process, use of operator primary advocates for each project, customer AMU liaisons with experience in both operations and research, flexibility in adapting the project plan based on lessons learned during the project, and incorporating training and other transition assistance into the project plans. Operator involvement has been critical to the AMU's remarkable success and many awards from NASA, the National Weather Association, and two citations from the Navy's Center of Excellence for Best Manufacturing Practices. This paper will present the AMU's proven methods and explain how they may be applied by other organizations to effectively transition research into operations.

Madura, John T.↗

Machine learning a molecular Hamiltonian for predicting electron dynamics

We develop a computational method to learn a molecular Hamiltonian matrix from matrix-valued time series of the electron density. As we demonstrate for three small molecules, the resulting Hamiltonians can be used for electron density evolution, producing highly accurate results even when propagating 1,000 time steps beyond the training data. As a more rigorous test, we use the learned Hamiltonians to simulate electron dynamics in the presence of an applied electric field, extrapolating to a problem that is beyond the field-free training data. We find that the resulting electron dynamics predicted by our learned Hamiltonian are in close quantitative agreement with the ground truth. Our method relies on combining a reduced-dimensional, linear statistical model of the Hamiltonian with a time-discretization of the quantum Liouville equation within time-dependent Hartree Fock theory. In conclusion, we train the model using a least-squares solver, avoiding numerous, CPU-intensive optimization steps. For both field-free and field-on problems, we quantify training and propagation errors, highlighting areas for future development.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A neural network for determination of latent dimensionality in Nonnegative Matrix Factorization

Non-negative Matrix Factorization (NMF) has proven to be a powerful unsupervised learning method for uncovering hidden features in complex and noisy datasets with applications in data mining, text recognition, dimension reduction, face recognition, anomaly detection, blind source separation, and many other fields. An important input for NMF is the latent dimensionality of the data, that is, the number of hidden features, K, present in the explored dataset. Unfortunately, and this quantity is rarely known a priori. The existing methods for determining latent dimensionality, such as Automatic Relevance Determination (ARD), are mostly heuristic and utilize different characteristics to estimate the number of hidden features. However, all of them require human presence to make a final determination of K. Here we utilize a supervised machine learning approach in combination with a recent method for model determination, called NMFk, to determine the number of hidden features automatically. NMFk performs a set of NMF simulations on an ensemble of matrices, obtained by bootstrapping the initial dataset, and estimates which K produces stable groups of latent features that reconstruct the initial dataset well. We then train a Multi-Layter Perceptron (MLP) classifier network to determine the correct number of latent features utilizing the statistics and characteristics of the NMF solution, obtained from NMFk. In order to train the MLP classifier, a training set of 58,660 matrices with predetermined latent features were factorized with NMFk. The MLP classifier in conjunction with NMFk maintains a greater than 95% success rate when applied to a held out test set. Additionally, when applied to two well-known benchmark datasets, the swimmer and MIT face data, NMFk/MLP correctly recovers the established number of hidden features. Finally, we compare the accuracy of our method to the ARD, AIC and Stability-based methods.

97 MATHEMATICS AND COMPUTING↗

Absolutely Stable Time Crystals at Finite Temperature

Here, we show that locally interacting, periodically driven (Floquet) Hamiltonian dynamics coupled to a Langevin bath support finite-temperature discrete time crystals (DTCs) with an infinite autocorrelation time. By contrast to both prethermal and many-body localized DTCs, the time crystalline order we uncover is stable to arbitrary perturbations, including those that break the time translation symmetry of the underlying drive. Our approach utilizes a general mapping from probabilistic cellular automata to open classical Floquet systems undergoing continuous-time Langevin dynamics. Applying this mapping to a variant of the Toom cellular automaton, which we dub the "π-Toom time crystal," leads to a 2D Floquet Hamiltonian with a finite-temperature DTC phase transition. We provide numerical evidence for the existence of this transition, and analyze the statistics of the finite temperature fluctuations. Finally, we discuss how general results from the field of probabilistic cellular automata imply the existence of discrete time crystals (with an infinite autocorrelation time) in all dimensions, d ≥ 1.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Effect of Air Plasma Sprayed Flash Bond Coatings on Furnace Cycle Lifetime of Disks and Rods

Air plasma sprayed (APS) flash coatings on high velocity oxygen fuel (HVOF) bond coatings are well known to extend the lifetime of thermal barrier coatings (TBCs). Recent work compared flash coatings of NiCoCrAlY and NiCoCrAlYHfSi applied to both rods and disk substrates of alloy 247. For rod specimens, 100 h cycles were used at 1100 °C in wet air. Both flash coatings significantly improved the lifetime compared to HVOF-only and vacuum plasma spray (VPS)-only MCrAlY bond coatings with no statistical difference between the two flash coatings. For disk specimens tested in 1 h cycles at 1100 °C in wet air, the NiCoCrAlY flash coating significantly outperformed an HVOF-only NiCoCrAlYHfSi bond coating and a NiCoCrAlYHfSi flash coating. The flash coatings formed a mixed oxide-metal zone that appeared to inhibit crack formation and therefore extend lifetime. In addition to the flash coating increasing the bond coating roughness, the underlying HVOF layer acted as a source of Al for this intermixed zone and prevented the oxide from penetrating deeper into the bond coating. Here, the lower Y+Hf content in the Y-only flash coating appeared to minimize oxidation in the flash layer, thereby increasing the benefit compared to a NiCoCrAlYHfSi flash coating.

42 ENGINEERING↗

End-Use Savings Shapes Measure Documentation: Variable Refrigerant Flow with Heat Recovery and Dedicated Outdoor Air System

Building on the successfully completed effort to calibrate and validate the U.S. Department of Energy's ResStock™ and ComStock™ models over the past 3 years, the objective of this work is to produce national data sets that empower analysts working for federal, state, utility, city, and manufacturer stakeholders to answer a broad range of analysis questions. The goal of this work is to develop energy efficiency, electrification, and demand flexibility end-use load shapes (electricity, gas, propane, or fuel oil) that cover a majority of the high-impact, market-ready (or nearly market-ready) measures. "Measures" refers to energy efficiency variables that can be applied to buildings during modeling. An end-use savings shape is the difference in energy consumption between a baseline building and a building with an energy efficiency, electrification, or demand flexibility measure applied. It results in a time-series profile that is broken down by end use and fuel (electricity or on-site gas, propane, or fuel oil use) at each time step. ComStock is a highly granular, bottom-up model that uses multiple data sources, statistical sampling methods, and advanced building energy simulations to estimate the annual subhourly energy consumption of the commercial building stock across the United States. The baseline model intends to represent the U.S. commercial building stock as it existed in 2018. The methodology and results of the baseline model are discussed in the final technical report of the End-Use Load Profiles project. This documentation focuses on a single heating, ventilation, and air-conditioning (HVAC) end-use savings shape measure - a variable refrigerant flow with heat recovery (VRF HR) heating and cooling system coupled with a dedicated outdoor air system (DOAS) for ventilation. This measure replaces existing multi-zone variable air volume (VAV) systems or single-zone rooftop units (RTU) with a VRF HR system coupled with a DOAS that includes an energy/heat recovery ventilator (E/HRV). The measure covers 53% of exisiting building stock's floor area and is not applicable to HVAC system types using district heating or cooling or buildings/spaces that include high-ventilation spaces such as kitchens where the amount of exhaust air is large. A DOAS with E/HRV is used to provide required outdoor ventilation air to spaces since ventilation air is generally not supplied by a VRF HR system. An exhaust air energy recovery ventilator (ERV ) with sensible and latent heat exchange is added to humid climate zones while a heat recovery ventilator (HRV ) with sensible only exchange is added to drier climate zones. The ERV is modeled as a fixed membrane plate counterflow heat exchanger, while the HRV is modeled as a sensible-only fixed aluminum plate counterflow heat exchanger. Both systems include a bypass (for temperature control and economizer lockout) and minimum exhaust temperature control for frost prevention.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

End-Use Savings Shapes Measure Documentation: Economizers

Building on the successfully completed effort to calibrate and validate the U.S. Department of Energy's ResStock™ and ComStock™ models over the past 3 years, the objective of this work is to produce national data sets that empower analysts working for federal, state, utility, city, and manufacturer stakeholders to answer a broad range of analysis questions. The goal of this work is to develop energy efficiency, electrification, and demand flexibility end-use load shapes (electricity, gas, propane, or fuel oil) that cover a majority of the high-impact, market-ready (or nearly market-ready) measures. "Measures" refers to energy efficiency variables that can be applied to buildings during modeling. An end-use savings shape is the difference in energy consumption between a baseline building and a building with an energy efficiency, electrification, or demand flexibility measure applied. It results in a time-series profile that is broken down by end use and fuel (electricity or on-site gas, propane, or fuel oil use) at each time step. ComStock is a highly granular, bottom-up model that uses multiple data sources, statistical sampling methods, and advanced building energy simulations to estimate the annual subhourly energy consumption of the commercial building stock across the United States. The baseline model intends to represent the U.S. commercial building stock as it existed in 2018. The methodology and results of the baseline model are discussed in the final technical report of the End-Use Load Profiles project.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

The Evolution of AGN Activity in Brightest Cluster Galaxies

We present the results of an analysis of Wide-field Infrared Survey Explorer (WISE) observations of the full 2500 deg 2 South Pole Telescope (SPT)-Sunyaev–Zel'dovich cluster sample. We describe a process for identifying active galactic nuclei (AGN) in brightest cluster galaxies (BCGs) based on WISE mid-IR color and redshift. Applying this technique to the BCGs of the SPT-SZ sample, we calculate the AGN-hosting BCG fraction, which is defined as the fraction of BCGs hosting bright central AGNs over all possible BCGs. Assuming an evolving single-burst stellar population model, we find statistically significant evidence (>99.9%) for a mid-IR excess at high redshift compared to low redshift, suggesting that the fraction of AGN-hosting BCGs increases with redshift over the range of 0 < z < 1.3. The best-fit redshift trend of the AGN-hosting BCG fraction has the form (1 + z) 4.1±1.0 . These results are consistent with previous studies in galaxy clusters as well as as in field galaxies. One way to explain this result is that member galaxies at high redshift tend to have more cold gas. While BCGs in nearby galaxy clusters grow mostly by dry mergers with cluster members, leading to no increase in AGN activity, BCGs at high redshift could primarily merge with gas-rich satellites, providing fuel for feeding AGNs. If this observed increase in AGN activity is linked to gas-rich mergers rather than ICM cooling, we would expect to see an increase in scatter in the P cav versus L cool relation at z > 1. Last, this work confirms that the runaway cooling phase, as predicted by the classical cooling-flow model, in the Phoenix cluster is extremely rare and most BCGs have low (relative to Eddington) black hole accretion rates.

79 ASTRONOMY AND ASTROPHYSICS↗