Search NASA⌕ Search

SEARCH · Search NASA

Results for “POINT MATCHING METHOD”

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 37 records · Page 2

Evaluating pulse-shaping capabilities of next-generation pulsed power architectures

This project evaluated the pulse shaping capabilities of next-generation pulsed power (NGPP) architectures. NGPP architectures share several common attributes including multiple independent pulse-generation lines, a radial water-insulated impedance transformer, and a central vacuum insulated load region. A multi-module circuit model was developed, incorporating independent pulse-generation lines and a 2-D transmission line mesh of the radial impedance transformer to assess the effects of azimuthal asymmetry in pulse-shaped experiments. Circuit model simulations demonstrated that NGPP architectures are able to produce the the desired current pulse shapes for exemplar NGPP experiments. Additionally, the project explored automated methods for experiment design, including derivative -ree optimization and machine learning. Pulse-shaped experiments require designers to determine machine parameters that reliably produce the desired current pulse at the load, a process that typically relies on expert knowledge and iterative adjustments using the Z circuit model. Given the increased complexity of NGPP systems, this manual approach may be impractical. While the evaluated methods do not eliminate the need for manual iteration, they can reduce the time required for experiment design. Derivative-free optimization automates much of the trial-and-error process, providing a close starting point for manual adjustments or making small modifications to near-final designs. Meanwhile, deep neural network methods can generate a good qualitative match to the desired current pulse in under one second without requiring circuit model simulations.

42 ENGINEERING↗

Suppressing the sample variance of DESI-like galaxy clustering with fast simulations

Ongoing and upcoming galaxy redshift surveys, such as the Dark Energy Spectroscopic Instrument (DESI) survey, will observe vast regions of sky and a wide range of redshifts. In order to model the observations and address various systematic uncertainties, N-body simulations are routinely adopted, however, the number of large simulations with sufficiently high mass resolution is usually limited by available computing time. Therefore, achieving a simulation volume with the effective statistical errors significantly smaller than those of the observations becomes prohibitively expensive. In this study, we apply the Convergence Acceleration by Regression and Pooling (CARPool) method to mitigate the sample variance of the DESI-like galaxy clustering in the AbacusSummit simulations, with the assistance of the quasi-N-body simulations FastPM. Based on the halo occupation distribution (HOD) models, we construct different FastPM galaxy catalogs, including the luminous red galaxies (LRGs), emission line galaxies (ELGs), and quasars, with their number densities and two-point clustering statistics well matched to those of AbacusSummit. We also employ the same initial conditions between AbacusSummit and FastPM to achieve high cross-correlation, as it is useful in effectively suppressing the variance. Our method of reducing noise in clustering is equivalent to performing a simulation with volume larger by a factor of 5 and 4 for LRGs and ELGs, respectively. We also mitigate the standard deviation of the LRG bispectrum with the triangular configurations k 2 = 2k 1 = 0.2 h Mpc -1 by a factor of 1.6. With smaller sample variance on galaxy clustering, we are able to constrain the baryon acoustic oscillations (BAO) scale parameters to higher precision. The CARPool method will be beneficial to better constrain the theoretical systematics of BAO, redshift space distortions (RSD) and primordial non-Gaussianity (NG).

79 ASTRONOMY AND ASTROPHYSICS↗

Score-based denoising for atomic structure identification

We propose an effective method for removing thermal vibrations that complicate the task of analyzing complex dynamics in atomistic simulation of condensed matter. Our method iteratively subtracts thermal noises or perturbations in atomic positions using a denoising score function trained on synthetically noised but otherwise perfect crystal lattices. The resulting denoised structures clearly reveal underlying crystal order while retaining disorder associated with crystal defects. Purely geometric, agnostic to interatomic potentials, and trained without inputs from explicit simulations, our denoiser can be applied to simulation data generated from vastly different interatomic interactions. The denoiser is shown to improve existing classification methods, such as common neighbor analysis and polyhedral template matching, reaching perfect classification accuracy on a recent benchmark dataset of thermally perturbed structures up to the melting point. Demonstrated here in a wide variety of atomistic simulation contexts, the denoiser is general, robust, and readily extendable to delineate order from disorder in structurally and chemically complex materials.

36 MATERIALS SCIENCE↗

