Search NASA⌕ Search

SEARCH · Search NASA

Results for “Fast Radiative Transfer Model”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
NASA NTRS2008

Mars Ozone Absorption Line Shapes from Infrared Heterodyne Spectra Applied to GCM-Predicted Ozone Profiles and to MEX/SPICAM Column Retrievals

We present the application of infrared heterodyne line shapes of ozone on Mars to those produced by radiative transfer modeling of ozone profiles predicted by general circulation models (GCM), and to contemporaneous column abundances measured by Mars Express SPICAM. Ozone is an important tracer of photochemistry Mars' atmosphere, serving as an observable with which to test predictions of photochemistry-coupled GCMs. Infrared heterodyne spectroscopy at 9.5 microns with spectral resolving power >1,000,000 is the only technique that can directly measure fully-resolved line shapes of Martian ozone features from the surface of the Earth. Measurements were made with Goddard Space Flight Center's Heterodyne instrument for Planetary Wind And Composition (HIPWAC) at the NASA Infrared Telescope Facility (IRTF) on Mauna Kea, Hawaii on February 21-24 2008 UT at Ls=35deg on or near the MEX orbital path. The HIPWAC observations were used to test GCM predictions. For example, a GCM-generated ozone profile for 60degN 112degW was scaled so that a radiative transfer calculation of its absorption line shape matched an observed HIPWAC absorption feature at the same areographic position, local time, and season. The RMS deviation of the model from the data was slightly smaller for the GCM-generated profile than for a line shape produced by a constant-with-height profile, even though the total column abundances were the same, showing potential for testing and constraining GCM ozone-profiles. The resulting ozone column abundance from matching the model to the HIPWAC line shape was 60% higher than that observed by SPICAM at the same areographic position one day earlier and 2.5 hours earlier in local time. This could be due to day-to-day, diurnal, or north polar region variability, or to measurement sensitivity to the ozone column and its distribution, and these possibilities will be explored. This work was supported by NASA's Planetary Astronomy Program.

Fast, Kelly E.↗
NASA NTRS2016

Fast Ionized X-Ray Absorbers in AGNs

We investigate the physics of the X-ray ionized absorbers often identified as warm absorbers (WAs) and ultra-fast outflows (UFOs) in Seyfert AGNs from spectroscopic studies in the context of magnetically-driven accretion-disk wind scenario. Launched and accelerated by the action of a global magnetic field anchored to an underlying accretion disk around a black hole, outflowing plasma is irradiated and ionized by an AGN radiation field characterized by its spectral energy density (SED). By numerically solving the Grad-Shafranov equation in the magnetohydrodynamic (MHD) framework, the physical property of the magnetized disk-wind is determined by a wind parameter set, which is then incorporated into radiative transfer calculations with xstar photoionization code under heating-cooling equilibrium state to compute the absorber's properties such as column density N(sub H), line-of-sight (LoS) velocity v, ionization parameter xi, among others. Assuming that the wind density scales as n varies as r(exp. -1), we calculate theoretical absorption measure distribution (AMD) for various ions seen in AGNs as well as line spectra especially for the Fe K alpha absorption feature by focusing on a bright quasar PG 1211+143 as a case study and show the model's plausibility. In this note we demonstrate that the proposed MHD-driven disk-wind scenario is not only consistent with the observed X-ray data, but also help better constrain the underlying nature of the AGN environment in a close proximity to a central engine.

accretion↗
DOE OSTI2024

SAM: A Modern System Code for Advanced Non-LWR Safety Analysis

The System Analysis Module (SAM), developed at Argonne National Laboratory and by collaborators at other organizations, is for advanced non–light water reactor safety analysis. SAM aims to provide fast-running, modest-fidelity, whole-plant transient analysis capabilities that are essential for fast-turnaround design scoping and engineering analyses of advanced reactor concepts. To facilitate code development, SAM utilizes the MOOSE object-oriented application framework, its underlying finite element library, and linear and nonlinear solvers to leverage modern advanced software environments and numerical methods. SAM aims to solve tightly coupled physical phenomena, including fission reaction, heat transfer, fluid dynamics, and thermal-mechanical responses in advanced reactor structures, systems, and components with high accuracy and efficiency. Finally, this paper gives an overview of the SAM code development, including goals and functional requirements, physical models, current capabilities, verification and validation, software quality assurance, and examples of simulations for advanced nuclear reactor applications.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS2008

