SEARCH · Engineering Papers

Results for “ENERGY DEPENDENCE”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 685 records · Page 38

All sources NASA NTRS DOE OSTI

DOE OSTI2024

New Nonreactive Force Field for Accurate Molecular Dynamics Simulations of TATB at Extreme Conditions

Insensitive high explosives based on TATB (1,3,5-triamino-2,4,6-trinitrobenzene) are needed in applications when safety is of paramount importance, but the basic material properties that give rise to its insensitivity are not fully understood. Molecular dynamics modeling using empirical force fields (FFs) has been the main route to characterize many complicated dynamical properties of TATB single crystal, but these FFs have not been comprehensively tested at extreme conditions typical of detonation. We collect a benchmark data set of (quasi)static TATB physical properties as determined by experiments and electronic structure calculations and apply this data set to validate four existing TATB FFs along with a new TATB FF that we develop here and denote as the CEA-LLNL-Missouri (CLM) FF. Benchmark data include vibrational spectra, the TATB crystal temperature–pressure–volume equation of state and lattice parameters, properties of TATB crystal polymorphs and transitions to the gaseous and liquid states, dimer energy landscapes, the pressure-dependent elastic tensor, and the energy landscape for inelastic deformation via sliding of TATB crystal layers. As a general assessment, we find that the two existing nonreactive FFs are more accurate in describing TATB’s physical properties compared to the two variants of the ReaxFF reactive FF considered. The new CLM FF is found to consistently yield similar or better agreement with experiments and electronic structure theory than any of the existing FF models, and it presents a distinct improvement in accurately modeling TATB elasticity and equation of state. So this work is expected to help improve the accuracy of FF-based modeling of complicated dynamic responses that ultimately govern the safety and performance characteristics of this material.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Interaction energy and isosteric heat of adsorption between hydrogen and magnesium diboride

Hydrogen storage materials form a crucial research topic for future energy utilization employing hydrogen and among those of interest magnesium diboride (MgB 2 ) has shown its prevalence. In this study, a first-principles analytical adsorption model of one hydrogen molecule in the vicinity of various magnesium diboride crystal surfaces was developed in order to obtain surface thermodynamic properties as a function of molecular and lattice properties. Henry's law constant (K H ) and isosteric heat of adsorption (ΔH ads ) indicators of the affinity between a gaseous molecule and a solid surface are thus calculated. The results in this paper not only address questions pertaining to the first stage of hydrogen storage processes but also advance the understanding of physisorption thermodynamics of a neutral molecule (H 2 ) coming in contact with a layered metallic-like surface (MgB 2 ). Although the model is built from a framework of classical calculations, quantum effects are incorporated as the fractional charge of the ions on the free surfaces, which is essential for the calculation of analytic thermodynamic values that approximate calculations from other methods. To benchmark our theoretical models, periodic density functional calculations were performed to determine the interactions between H 2 and different MgB 2 surfaces from first-principles. By considering both the top and sublayers of MgB 2 in calculating interaction energy, we have analytically and computationally calculated the interaction energies of H 2 molecules and MgB 2 's terminated planes, and witnessed the strong dependence of interaction energies on surface charges. We have also observed a dipole flipping phenomenon which explains the discontinuity seen in the interaction energy graph of Mg(0001). Both analytical and computational results showed heat of adsorption at zero coverage varying at a very low range (<7 kJ mol -1 ).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Neutron-Neutron Correlations in the Photofission of 238U

Background: In the fission of actinides, the nearly back-to-back motion of the fission fragments has a strong effect on the kinematics of fission neutrons. This leads to a favoring of opening angles near 0? and 180? in the neutron-neutron (n-n) opening angle distributions of correlated neutron pairs from the same fission event. Purpose: To measure the n-n opening angle and energy correlations in the photofission of 238U. As of this writing, measurements of correlated n-n opening angle distributions have been reported only for the spontaneous and neutron-induced fission of actinides. This work is the first to report such a measurement using photofission and will provide useful experimental input for photofission models used in codes such as MCNP and FREYA. Method: Fission is induced using bremsstrahlung photons produced via a low duty factor, pulsed, linear electron accelerator. The bremsstrahlung photon beam impinges upon a 238U target that is surrounded by a large neutron scintillation detection system capable of measuring particle position and time of flight, from which n-n opening angle and energy are measured. Neutron-neutron angular correlations are determined by taking the ratio between a correlated neutron distribution and an uncorrelated neutron distribution formed by the pairing of neutrons produced during different beam pulses. This analysis technique greatly diminishes effects due to detector efficiencies, acceptance, and experimental drifts. Results: The angular correlation of neutrons from the photofission of 238U shows a high dependence on neutron energy as well as a dependence on the angle of the emitted neutrons with respect to the incoming photon beam. Angular correlations were also measured using neutrons from the spontaneous fission of 252Cf, showing good agreement with past measurements. Conclusions: The measured angular correlations reflect the underlying back-to-back nature of the fission fragments. An anomalous decline in n-n yield was observed for opening angles near 180? for 238U.