Technical and Economic Evaluation of the First Ever Polymer Flood Field Pilot to Enhance the Recovery of Heavy Oils on Alaska's North Slope via Machine Assisted History Matching

Polymer flooding has become globally established as a potential enhanced oil recovery method for heavy oils. To determine whether this technology may be useful in developing the substantial heavy oil resources on the Alaska North Slope, a polymer flood field pilot commenced at the Milne Point Unit in August 2018. This study seeks to evaluate the results of the field pilot on a technical and economic basis. A reservoir simulation model is constructed and calibrated to predict the oil recovery performance of the pilot through machine-assisted reservoir simulation techniques. To replicate the early water breakthrough observed during waterflooding, transmissibility contrasts are introduced into the simulation model, forcing viscous fingering effects. In the ensuing polymer flood, these transmissibility contrasts are reduced to replicate the restoration of injection conformance during polymer flooding. Transmissibility contrasts are later reinstated to replicate fracture overextension interpreted in one of the producing wells. The calibrated simulation models produced at each stage of the history matching process are used to forecast oil recovery. These forecasts are used as input for economic analysis, incremental to waterflooding expectations. The simulation forecasts indicate that polymer flooding significantly increases the heavy oil production for this field pilot compared to waterflooding alone, yielding attractive project economics. However, meaningful variations between simulation scenarios demonstrate that a simulation model is only valid for prediction if flow behavior in the reservoir remains consistent with that observed during the history matched period. Critically, this means that a simulation model calibrated for waterflooding may not fully capture the technical and economic benefits of an enhanced oil recovery process such as polymer flooding. Subsequently, the simulation model and economic model are used in conjunction to conduct a sensitivity analysis for polymer flood design parameters, from which recommendations are provided for both the continued operation of the current field pilot and future polymer flood designs. The results demonstrate that a higher polymer concentration can be injected due to the development of fractures in the reservoir. The throughput rate should remain high without exceeding operating constraints. A calculated point-forward polymer utilization parameter demonstrates the decreasing efficiency of the polymer flood at later times in the pattern life. Future projects will benefit from starting polymer injection earlier in the pattern life. A pattern with tighter horizontal well spacing will observe a greater incremental benefit from polymer flooding.

Keith, Cody↗

Design of protein-binding proteins from the target structure alone

The design of proteins that bind to a specific site on the surface of a target protein using no information other than the three-dimensional structure of the target remains a challenge. Here we describe a general solution to this problem that starts with a broad exploration of the vast space of possible binding modes to a selected region of a protein surface, and then intensifies the search in the vicinity of the most promising binding modes. We demonstrate the broad applicability of this approach through the de novo design of binding proteins to 12 diverse protein targets with different shapes and surface properties. Biophysical characterization shows that the binders, which are all smaller than 65 amino acids, are hyperstable and, following experimental optimization, bind their targets with nanomolar to picomolar affinities. We succeeded in solving crystal structures of five of the binder–target complexes, and all five closely match the corresponding computational design models. Experimental data on nearly half a million computational designs and hundreds of thousands of point mutants provide detailed feedback on the strengths and limitations of the method and of our current understanding of protein–protein interactions, and should guide improvements of both. Our approach enables the targeted design of binders to sites of interest on a wide variety of proteins for therapeutic and diagnostic applications.

59 BASIC BIOLOGICAL SCIENCES↗

Target Detection via Cognitive Radars Using Change-Point Detection, Learning, and Adaptation

Many radar detection algorithms that assume a stationary environment (clutter) have been proposed and analyzed over the years. However, in practice, changes in the nonstationary environment can perturb the parameters of the clutter distribution, or even alter the clutter distribution family, which can greatly deteriorate the target detection capability. To avoid such potential performance degradation, cognitive radar systems are envisioned which are required to rapidly realize the nonstationarity, accurately learn the new characteristics of the environments, and adaptively update the detector. In this paper, aiming to develop a fully cognitive radar for target detection in nonstationary environments, we propose a unifying framework that integrates (i) change-point detection of clutter distributions by using a data-driven cumulative sum (CUSUM) algorithm and its extended version, (ii) learning/identification of clutter distribution by applying sparse theory and kernel density estimation methods, and (iii) adaptive target detection by automatically modifying the likelihood-ratio test and corresponding detection threshold. Further, with extensive numerical examples, we demonstrate the achieved improvements in detection performance due to the proposed framework in comparison to a nonadaptive case, an adaptive matched filter (AMF) method, and the clairvoyant case. Herein, we also use Wilcoxon rank-sum tests to evaluate the statistical significance of the performance improvements

