Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Radiative transfer simulations”

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,027 records · Page 57

Electronic effects in radiation-induced collision cascades in nickel

The accurate treatment of electronic effects in multi-million-atom simulations of radiation-induced collision cascades is crucial for reliable predictions of primary radiation damage. In this work, we explore the fidelity of a recently developed two-temperature molecular dynamics model implementing an electron density-dependent coupling of electronic and atomic subsystems for cascade simulations in nickel. We show that the parameter-free model realistically captures the instantaneous energy losses during all stages of the highly nonequilibrium cascade process. Our simulations predict two distinct coupling regimes, corresponding to the rapid energy losses through electronic stopping in the early stages of the cascade and to the slow equilibration through the electron-phonon coupling mechanism in the later stages, without the use of separate models or coupling terms. The intermediate stage of the cascade dynamics displays a complex energy transfer between the subsystems, which cannot be interpreted by comparison to either electronic stopping or electron-phonon coupling theories. We therefore compare the predicted atomic mixing, which is sensitive to the energy losses during the intermediate cascade stage, with experimental ion beam mixing measurements. We find good agreement with the experiments, validating the coupling model for the intermediate stage of the cascade. Predictions of final defect numbers and cluster sizes are found in line with the results from conventional electronic stopping-based methods, while significantly reducing the theoretical uncertainty in the outcomes of conventional models stemming from arbitrary choices of thresholds for different coupling terms. Our results represent a notable improvement in cascade damage predictions in nickel, providing validation of the electron density-dependent coupling model for radiation damage simulations in general. The results lead us to propose an interpretation of the electronic energy losses in the intermediate regime of velocities, where we find an effectively nonlinear dissipation.

Crystal defects↗

Global 3D Radiation Magnetohydrodynamic Simulations of Accretion onto a Stellar-mass Black Hole at Sub- and Near-critical Accretion Rates

We present global 3D radiation magnetohydrodynamic simulations of accretion onto a 6.62 solar-mass black hole, with quasi-steady-state accretion rates reaching 0.016–0.9 times the critical accretion rate, which is defined as the accretion rate for powering the Eddington luminosity, assuming a 10% radiative efficiency, in three different runs. The simulations show no sign of thermal instability over hundreds of thermal timescales at 10 r g . The energy dissipation occurs close to the mid-plane in the near-critical runs and near the disk surface in the low–accretion rate run. The total radiative luminosity inside ~20 r g is about 1%–30% of the Eddington limit, with radiative efficiencies of about 6% and 3%, respectively, in the sub- and near-critical accretion regimes. In both cases, self-consistent turbulence generated by the magnetorotational instability leads to angular momentum transfer, and the disk is supported by magnetic pressure. Outflows from the central low-density funnel, with a terminal velocity of ~0.1c, are seen only in the near-critical runs. We conclude that these magnetic pressure–dominated disks are thermally stable and thicker than the α disk, and that the effective temperature profiles are much flatter than those in the α disks. The magnetic pressures of these disks are comparable within an order of magnitude to the previous analytical magnetic pressure–dominated disk model.

79 ASTRONOMY AND ASTROPHYSICS↗

Investigation of Best-Practices and Computationally Inexpensive Radiative Exchange Models for Discrete Element Method Modeling of Aluminosilicate Particles in Concentrating Solar Power Environments