Probing the Distribution of Ozone on Mars

We present the application of infrared heterodyne line shapes of ozone on Mars to those produced by radiative transfer modeling of ozone profiles predicted by photochemistry-coupled general circulation models (GCM), and to contemporaneous column abundances measured by Mars Express SPICAM. Ozone is an important tracer of photochemistry in Mars' atmosphere, serving as an observable with which to test predictions of photochemical models. Infrared heterodyne measurements of ozone absorption features on Mars have been obtained at various Martian seasons from 1988 until present at the NASA Infrared Telescope Facility (IRTF) on Mauna Kea, Hawai'i [I]. The NASAiGoddard Space Flight Center spectrometers used were the Infrared Heterodyne Spectrometer (IRHS) [2, 3] and, since 2003, the Heterodyne Instrument for Planetary Wind and Composition (HIPWAC) [4]. A description the infrared heterodyne technique applied to ground-base observations of Martian ozone can be found in [I]. The most recent measurements on February 21-24 2008 UT at Ls=35deg were made by HIPWAC on or near the Mars Express orbital path with the goal of acquiring spectra that can be directly compared to nadir observations by SPICAM.

Fast, K. E.↗
DOE OSTI2026

Dynamics and observational signatures of core-collapse supernovae with central engines: hydrodynamics simulations with Monte Carlo post-processing

A long-lived central engine embedded in expanding supernova ejecta can alter the dynamics and observational signatures of the event, producing an unusually luminous, energetic, and/or rapidly evolving transient. We use 2D hydrodynamics simulations to study the effect of a central energy source, varying the amount, rate, and isotropy of the energy deposition. We post-process the results with a time-dependent Monte Carlo radiation transport code to extract observational signatures. The engine excavates a bubble at the centre of the ejecta, which becomes Rayleigh–Taylor unstable. Sufficiently powerful engines are able to break through the edge of the bubble and accelerate, shred, and compositionally mix the entire ejecta. The breakout of the engine-driven wind occurs at distinct rupture points, and the outflowing high-velocity gas may eventually give rise to radio emission. The dynamical impact of the engine leads to faster rising optical light curves, with photon escape facilitated by the faster expansion of the ejecta and the opening of low-density channels. For models with strong engines, the spectra are initially hot and featureless, but later evolve to resemble those of broad-line Ic supernovae. Under certain conditions, line emission from ionized, low-velocity material near the centre of the ejecta may be able to escape and produce narrow emission similar to that seen in interacting supernovae. We discuss how variability in the engine energy reservoir and injection rate could give rise to a heterogeneous set of events spanning multiple observational classes, including the fast blue optical transients, broad-line Ic supernovae, and superluminous supernovae.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1989

Algorithm for Atmospheric Corrections of Aircraft and Satellite Imagery

A simple and fast atmospheric correction algorithm is described which is used to correct radiances of scattered sunlight measured by aircraft and/or satellite above a uniform surface. The atmospheric effect, the basic equations, a description of the computational procedure, and a sensitivity study are discussed. The program is designed to take the measured radiances, view and illumination directions, and the aerosol and gaseous absorption optical thickness to compute the radiance just above the surface, the irradiance on the surface, and surface reflectance. Alternatively, the program will compute the upward radiance at a specific altitude for a given surface reflectance, view and illumination directions, and aerosol and gaseous absorption optical thickness. The algorithm can be applied for any view and illumination directions and any wavelength in the range 0.48 micron to 2.2 micron. The relation between the measured radiance and surface reflectance, which is expressed as a function of atmospheric properties and measurement geometry, is computed using a radiative transfer routine. The results of the computations are presented in a table which forms the basis of the correction algorithm. The algorithm can be used for atmospheric corrections in the presence of a rural aerosol. The sensitivity of the derived surface reflectance to uncertainties in the model and input data is discussed.