42 ENGINEERING↗

Renormalization group improvement of the effective potential: an EFT approach

We apply effective field theory (EFT) methods to compute the renormalization group improved effective potential for theories with a large mass hierarchy. Our method allows one to compute the effective potential in a systematic expansion in powers of the mass ratio, as well as to sum large logarithms of mass ratios using renormalization group evolution. The effective potential is the sum of one-particle irreducible diagrams (1PI) but information about which diagrams are 1PI is lost after matching to the EFT, since heavy lines get shrunk to a point. We therefore introduce a tadpole condition in place of the 1PI condition, and use the renormalization group improved value of the tadpole in computing the effective potential. We explain why the effective potential computed using an EFT is not the same as the effective potential of the EFT. We illustrate our method using the O(N) model, a theory of two scalars in the unbroken and broken phases, and the Higgs-Yukawa model. Our leading-log result, obtained by integrating the one-loop β-functions, correctly reproduces the log-squared term in explicit two-loop calculations. Our method does not have a Goldstone boson infrared divergence problem.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Sensors and methods for monitoring flying objects

Described herein are sensing methods, sensor systems, and non-transitory, computer-readable, storage media having programs for long-duration, continuous monitoring of flying objects during the day or the night and regardless of weather conditions. The methods and systems are computationally efficient and can provide compact, three-dimensional representations of motion from the observed object. A 3D track of the flying object can be generated from a point-matched pair of stereo composite motion track images and not directly from the videos, wherein each composite motion track image is based on a composite of a plurality of video frames composited in part according to video frame numbers.

Matzner, Shari↗

Challenge Solutions for a Significant Water Cut Reduction History Match on Heavy Oil Polymer EOR

Milne Point Field is conducting the first-ever active polymer flood pilot program designed to enhance oil recovery (EOR) of heavy oils on the Alaska North Slope (ANS). The experimental field program consists of two pairs of horizontal injection and production wells deployed in an isolated fault block in the Schrader Bluff reservoir. The EOR method combines polymer flooding and low salinity water flooding to improve recovery. Prior average water cut reduction in polymer EOR fields ordinarily ranges between 10% to 15% for most light oil and heavy oil reservoirs. When comparing other field applications, the long-lasting (24 months) and much lower water cut poses a challenge for numerical history matching when applying typical simulation methods. This abstract discusses two practical new approaches used for history matching of the long-lasting ultralow water cut. These new applications include 1) multiple relative permeability curves with J-function incorporation with low endpoints of relative permeability of water (Krw). This approach integrated the permeability differential adjustment in the geometries of horizontal and vertical directions based on geological review. 2) A specific permeability channel strips along with the horizontal or vertical planar assumption. The current history matching results indicate success in achieving an excellent agreement with the field observed low water-cut.

Namie, S.↗

Total metals & anion concentration data; Slate River floodplain, Crested Butte, CO; May 2020-September 2020

