Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Sensitivity Simulation”

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,063 records · Page 59

Micromechanical Surrogate Machine Learning Model for Creep Deformation Modeling

Process variability during the manufacture of gas turbine engine hot section components can significantly affect the material’s resulting microstructure. In casting, for instance, geometric variation within a component (thin sections versus thick sections, radial location) influences cooling rates and the resulting grain size. The high temperature creep response is known to be sensitive to grain size owing to a diffusional creep mechanism which occurs more readily along grain boundaries. Microstructural variation correspondingly drives mechanical behavior which propagates into component scale performance uncertainty. These factors are essential when planning inspection, maintenance, and repair strategies within a reliability framework. These benefits provide opportunities to increase overall energy efficiency through refined margins. Critically, there is an opportunity to bolster existing data-driven reliability models using physics-driven process-structure-property relations. Here we present recent work establishing a framework for evaluating the probabilistic creep performance of high-temperature materials. A novel microstructure-sensitive crystal plasticity finite element model is established that captures both grain boundary and crystallographic deformation effects. The computationally expensive physics model is calibrated using a statistical approach and this high-fidelity model is subsequently used to train a computationally efficient machine learning surrogate model. The surrogate model is essential for sampling a large ensemble of simulated structure-property pair results. The ensemble data are then mined to extract salient trends to be incorporated into a microstructure-sensitive reliability model. The proposed approach represents a novel way to capture microstructure-sensitive trends from physics-based models within a modern reliability framework.

Fernandez-Zelaia, Patxi [ORNL]↗

Uncertainty quantification and sensitivity analysis for SPERT III E-core reactivity measurement benchmarking

The Special Power Excursion Reactor Test (SPERT) III E-core experiment is important because it provides critical data on reactor behavior under significant reactivity insertions, which is essential for validating computer simulations and ensuring the safety of modern light water reactors. Its design similarities to contemporary reactors make it a valuable resource for understanding and mitigating extreme hazards in nuclear operations. The current study details the application of formal parametric uncertainty quantification and sensitivity analysis to a model of the SPERT-III E-core model for zero power reactivity benchmarking. Additionally, the reactivity impact from various modeling assumptions is quantified. Overall, an conservative estimate for an uncertainty in k$_{\text{eff}}$ of $\pm$1257 was observed. A less conservative, more realistic, uncertainty estimate of $\pm$1096 pcm can be justified by the potential for various parametric uncertainties to become negligible when sampled independently across the ~1400 pins in the core. The experimental results fall within both of these uncertainty bounds. Standardized regression coefficient as well as Sobol indices are used to identify the guide tube thicknesses as the primary contributors to the uncertainty in k$_{\text{eff}}$. Overall, this study provides information on how the uncertainties in input parameters and modeling methods impact simulated k$_{\text{eff}}$ values and can be used to aid model building efforts for future code validation with the SPERT-III E-core experiment.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗

Applicability of 1D site response analysis to shallow sedimentary basins: A critical evaluation through physics-based 3D ground motion simulations