Burggraf, Jeff↗

DOE OSTI2022

Role of ripples in altering the electronic and chemical properties of graphene

Ripples of graphene are known to manipulate electronic and hydrogenation properties of graphitic materials. More detailed work is needed to elucidate the structure–property relationship of these systems. In this work, the density functional theory is used to compute the energy and electronic structure of the graphene models with respect to variable curvatures and hydrogen adsorption sites. The magnitude of finite bandgap opening depends on the orientation of ripples, and the hydrogen adsorption energy depends on the local curvature of graphene. An adsorbed hydrogen alters the local curvature, resulting in relatively weakened adsorption on the neighboring three sites, which gives a rationale to experimentally observed dynamic equilibrium stoichiometry (H:C = 1:4) of hydrogenated graphene. The surface diffusion transition state energy of adsorbed hydrogen is computed, which indicates that the Eley–Rideal surface recombination mechanism may be important to establish the dynamic equilibrium, instead of the commonly assumed Langmuir–Hinshelwood mechanism.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Analysis and visualization of energy densities. I. Insights from real-time time-dependent density functional theory simulations

In this article, we report a scheme to analyze and visualize the energy density fluctuations during the real-time time-dependent density functional theory (RT-TDDFT) simulations. Using Ag 4 –N 2 complexes as examples, here we show that the grid-based Kohn–Sham energy density can be computed at each time step using a procedure from Nakai and coworkers. Then the instantaneous energy of each molecular fragment (such as Ag 4 and N 2 ) can be obtained by partitioning the Kohn–Sham energy densities using Becke or fragment-based Hirshfeld (FBH) scheme. A strong orientation-dependence is observed for the energy flow between the Ag 4 cluster and a nearby N 2 molecule in the RT-TDDFT simulations. Future applications of such an energy density analysis in electron dynamics simulations are discussed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Switching of the electron-phonon interaction in 1 T –VSe 2 assisted by hot carriers

We apply an intense infrared laser pulse in order to perturb the electronic and vibrational states in the three-dimensional charge density wave material 1T-VSe 2 . Ultrafast snapshots of the light-induced hot carrier dynamics and nonequilibrium quasiparticle spectral function are collected using time- and angle-resolved photoemission spectroscopy. Furthermore, the hot carrier temperature and time-dependent electronic self-energy are extracted from the time-dependent spectral function, revealing that incoherent electron-phonon interactions heat the lattice above the charge density wave critical temperature on a timescale of (200 ± 40) fs. Density functional perturbation theory calculations establish that the presence of hot carriers alters the overall phonon dispersion and quenches efficient low-energy acoustic phonon scattering channels, which results in a new quasiequilibrium state that is experimentally observed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2024

Environmental life-cycle analysis of hydrogen technology pathways in the United States

Hydrogen is a zero-carbon energy carrier with potential to decarbonize industrial and transportation sectors, but its life-cycle greenhouse gas (GHG) emissions depend on its energy supply chain and carbon management measures (e.g., carbon capture and storage). Global support for clean hydrogen production and use has recently intensified. In the United States, Congress passed several laws that incentivize the production and use of renewable and low-carbon hydrogen, such as the Bipartisan Infrastructure Law (BIL) in 2021 and the Inflation Reduction Act (IRA) in 2022, which provides tax credits of up to $3/kg depending on the carbon intensity of the produced hydrogen. A comprehensive life-cycle accounting of GHG emissions associated with hydrogen production is needed to determine the carbon intensity of hydrogen throughout its value chain. In the United States, Argonne’s R&D GREET ® (Greenhouse Gases, Regulated emissions, and Energy use in Technologies) model has been widely used for hydrogen carbon intensity calculations. This paper describes the major hydrogen technology pathways considered in the United States and provides data sources and carbon intensity results for each of the hydrogen production and delivery pathways using consistent system boundaries and most recent technology performance and supply chain data.