Chemically inert, aluminosilicate based particles have been investigated as both a thermal transport and sensible energy storage medium for concentrating solar power facilities. These particles will experience a wide range of operating temperatures (300-1000 K) and handling conditions (dense to dilute falling particle curtains, dense granular flows, or dense structures), requiring specially-designed and optimized infrastructures. The relative influence of collisional and frictional interactions between particles varies based on temperature-dependent particulate properties and greatly impacts the bulk, granular flow behavior. These underlying physics are captured using discrete element method modeling tools. However, this modeling method is computationally expensive as each particle position and interaction is tracked during the simulation. These modeling methods are further complicated by introducing temperature-dependent particle properties, high-temperature radiative exchange, and directional irradiation sources experienced by granular flows in concentrating solar power environments. In this study, coupled experimental and numerical slump testing of aluminosilicate particles was performed and computationally efficient radiative exchange models were evaluated to establish best-practices for discrete element method models for concentrating solar power environments. The three particle types investigated included Carbobead HSP 30 /60, Carbobead CP 30/60, and Granusil 4030. Existing modeling limitations and computationally-efficient multi-modal heat transfer models were evaluated using Aspherix®, a commercial discrete element method software. High-temperature (< 1073 K) slump testing of aluminosilicate particles was performed to investigate the deviation between experimentally-observed and numerically-predicted angles of repose introduced by computation-time reduction practices including the relaxation of the particle elastic modulus and coarse-graining. Coarse-graining is used to use a single modeled particle that is representative of a collection of smaller particles, decreasing the computational cost at the expense of geometric accuracy. Additionally, relaxation of the elastic modulus is used to reduce computational time at the expense of an increased, modeled particle overlap. Prior studies have determined that aluminosilicate particles retain a high elastic modulus at high temperatures (< 1073 K), requiring small simulation timesteps to ensure resolved contact forces resemble appropriate solid mechanics. A parametric study was performed to evaluate the influence of computation time improvements on the deviation between experimental and modeled angle of repose across high temperatures < 1073 K. Additionally, numerical case studies were performed on candidate particle systems at varying porosities and temperatures. These studies were performed to investigate the influence of computationally-efficient radiative-exchange modeling methods coupled to Aspherix® on modeled accuracy and computation time. The recently-developed distance-based approximation was evaluated in estimating radiative exchange between particles and participating surfaces located in close proximity. The distance based approximation was developed to use tabulated estimates of the radiative distribution factor between individual particles and surfaces in close proximity (< 40 particle radii). These methods were expanded to the aluminosilicate particles of interest, including the influence of particle size distributions. To capture radiative exchange between particles and surfaces not in close proximity (> 40 particle radii) and to capture the absorption of directional irradiation from concentrating solar resources, a volumetrically-averaged radiative distribution factor was calculated between the modeled granular flow and surfaces using Monte Carlo ray-tracing for participating media. Volume-averaged absorption and scattering coefficients were predicted using a volumetric discretization of the modeled domain with monodisperse approximations based on geometric optics and experimentally-determined scattering phase functions for aluminosilicate particles.

14 SOLAR ENERGY↗

Continued Verification of MOOSE Structural Mechanics Tools for Modeling Core Bowing Phenomena in Fast Reactors

Under the U.S. Department of Energy Office of Nuclear Energy’s Advanced Modeling and Simulation (NEAMS) Program, an integrated multiphysics approach is being developed to model the core bowing phenomena important to liquid metal-cooled fast reactors. Core bowing is an important passive safety mechanism whereby increased power (which leads to temperature and flux gradients) influences the core to bow into less reactive configurations when the restraint system is properly designed. The phenomenon includes a complex interplay of radiation transport, duct temperature calculations involving fluid flow and heat transfer, and thermo-mechanical responses to the induced temperature and flux gradients. Structural material properties are also important to determining inelastic response to longer term flux gradients which cause irradiation creep and swelling. While core bowing provides a strong negative reactivity feedback when the restraint system is designed properly, it also results in additional forces between assemblies which increase the loads required to extricate them during refueling or control rod movement. Therefore, the restraint system must be designed with these tradeoffs in mind. The first stage of the work, which commenced in FY21 and continues through FY22, assesses thermo-mechanical modeling tools for producing core bowing predictions consistent with conventional tools. The Multiphysics Object Oriented Simulation Environment (MOOSE) Tensor Mechanics and Contact Modules are employed. This status report describes work on additional thermo-mechanical benchmark verification problems with increased complexity from the examples demonstrated in FY21. Several benchmark verification examples were selected from the IAEA verification and validation report. These examples involve clusters of ducts representative of a sector of a hexagonal reactor core which bow into each other and cause contact and load pad elevations, as well as single ducts subjected to irradiation fields undergoing swelling and subsequent bowing. The MOOSE-based results were compared to both IAEA benchmark participants’ results, analytic equations as available, and NUBOW-3D, a beam model code developed by Argonne National Laboratory. In every case, the MOOSE results agreed with other simulations results, providing additional verification basis of the tools for this particular physics application.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

FY23 Status Report on MOOSE-Based Approaches to Modeling Core Bowing in Fast Reactors

