Search NASA⌕ Search

SEARCH · Search NASA

Results for “Model fitting”

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 649 records · Page 36

The clustering of DESI-like luminous red galaxies using photometric redshifts

ABSTRACT We present measurements of the redshift-dependent clustering of a DESI-like luminous red galaxy (LRG) sample selected from the Legacy Survey imaging data set, and use the halo occupation distribution (HOD) framework to fit the clustering signal. The photometric LRG sample in this study contains 2.7 million objects over the redshift range of 0.4 < z < 0.9 over 5655 deg2. We have developed new photometric redshift (photo-z) estimates using the Legacy Survey DECam and WISE photometry, with σNMAD = 0.02 precision for LRGs. We compute the projected correlation function using new methods that maximize signal-to-noise ratio while incorporating redshift uncertainties. We present a novel algorithm for dividing irregular survey geometries into equal-area patches for jackknife resampling. For a five-parameter HOD model fit using the MultiDark halo catalogue, we find that there is little evolution in HOD parameters except at the highest redshifts. The inferred large-scale structure bias is largely consistent with constant clustering amplitude over time. In an appendix, we explore limitations of Markov chain Monte Carlo fitting using stochastic likelihood estimates resulting from applying HOD methods to N-body catalogues, and present a new technique for finding best-fitting parameters in this situation. Accompanying this paper, we have released the Photometric Redshifts for the Legacy Surveys catalogue of photo-z’s obtained by applying the methods used in this work to the full Legacy Survey Data Release 8 data set. This catalogue provides accurate photometric redshifts for objects with z < 21 over more than 16 000 deg2 of sky.

79 ASTRONOMY AND ASTROPHYSICS↗

Effects of Catalytic Action and Ligand Binding on Conformational Ensembles of Adenylate Kinase

Crystal structures of adenylate kinase (AdK) from Escherichia coli capture two states: an “open” conformation (apo) obtained in the absence of ligands and a “closed” conformation in which ligands are bound. Other AdK crystal structures suggest intermediate conformations that may lie on the transition pathway between these two states. To characterize the transition from open to closed states in solution, X-ray solution scattering data were collected from AdK in the apo form and with progressively increasing concentrations of five different ligands. Scattering data from apo AdK are consistent with scattering predicted from the crystal structure of AdK in the open conformation. In contrast, data from AdK samples saturated with Ap5A do not agree with that calculated from AdK in the closed conformation. Using cluster analysis of available structures, we selected representative structures in five conformational states: open, partially open, intermediate, partially closed, and closed. We used these structures to estimate the relative abundances of these states for each experimental condition. X-ray solution scattering data obtained from AdK with AMP are dominated by scattering from AdK in the open conformation. For AdK in the presence of high concentrations of ATP and ADP, the conformational ensemble shifts to a mixture of partially open and closed states. Even when AdK is saturated with Ap5A, a significant proportion of AdK remains in a partially open conformation. Here, these results are consistent with an induced-fit model in which the transition of AdK from an open state to a closed state is initiated by ATP binding.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Evaluation of CMIP6 Streamflow in the Arctic

Earth system models (ESMs) serve as the primary basis for projecting future streamflow changes, but they are biased in terms of the their ability to reproduce historical observations of streamflow. Rigorous downscaling and bias-correction procedures, which are time-consuming and introduce uncertainties, are commonly used to address these biases in streamflow or runoff. However, the applicable limits of streamflow projections remain largely untested through direct comparisons with observations. In this work, we compare historical streamflow time series observations with a suite of Coupled Model Intercomparison Project phase 6 (CMIP6) ESMs to assess the space and time scales at which each ESM represents historic streamflow observations. We focus our analysis on streamflow across the Arctic given the higher degree of warming experienced over this region, underscoring its importance for understanding the future changes. A series of metrics—volume, seasonality, extreme capturing streamflow variability, drought, and flooding is tested to assess ESM simulation skill for total events and overall distributions. Our results indicate that although improvements are necessary in ESM runoff and streamflow projections, the use of temporal averaging and targeting the best-fit model for comparison to observations when a suite of ESMs is considered can provide robust streamflow projections. Moreover, this work provides a basis for contextualizing the limits of ESM applicability for studies aiming to assess streamflow in the future.