Elgowainy, Amgad↗

DOE OSTI2024

Finite element procedure for thermomechanical and structural integrity analysis of beam intercepting devices subjected to free electron laser

High-energy particles, including photons (x-ray, γ-ray, bremsstrahlung), electrons, and protons, possess the capability to penetrate materials and deposit energy within them. The degree of absorption depends on both the energy and type of particles, as well as the properties of the materials with which they interact. This energy deposition can manifest either at the material's surface or throughout its volume, potentially resulting in various failure modes. The primary aim of this paper is to establish a structured analysis methodology for evaluating the structural integrity of beam-intercepting devices when subjected to high-energy particles. Here, the paper also reviews some of the underlying physics, pertinent to the scope of the thermomechanical analysis, potential failure modes, and introduces verification and validation methodologies. Engineers and researchers can utilize the guidelines presented in this paper to effectively plan the development of beam intercepting devices, thereby ensuring their reliability and performance in the presence of high-energy particle exposure.

Finite element analysis↗

DOE OSTI2021

Improving Data Center Energy Efficiency Through End-to-End Cooling Modeling and Optimization. Final Report

This report summarizes the development and demonstration of the three features developed for the DOE project “Improving Data Center Energy Efficiency through End-to-End Cooling Modeling and Optimization” (Award No. DE-EE0007688). The development of the three software features, Feature 1: Optimization of Cooling System Operation, Feature 2: Optimization of Airflow Management System, and Feature 3: Simultaneous Optimization of Airflow Management and Cooling Systems, are described first. Feature 1 is implemented in Phase 1 of the project and demonstrates 15% energy savings from cooling system optimization by applying the proposed retrofit solutions. Feature 2 is implemented in Phase 2 and the assessments show cooling energy-saving potential between 19.3% and 54.0%, depending on the energy-saving measures employed. In addition, investigating the energy saving potential for a hypothetical simple data center shows about 6.8% energy savings by optimizing the server load placement. Feature 3 is implemented in Phase 3 and the results show energy savings of 53.4% and 74.6% for the two partner data centers, respectively, compared to the baselines defined in Phase 1. Lastly, the project outcomes including software releases and publications are discussed.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2026

First inclusive triple-differential measurement of the muon-antineutrino charged-current cross section using the NOvA Near Detector

We present the first measurement of the triple-differential muon antineutrino charged-current inclusive cross section, using the NOvA Near Detector and $12.5 \times 10^{20}$ protons on target in the NuMI beam. This sample of muon antineutrino interactions is the largest ever published, with approximately 1 million selected muon antineutrino events. The triple-differential cross section is measured in the final-state kinetic energy, the scattering angle, and the available energy of the interaction. The measurement enables phase-space regions populated by differing neutrino reaction processes to be isolated and the transition regions between them to be defined. The results are compared with the predictions of the main event generators used in the neutrino community and we observe energy- and angle-dependent discrepancies across a broad range of energies and interaction types

Abubakar, S. [Erciyes U.]↗

DOE OSTI2025

Be₂C as a neutron moderator for molten salt reactors: A first-principles study of structural, electronic, and defect properties

Beryllium carbide (Be 2 C), valued for its high neutron moderation efficiency and low absorption cross section, is a promising high-temperature neutron moderator for molten salt reactors. Its practical adoption, however, demands significant technological maturation, requiring comprehensive theoretical and experimental studies of its response to different conditions, including high temperature and irradiation. Here, we report initial results on the fundamental properties and radiation-induced defects of Be 2 C, focusing on antisites, vacancies, interstitial atoms, and Frenkel pairs in the Be and C sublattices. Using density functional theory (DFT) and ab initio molecular dynamics (AIMD), we calculate the defects formation and binding energies, evaluating their dependence on the supercell size, charge states, and chemical environment. In general, carbon defects exhibit higher formation energies, greater sensitivity to cell size, and stronger impacts on the density of states compared to beryllium defects, with charged state the effects being more pronounced. Static DFT reveals multiple metastable interstitial configurations, while AIMD identifies ground states as C-C <100> dumbbells and octahedral Be interstitials. In conclusion, the diversity of metastable configurations and defect states complicates the diffusion mechanisms, requiring further molecular dynamics analysis to elucidate the mechanisms and rates of radiation-induced atomic transport, as well as the structural stability of Be2C.