Under the U.S. Department of Energy Office of Nuclear Energy’s Advanced Modeling and Simulation (NEAMS) Program, an integrated multiphysics approach is being developed to model the core bowing phenomena important to liquid metal-cooled fast reactors. Core bowing is an important passive safety mechanism whereby increased power (which leads to temperature and flux gradients) influences the core to bow into less reactive configurations when the restraint system is properly designed. The phenomenon includes a complex interplay of radiation transport, duct temperature calculations involving fluid flow and heat transfer, and thermo-mechanical responses to the induced temperature and flux gradients. Structural material properties are also important to determining inelastic response to longer term flux gradients which cause irradiation creep and swelling. While core bowing provides a strong negative reactivity feedback when the restraint system is designed properly, it also results in additional forces between assemblies which increase the loads required to extricate them during refueling or control rod movement. Therefore, the restraint system must be designed with these tradeoffs in mind. The first stage of the work, which commenced in FY21 and continues through FY23, assesses thermo-mechanical modeling tools for producing core bowing predictions consistent with conventional tools. The Multiphysics Object Oriented Simulation Environment (MOOSE) Tensor Mechanics and Contact Modules are employed. This status report describes work on additional thermo-mechanical benchmark verification problems with increased complexity from the examples demonstrated in FY21 and FY22. Several benchmark verification examples were selected from the IAEA verification and validation report with increased number of ducts and more complex contact interaction behavior. These examples involve a full symmetric sector with restraint rings at multiple load pad locations to simulate a limited free-bow restraint system concept, as well as irradiation induced swelling and creep effects in a sector. In addition, improvements to the contact module sideset assignment were assessed and compared with previous MOOSE results to verify the contact behavior. The MOOSE-based results were compared to IAEA benchmark participants’ results. In every case, the MOOSE results agreed with other simulations results for estimating bowing behavior, providing additional verification basis of the tools for this particular physics application. Estimation of contact forces were mostly in agreement, with a few outlier results. A plan was suggested for dealing with the discrepancies with estimating contact force values. In addition, a 1-way coupling demonstration was performed using subchannel analysis code Pronghorn-SC and MOOSE on an ABR-1000-design sodium-cooled fast reactor assembly to evaluate coolant and duct wall temperatures.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Galactic cosmic ray simulation at the NASA space radiation laboratory – Progress, challenges and recommendations on mixed-field effects

For missions beyond low Earth orbit to the moon or Mars, space explorers will encounter a complex radiation field composed of various ion species with a broad range of energies. Such missions pose significant radiation protection challenges that need to be solved in order to minimize exposures and associated health risks. An innovative galactic cosmic ray simulator (GCRsim) was recently developed at the NASA Space Radiation Laboratory (NSRL) at Brookhaven National Laboratory (BNL). The GCRsim technology is intended to represent major components of the space radiation environment in a ground analog laboratory setting where it can be used to improve understanding of biological risks and serve as a testbed for countermeasure development and validation. The current GCRsim consists of 33 energetic ion beams that collectively simulate the primary and secondary GCR field encountered by humans in space over the broad range of particle types, energies, and linear energy transfer (LET) of interest to health effects. A virtual workshop was held in December 2020 to assess the status of the NASA baseline GCRsim. Workshop attendees examined various aspects of simulator design, with a particular emphasis on beam selection strategies. Experimental results, modeling approaches, areas of consensus, and questions of concern were also discussed in detail. This report includes a summary of the GCRsim workshop and a description of the current status of the GCRsim. This information is important for future advancements and applications in space radiobiology.

79 ASTRONOMY AND ASTROPHYSICS↗

Design Considerations and Operational Features of the Collimators for the Fermilab Main Injector and Recycler

The Fermilab Main Injector system delivers 700 kW of 120 GeV Proton beam for neutrino experiments. Since 2013 this has been achieved using slip stacking accumulation in the Recycler with up to 12 batches from the Fermilab Booster per Main Injector Ramp Cycle. To control activation from beam loss, collimation systems in the Booster to Recycler transfer line, in the Recycler and in the Main Injector are employed. Residual radiation measurements around the ring with detailed studies at the collimators are required to maintain adequate loss control. We will review design considerations, operational parameters and activation results for more than ten years of operation. Simulations with MARS15 are used to explore the activation rates and the isotopic composition of the resulting activation.

43 PARTICLE ACCELERATORS↗

Evidence of Recent Material Transport within a Binary Asteroid System