54 ENVIRONMENTAL SCIENCES↗

A linear, separable two-parameter model for dual energy CT imaging of proton stopping power computation

Purpose: To evaluate the accuracy and robustness of a simple, linear, separable, two-parameter model (basis vector model, BVM) in mapping proton stopping powers via dual energy computed tomography (DECT) imaging. Methods: The BVM assumes that photon cross sections (attenuation coefficients) of unknown materials are linear combinations of the corresponding radiological quantities of dissimilar basis substances (i.e., polystyrene, CaCl{sub 2} aqueous solution, and water). The authors have extended this approach to the estimation of electron density and mean excitation energy, which are required parameters for computing proton stopping powers via the Bethe–Bloch equation. The authors compared the stopping power estimation accuracy of the BVM with that of a nonlinear, nonseparable photon cross section Torikoshi parametric fit model (VCU tPFM) as implemented by the authors and by Yang et al. [“Theoretical variance analysis of single- and dual-energy computed tomography methods for calculating proton stopping power ratios of biological tissues,” Phys. Med. Biol. 55, 1343–1362 (2010)]. Using an idealized monoenergetic DECT imaging model, proton ranges estimated by the BVM, VCU tPFM, and Yang tPFM were compared to International Commission on Radiation Units and Measurements (ICRU) published reference values. The robustness of the stopping power prediction accuracy of tissue composition variations was assessed for both of the BVM and VCU tPFM. The sensitivity of accuracy to CT image uncertainty was also evaluated. Results: Based on the authors’ idealized, error-free DECT imaging model, the root-mean-square error of BVM proton stopping power estimation for 175 MeV protons relative to ICRU reference values for 34 ICRU standard tissues is 0.20%, compared to 0.23% and 0.68% for the Yang and VCU tPFM models, respectively. The range estimation errors were less than 1 mm for the BVM and Yang tPFM models, respectively. The BVM estimation accuracy is not dependent on tissue type and proton energy range. The BVM is slightly more vulnerable to CT image intensity uncertainties than the tPFM models. Both the BVM and tPFM prediction accuracies were robust to uncertainties of tissue composition and independent of the choice of reference values. This reported accuracy does not include the impacts of I-value uncertainties and imaging artifacts and may not be achievable on current clinical CT scanners. Conclusions: The proton stopping power estimation accuracy of the proposed linear, separable BVM model is comparable to or better than that of the nonseparable tPFM models proposed by other groups. In contrast to the tPFM, the BVM does not require an iterative solving for effective atomic number and electron density at every voxel; this improves the computational efficiency of DECT imaging when iterative, model-based image reconstruction algorithms are used to minimize noise and systematic imaging artifacts of CT images.

61 RADIATION PROTECTION AND DOSIMETRY↗

The Swift Burst Alert Telescope Detected Seyfert 1 Galaxies: X-Ray Broadband Properties and Warm Absorbers

We present results from an analysis of the broadband, 0.3-195 keV, X-ray spectra of 48 Seyfert 1-1.5 sources detected in the very hard X-rays with the Swift Burst Alert Telescope (BAT). This sample is selected in an all-sky survey conducted in the 14-195 keV band. Therefore, our sources are largely unbiased toward both obscuration and host galaxy properties. Our detailed and uniform model fits to Suzaku/BAT and XMM-Newton/BAT spectra include the neutral absorption, direct power-law, reflected emission, soft excess, warm absorption, and narrow Fe I K[alpha] emission properties for the entire sample. We significantly detect O VII and O VIII edges in 52% of our sample. The strength of these detections is strongly correlated with the neutral column density measured in the spectrum. Among the strongest detections, X-ray grating and UV observations, where available, indicate outflowing material. The ionized column densities of sources with O VII and O VIII detections are clustered in a narrow range with Nwarm [approx] 1021 cm-2, while sources without strong detections have column densities of ionized gas an order of magnitude lower. Therefore, we note that sources without strong detections likely have warm ionized outflows present but at low column densities that are not easily probed with current X-ray observations. Sources with strong complex absorption have a strong soft excess, which may or may not be due to difficulties in modeling the complex spectra of these sources. Still, the detection of a flat [Gamma] [approx] 1 and a strong soft excess may allow us to infer the presence of strong absorption in low signal-to-noise active galactic nucleus spectra. Additionally, we include a useful correction from the Swift BAT luminosity to bolometric luminosity, based on a comparison of our spectral fitting results with published spectral energy distribution fits from 33 of our sources.