This data package includes processed and undiluted measurements for metal and anion concentrations from pore water (groundwater) samples from the Slate River floodplain of Crested Butte, CO, a focus field site for the SLAC Floodplain Hydro-Biogeochemistry SFA. The data was generated as part of the work targeting the overarching research question for the SLAC SFA: How do ubiquitous subsurface interfaces mediate molecular-scale biogeochemical processes and groundwater quality in floodplains and watersheds? Samples were collected between May and September of 2020. These measurements were all recorded at the Arizona Laboratory for Emerging Contaminants (ALEC) at the University of Arizona located in Tucson, AZ. Groundwater samples were extracted from a network of installed rhizon (Rhizosphere Research Products, part no. 19.60.21F, 0.6 micrometer mesh size) and piezometer wells within the river floodplain. All water samples were shaded from sun exposure during extraction from the subsurface and preserved at 4C until measured at ALEC.Analysis by ICP-MS:Measurements for total metals were made on the Agilent 7700x ICP-MS (for total metals) – Agilent Technologies, Santa Clara, CA.The analytical QA/QC protocol was adapted from US EPA Method 200.8 for analysis by ICP-MS. Calibration standards were prepared from multi-element stock solution (Sigma-Aldrich Multielement standard solution for ICP, St. Louis, MO) using matrix matched to sample solutions (either 2% HCl or HNO3 from AriStar Plus,grade acids from VWR Scientific). Calibration curves include at least 7 points with correlation coefficients > 0.995. The QC protocol includes a continuing calibration blank (CCB), a continuing calibration verification (CCV) solution and at least one quality control sample (QCS) to be analyzed just after calibration and again after every 12 samples and at the completion of the run. The QCS solutions are from an independent source, such as NIST SRM 1643e - Trace Elements in Water, or QCS solutions from High Purity Standards (Charleston, SC). Acceptable QC responses must be between 90 and 110% of the certified value. An internal standard (Rh) is added via on-line addition into the sample line using a mixing tee.Analysis by Ion Chromatography (Anions):The protocol follows Method 4110 in Standard Methods for Examination of Water and Wastewater.The instrument used is the Thermo Scientific Dionex ICS-6000 using AS+AG22 column set for anion analysis with isocratic method using sodium carbonate eluent. Detection is by chemical suppression of eluent conductivity. Quality control solutions and mixed analyte standards purchased from Inorganic Ventures, Christiansburg, VA.All files are in csv format.

54 ENVIRONMENTAL SCIENCES↗

Methods to Evaluate Subcolumn Profiles Based on Two-Point Diagnostics

In atmospheric models, stochastic generation of subgrid-scale profiles or “subcolumns” has been used for a variety of purposes. Such subcolumns can be generated from subgrid probability density functions (PDFs) at different vertical levels, when such PDFs are available. To do so, the generator needs to decide how strongly points should be correlated in the vertical, that is, how much the values should be overlapped. This is sometimes called “PDF overlap.” To assess vertical correlation in a simplified, observable setting, here the vertical correlation of vertical velocity in subcloud layers is examined. Doppler lidar is used to evaluate the vertical profiles of vertical velocity produced by a large-eddy simulation (LES) model and the Subgrid Importance Latin Hypercube Sampler (SILHS) subcolumn generator. In order to diagnose unrealistic features in subcolumn profiles, various statistical diagnostics are examined here, including the bivariate PDF of vertical velocity at two separated points (i.e., altitudes), the two-point velocity correlation, the integral correlation length, the PDF of two-point velocity differences, and the skewness and kurtosis of two-point velocity differences. The profiles produced by LES match lidar well, except that they are too smooth at small scales. The profiles produced by SILHS exhibit sharp jumps from updraft to downdraft that are not observed in the lidar data. To reduce the generation of these unrealistically sharp jumps, the SILHS sampling method is revised. The diagnostics confirm that the revised sampling method reduces the overprediction of sharp jumps.

54 ENVIRONMENTAL SCIENCES↗

Energetic and Spectroscopic Properties of Astrophysically Relevant MgC 4 H Radicals Using High-Level Ab Initio Calculations

Considering the importance of magnesium-bearing hydrocarbon molecules (MgC n H; n = 2, 4, and 6) in the carbon-rich circumstellar envelopes (e.g., IRC+10216), a total of 28 constitutional isomers of MgC 4 H have been theoretically investigated using density functional theory (DFT) and coupled-cluster methods. The zero-point vibrational energy corrected relative energies at the ROCCSD(T)/cc-pCVTZ level of theory reveal that the linear isomer, 1-magnesapent-2,4-diyn-1-yl (1, 2 Σ + ), is the global minimum geometry on the MgC 4 H potential energy surface. The latter has been detected both in the laboratory and in the evolved carbon star, IRC+10216. The calculated spectroscopic data for 1 match well with the experimental observations (error ~ 0.78%) which validates our theoretical methodology. Plausible isomerization processes happening among different isomers are examined using DFT and coupled-cluster methods. CASPT2 calculations have been performed for a few isomers exhibiting multireference characteristics. The second most stable isomer, 1-ethynyl-1λ 3 -magnesacycloprop-2-ene-2,3-diyl (2, 2 A 1 , μ = 2.54 D), is 146 kJ mol –1 higher in energy than 1 and possibly the next promising candidate to be detected in the laboratory or in the interstellar medium in future.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Millimeter-wave Point Sources from the 2500 Square Degree SPT-SZ Survey: Catalog and Population Statistics