One-dimensional site response analysis (1D SRA) remains the standard practice in considering the effect of local soil deposits and predicting site-specific ground motions, although its range of applicability to realistic seismic wavefields is still in question. In this 1D approach, horizontal and vertical ground shaking are assumed to be induced by vertically propagating shear and compressional waves, respectively. A recent study based on analytical two-dimensional (2D) plane waves and simple point source earthquake simulations has shown two mechanistic limitations in this 1D modelling technique for general inclined seismic waves, that is, systematic over-prediction of the vertical motion and wave trapping in the 1D soil column. In this article, we evaluate in detail the applicability of this 1D modelling approach to realistic three-dimensional (3D) simulated seismic wavefields in shallow sedimentary basins. Linear-viscoelastic 1D SRA predictions using two types of input motions that are commonly used in practice—rock outcrop and in-column motions, are compared with the reference true site response results from 3D earthquake simulations in terms of various measures in the frequency and time domain. It is shown that the horizontal motion in the 3D seismic wavefield exhibits dominant shear wave propagation phenomenon, while the vertical motion is a combined effect of compressional and shear waves and can be over-predicted by the 1D approach when the incident seismic waves are inclined. Direct evidence of the wave refraction process that leads to the vertical motion over-prediction is provided. 1D SRA with in-column inputs can yield motions that have significantly longer duration compared to the true 3D site response solution due to trapped waves, casting in doubt the frequent need for increased soil damping in existing site studies to compensate for wave attenuation due to scattering alone. Sensitivity investigation on the increase of soil profile damping by a multiplier D mul shows D mul values compatible with those found in the literature for both horizontal and vertical motions. It is shown that the level of D mul optimized for a best match of the spectral acceleration is dependent on the characteristic of the input motion and a larger D mul is typically required for the vertical component. In contrast, 1D SRA with outcrop motions predicts motions with shorter significant duration due to its inability to capture the basin-edge generated surface waves. A suite of ground motion simulations was performed to assess the sensitivity of the observations to the basin geologic structure including the velocity gradient, rock-basin impedance contrast and basin depth. The analysis results show that the accuracy of the simplified 1D procedure is dependent on the wavefield composition of both the input motions and the true 3D site response solution. While the horizontal motions in shallow sedimentary basins can, to the first order, be reasonably captured by the simplified 1D approach, 1D SRA for the vertical component is in general not reliable and contributions from inclined shear waves should be accounted for in site-specific evaluation of the vertical design ground motion.

1D assumption↗

Model Formulations: Integrating Distributed Energy Resources (DER) using Advanced Unit Commitment Models and DER Aggregation Methodologies

A distribution energy resource aggregator (DERA) constitutes a group of distribution energy resources with small generation capacities which meet the threshold to participate in the electricity wholesale market. This document provides the proposed DERA model formulation that will be implemented in the SCUC simulation’s architecture for the SCUC-DER project. Different economical assessment methodologies have been developed to incorporated bids for individual distributed resources, which include solar cost dispatch and cost model, BESS opportunity cost offer algorithm, and price sensitive demand response model. Detailed methods are proposed to aggregate individual cost offers to a DERA cost curve to bid in SCUC market while three methods are proposed to simulate DER actual dispatch. Based on the DERA models in this document, the SCUC-DER project will be able to assess the impacts of DERA on the distribution system’s operation and reliability.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Preliminary Sensitivity Analysis for Sensors Impacts on Building Control Performance

This report describes the preliminary sensitivity analysis for sensor impacts on building control performance through the US Department of Energy’s Oak Ridge National Laboratory’s Flexible Research Platform (FRP-2) building. The rooftop unit system provides cooling and heating to the building. The main heating coil is a gas heating coil. Each zone is served by a variable air volume box with an electricity reheat coil. The rooftop unit and variable air volume box controls adopted the practical control sequences from ASHRAE Guideline 36-2018: High-Performance Sequences of Operation. For sensors, the incipient (time-changing) sensor errors, including bias sensor error and precision sensor error, are the inputs of interest. The outputs are energy consumption and thermal comfort (e.g., the predicted percentage of dissatisfied occupants). The large-scale simulation (3,600 cases) was conducted on a cloud platform by integrating sensor errors and ASHRAE Guideline 36 control sequences into an emulator based on the EnergyPlus simulation program with Python energy management system feature. The surrogate models were developed based on cloud simulation results. The uncertainty analysis showed that the sensor errors substantially affect building energy consumption and thermal comfort. The sensitivity analysis shows a ranking of sensor error impacts for each interested output item (e.g., cooling energy, reheat coil heating energy, predicted percentage of dissatisfied occupants). In FY 2022, sensor locations, types, and costs will be evaluated. The field test in Oak Ridge National Laboratory’s Flexible Research Platform building regarding sensor impacts will also be performed. Finally, a comparative analysis will be conducted based on the field test results and emulator results.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Cosmogenic 137 Xe Analysis of Alternatives for nEXO