Winter, Lisa M.↗

OBSERVATIONS OF INFALLING AND ROTATIONAL MOTIONS ON A 1000 AU SCALE AROUND 17 CLASS 0 AND 0/I PROTOSTARS: HINTS OF DISK GROWTH AND MAGNETIC BRAKING?

We perform imaging and analyses of SMA 1.3 mm continuum, C{sup 18}O (2-1) and {sup 12}CO (2-1) line data of 17 Class 0 and 0/I protostars to study their gas kinematics on a 1000 AU scale. Continuum and C{sup 18}O (2-1) emission are detected toward all the sample sources and show central primary components with sizes of ∼600-1500 AU associated with protostars. The velocity gradients in C{sup 18}O (2-1) have wide ranges of orientations from parallel to perpendicular to the outflows, with magnitudes from ∼1 to ∼530 km s{sup –1} pc{sup –1}. We construct a simple kinematic model to reproduce the observed velocity gradients, estimate the infalling and rotational velocities, and infer the disk radii and the protostellar masses. The inferred disk radii range from <5 AU to >500 AU with estimated protostellar masses from <0.1 M {sub ☉} to >1 M {sub ☉}. Our results hint that both large and small disks are possibly present around Class 0 protostars, which could be a sign of disk growth at the Class 0 stage. In addition, the directions of the overall velocity gradients in 7 out of the 17 sources are close to perpendicular to their outflow axes (Δθ > 65°), which is a signature of significant rotational motions. From our model fitting, the specific angular momenta in these sources are estimated to be >2 × 10{sup –4} km s{sup –1} pc, suggesting that magnetic braking is unlikely efficient on a 1000 AU scale in these Class 0 and 0/I sources. In a sub-sample with observed magnetic field orientations, we find no source with large specific angular momenta together with closely aligned magnetic field and outflow axes. This possibly hints that the magnetic field, if originally aligned with the rotational axis, can play a role in removing angular momentum from infalling material at the Class 0 stage. We discuss our results in comparison with theoretical models of collapsing dense cores with and without magnetic fields in the context of disk formation.

79 ASTRONOMY AND ASTROPHYSICS↗

The European Southern Observatory-MIDAS table file system

The new and substantially upgraded version of the Table File System in MIDAS is presented as a scientific database system. MIDAS applications for performing database operations on tables are discussed, for instance, the exchange of the data to and from the TFS, the selection of objects, the uncertainty joins across tables, and the graphical representation of data. This upgraded version of the TFS is a full implementation of the binary table extension of the FITS format; in addition, it also supports arrays of strings. Different storage strategies for optimal access of very large data sets are implemented and are addressed in detail. As a simple relational database, the TFS may be used for the management of personal data files. This opens the way to intelligent pipeline processing of large amounts of data. One of the key features of the Table File System is to provide also an extensive set of tools for the analysis of the final results of a reduction process. Column operations using standard and special mathematical functions as well as statistical distributions can be carried out; commands for linear regression and model fitting using nonlinear least square methods and user-defined functions are available. Finally, statistical tests of hypothesis and multivariate methods can also operate on tables.

Peron, M.↗

Statistical studies of impulsive events at high latitudes