In this work, we present a catalog of emissive point sources detected in the SPT-SZ survey, a contiguous 2530- square-degree area surveyed with the South Pole Telescope (SPT) from 2008 - 2011 in three bands centered at 95, 150, and 220GHz. The catalog contains 4845 sources measured at a significance of 4.5σ or greater in at least one band, corresponding to detections above approximately 9.8, 5.8, and 20.4 mJy in 95, 150, and 220 GHz, respectively. Spectral behavior in the SPT bands is used for source classification into two populations based on the underlying physical mechanisms of compact, emissive sources that are bright at millimeter wavelengths: synchrotron radiation from active galactic nuclei and thermal emission from dust. The latter population includes a component of high-redshift sources often referred to as submillimeter galaxies (SMGs). In the relatively bright flux ranges probed by the survey, these sources are expected to be magnified by strong gravitational lensing. The survey also contains sources consistent with protoclusters, groups of dusty galaxies at high redshift undergoing collapse. We cross-match the SPT-SZ catalog with external catalogs at radio, infrared, and X-ray wavelengths and identify available redshift information. The catalog splits into 3980 synchrotron-dominated and 865 dust-dominated sources and we determine a list of 506 SMGs. Ten sources in the catalog are identified as stars. We calculate number counts for the full catalog, and synchrotron and dusty components, using a bootstrap method and compare our measured counts with models. This paper represents the third and final catalog of point sources in the SPT-SZ survey.

79 ASTRONOMY AND ASTROPHYSICS↗

Modelling a Heaving Point-Absorber with a Closed-Loop Control System Using the DualSPHysics Code

The present work addresses the need for an efficient, versatile, accurate and open-source numerical tool to be used during the design stage of wave energy converters (WECs). The device considered here is the heaving point-absorber developed and tested by Sandia National Laboratories. The smoothed particle hydrodynamics (SPH) method, as implemented in DualSPHysics, is proposed since its meshless approach presents some important advantages when simulating floating devices. The dynamics of the power take-off system are also modelled by coupling DualSPHysics with the multi-physics library Project Chrono. A satisfactory matching between experimental and numerical results is obtained for: (i) the heave response of the device when forced via its actuator; (ii) the vertical forces acting on the fixed device under regular waves and; (iii) the heave response of the WEC under the action of both regular waves and the actuator force. This proves the ability of the numerical approach proposed to simulate accurately the fluid–structure interaction along with the WEC’s closed-loop control system. In addition, radiation models built from the experimental and WAMIT results are compared with DualSPHysics by plotting the intrinsic impedance in the frequency domain, showing that the SPH method can be also employed for system identification.

42 ENGINEERING↗

Survey of Surveys: I. The largest compilation of radial velocities for the Galaxy

In the present-day panorama of large spectroscopic surveys, the amount, diversity, and complexity of the available data continuously increase. The overarching goal of studying the formation and evolution of our Galaxy is hampered by the heterogeneity of instruments, selection functions, analysis methods, and measured quantities. Here, we present a comprehensive catalogue, the Survey of Surveys (SoS), built by homogeneously merging the radial velocity (RV) determinations of the largest ground-based spectroscopic surveys to date, such as APOGEE, GALAH, Gaia-ESO, RAVE, and LAMOST, using Gaia as a reference. This pilot study serves to prove the concept and to test the methodology that we plan to apply in the future to the stellar parameters and abundance ratios as well. We have devised a multi-staged procedure that includes: (i) the cross match between Gaia and the spectroscopic surveys using the official Gaia cross-match algorithm, (ii) the normalisation of uncertainties using repeated measurements or the three-cornered hat method, (iii) the cross calibration of the RVs as a function of the main parameters on which depend (magnitude, effective temperature, surface gravity, metallicity, and signal-to-noise ratio) to remove trends and zero point offsets, and (iv) the comparison with external high-resolution samples, such as the Gaia RV standards and the Geneva-Copenhagen survey, to validate the homogenisation procedure and to calibrate the RV zero-point of the SoS catalogue. We provide the largest homogenised RV catalogue to date, containing almost 11 million stars, of which about half come exclusively from Gaia and half in combination with the ground-based surveys. We estimate the accuracy of the RV zero-point to be about 0.16–0.31 km s –1 and the RV precision to be in the range 0.05–1.50 km s –1 depending on the type of star and on its survey provenance. We validate the SoS RVs with open clusters from a high resolution homogeneous samples and provide the systemic velocity of 55 individual open clusters. Additionally, we provide median RVs for 532 clusters recently discovered by Gaia data. The SoS is publicly available and ready to be applied to various research projects, such as the study of star clusters, Galactic archaeology, stellar streams, or the characterisation of planet-hosting stars, to name a few. We also plan to include survey updates and more data sources in future versions of the SoS.