Neutrinoless double-beta decay is a hypothetical nuclear decay in which two matter particles are produced without corresponding antimatter particles. Observation of neutrinoless double-beta decay would shed light on a potential explanation of matter-antimatter observed in the universe and also indicate the Majorana nature of neutrinos, meaning that the neutrino can act as its own antiparticle. The nEXO experiment is a proposed 5-ton liquid xenon time projection chamber detector with a projected 1.35 × 10 28 yr half-life sensitivity to neutrinoless double-beta decay in 136 Xe [J. Phys. G 49, 015104 (2022)]. While nEXO has a baseline location of SNOLAB, there are several other underground facilities in the world capable of hosting the experiment. The primary change at each site will be cosmogenic 137 Xe backgrounds, caused by varied cosmogenic muon activity at each site. This study finds that while there is a negligible increase in nEXO’s half-life sensitivity at 90% C.L. from improving the cosmogenic 137 Xe veto rejection procedure for SNOLAB, much higher daily muon rates at three other sites will have an adversely negative effect on nEXO’s sensitivity and would likely require further modifications to nEXO’s baseline design to reach nEXO’s scientific goals. Monte Carlo simulations of muons inducing 137 Xe backgrounds were conducted at four underground sites, SNOLAB, SURF, LNGS, and WIPP, each with roughly a factor of 10 times higher daily muon activity than the previous site. There were two primary inputs to these simulations. The first was sample data of cosmogenic muon flight through nEXO’s outer detector, a 12.8m tall and 12.3m wide cylindrical water tank lined with photomultiplier tubes to detect Cerenkov radiation, and all of the secondary interactions produced in the detector. The second input was a parametrization of 137 Xe* production tagging when 137 Xe de-excites and release prompt gammas, as well as accidental tagging caused by neutron capture on 63 Cu, 65 Cu, 19 F, and 1 H throughout the detector. Together, the number of backgrounds relative to the previous 2021 Sensitivity Paper [J. Phys. G 49, 015104 (2022)] was found at each of these four sites across three energy windows, from 400, 1400, and 2600 keV all up to 4600 keV, just above the end point of the prompt gamma spectrum for 137 Xe*. Using level curves for the sensitivity as 137 Xe backgrounds scale relative to baseline calculations, the half-life sensitivity at 90% C.L. for each of trial was found. Even with reduced 137 Xe backgrounds, the results for SNOLAB yielded neglible improvements in nEXO’s predicted sensitivity. On the contrary, the increased backgrounds for the three other considered sites yielded sensitivity estimates below nEXO’s previous projections. This study demonstrates that the shielding from cosmogenic muons at SNOLAB makes it highly preferred even amongst other world leading underground laboratories.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

The Global‐Mean Precipitation Response to CO 2 ‐Induced Warming in CMIP6 Models

Abstract We examine the response of globally averaged precipitation to global warming—the hydrologic sensitivity (HS)—in the Coupled Model Intercomparison Project phase 6 (CMIP6) multi‐model ensemble. Multi‐model mean HS is 2.5% K −1 (ranging from 2.1–3.1% K −1 across models), a modest decrease compared to CMIP5 (where it was 2.6% K −1 ). This new set of simulations is used as an out‐of‐sample test for observational constraints on HS proposed based on CMIP5. The constraint based on clear‐sky shortwave absorption sensitivity to water vapor has weakened, and it is argued that a proposed constraint based on surface low cloud longwave radiative effects does not apply to HS. Finally, while a previously proposed mechanism connecting HS and climate sensitivity via low clouds is present in the CMIP6 ensemble, it is not an important factor for variations in HS. This explains why HS is uncorrelated with climate sensitivity across the CMIP5 and CMIP6 ensembles.

Pendergrass, A. G.↗

Hydrogen Burning of 29 Si and Its Impact on Presolar Stardust Grains from Classical Novae