A statistical study has been made of the high-latitude impulsive events that were observed during the 1985-1986 South Pole Balloon Campaign. The events were selected by searching for unipolar pulses greater than or equal to 10 nT above background in the vertical component of the magnetic field on the ground and/or pedestal or 'W' shaped horizontal electric field perturbations greater than or equal to 10 mV/m in amplitude and accompanied by perturbations in the vertical electric field at balloon altitude. A main event list comprising 112 events was compiled from the 468 hours of data available. Three aspects of the events were examined: the solar wind conditions prior to the event, local time of observation, and intrinsic properties of the events. The local time distribution was obtained from the 112 entry main event list and was found to be nearly uniform across the dayside, with no midday gap. The event rate found using our low-amplitude selection criteria was 0.7 event/hr, comparable to expectations based on in situ studies of the magnetopause. A total of 42 events were found for which data were available from Interplanetary Monitoring Platform (IMP) 8. Of these events, 12 occurred when the Z(sub GSM) component (B(sub Z)) of the interplanetary magnetic field (IMF) was northward and 30 occurred when B(sub Z) was southward or fluctuating. Only three of the B(sub Z) northward cases and only five of the B(sub Z) southward cases were preceded by pressure pulses greater than 0.4 nPa in amplitude. Ten of the events were studied in detail by means of a model-fitting method discussed elsewhere. This method infers values of several parameters, including the total current flowing in a coaxial or monopole system and a two-dimensional dipole system. The intrinsic properties of the events showed that only approximately 10% of the total current contributed to momentum transfer to the high-latitude ionosphere, that the direction of the motion depended more on local time of observation than IMF B(sub y), and that events were usually several hundred kilometers in size. The observed B(sub z) control found in the 42 event list and the prevalence of coaxial current dominated events are inconsistent with the predictions of the pressure pulse model.

Lin, Z. M.↗

Interferometry of Saturn and its rings at 1.30-cm wavelength

Interferometric observations of Saturn and its rings at a wavelength of 1.30 cm are presented in an attempt to place constraints upon the amount of thermal radiation emitted by the ring particles. Model-fitting and aperture synthesis techniques were used to analyze the data obtained on nine baselines at a frequency of 23 GHz. Ring optical depth is found to be close to that observed at visible wavelengths, while ring brightness temperature is only 7 + or - 1 K, requiring the ring particles to be nearly conservative scatterers at this wavelength and implying an upper limit of 2.4 m to the radius of a typical ring particle with a lower limit of 0.95 to its single scattering albedo. An observed difference between planetary radii observed at 1.30 and 3.71 cm is interpreted in terms of limb darkening and found to be marginally different from the predictions of atmospheric models in which NH3 is the principal source of microwave opacity.

Schloerb, F. P.↗

Cosmological inference from an emulator based halo model. II. Joint analysis of galaxy-galaxy weak lensing and galaxy clustering from HSC-Y1 and SDSS

Here, we present high-fidelity cosmology results from a blinded joint analysis of galaxy-galaxy weak lensing (ΔΣ) and projected galaxy clustering (w p ) measured from the Hyper Suprime-Cam Year-1 (HSC-Y1) data and spectroscopic Sloan Digital Sky Survey (SDSS) galaxy catalogs in the redshift range 0.15 < z < 0.7. We define luminosity-limited samples of SDSS galaxies to serve as the tracers of w p in three spectroscopic redshift bins, and as the lens samples for ΔΣ. For the ΔΣ measurements, we select a single sample of 4×10 6 source galaxies over 140 deg 2 from HSC-Y1 with photometric redshifts (photo z) greater than 0.75, enabling a better handle of photo- z errors by comparing the ΔΣ amplitudes for the three lens redshift bins. The deep, high-quality HSC-Y1 data enable significant detections of the ΔΣ signals, with integrated signal-to-noise ratio S/N ~ 15 in the range 3 ≤ R/[h –1 Mpc] ≤ 30 for the three lens samples, despite the small area coverage. For cosmological parameter inference, we use an input galaxy-halo connection model built on the dark emulator package (which uses an ensemble set of high-resolution N-body simulations and enables fast, accurate computation of the clustering observables) with a halo occupation distribution that includes nuisance parameters to marginalize over modeling uncertainties. We model the ΔΣ and wp measurements on scales from R≃3 and 2h –1 Mpc , respectively, up to 30 h –1 Mpc (therefore excluding the baryon acoustic oscillations information) assuming a flat Λ CDM cosmology, marginalizing over about 20 nuisance parameters and demonstrating the robustness of our results to them. With various tests using mock catalogs described in Miyatake et al. [preceding paper, Phys. Rev. D 106, 083519 (2022)], we show that any bias in the clustering amplitude S 8 ≡ σ 8 (Ω m /0.3) 0.5 due to uncertainties in the galaxy-halo connection is less than ~ 50 % of the statistical uncertainty of S 8 , unless the assembly biaseffect is unexpectedly large. Our best-fit models have S 8 = 0.795$_{-0.042}^{+0.049}$ (mode and 68% credible interval) for the flat Λ CDM model; we find tighter constraints on the quantity S 8 (α = 0.17)≡σ 8 (Ω m /0.3) 0.17 = 0.745$_{-0.031}^{+0.039}$.