Ab initio modeling↗

DOE OSTI2021

Hyperglycemia alters mitochondrial respiration efficiency and mitophagy in human podocytes

Highlights: • Impairment of mitochondrial dynamics underlies the pathogenesis of type 2 diabetes. • Mitochondrial respiration is altered in glomeruli isolated from diabetic STZ-rats. • High glucose induces mitochondrial fission and fragmentation in human podocytes. • Mitochondrial fusion, biogenesis and mitophagy are reduced in HG-treated podocytes. • HG can injure the kidney by alteration of podocytes and glomerular bioenergetics. Podocytes constitute the outer layer of the renal glomerular filtration barrier. Their energy requirements strongly depend on efficient oxidative respiration, which is tightly connected with mitochondrial dynamics. We hypothesized that hyperglycemia modulates energy metabolism in glomeruli and podocytes and contributes to the development of diabetic kidney disease. We found that oxygen consumption rates were severely reduced in glomeruli from diabetic rats and in human podocytes that were cultured in high glucose concentration (30 mM; HG). In these models, all of the mitochondrial respiratory parameters, including basal and maximal respiration, ATP production, and spare respiratory capacity, were significantly decreased. Podocytes that were treated with HG showed a fragmented mitochondrial network, together with a decrease in expression of the mitochondrial fusion markers MFN1, MFN2, and OPA1, and an increase in the activity of the fission marker DRP1. We showed that markers of mitochondrial biogenesis, such as PGC-1α and TFAM, decreased in HG-treated podocytes. Moreover, PINK1/parkin-dependent mitophagy was inhibited in these cells. These results provide evidence that hyperglycemia impairs mitochondrial dynamics and turnover, which may underlie the remarkable deterioration of mitochondrial respiration parameters in glomeruli and podocytes.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2024

Retrofitting Holcim Ste. Genevieve Cement Plant with CO2 Capture Plant Using Air Liquide Cryocap™ FG Technology