Abstract Presolar stardust grains found in primitive meteorites are believed to retain the isotopic composition of stellar outflows at the time of grain condensation. Therefore, laboratory measurements of their isotopic ratios represent sensitive probes for investigating open questions related to stellar evolution, stellar explosions, nucleosynthesis, mixing mechanisms, dust formation, and galactic chemical evolution. For a few selected presolar grains, classical novae have been discussed as a potential source. For SiC, silicate, and graphite presolar grains, the association is based on the observation of small N ( 12 C)/ N ( 13 C) and N ( 14 N)/ N ( 15 N) number abundance ratios compared to solar values, and abundance excesses in 30 Si relative to 29 Si, as previously predicted by models of classical novae. We report on a direct measurement of the 29 Si(p, γ ) 30 P reaction, which strongly impacts simulated δ 29 Si values from classical novae. Our new experimental 29 Si(p, γ ) 30 P thermonuclear reaction rate differs from previous results by up to 50% in the classical nova temperature range ( T = 100–400 MK), while the rate uncertainty is reduced by up to a factor of 3. Using our new reaction rate in Monte Carlo reaction network and hydrodynamic simulations of classical novae, we estimate δ 29 Si values with much reduced uncertainties. Our results establish δ 29 Si values measured in presolar grains as a sensitive probe for assessing their classical nova paternity. We also demonstrate that δ 30 Si values from nova simulations are currently not a useful diagnostic tool unless the large uncertainty of the 30 P(p, γ ) 31 S reaction rate can be significantly reduced.

79 ASTRONOMY AND ASTROPHYSICS↗

Randomized low-rank decompositions of nuclear three-body interactions

First-principles simulations of many-fermion systems are commonly limited by the computational requirements of processing large data objects. As a remedy, we propose the use of low-rank approximations of three-body interactions, which are the dominant such limitation in nuclear physics. We introduce a randomized decomposition technique to handle the excessively large matrix dimensions and study the sensitivity of low-rank properties to interaction details. The developed low-rank three-nucleon interactions are benchmarked in ab initio simulations of few- and many-body systems. Exploiting low-rank properties provides a promising route to extend the microscopic description of atomic nuclei to large systems where storage requirements exceed the computational capacities of the most advanced high-performance computing facilities.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Intermodel Spread in the Pattern Effect and Its Contribution to Climate Sensitivity in CMIP5 and CMIP6 Models

Radiative feedbacks depend on the spatial patterns of sea surface temperature (SST) and thus can change over time as SST patterns evolve—the so-called pattern effect. In this study, we investigate intermodel differences in the magnitude of the pattern effect and how these differences contribute to the spread in effective equilibrium climate sensitivity (ECS) within CMIP5 and CMIP6 models. Effective ECS in CMIP5 estimated from 150-yr-long abrupt4×CO2 simulations is on average 10% higher than that estimated from the early portion (first 50 years) of those simulations, which serves as an analog for historical warming; this difference is reduced to 7% on average in CMIP6. The (negative) net radiative feedback weakens over the course of the abrupt4×CO2 simulations in the vast majority of CMIP5 and CMIP6 models, but this weakening is less dramatic on average in CMIP6. For both ensembles, the total variance in the effective ECS is found to be dominated by the spread in radiative response on fast time scales, rather than the spread in feedback changes. Using Green’s functions derived from two AGCMs shows that the spread in feedbacks on fast time scales may be primarily due to differences in atmospheric model physics, whereas the spread in feedback evolution is primarily governed by differences in SST patterns. Intermodel spread in feedback evolution is well explained by differences in the relative warming in the West Pacific warm-pool regions for the CMIP5 models, but this relation fails to explain differences across the CMIP6 models, suggesting that a stronger sensitivity of extratropical clouds to surface warming may also contribute to feedback changes in CMIP6.

54 ENVIRONMENTAL SCIENCES↗

3D simulations of inertial confinement fusion implosions part 1: inline modeling of polarized cross beam energy transfer and subsequent drive anomalies on OMEGA and NIF