79 ASTRONOMY AND ASTROPHYSICS↗

"Efficiency Space" - A Framework for Evaluating Joint Evaporation and Runoff Behavior

At the land surface, higher soil moisture levels generally lead to both increased evaporation for a given amount of incoming radiation (increased evaporation efficiency) and increased runoff for a given amount of precipitation (increased runoff efficiency). Evaporation efficiency and runoff efficiency can thus be said to vary with each other, motivating the development of a unique hydroclimatic analysis framework. Using a simple water balance model fitted, in different experiments, with a wide variety of functional forms for evaporation and runoff efficiency, we transform net radiation and precipitation fields into fields of streamflow that can be directly evaluated against observations. The optimal combination of the functional forms the combination that produces the most skillful stream-flow simulations provides an indication for how evaporation and runoff efficiencies vary with each other in nature, a relationship that can be said to define the overall character of land surface hydrological processes, at least to first order. The inferred optimal relationship is represented herein as a curve in efficiency space and should be valuable for the evaluation and development of GCM-based land surface models, which by this measure are often found to be suboptimal.

GCM↗

Multistage Core Formation in Planetesimals Revealed by Numerical Modeling and Hf-W Chronometry of Iron Meteorites

Iron meteorites provide some of the most direct insights into the processes and timescales of core formation in planetesimals. Of these, group IVB irons stand out by having one of the youngest 182 Hf- 182 W model ages for metal segregation (2.9 ± 0.6 Ma after solar system formation), as well as the lowest bulk sulfur content and hence highest liquidus temperature. Here in this paper, using a new model for the internal evolution of the IVB parent body, we show that a single stage of metal-silicate separation cannot account for the complete melting of pure Fe metal at the relatively late time given by the Hf-W model age. Instead, a complex metal-silicate separation scenario is required that includes migration of partial silicate melts, formation of a shallow magma ocean, and core formation in two distinct stages of metal segregation. In the first stage, a protocore formed at ≈1.5 Ma via settling of metal particles in a mantle magma ocean, followed by metal segregation from a shallow magma ocean at ≈5.4 Ma. As these stages of metal segregation occurred at different times, the two metal fractions had different 182 W compositions. Consequently, the final 182W composition of the IVB core does not correspond to a single differentiation event, but represents the average composition of early- and late-segregated core fractions. Our best fit model indicates an ≈100 km radius for the IVB parent body and provides an accretion age of ≈0.1–0.5 Ma after solar system formation. The computed solidification time is, furthermore, consistent with the Re-Os age for crystallization of the IVB core.

58 GEOSCIENCES↗

A Coupled Mathematical Model of the Intracellular Replication of Dengue Virus and the Host Cell Immune Response to Infection