Images acquired by the Double Asteroid Redirection Test (DART) before its impact into Dimorphos—the moon of the Didymos binary asteroid system—unexpectedly reveal a set of discontinuous rays spanning the imaged surface. These albedo features are visible only after correcting for lighting differences, which normalize the effects of the pervasive boulders that cover the surface. Comparison to physical experiments and numerical simulations shows that these bright, fan-shaped albedo features are consistent with the aftermath of a low-velocity impact into a field of boulders and provide our first direct evidence for the transfer of material in a binary asteroid system due to the shedding material from radiative torques. Parts of these and possible additional albedo markings from the reentry of boulders after the DART impact may be visible to the Hera mission when it explores the post-DART Didymos system.

Astronomy and AstroPhysics↗

Evaluation of spatial, radiometric and spectral Thematic Mapper performance for coastal studies

A series of experiments were initiated to determine the feasibility of using thematic mapper spectral data to estimate wetlands biomass. The experiments were conducted using hand-held radiometers simulating thematic mapper wavebands 3, 4 and 5. Spectral radiance data were collected from the ground and from a low altitude aircraft in an attempt to gain some insight into the potential utility of actual thematic mapper data for biomass estimation in wetland plant communities. In addition, radiative transfer models describing volume reflectance of eight water column containing submerged aquatic vegetation were refined.

Klemas, V.↗

The response of cloud tops to radiative forcing

A longwave radiative transfer model was developed to study the effects on clouds of radiative loss of thermal energy from boundary layer cloud tops. Spectral properties of the drop size distributions are considered with Mie theory and LOWTRAN 6 gas transmission, which is combined with pathlength distributions from a Monte Carlo simulation, are used to obtain a 20/cm spectral resolution. Results are provided from a simulation of radiative cooling, at various wavelengths, of a 500 m thick cloud with a top at 2 km altitude. The very high values generated for the rates of cooling, of change of the liquid water content and of droplet growth rates indicate that cloud tops respond dynamically to strong thermal forcing.

Alves, A. R.↗

Radiative transfer to space through a precipitating cloud at multiple microwave frequencies. III - Influence of large ice particles

The time-dependent influence of large ice particles on the brightness temperatures, T(B)s, that would be observed from space above severe thunderstorms is investigated using a vertically and angularly detailed plane-parallel microwave radiative transfer model. The impact of large ice particles on passive microwave T(B)s over an evolving model rain cloud was examined at ten separate frequencies in the EHF/SHF spectrum. It is shown that the frequency-dependent vertically distributed generalized emission/scattering weighting functions, introduced in this simulation to vertically resolve the contributions by individual cloud and precipitation layers to the T(B)s, can be used to identify the specific layers responsible for regulating the magnitude of top-of-atmoshere T(B)s.

Smith, Eric A.↗

The Effect of Cumulus Cloud Field Anisotropy on Domain-Averaged Solar Fluxes and Atmospheric Heating Rates

Cumulus clouds can become tilted or elongated in the presence of wind shear. Nevertheless, most studies of the interaction of cumulus clouds and radiation have assumed these clouds to be isotropic. This paper describes an investigation of the effect of fair-weather cumulus cloud field anisotropy on domain-averaged solar fluxes and atmospheric heating rate profiles. A stochastic field generation algorithm was used to produce twenty three-dimensional liquid water content fields based on the statistical properties of cloud scenes from a large eddy simulation. Progressively greater degrees of x-z plane tilting and horizontal stretching were imposed on each of these scenes, so that an ensemble of scenes was produced for each level of distortion. The resulting scenes were used as input to a three-dimensional Monte Carlo radiative transfer model. Domain-average transmission, reflection, and absorption of broadband solar radiation were computed for each scene along with the average heating rate profile. Both tilt and horizontal stretching were found to significantly affect calculated fluxes, with the amount and sign of flux differences depending strongly on sun position relative to cloud distortion geometry. The mechanisms by which anisotropy interacts with solar fluxes were investigated by comparisons to independent pixel approximation and tilted independent pixel approximation computations for the same scenes. Cumulus anisotropy was found to most strongly impact solar radiative transfer by changing the effective cloud fraction, i.e., the cloud fraction when the field is projected on a surface perpendicular to the direction of the incident solar beam.

Hinkelman, Laura M.↗