Inertial confinement fusion experiments are sensitive to cross-beam energy transfer (CBET), a laser-plasma instability that redistributes laser energy in the coronal plasma through self-generated ion acoustic wave (IAW) gratings. The detailed CBET coupling depends on the polarization state of the crossing wavefields. CBET itself can also scramble the beam polarizations by inducing ellipticity through the IAW grating, and rotating the seed polarization toward that of the pump. We develop a ray-based model that describes the polarized CBET coupling and that is compatible with the framework of 3D inline radiative hydrodynamics simulations. The model is implemented in the ASTER/IFRIIT code and verified against an academic test case and an offline polarized CBET post-processor. It is then applied to the detailed configuration of the distributed polarization rotator system on OMEGA, where results highlight how polarized CBET induces significant low modes in the collisional absorption source term. Finally, the modeling is applied to a simple indirect-drive configuration, comparing CBET calculations with 96 unpolarized or polarized beams with 24 unpolarized quads. It is shown that these cases produce similar power amplification per cone of beams grouped with similar polar angles. However, the 96 beam geometry itself is found to reduce azimuthal variations in quad power after the interaction and favors beams with larger polar angles within the cones, an effect that is amplified by the polarized CBET. Application of the model to inline calculations of OMEGA implosions are presented in a companion paper.

3D inline modeling↗

Bimolecular Peroxy Radical (RO 2 ) Reactions and Their Relevance in Radical Initiated Oxidation of Hydrocarbons

The kinetics of peroxy radical (RO 2 ) reactions have been of long-standing interest in atmospheric and combustion chemistry. Nevertheless, the lack of kinetic studies at higher temperatures for their reactions with other radicals such as OH has precluded the inclusion of this class of reactions in detailed kinetics models developed for combustion applications. In this work, guided by the limited room temperature experimental studies on selected alkyl-peroxy radicals and literature theoretical kinetics on the prototypical CH 3 O 2 + OH system, we have performed parametric studies on the effect of uncertainties in the rate coefficients and branching ratios to potential product channels for RO 2 + OH reactions at higher temperatures. Literature kinetics models were used to simulate autoignition delays, laminar flame speeds, and speciation profiles in flow and stirred reactors for a variety of common combustion-relevant fuels. Inclusion of RO 2 + OH reactions was found to retard autoignition in fuel-lean (φ = 0.5) mixtures of ethane and dimethyl ether in air. The observed effects were noticeably more pronounced in ozone-enriched combustion of ethane and dimethyl ether. The simulations also examined the influence of ozone doping levels, pressures, and equivalence ratios for both ethane and dimethyl ether oxidation. Sensitivity and flux analyses revealed that the RO 2 + OH reaction is a significant sink of RO 2 radicals at the early stage of autoignition, affecting fuel oxidation through RO 2 ↔ QOOH, RO 2 ↔ alkene + HO 2 , or RO 2 + HO 2 ↔ ROOH + O 2 . Additionally, the kinetic stability of the trioxide formed from RO 2 + OH reactions was investigated using master equation analyses. Last, we discuss other bimolecular reactions that are missing in literature kinetics models but are relevant to hydrocarbon oxidation initiated by external radical sources (plasma-enhanced, ozone-enriched combustion, etc.). In conclusion, the present simulations provide a strong motivation for better characterizing the bimolecular kinetics of peroxy radicals.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Hadronization of heavy quarks

Heavy-flavor hadrons produced in ultrarelativistic heavy-ion collisions are a sensitive probe for studying hadronization mechanisms of the quark-gluon-plasma. In this paper, we survey how different transport models for the simulation of heavy-quark diffusion through a quark-gluon plasma in heavy-ion collisions implement hadronization and how this affects final state observables. Utilizing the same input charm-quark distribution in all models at the hadronization transition, we find that the transverse-momentum dependence of the nuclear modification factor of various charm hadron species has significant sensitivity to the hadronization scheme. In addition, the charm-hadron elliptic flow exhibits a nontrivial dependence on the elliptic flow of the hadronizing partonic medium.

Atomic, Nuclear and Particle Physics↗