The global cement manufacturing industry is a major contributor to carbon dioxide emissions. The International Energy Agency's "Net Zero Emissions by 2050 Scenario" identifies CCS as a major strategy for meeting that goal. This project is among the first attempts to transfer capture technology developed at coal-fired power plants to the cement industry. The main objective of the project is to execute and complete a front-end engineering and design (FEED) studies for commercial-scale, carbon capture projects that separates 95% of the total CO2 emissions at the Holcim (US) Ste. Genevieve cement manufacturing facility using Air Liquide’s Pressure Swing Adsorption system (PSA) assisted Cryocap™ technology. The Holcim Ste. Genevieve cement plant in Missouri, US, boasts one of the largest single cement production lines in the world, with a capacity of approximately 12,000 t/day. The plant currently uses traditional fuels, namely coal and petcoke. The captured CO2 will be pipeline and geological storage grade. The industrial host site emits approximately 3.0 million tonne CO2/yr. Air Liquide’s Cryocap™ technology has been developed over the last 18+ years for CO2 capture applications. It has been shown to be applicable to a variety of industrial applications (e.g., steel, cement, SMR, Fluidized Catalytic Crackers (FCCs)). Cryocap™ FG consists of a Pressure Swing Adsorption (PSA) unit coupled with a Cryogenic System. The PSA pre-concentrates the CO2 from the flue gas, while the cryogenic unit enables the CO2 purity to be increased to the desired level. The project team is led by the Prairie Research Institute at the University of Illinois at Urbana-Champaign. The tasks include: complete FEED study for retrofitting the industrial facility with a carbon capture system to support developing a detailed cost estimate; business case analysis outlining the anticipated revenue and credits if projects was built and operated; technoeconomic analysis (TEA) outlining how capture system achieves DOE capture goals; and life cycle (LCA) analysis demonstrating zero net carbon emissions. The FEED study was successfully completed. This includes completing the process basis of design; preliminary engineering; outside battery limits (OSBL) detailed engineering including a Zero Liquid Discharge (ZLD) wastewater treatment system; inside battery limits (ISBL) detailed engineering [1]. An overall project capital cost estimate within a -20%/+30% accuracy was developed. The major contributors to the Total Plant Cost (TPC), by system, are the costs associated with the Outside Battery Limit (OSBL) section of the plant which includes a new river water intake structure and a Zero Liquid Discharge (ZLD) system. By cost category, the major contributors to the TPC are equipment and subcontractor costs, followed closely by engineering, construction management, home office and contractor fees. The TEA has been created to reflect the findings of the project. It analyzes the economic performance of the Cryocap™ technology by reviewing the estimated capital costs, operating cost, and revenue. The Cost of Capture (COC) associated with the Cryocap™ technology for 95% CO2 capture, when considering NETL 2018 economic assumptions (42/58 debt/equity ratio, 5.15% interest on debt and 1.42% return on equity in real dollars) and 2022 economic assumptions (42/58 debt/equity ratio, 8.82% interest on debt and 4.90% return on equity in real dollars) was found to be much lower than that for the DOE-NETL’s base-line cases. The highest contributors to the COC are annualized capital expenditures (CAPEX) and electricity consumption which can be offset by using lower cost renewable sources. The LCA was conducted using OpenLCA which is an open-source software that is recommended by NETL. The database utilized for this study was a modified version of TRACI 2.1 (developed by the US. Environmental Protection Agency’s National Risk Management Research Laboratory and modified by NETL). The Cryocap™ FG technology does not consume fuels in significant quantities and does not utilize specialized chemical solvents subject to decomposition, such as those utilized in amine-based carbon capture systems. The Cryocap™ FG technology mainly utilizes electricity as its energy input; hence, its calculated emissions are mainly associated with the generation of electricity offsite and are dependent on the energy matrix of the grid at the time of project implementation. The water consumption impact of the Cryocap™ FG is mostly for makeup of the water lost by evaporation in the cooling tower; however, the carbon capture plant will be equipped with a ZLD system to avoid effluent streams and minimize water consumption. The successful construction and operation of this plant based on this study results will provide a means to demonstrate an economically attractive and transformational capture technology that can be used to retrofit existing plants and be deployed at new plants.

01 COAL, LIGNITE, AND PEAT↗

DOE OSTI2020

Characterization of light-induced potentials in the strong-field dissociation of O 2 +

In this work, we investigate theoretically the imprints of light-induced potentials (LIPs) on the dissociation dynamics of O 2 + molecular ions, as observable in angle-resolved fragment kinetic-energy-release (KER) spectra. Following the vibrational and rotational dynamics of the initial pump-laser-excited cationic nuclear wave packet, while accounting for the dipole coupling between the O 2 + ( a 4 Π u ) and O 2 + ( f 4 Π g ) electronic states in 800-nm 40-fs probe-laser pulses with peak intensities between 10 13 and 10 14 W / cm 2 , we calculate angle-resolved KER spectra which reveal characteristic energy- and angle-dependent fringe structures. These fringes shift downward in energy as the molecular alignment angle θ relative to the probe-pulse polarization direction increases from 0 to π / 2 . The angle-dependent shifts in the KER fringes increase for larger probe-pulse peak intensities and follow the angle and light-wave-intensity dependence of the vibrational spectrum in the associated Floquet bond-hardening well, which is a manifestation of transient O 2 + nuclear-probability trapping in the LIP during dissociation. By examining the rovibrational dynamics of the dissociating molecular cation near the light-induced conical intersection (LICI) in the cationic LIP surface at θ = π / 2 , we identify related angle-dependent structures in the KER spectra, suggesting a means for assessing the significance of LICIs in molecular dissociation pathways.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2022

Single-nucleon energies changing with nucleon number

The broad range of accumulated experimental data on the binding energies for single-particle states in nuclei is examined as a function of the constituent number of neutrons and protons and an unexpectedly simple pattern emerges. The dependence of the energies of neutron states on the number of constituent protons, or of proton states on the number of neutrons, are very similar to each other and the sign reflects the well-known strong attraction. For the same kind of nucleons changing as in the state – energies for neutron states with neutron number changing or proton states with protons – the dependence is at least a factor of four weaker in magnitude and slightly repulsive, except when the changing nucleons are only within the same orbit as the state. The systematics of the accumulated data are presented with a minimum of use made of model assumptions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2025