Fraser, Robert S.↗
NASA NTRS1992

Algorithm for atmospheric corrections of aircraft and satellite imagery

A simple and fast atmospheric correction algorithm is described which is used to correct radiances of scattered sunlight measured by aircraft and/or satellite above a uniform surface. The atmospheric effect, the basic equations, a description of the computational procedure, and a sensitivity study are discussed. The program is designed to take the measured radiances, view and illumination directions, and the aerosol and gaseous absorption optical thickness to compute the radiance just above the surface, the irradiance on the surface, and surface reflectance. Alternatively, the program will compute the upward radiance at a specific altitude for a given surface reflectance, view and illumination directions, and aerosol and gaseous absorption optical thickness. The algorithm can be applied for any view and illumination directions and any wavelength in the range 0.48 micron to 2.2 microns. The relation between the measured radiance and surface reflectance, which is expressed as a function of atmospheric properties and measurement geometry, is computed using a radiative transfer routine. The results of the computations are presented in a table which forms the basis of the correction algorithm. The algorithm can be used for atmospheric corrections in the presence of a rural aerosol. The sensitivity of the derived surface reflectance to uncertainties in the model and input data is discussed.

Fraser, R. S.↗
NASA NTRS2007

A Simple Stochastic Model for Generating Broken Cloud Optical Depth and Top Height Fields

A simple and fast algorithm for generating two correlated stochastic twodimensional (2D) cloud fields is described. The algorithm is illustrated with two broken cumulus cloud fields: cloud optical depth and cloud top height retrieved from Moderate Resolution Imaging Spectrometer (MODIS). Only two 2D fields are required as an input. The algorithm output is statistical realizations of these two fields with approximately the same correlation and joint distribution functions as the original ones. The major assumption of the algorithm is statistical isotropy of the fields. In contrast to fractals and the Fourier filtering methods frequently used for stochastic cloud modeling, the proposed method is based on spectral models of homogeneous random fields. For keeping the same probability density function as the (first) original field, the method of inverse distribution function is used. When the spatial distribution of the first field has been generated, a realization of the correlated second field is simulated using a conditional distribution matrix. This paper is served as a theoretical justification to the publicly available software that has been recently released by the authors and can be freely downloaded from http://i3rc.gsfc.nasa.gov/Public codes clouds.htm. Though 2D rather than full 3D, stochastic realizations of two correlated cloud fields that mimic statistics of given fields have proved to be very useful to study 3D radiative transfer features of broken cumulus clouds for better understanding of shortwave radiation and interpretation of the remote sensing retrievals.

Prigarin, Sergei M.↗
DOE OSTI2020

Remote Radiation Sensing Using Aerial and Ground Platforms

Remote sensing of ionizing radiation has a significant role in waste management, nuclear material management and nonproliferation, and radiation safety. Robotic platforms can surpass the number of tasks that are achieved by humans. With this technique, the operator's radiation exposure can be decreased. Remote sensing allows for the evaluation and monitoring of radiological contamination. Gamma-ray and neutron sensors were integrated onto the robotic platforms. This approach allows for the radiation sensor data to be dynamically tracked and mapped thus enabling further analysis of the radiation flux in temporal and spatial domains. The goal is to complete scheduled tasks while the robot is being irradiated. To achieve this, electronic components must be shielded and radiation hardened. CZT Detector: Cadmium Zinc Telluride (CZT) detector technology has been a promising solution for gamma-ray and x-ray measurements. Detector data is transferred to the Odroid minicomputer that controls and powers the module via the USB. Robot Operating System (ROS) was utilized for data acquisition and data fusion. The Mariscotti method was employed for the spectrum analysis. A function was programmed in ROS for the automatic identification of photopeaks. CLYC Detector: A Cs{sub 2}LiYCl{sub 6}:Ce{sup 3+} (CLYC) detector was used for simultaneous medium-resolution gamma-ray measurements and neutron counting. A 2.54 cm diameter photomultiplier tube (PMT) was equipped with a high voltage supply and a miniature digitizer. Gamma-ray excitation: fast core-to-valence luminescence (CVL) with 1 ns decay constant, and prompt Ce{sup 3+} emission with 50 ns decay constant. Neutron excitation: slow cerium self-trapped excitation (Ce{sup 3+} STE), 1000 ns decay constant. Radiation Source Localization: Maximum Likelihood Estimation (MLE) and gradient-based methods were used to locate the position of a radiation source based on measured radiation intensities. Multi-Particle Transport Code FLUKA: Estimation of radiation damage of the electronic components is important in order to optimize the robot's operational time while it is irradiated. Displacement per atom (DPA) represents the radiation damage in materials exposed to the ionizing radiation. Various shielding layers of different thickness t were analyzed (< 5% statistical error). The model of the controller of the UAS was designed in FLUKA. Conclusion: CZT and CLYC detectors were integrated onto the robotic platforms. Radiation source localization and contour mapping using robotic platforms were studied. Functions for data analysis and fusion were developed in ROS. FLUKA code was utilized to analyze DPA values. Layers of low-density and high-density materials were used to shield the UAS electronics.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2024