Impact of Simulated Galactic Cosmic Radiation and Hindlimb Unloading on the Mouse Adrenal Gland Morphology and Histology

Stress related impact of Galactic Cosmic Radiation (GCR) combined with microgravity in space are concerns for astronauts going beyond the Van Allen belts on their way to the Moon and Mars. The adrenal gland is the central endocrine tissue involved in stress management, adrenal hormones regulate the immune system, cellular and tissue metabolism, and glucose and sodium balances. Adult (6 m old) female (n=18) and male (n=18) C57BL/6J mice were divided evenly into 6 treatment groups (control, hindlimb unloading (HU), irradiated (IR), HU+IR, group house control, and group house IR). Adrenals were dissected and fixed for 24 h in 4% paraformaldehyde, followed by transfer to 70% EtOH, and shipment to KUMC. Adrenals were then embedded in paraffin and serial sectioned and 4 central sections from each were stained with H&E, prior to imaging on a Nikon 80i microscope. Whole adrenal cross section area, medullary area, and cortex area were determined. General adrenal gland morphology and histology did not appear to be affected by hindlimb unloading or GCR exposure. Total cross-sectional area, medullary region, and cortex of the adrenal were not different across treatment groups. Male adrenal cross-sectional areas are half the size of the female adrenals, due to a difference in cortex area. In conclusion, simulated microgravity and GCR did not dramatically affect the adrenal gland histology or size. Support from Accelerate Cancer Education Internship Program (KUCC) and KINBRE Summer Scholars Program (KUMC) and NASA Human Research Program (HRP) Human Factors Behavioral Performance Element Grant 18 18FLAG 2 0028

radiation↗

A numerical study of the spectral radiative properties of packed bed with mixed bauxite and silica spheres

Bauxite and silica particles have gained increasing attention for applications in the field of concentrated solar power. In this work, a Monte Carlo ray-tracing simulation is performed to predict the radiative properties (absorptance, reflectance, and transmittance) of packed beds with mixed bauxite and silica spherical particles at wavelengths of 0.5, 2.6, 9.0, and 9.35 µm. These wavelengths are representative for the visible, near-infrared, and the mid-infrared regions that are important for solar and thermal radiation, respectively. A repeating unit-column approach is used to mathematically represent the particle bed. The effects of particle mixing ratio, volume fraction, and wavelength on the predicted radiative properties are examined. The obtained radiative properties are inputted to an inverse method to retrieve the effective absorption and scattering coefficients as well as the scattering albedo, which may be later used in a continuous-scale radiative heat transfer analysis. Furthermore, the independent scattering model is used to obtain the absorption coefficient and scattering albedo based on the absorption and scattering cross sections predicted by a Monte Carlo algorithm for a single particle. It is shown that the independent scattering model underpredicts the scattering coefficient for opaque particles but overpredicts the scattering coefficient for semitransparent particles for sufficiently high particle volume fractions. Here, the radiative properties calculated from the independent scattering model are compared to the full Monte Carlo simulation of the particle bed to examine the influence of particle mixing on dependent scattering.

42 ENGINEERING↗

Advancements in NEAMS Tool Capabilities for Multiphysics Simulation of Fast Reactor Core Bowing and Identification of Validation Test Data

Under the U.S. Department of Energy Office of Nuclear Energy Advanced Modeling and Simulation (NEAMS) Program, an integrated multiphysics approach is being developed to model the core bowing phenomena important to liquid metal-cooled fast reactors. Core bowing is an important passive safety mechanism in liquid metal-cooled fast reactors and involves Multiphysics effects including radiation transport, fluid flow, heat transfer, and mechanical response to temperature and flux gradients. This report summarizes recent progress on developing a multiphysics, MOOSE-based workflow to predict core bowing and associated reactivity feedback. Last year, thermal fluids and mechanics were coupled on a multi-assembly benchmark problem based on ABR-1000 design. This year, the reactor physics code Griffin was assessed for readiness of core bowing calculations. Preliminary integration of Griffin’s ring-heterogeneous model with thermal fluids and thermal mechanics solvers was performed. Specifically, thermal-mechanics and reactor physics were coupled for single- and multi-assembly problems, and reactor physics and subchannel methods were coupled for a single assembly model. Finally, the workflow of all three physics was preliminarily demonstrated on a single assembly model. Caveats and future development needed have been identified. To supplement the multiphysics demonstration, verification and assessment efforts of thermos-mechanical capabilities for modeling thermo-mechanical core bowing behavior were continued by analyzing IAEA Verification Problem 5 which includes radiation swelling and creep. Additionally, a small core reactor physics benchmark defined by Japan Atomic Energy Agency (JAEA) was performed to assess neutronics models for estimating reactivity feedback. Finally, Fast Flux Test Facility (FFTF) validation test data for core bowing phenomena has been identified and summarized, with a recommended path forward for validation once this capability is mature.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Technical feasibility of a modular dish solar electric system