79 ASTRONOMY AND ASTROPHYSICS↗

Accelerating crystal structure determination with iterative AlphaFold prediction

Experimental structure determination can be accelerated with artificial intelligence (AI)-based structure-prediction methods such as AlphaFold . Here, an automatic procedure requiring only sequence information and crystallographic data is presented that uses AlphaFold predictions to produce an electron-density map and a structural model. Iterating through cycles of structure prediction is a key element of this procedure: a predicted model rebuilt in one cycle is used as a template for prediction in the next cycle. This procedure was applied to X-ray data for 215 structures released by the Protein Data Bank in a recent six-month period. In 87% of cases our procedure yielded a model with at least 50% of C α atoms matching those in the deposited models within 2 Å. Predictions from the iterative template-guided prediction procedure were more accurate than those obtained without templates. It is concluded that AlphaFold predictions obtained based on sequence information alone are usually accurate enough to solve the crystallographic phase problem with molecular replacement, and a general strategy for macromolecular structure determination that includes AI-based prediction both as a starting point and as a method of model optimization is suggested.

59 BASIC BIOLOGICAL SCIENCES↗

The DESI One-Percent Survey: Modelling the clustering and halo occupation of all four DESI tracers with U CHUU

We present results from a set of mock lightcones for the DESI One-Percent Survey, created from the UCHUU simulation. This 8 h −3 Gpc 3 N-body simulation comprises 2.1 trillion particles and provides high-resolution dark matter (sub)haloes in the framework of the Planck-based ΛCDM cosmology. Employing the subhalo abundance matching (SHAM) technique, we populated the UCHUU (sub)haloes with all four DESI tracers – Bright Galaxy Survey (BGS), luminous red galaxies (LRGs), emission line galaxies (ELGs), and quasars (QSOs) – to z = 2.1. Our method accounts for redshift evolution as well as the clustering dependence on luminosity and stellar mass. The two-point clustering statistics of the DESI One-Percent Survey generally agree with predictions from UCHUU across scales ranging from 0.3 h −1 Mpc to 100 h −1 Mpc for the BGS and across scales ranging from 5 h −1 Mpc to 100 h −1 Mpc for the other tracers. We observed some differences in clustering statistics that can be attributed to incompleteness of the massive end of the stellar mass function of LRGs, our use of a simplified galaxy-halo connection model for ELGs and QSOs, and cosmic variance. We find that at the high precision of UCHUU, the shape of the halo occupation distribution (HOD) of the BGS and LRG samples is smaller bias values, likely due to cosmic variance. The bias dependence on absolute magnitude, stellar mass, and redshift aligns with that of previous surveys. These results provide DESI with tools to generate high-fidelity lightcones for the remainder of the survey and enhance our understanding of the galaxy-halo connection.

cosmology↗

Path integral factorization and the gravitational effective action

Abstract We discuss the factorization and continuity properties of fields in the Euclidean gravitational path integral with higher dimension operators constructed from powers of the Riemann tensor. We construct the boundary terms corresponding to the microcanonical ensemble and show that the saddle point approximation to the path integral with a quasilocal energy constraint generally yields a saddle point with discontinuous temperature. This extends a previous result for the Euclidean Schwarzschild–de Sitter geometry in Einstein gravity and shows that it is robust against at least some types of quantum corrections from heavy fields. As an application, we compute the entropy of SdS in D = 4 using the BTZ method. Our result matches the entropy calculated using Wald’s formula.

Astronomy & Astrophysics↗