Methodological guidelines on concrete degradation based on predictive models and the release of MOSAIC for industry use

In light-water reactors (LWRs), the concrete biological shield (CBS) designates the concrete structure directly facing the reactor pressure vessel (RPV). The primary function of the CBS is to protect equipment and personnel from the neutron and gamma radiation exiting the RPV. Most CBSs in operation in the United States also provide structural support to the RPV. Based on a long-term structural analysis, the structural integrity of the CBS to transfer the in-service passive load to the foundation, dynamic loading during an earthquake, and thermal loading during a loss-of-coolant accident (LOCA) must be maintained. Based on test reactor data for accelerated conditions, the mechanical properties of irradiated concrete are affected when exposed to neutron fluence higher than ≈ 10 19 n.cm -2 (E > 0.1 MeV). At 80 years of operation, the fluence estimates at the surface of the CBS range between 1×10 19 and 7×10 19 n.cm -2 (E > 0.1 MeV): the fluence value depends mostly on the design. Although fluence attenuates through the CBS by virtue of the shielding properties of concrete, about 5% to 10% of the CBS wall depth is affected by high fluence. The main degradation mechanism is the so-called radiation-induced volumetric expansion (RIVE) caused by the amorphization of aggregate-forming minerals. RIVE causes large deformation of up to 18% in pure quartz, and it also reduces the structural properties of concrete. The structural effects of long-term exposure to fast neutron irradiation in LWRs was studied in the report entitled “Assessment of the Effect of the Irradiation-Induced Degradation on the In-Service Structural Performance of the Concrete Biological Shields” (M3LW-23OR0403013). This report illustrates that (1) irradiation-induced damage expands further than the region subject to a cumulative fluence of 10 19 n.cm -2 (E > 0.1 MeV) to reach the steel reinforcement location, and (2) the irradiation-induced damage region shows no substantial residual bearing capacity and should not be accounted for in the integrity assessment of the CBS.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
NASA NTRS2017

Treatment of Overlapping Gaseous Absorption with the Correlated-k Method in Hot Jupiter and Brown Dwarf Atmosphere Models

The correlated-k method is frequently used to speed up radiation calculations in both one-dimensional and three-dimensional atmosphere models. An inherent difficulty with this method is how to treat overlapping absorption, i.e. absorption by more than one gas in a given spectral region. We have evaluated the applicability of three different methods in hot Jupiter and brown dwarf atmosphere models, all of which have been previously applied within models in the literature: (i) Random overlap, both with and without resorting and rebinning, (ii) equivalent extinction and (iii) pre-mixing of opacities, where (i) and (ii) combine k-coefficients for different gases to obtain k-coefficients for a mixture of gases, while (iii) calculates k-coefficients for a given mixture from the corresponding mixed line-by-line opacities. We find that the random overlap method is the most accurate and flexible of these treatments, and is fast enough to be used in one-dimensional models with resorting and rebinning. In three-dimensional models such as GCMs it is too slow, however, and equivalent extinction can provide a speed-up of at least a factor of three with only a minor loss of accuracy while at the same time retaining the flexibility gained by combining k-coefficients computed for each gas individually. Pre-mixed opacities are significantly less flexible, and we also find that particular care must be taken when using this method in order to to adequately resolve steep variations in composition at important chemical equilibrium boundaries. We use the random overlap method with resorting and rebinning in our one-dimensional atmosphere model and equivalent extinction in our GCM, which allows us to e.g. consistently treat the feedback of non-equilibrium chemistry on the total opacity and therefore the calculated P-T profiles in our models.