Combining Reactive Quantum-Mechanical Molecular-Dynamics Simulations with Mutagenesis, Crystallography, and Enzyme Kinetics to Reveal Plausible Steps of Isocyanide Hydratase Catalysis

A complete understanding of enzyme mechanisms requires atomistic details of chemical reactions. Quantum-based molecular dynamics simulations (QMD) are a potential source of this information, but trade-offs between accuracy and computational cost have limited their use. We previously developed extended Lagrangian Born–Oppenheimer molecular dynamics (XL-BOMD) methods that leverage a negligible compromise in accuracy to substantially decrease the cost of QMD simulations. Here, we develop a reactive QMD approach using the latest XL-BOMD formulation, which enables efficient simulations of highly reactive systems, and use it to investigate mechanisms of intermediate formation in isocyanide hydratase (ICH) catalysis. In QMD simulations, molecular analogs of ICH active site residues reacted with para-nitrophenyl isocyanide, forming a thioimidate. Analysis of simulated atomic configurational and charge dynamics revealed a pathway where protonation of the isocyanide carbon occurs prior to thioimidate formation and suggested a possible role of Asp17 as a proton donor in the early phase of ICH catalysis. To test whether the pathway seen using the reactive QMD approach might be relevant to ICH catalysis, we performed X-ray crystallography and pre-steady-state enzyme kinetics studies of wild-type and D17N mutant ICH. Both the structure and kinetics are sensitive to the D17N mutation in a manner that is consistent with the order of the reaction steps seen in the simulations. Mobile protons play essential roles in many enzymes, yet they are difficult to observe experimentally, making the ordering of proton-dependent steps ambiguous in many enzyme mechanisms. The ability to directly simulate model reactions for the design of experiments that provide information about enzyme mechanisms involving mobile protons demonstrates the significance of our reactive QMD approach and motivates further biological applications.

36 MATERIALS SCIENCE↗

A molecular dynamics simulation study of thermal conductivity anisotropy in β-octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine (β-HMX)