The paper presents the technical feasibility of generating electricity by using a gas turbine in conjunction with a paraboloid of revolution (dish) solar concentrator. A conceptual design of a dish concentrator is obtained by parametrically examining the significant optical parameters. The optical performance analysis was conducted using a Monte Carlo ray-trace simulation program. The performance of four candidate thermodynamic cycles utilizing two working fluids was analyzed; the cycles were the regenerative and nonregenerative Brayton cycles of both open and closed type. Air and helium were the working fluids for the open and closed type respectively. Heat transfer from the solar radiation to the gas flowing in metal tubes at the receiver was also analyzed. A paraboloid of revolution dish with a cavity receiver using an open air regenerative Brayton cycle turbine emerged as a technically feasible concept in the power range from 30 to 100 KWth per module.

Gupta, B. P.↗

Microphysical Retrievals Over Stratiform Rain Using Measurements from an Airborne Dual-Wavelength Radar-Radiometer

The need to understand the complementarity of the radar and radiometer is important not only to the Tropical Rain Measuring Mission (TRMM) program but to a growing number of multi-instrumented airborne experiment that combine single or dual-frequency radars with multichannel radiometers. The method of analysis used in this study begins with the derivation of dual-wavelength radar equations for the estimation of a two-parameter drop size distribution (DSD). Defining a "storm model" as the set of parameters that characterize snow density, cloud water, water vapor, and features of the melting layer, then to each storm model there will usually correspond a set of range-profiled drop size distributions that are approximate solutions of the radar equations. To test these solutions, a radiative transfer model is used to compute the brightness temperatures for the radiometric frequencies of interest. A storm model or class of storm models is considered optimum if it provides the best reproduction of the radar and radiometer measurements. Tests of the method are made for stratiform rain using simulated storm models as well as measured airborne data. Preliminary results show that the best correspondence between the measured and estimated radar profiles usually can be obtained by using a moderate snow density (0.1-0.2 g/cu cm), the Maxwell-Garnett mixing formula for partially melted hydrometeors (water matrix with snow inclusions), and low to moderate values of the integrated cloud liquid water (less than 1 kg/sq m). The storm-model parameters that yield the best reproductions of the measured radar reflectivity factors also provide brightness temperatures at 10 GHz that agree well with the measurements. On the other hand, the correspondence between the measured and modeled values usually worsens in going to the higher frequency channels at 19 and 34 GHz. In searching for possible reasons for the discrepancies, It is found that changes in the DSD parameter Mu, the radar constants, or the path-integrated attenuation can affect the high frequency channels significantly. In particular, parameters that cause only modest increases in the median mass diameter of the snow, and which have a minor effect on the radar returns or the low frequency brightness temperature, can produce a strong cooling of the 34 GHz brightness temperature.

Meneghini, Robert↗

Electronic Excitation Response of DNA to High-Energy Proton Radiation in Water

The lack of molecular-level understanding for the electronic excitation response of DNA to charged particle radiation, such as high-energy protons, remains a fundamental scientific bottleneck in advancing proton and other ion beam cancer therapies. In particular, the dependence of different types of DNA damage on high-energy protons represents a significant knowledge void. Here, in this work, we employ first-principles real-time time-dependent density functional theory simulation, using a massively parallel supercomputer, to unravel the quantum-mechanical details of the energy transfer from high-energy protons to DNA in water. The calculations reveal that protons deposit significantly more energy onto the DNA sugar-phosphate side chains than onto the nucleobases, and greater energy transfer is expected onto the DNA side chains than onto water. As a result of this electronic stopping process, highly energetic holes are generated on the DNA side chains as a source of oxidative damage.

74 ATOMIC AND MOLECULAR PHYSICS↗