atmospheres↗
NASA NTRS2018

Pilot Study on the Investigation of Tear Fluid Biomarkers as an Indicator of Ocular, Neurological, and Immunological Health in Astronauts

The purpose of this pilot study is to investigate the collection, preparation, and analysis of tear biomarkers as a means of assessing ocular, neurological, and immunological health. At present, no published data exists on the cytokine profiles of tears from astronauts exposed to long periods of microgravity and space irradiations. In addition, no published data exist on cytokine (biomarker) profiles of tears that have been collected from irradiated non-human biological systems (primates and other animal models). A goal for the proposed pilot study is to discover novel tear biomarkers which can help inform researchers, clinicians, epidemiologist and healthcare providers about the health status of a living biological system, as well as informing them when a disease state is triggered. This would be done via analysis of the onset of expression of pro-inflammatory cytokines, leading up to the full progression of a disease (i.e. cancer, loss of vision, radiation-induced oxidative stress, cardiovascular disorders, fibrosis in major organs, bone loss). Another goal of this pilot study is to investigate the state of disease against proposed medical countermeasures, in order to determine whether the countermeasures are efficacious in preventing or mitigating these injuries. An example of an up and coming tear biomarker technology, Ascendant Dx, a clinical stage diagnostic company, is developing a screening test to detect breast cancer using proteins from tears. The team utilized Liquid Chromatography -Mass Spectrometry with Mass analysis (LC MS/MS) as a discovery platform followed by validation with ELISA to come up with a panel of protein biomarkers that can differentiate breast cancer samples from control ("cancer free") samples with results far surpassing the results of imaging techniques in use today. Continued research into additional proteins is underway to increase the sensitivity and specificity of the test and development efforts are on the way to transfer the test onto a fast, accurate and inexpensive point of care platform. In conclusion, the expected results from this proposed pilot study are to: a) establish an SOP for retrieving/storing/transporting tear fluid samples from multicentre sites b) establish a normal range for relevant biomarkers in tears; and c) establish a database (biobank) of tears of space naïve versus veteran astronauts, to establish a personal baseline for long-term ocular health monitoring

Morton, Stephen↗
DOE OSTI2023

Development of HELIOS/BIPR/PARCS/MCNP6 Computation Route for WWER RPV Neutron Fluence Analysis and Validation Against Ex-Vessel Detector Measurement Data