Molecular dynamics (MD) simulations were used to predict the thermal conductivity of β-octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine (β-HMX) along directions normal to the (011), (110), and (010) crystal planes. These directions were selected based on the measured morphological importance of the corresponding crystal surfaces. A reverse non-equilibrium MD approach was used wherein a constant heat flux is imposed along a prescribed direction and the resulting steady-state temperature gradient determined. The coefficient of thermal conductivity λ is the quotient of heat flux and temperature gradient (i.e. Fourier's law). Finite-size effects and sensitivity to imposed heat flux were investigated. The results reveal a modest dependence of the conductivity on crystal orientation, significant finite-size effects, and low sensitivity to imposed flux so long as the Fourier's law analysis is limited to the spatial interval in the simulation cell for which the temperature gradient is constant. Infinite-length thermal conductivities were estimated by linear regression of λ -1 (L) versus reciprocal cell length L -1 for each direction. The predicted values are systematically larger, but within a factor of two, than most published experimental determinations, the latter of which were obtained for pressed-powder or composite samples rather than oriented single crystals.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Time-Resolved Optical Pump-Resonant X-ray Probe Spectroscopy of 4-Thiouracil: A Simulation Study

We theoretically monitor the photoinduced ππ* → nπ* internal conversion process in 4-thiouracil (4TU), triggered by an optical pump. The element-sensitive spectroscopic signatures are recorded by a resonant X-ray probe tuned to the sulfur, oxygen, or nitrogen K-edge. We employ high-level electronic structure methods optimized for core-excited electronic structure calculation combined with quantum nuclear wavepacket dynamics computed on two relevant nuclear modes, fully accounting for their quantum nature of nuclear motions. We critically discuss the capabilities and limitations of the resonant technique. For sulfur and nitrogen, we document a pre-edge spectral window free from ground-state background and rich with ππ* and nπ* absorption features. The lowest sulfur K-edge shows strong absorption for both ππ* and nπ*. In the lowest nitrogen K-edge window, we resolve a state-specific fingerprint of the ππ* and an approximate timing of the conical intersection via its depletion. A spectral signature of the nπ* transition, not accessible by UV–vis spectroscopy, is identified. The oxygen K-edge is not sensitive to molecular deformations and gives steady transient absorption features without spectral dynamics. Here, the ππ*/nπ* coherence information is masked by more intense contributions from populations. Altogether, element-specific time-resolved resonant X-ray spectroscopy provides a detailed picture of the electronic excited-state dynamics and therefore a sensitive window into the photophysics of thiobases.

97 MATHEMATICS AND COMPUTING↗

Redshift inference from the combination of galaxy colours and clustering in a hierarchical Bayesian model – Application to realistic N -body simulations

ABSTRACT Photometric galaxy surveys constitute a powerful cosmological probe but rely on the accurate characterization of their redshift distributions using only broad-band imaging, and can be very sensitive to incomplete or biased priors used for redshift calibration. A hierarchical Bayesian model has recently been developed to estimate those from the robust combination of prior information, photometry of single galaxies, and the information contained in the galaxy clustering against a well-characterized tracer population. In this work, we extend the method so that it can be applied to real data, developing some necessary new extensions to it, especially in the treatment of galaxy clustering information, and we test it on realistic simulations. After marginalizing over the mapping between the clustering estimator and the actual density distribution of the sample galaxies, and using prior information from a small patch of the survey, we find the incorporation of clustering information with photo-z’s tightens the redshift posteriors and overcomes biases in the prior that mimic those happening in spectroscopic samples. The method presented here uses all the information at hand to reduce prior biases and incompleteness. Even in cases where we artificially bias the spectroscopic sample to induce a shift in mean redshift of $\Delta \bar{z} \approx 0.05,$ the final biases in the posterior are $\Delta \bar{z} \lesssim 0.003.$ This robustness to flaws in the redshift prior or training samples would constitute a milestone for the control of redshift systematic uncertainties in future weak lensing analyses.

79 ASTRONOMY AND ASTROPHYSICS↗

Understanding the origin of structure sensitivity in hydrodechlorination of trichloroethylene on a palladium catalyst

Periodic density functional theory (DFT) calculations are employed to understand the origin of structure sensitivity in hydrodechlorination (HDC) of trichloroethylene (TCE) over different facets of a palladium catalyst. The HDC reaction is simulated on the terrace (Pd (111) and Pd (100)) and undercoordinated (Pd (211) and Pd (110)) sites of the Pd catalyst. The most stable binding configuration of TCE on the Pd surfaces is observed to be through the di-σ mode of binding, wherein each carbon atom of the TCE molecule is adsorbed atop of the Pd atom. On comparing TCE adsorption over different facets of Pd, a maximum binding energy of –178 kJ mol –1 is calculated over the Pd (110) surface. TCE, upon adsorption on Pd catalyst, undergoes dechlorination followed by hydrogenation of the hydrocarbon intermediates. The activation energies for C–Cl bond dissociation steps are significantly low when compared to the hydrogenation steps. The chlorine released from dechlorination tends to block the active sites, thereby poisoning the surface with high binding energies (B.E > –160 kJ mol –1 ) on all the surfaces. The trend in chlorine binding energies on Pd facets follows: Pd (110) > Pd (211) > Pd (100) > Pd (111). The removal of surface chlorine is facilitated by its reaction with surface hydrogen to form hydrogen chloride. The activation energy for hydrogen chloride formation is calculated to be 90 kJ mol –1 and 88 kJ mol –1 on Pd (111) and Pd (100) terrace sites, respectively and 109 kJ mol –1 and 118 kJ mol –1 on the step Pd (211) and corrugated Pd (110) facets, respectively. This suggests the ease of removal of Cl as HCl from the terrace sites as compared to the step and corrugated sites. The structure sensitivity in the TCE HDC reaction could possibly arise due to the differences in the energetics of Cl removal on different Pd facets. As a result, this mechanistic understanding could provide a rationale for designing suitable catalysts for the HDC of TCE.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