Dengue virus (DV) is a positive-strand RNA virus of the Flavivirus genus. It is one of the most prevalent mosquito-borne viruses, infecting globally 390 million individuals per year. The clinical spectrum of DV infection ranges from an asymptomatic course to severe complications such as dengue hemorrhagic fever (DHF) and dengue shock syndrome (DSS), the latter because of severe plasma leakage. Given that the outcome of infection is likely determined by the kinetics of viral replication and the antiviral host cell immune response (HIR) it is of importance to understand the interaction between these two parameters. In this study, we use mathematical modeling to characterize and understand the complex interplay between intracellular DV replication and the host cells' defense mechanisms. We first measured viral RNA, viral protein, and virus particle production in Huh7 cells, which exhibit a notoriously weak intrinsic antiviral response. Based on these measurements, we developed a detailed intracellular DV replication model. We then measured replication in IFN competent A549 cells and used this data to couple the replication model with a model describing IFN activation and production of IFN stimulated genes (ISGs), as well as their interplay with DV replication. By comparing the cell line specific DV replication, we found that host factors involved in replication complex formation and virus particle production are crucial for replication efficiency. Regarding possible modes of action of the HIR, our model fits suggest that the HIR mainly affects DV RNA translation initiation, cytosolic DV RNA degradation, and naïve cell infection. We further analyzed the potential of direct acting antiviral drugs targeting different processes of the DV lifecycle in silico and found that targeting RNA synthesis and virus assembly and release are the most promising anti-DV drug targets.

59 BASIC BIOLOGICAL SCIENCES↗

Data-guided Multi-Map variables for ensemble refinement of molecular movies

Driving molecular dynamics simulations with data-guided collective variables offer a promising strategy to recover thermodynamic information from structure-centric experiments. In this study, the three-dimensional electron density of a protein, as it would be determined by cryo-EM or x-ray crystallography, is used to achieve simultaneously free-energy costs of conformational transitions and refined atomic structures. Unlike previous density-driven molecular dynamics methodologies that determine only the best map-model fits, our work employs the recently developed Multi-Map methodology to monitor concerted movements within equilibrium, non-equilibrium, and enhanced sampling simulations. Construction of all-atom ensembles along the chosen values of the Multi-Map variable enables simultaneous estimation of average properties, as well as real-space refinement of the structures contributing to such averages. Using three proteins of increasing size, we demonstrate that biased simulation along the reaction coordinates derived from electron densities can capture conformational transitions between known intermediates. The simulated pathways appear reversible with minimal hysteresis and require only low-resolution density information to guide the transition. The induced transitions also produce estimates for free energy differences that can be directly compared to experimental observables and population distributions. The refined model quality is superior compared to those found in the Protein Data Bank. We find that the best quantitative agreement with experimental free-energy differences is obtained using medium resolution density information coupled to comparatively large structural transitions. Practical considerations for probing the transitions between multiple intermediate density states are also discussed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Constraints on the Emission and Viewing Geometry of the Transient Anomalous X-ray Pulsar XTE J1810-197

The temporal decay of the flux components of the transient anomalous X-ray pulsar XTE J1 810-197 following its 2002 outburst presents a unique opportunity to probe the emission geometry of a magnetar. Toward this goal, we model the magnitude of the pulsar's modulation in narrow spectral bands over time. Following previous work, we assume that the postoutburst flux is produced in two distinct thermal components arising from a hot spot and a warm concentric ring. We include general relativistic effects on the blackbody spectra due to gravitational redshift and light bending near the stellar surface, which strongly depend on radius. This affects the model fits for the temperature and size of the emission regions. For the hot spot, the observed temporal and energy-dependent pulse modulation is found to require an anisotropic, pencil-beamed radiation pattern. We are able to constrain an allowed range for the angles that the line of sight (psi) and the hot spot pole (xi) make with respect to the spin axis. Within errors, this is defined by the locus of points in the xi-psi plane that lie along the line [xi + beta(R)] [psi + [beta(R)] = const, where beta(R) is a function of the radius R of the star. For a canonical value of R = 12 km, the viewing parameters range from psi = xi = 37deg to (psi, xi) = (85deg, 15deg). We discuss our results in the context of magnetar emission models.

Perna, Rosalba↗