Here, a computational route was developed for precise calculation of fast neutron fluence on a WWER-type reactor pressure vessel (RPV). The method is based on the transfer of neutronics data from HELIOS-2 lattice calculations and nodal diffusion neutronics data (power, density, and temperature) from BIPR7.1 and PARCS 3.36/PATHS core calculations into a three-dimensional (pinwise axially distributed) fixed neutron source for modeling of transport of fast neutrons from the reactor core to the outer surface of the RPV using MCNP6.2. Validation of the proposed computational method was carried out based on comparative analysis of MCNP6.2-predicted and neutron dosimetry–measured reaction rates [ 54 Fe(n,p) 54 Мn, 93 Nb(n,n') 93 mNb, and 58 Ni(n,р) 58 Со] on the outer surface of the Armenian Nuclear Power Plant (ANPP) Unit 2 RPV. Validation revealed that the MCNP6.2-predicted fast neutron fluence results are very sensitive to the ENDF-B neutron data. Particularly, MCNP6.2 with ENDF/B-VIII.0 significantly underpredicts (20% to 30%) fast neutron fluence while using ENDF/B-VII.1 data overpredicts it. Adding revised beta-released evaluations of 54 Fe, 56 Fe, 57 Fe, and 16 O from the International Nuclear Data Evaluation Network (INDEN) to ENDF/B-VIII.0 allows one to obtain reasonable agreement with measurement results for all types of measured reaction rates.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Strong Proton‐Phonon Coupling Drives Fast Ion Transport in Perovskites

Conduction of protons in solids is a cooperative process propelled by phonons, with molecular details obscured by the irregular movements in the thermal bath. It is shown that substitution with Y forms an imaginary phonon mode, instrumental for the function as proton conductor and effectively lowering the activation barrier for proton transport. To untangle the interplay in the exemplary proton conductor BaSn 0.9 Y 0.1 O 3 , its crystallographic structure is determined with high resolution neutron diffractometry and its phonon density of states with density functional theory calculations, experimentally validated by element specific nuclear resonant vibration spectroscopy. Based on phonon analysis, a quantitative transport model is present, which predicts the activation energy and performance by the ratio of ionic radii. Rather than individual vibrational modes, it is the oxygen sub-lattice which exerts its momentum on the protons. The extent of this momentum transfer is governed by the ratio of ionic radii. This model extends the transition state theory by the phonon-phonon interaction and complements the previously proposed idea that lattice dynamics is decisive for proton transport and specifies which properties of the material exactly define the vibration properties.

activation barriers↗
DOE OSTI2020

Sensitivity of aerosol optical properties to the aerosol size distribution over central Europe and the Mediterranean Basin using the WRF-Chem v.3.9.1.1 coupled model

The size distribution of atmospheric aerosols plays a key role for understanding and quantifying the uncertainties related to aerosol–radiation and aerosol–cloud interactions. These interactions ultimately depend on the size distribution through optical properties (such as aerosol optical depth, AOD) or cloud microphysical properties. Hence, the main objective of this contribution is to disentangle the impact of the representation of aerosol size distribution on aerosol optical properties over central Europe, particularly over the Mediterranean Basin, during a summertime aerosol episode. To fulfill this objective, a sensitivity test has been conducted using the coupled chemistry–meteorology model WRF-Chem (Weather Research Forecast model coupled with Chemistry). The test modified the parameters defining a lognormal size distribution (geometric diameter and standard deviation) by 10 %, 20 %, and 50 %. Results reveal that the reduction in the standard deviation of the accumulation mode leads to the largest impacts on AOD due to a transfer of particles from the accumulation mode to the coarse mode. A reduction in the geometric diameter of the accumulation mode also has an influence on AOD representation since particles in this mode are assumed to be smaller. In addition, an increase in the geometric diameter of the coarse mode produces a redistribution through the total size distribution by relocating particles from the finer modes to the coarse.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Effects of Quantum Dot Loading on the Radioluminescence Efficiency in Quantum-Dot-Embedded Composites

Nanoparticle-embedded plastic scintillators are an emerging technology for fast, large-area, high-resolution radiation detection and imaging. Here, this study investigates the properties of such composites, focusing on the effects of the quantum dot (QD) concentration on the radioluminescence (RL) intensity, spectra, and dynamics. Experiments using CdSe/CdS QDs in a polymer reveal a superlinear increase in RL with the QD concentration despite optical losses from inner filtering and interparticle interactions. When corrected for inner filtering, RL shows a quadratic concentration dependence, consistent with simple analytical models of improving the secondary electron capture. Practically, the benefits of high QD concentrations are muted by optical losses, but the findings apply to other systems with insulating hosts. In addition to manipulating emission for large effective Stokes shifts, future improvements may come from hosts with higher stopping power and better charge transport, which enable more effective funneling of excitations but without concomitant optical losses associated with high nanoparticle concentrations.

Auger recombination↗
DOE OSTI2020

Quantal diffusion approach for multinucleon transfer processes in the 58,64 Ni + 208 Pb reactions: Toward the production of unknown neutron-rich nuclei

In recent years, substantial efforts have been made for the study of multinucleon transfer reactions at energies around the Coulomb barrier both experimentally and theoretically, aiming at the production of unknown neutron-rich heavy nuclei. It is crucial to provide reliable theoretical predictions based on microscopic theories with sufficient predictive power. Purpose: This article aims to clarify the applicability of the quantal diffusion approach based on the stochastic mean-field (SMF) theory for multinucleon transfer processes. Isotope production cross sections are evaluated for the reactions of 64 Ni + 208 Pb at E c.m. = 268 MeV and 58 Ni + 208 Pb at E c.m. = 270 MeV and are compared with available experimental data. Methods: Three-dimensional time-dependent Hartree-Fock (TDHF) calculations are carried out for a range of initial orbital angular momenta with Skyrme SLy4d functional. Quantal diffusion equations, derived based on the SMF theory, for variances and covariance of neutron and proton numbers of reaction products are solved, with microscopic drift and diffusion coefficients obtained from time evolution of occupied single-particle orbitals in TDHF. Secondary de-excitation processes, both particle evaporation and fission, are simulated by a statistical compound-nucleus de-excitation model, gemini++. Dynamics of a fast isospin equilibration process followed by a slow drift toward the mass symmetry are commonly observed, as expected. Various reaction outcomes are evaluated, including average mass and charge numbers of reaction products, total kinetic energy loss (TKEL), scattering angle, contact time, and production cross sections for primary and secondary products. By comparing with the experimental data, we find that SMF and TDHF quantitatively reproduce experimental data for few-nucleon-transfer channels around the average values. In contrast, for many-nucleon-transfer channels, we find that the SMF approach provides much better description of the experimentally measured isotopic distributions. The results underline the importance of beyond-mean-field effects, especially one-body (mean-field) fluctuations and correlations, in describing multinucleon transfer processes. Moreover, through a combined analysis of SMF with a statistical model, gemini++, we find a significant contribution of transfer-induced fission, which is consistent with the experimental observation. In some cases, the SMF approach overestimates the isotopic width, requiring further improvements of the theoretical description. Possible ways to improve the description are discussed. The SMF approach is designed to describe the quantum many-body problem according to an ensemble of mean-field trajectories, taking into account part of many-body correlations in the description. As it requires feasible computational costs comparable to the ordinary TDHF approach, together with further model improvements, it will be a promising tool in the search for optimal reaction conditions to produce yet-unknown neutron-rich heavy nuclei through the multinucleon transfer reaction.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS2011

Numerical Solution of the Radiative Transfer Equation: X-Ray Spectral Formation from Cylindrical Accretion onto a Magnetized Neutron Star

Predicting the emerging X-ray spectra in several astrophysical objects is of great importance, in particular when the observational data are compared with theoretical models. This requires developing numerical routines for the solution of the radiative transfer equation according to the expected physical conditions of the systems under study. Aims. We have developed an algorithm solving the radiative transfer equation in the Fokker-Planck approximation when both thermal and bulk Comptonization take place. The algorithm is essentially a relaxation method, where stable solutions are obtained when the system has reached its steady-state equilibrium. Methods. We obtained the solution of the radiative transfer equation in the two-dimensional domain defined by the photon energy E and optical depth of the system pi using finite-differences for the partial derivatives, and imposing specific boundary conditions for the solutions. We treated the case of cylindrical accretion onto a magnetized neutron star. Results. We considered a blackbody seed spectrum of photons with exponential distribution across the accretion column and for an accretion where the velocity reaches its maximum at the stellar surface and at the top of the accretion column, respectively. In both cases higher values of the electron temperature and of the optical depth pi produce flatter and harder spectra. Other parameters contributing to the spectral formation are the steepness of the vertical velocity profile, the albedo at the star surface, and the radius of the accretion column. The latter parameter modifies the emerging spectra in a specular way for the two assumed accretion profiles. Conclusions. The algorithm has been implemented in the XPEC package for X-ray fitting and is specifically dedicated to the physical framework of accretion at the polar cap of a neutron star with a high magnetic field (approx > 10(exp 12) G). This latter case is expected to be of typical accreting systems such as X-ray pulsars and supergiant fast X ray transients.

Fairnelli, R.↗