Evaluating the Structural Response of Amphiphilic Monolayers to Environmental Stimuli

Amphiphilic monolayers composed of end groups with distinct polar and nonpolar functional groups offer rapid and reversible interfacial adaptation in response to environmental stimuli such as a change in interfacial medium polarity. We have synthesized and characterized a suite of monolayers with functional groups of competing polarity designed to reconfigure their interfacial chemical composition in response to solvent polarity. In these films, the end group is designed to be able to reorient and expose the functional groups that minimize the interfacial free energy between the film and the environment. Using a combination of spectroscopic, computational, and wettability studies, we have investigated the responsive interfacial behavior of different end groups upon exposure to environments with varying polarities. Contact angle measurements across a series of polar and dispersive probe liquids reveal trends that reflect the underlying molecular flexibility and composition. Vibrational sum frequency generation (SFG) spectroscopy and atomistic molecular dynamics (MD) simulations confirm solvent-driven reorientation of the end groups, with restructuring observed at the interface. To quantify these effects, we have developed a surface energy calculation model that incorporates solvent-induced surface rearrangements into the estimations. Our findings reveal a strong dependence of surface energy and switching behavior on the length and flexibility of the functionalities in the end group, which affects the exposure of certain interfacial compositions under different solvents. These results offer new insights into the design of adaptive monolayers and provide a framework for evaluating solvent-responsive surfaces.

functional groups↗

DOE OSTI2021

Hot-carrier dynamics in InAs/AlAsSb multiple-quantum wells

A type-II InAs/AlAs 0.16 Sb 0.84 multiple-quantum well sample is investigated for the photoexcited carrier dynamics as a function of excitation photon energy and lattice temperature. Time-resolved measurements are performed using a near-infrared pump pulse, with photon energies near to and above the band gap, probed with a terahertz probe pulse. The transient terahertz absorption is characterized by a multi-rise, multi-decay function that captures long-lived decay times and a metastable state for an excess-photon energy of >100 meV. For sufficient excess-photon energy, excitation of the metastable state is followed by a transition to the long-lived states. Excitation dependence of the long-lived states map onto a nearly-direct band gap (E g ) density of states with an Urbach tail below E g . As temperature increases, the long-lived decay times increase < E g , due to the increased phonon interaction of the unintentional defect states, and by phonon stabilization of the hot carriers > E g . Additionally, Auger (and/or trap-assisted Auger) scattering above the onset of the plateau may also contribute to longer hot-carrier lifetimes. Meanwhile, the initial decay component shows strong dependence on excitation energy and temperature, reflecting the complicated initial transfer of energy between valence-band and defect states, indicating methods to further prolong hot carriers for technological applications.

14 SOLAR ENERGY↗

DOE OSTI2020

Neutrino Transport with Monte Carlo Method. I. Toward Fully Consistent Implementation of Nucleon Recoils in Core-collapse Supernova Simulations

The small energy exchange via nucleon recoils in neutrino–nucleon scattering is now supposed to be one of the important factors for successful explosions of core-collapse supernovae (CCSNe), as they can change neutrino spectra through accumulation of a large number of scatterings. In deterministic methods employed for neutrino transport in CCSN simulations, we normally cannot afford to deploy a large enough number of energy bins needed to resolve this small energy exchange, and subgrid techniques are employed one way or another. In this paper, we study quantitatively with the Monte Carlo (MC) method how well such a treatment performs. We first investigate the effects of nucleon recoils on the neutrino spectra and confirm that the average energy is reduced by ~15% for heavy-lepton neutrinos and much smaller amounts for other types of neutrinos in a typical postbounce situation. It is also observed that the nucleon scattering dominates the electron scattering in the thermalization of neutrino spectra in all flavors. We then study possible artifacts that the coarse energy grid may produce in the deterministic methods. In order to mimic the latter calculation, we redistribute MC particles in each energy bin after a certain interval in a couple of ways and study how the results are affected, depending on the energy resolution. We also discuss the possible implications of our results for the deterministic methods.

79 ASTRONOMY AND ASTROPHYSICS↗