Constraining the properties of gaseous halos via cross-correlations of upcoming galaxy surveys and thermal Sunyaev-Zel’dovich maps

The thermal Sunyaev-Zel’dovich (tSZ) effect induces a Compton-$\mathit{y}$ distortion in cosmic microwave background (CMB) temperature maps that is sensitive to a line-of-sight integral of the ionized gas pressure. By correlating the positions of galaxies with maps of the Compton-$\mathit{y}$ distortion, one can probe baryonic feedback processes and study the thermodynamic properties of a significant fraction of the gas in the Universe. Using a model fitting approach, we forecast how well future galaxy and CMB surveys will be able to measure these correlations, and show that powerful constraints on halo pressure profiles can be obtained. In this work, our forecasts are focused on correlations between galaxies and halos identified by the upcoming Dark Energy Spectroscopic Instrument survey and tSZ maps from the Simons Observatory and CMB-S4 experiments, but have general applicability to other surveys, such as the Large Synoptic Survey Telescope. We include prescriptions for observational systematics, such as halo miscentering and halo mass bias, demonstrating several important degeneracies with pressure profile parameters. Assuming modest priors on these systematics, we find that measurements of halo-$\mathit{y}$ and galaxy-$\mathit{y}$ correlations with future surveys yield tight constraints on the pressure profiles of group-scale dark matter halos, and enable current feedback models to either be confirmed or ruled out.

79 ASTRONOMY AND ASTROPHYSICS↗

A Qualitative Strategy for Fusion of Physics into Empirical Models for Process Anomaly Detection

To facilitate the automated online monitoring of power plants, a systematic and qualitative strategy for anomaly detection is presented. This strategy is essential to provide credible reasoning on why and when an empirical versus hybrid (i.e., physics-supported) approach should be used and to determine the ideal mix of these two approaches for a defined anomaly detection scope. Empirical methods are usually based on pattern, statistical, and causal inference. Hybrid methods include the use of physics models to train and test data methods, reduce data dimensionality, reduce data-model complexity, augment data, and reduce empirical uncertainty; hybrid methods also include the use of data to tune physics models. The presented strategy is driven by key decision points related to data relevance, simple modeling feasibility, data inference, physics-modeling value, data dimensionality, physics knowledge, method of validation, performance, data availability, and suitability for training and testing, cause-effect, entropy inference, and model fitting. The strategy is demonstrated through a pilot use case for the application of anomaly detection to capture a valve packing leak at the high-pressure coolant injection system of a nuclear power plant.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

A Weakening Compton Hump and Soft X-ray Excess Detected in the Seyfert-1 Galaxy MCG –02–58–22

We have carried out an extensive X-ray spectral study of the bare Seyfert-1 galaxy MCG –02–58–22to ascertain the nature of the X-ray reprocessing media, using observations from Suzaku(2009) and simultaneous observations from XMM-Newton and NuSTAR(2016) . The most significant results of our investigation are: 1. The primary X-ray emission from the corona is constant in these observations, both in terms of the power law slope (Γ = 1.80) and luminosity (L2−10 keV= 2.55×1044erg s−1). 2.The soft excess flux decreased by a factor of two in 2016, the Compton hump weakened/vanished in2016, and the narrow FeKαemission line became marginally broad (σ= 0.35±0.08 keV) and its flux doubled in 2016. 3. From physical model fits we find that the normalization of the narrow component of the FeKαline does not change in the two epochs, although the Compton hump vanishes in the same time span. Since the primary X-ray continuum does not change, we presume that any changes in there processed emission must arise due to changes in the reprocessing media. Our primary conclusions are: A. The vanishing of the Compton hump in 2016 can probably be explained by a dynamic clumpytorus which is infalling/outflowing, or by a polar torus wind. B. The torus in this AGN possibly has two structures: an equatorial toroidal disk (producing the narrow FeKαemission) and a polar component(producing the variable Compton hump), C. The reduction of the soft-excess flux by half and increase in the FeKαflux by a factor of two in the same period cannot be adequately explained by ionized disk reflection model alone.

Sibasish Laha↗