Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “RANDOM PHASE APPROXIMATION”

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 325 records · Page 18

Model Parameter Development for Complex Materials: Species-Specific Diffusion Barriers in 316 Stainless Steel from Systematic DFT Calculations

Vacancy-mediated diffusion barriers in 316 stainless steel have been systematically calculated using density functional theory to provide essential parameters for mesoscale microstructure evolution models. A statistical sampling approach employing 210 nudged elastic band calculations across multiple special quasi-random structures captures the effects of local chemical environments in this concentrated alloy. The computational methodology addresses challenges specific to chemically disordered systems, including proper magnetic treatment throughout multi-step calculations and validation against experimental structural properties. The calculated activation barriers reveal clear species-dependent diffusion behavior with the hierarchy Ni >> Fe ˜ Cr >> Mo. Nickel exhibits the highest barriers (0.74–1.31 eV, mean 1.045 eV), confirming its role as the slowest-diffusing major component. Iron and chromium show similar moderate barriers averaging 0.587 eV and 0.522 eV, respectively. Remarkably, molybdenum demonstrates exceptionally low barriers (0.12–0.28 eV, mean 0.194 eV), suggesting much higher mobility than previously recognized and potentially significant implications for precipitation kinetics and microstructure evolution. The barrier ranges remain consistent across different 316 SS compositions, supporting parameter transferability for modeling applications. The overall mean barrier of 0.64 eV provides a practical approximation for phase field simulations, while species-specific values enable detailed treatments of diffusion-controlled processes. This systematic approach establishes a validated framework for generating diffusion parameters in other concentrated alloys where experimental data are limited, while providing the first systematic set of species-specific barriers for predictive modeling of 316 stainless steel microstructure evolution.

36 MATERIALS SCIENCE↗

Hydrogen Adsorption on Ordered and Disordered Pt-Ni Alloys

The bulk properties and chemical reactivity of disordered Pt-Ni alloys in the A1 ( fcc ) structure are investigated herein using different methods: Virtual Crystal Approximation (VCA), Korringa–Kohn–Rostoker Coherent Potential Approximation (KKR-CPA), and large explicit supercells generated using Super-Cell Random Approximates (SCRAPs). While VCA predicts lattice constants that closely follow Vegard’s law, the large supercells and KKR-CPA predict lattice constants that are consistently larger than Vegard’s law. KKR-CPA results closely agree with those from the large supercells for the disordered alloys, producing similar projected density of states and magnetic moment across the composition range. For instance, while VCA predicts the disordered alloys to be non-magnetic at a Pt concentration ( x Pt ) ≥ 0.5, KKR-CPA and SCRAPs predict the disordered alloys to remain ferromagnetic to higher Pt concentrations. As x Pt decreases, the adsorption of H becomes more exothermic on bulk-terminated (111) surfaces but less exothermic on Pt monolayer-terminated (111) surfaces due largely to strain effects. (111) surfaces cut from the large supercells predict average H adsorption energies on the disordered alloys similar to those on the ordered phases of the same compositions, while VCA predicts H adsorption to be more exothermic.

36 MATERIALS SCIENCE↗

Analysis of Phase-Type Stochastic Petri Nets With Discrete and Continuous Timing

The Petri net formalism is useful in studying many discrete-state, discrete-event systems exhibiting concurrency, synchronization, and other complex behavior. As a bipartite graph, the net can conveniently capture salient aspects of the system. As a mathematical tool, the net can specify an analyzable state space. Indeed, one can reason about certain qualitative properties (from state occupancies) and how they arise (the sequence of events leading there). By introducing deterministic or random delays, the model is forced to sojourn in states some amount of time, giving rise to an underlying stochastic process, one that can be specified in a compact way and capable of providing quantitative, probabilistic measures. We formalize a new non-Markovian extension to the Petri net that captures both discrete and continuous timing in the same model. The approach affords efficient, stationary analysis in most cases and efficient transient analysis under certain restrictions. Moreover, this new formalism has the added benefit in modeling fidelity stemming from the simultaneous capture of discrete- and continuous-time events (as opposed to capturing only one and approximating the other). We show how the underlying stochastic process, which is non-Markovian, can be resolved into simpler Markovian problems that enjoy efficient solutions. Solution algorithms are provided that can be easily programmed.

Jones, Robert L.↗

Intermittent turbulence in solar wind from the south polar hole

The magnetic fields measured by the Ulysses spacecraft are used to study solar wind turbulence in the fast solar wind from the south polar hole. The spacecraft was at about 46 deg south latitude and 3.9 AU. For a magnetic field with a Gaussian distribution the power spectrum (second-order structure function) is sufficient to completely characterize the turbulence. However, the actual distribution is non-Gaussian so that the effects of intermittency must be taken into account. The observed spectral exponents include effects of intermittency and cannot be directly compared with the standard second-order spectral theories such as the Kolmogorov and Kraichnan theories. To permit a better comparison of the observations with the theoretical models, we study the structure characteristics of the data. We find the exponents of the second-order structure functions (power spectra) and the higher-order normalized structure functions for the components of the magnetic fields. We show that these sets of exponents can be approximately described by two basic numbers: the spectral exponent and the intermittency exponent. The intermittency exponent characterizes correlation properties of the energy cascade from large to small scales. Before comparing the observations to the theoretically expected values, a reduction must be made to the observed spectral exponent. The amount of the reduction depends on both the intermittency exponent and the model of the energy cascade assumed in the turbulence theory. We reduce the measured spectral indices according to a simple model for Alfven turbulence that is described here. We then compare our reduced spectral indices with second-order spectral theory. The reduced spectral indices for the period range of 1 min to about a half hour are remarkably constant and in good agreement with the value of 3/2. Thus our treatment is self-consistent. Our tentative conclusion is that the high-frequency turbulence appears to agree with the model of random-phased Alfven waves. This tentative conclusion must be tested by further theoretical and observational work.

Ruzmaikin, A. A.↗

Anomalous elastic softening in ferroelectric hafnia under pressure

his study employs first-principles density-functional theory (DFT) calculations to explore the elastic and mechanical properties of ferroelectric hafnia (HfO 2 ) in its polar orthorhombic 𝑃⁢𝑐⁢𝑎⁢2 1 phase under varying hydrostatic pressure conditions up to 30 GPa. Utilizing a plane-wave basis set and Perdew-Burke-Ernzerhof generalized-gradient approximation for solids in our DFT calculations, we investigate both pure and yttrium-substituted HfO 2 . Our findings reveal an anomalous reduction in the 𝐶 33 component of the elastic tensor with increasing pressure, which becomes significant above 15 GPa and signals a potential pressure-driven structural phase transition at higher pressure. The analysis of atomic displacements under pressure sheds light on the unusual mechanical behavior and phase stability of this material. Additionally, we observe a transition from an indirect band gap to a direct band gap with increasing pressure, which could have significant implications for optical applications. Here, the effects of yttrium substitution on the mechanical and electronic properties are further examined, revealing that yttrium substitution softens the elastic response of this material and reduces the electronic band gap. These results enhance our understanding of elastic and mechanical responses of ferroelectric hafnia and its potential for applications in microelectronics, piezoelectric devices, and nonvolatile ferroelectric random-access memories. Further experimental validation is recommended to confirm our predictions and explore the practical implications of the observed phase transitions and electronic behavior of the ferroelectric hafnia under high-pressure conditions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

An experimental study of light scattering by large, irregular particles

The intensity and polarization of light scattered by a variety of types of artificial partices large compared to the wavelength were measured as a function of phase angle. Shape, surface roughness, absorption coefficient, and internal scattering coefficient were varied systematically and their effects studied. Scattering by clear, smooth-surfaced spheres is in quantitative agreement with the predictions of the geometrical optics (ray theory) approximation to physical optics (Mie theory). The phase functions of almost all of the particles measured have both forward and backward scattering lobes. A two-parameter, double Henyey-Greenstein function generally provides reasonably good descriptions of the data, while keeping the number of free parameters to the minimum necessary. On a double Henyey- Greenstein parameter plot all of the particles fall into an L-shaped area of restricted size in which the location is characteristic of the particle type. Formalisms based on the equivalent slab model are also given for estimating the scattering efficiency of a large, irregular particle. For most dielectric particles the transmitted, forward scattered light is partially negatively polarized. It is this component that is respopnsible for the well-known maximum in the polarization curves of planetary regoliths at phase angles around 100 deg. For phase angles between about 30 deg and 70 deg the internally scattered light is found to be randomly polarized in the particles studied here, so that the only contribution to the second component of the Stokes vector is by Fresnel reflection from the particle surface. If this empirical result is general, measurement of the second Stokes vector of the light scattered from a regolith at these angles may provide a method of remotely measuring the mean refractive index.

Mcguire, Audrey F.↗

N-particle irreducible actions for stochastic fluids

We construct one- and two-particle irreducible (1PI and 2PI) effective actions for the stochastic fluid dynamics of a conserved density undergoing diffusive motion. We compute the 1PI action in one-loop order and the 2PI action in two-loop approximation. We derive a set of Schwinger-Dyson equations and regularize the resulting equations using Pauli-Villars fields. We numerically solve the Schwinger-Dyson equations for a non-critical fluid. We find that higher-loop effects summed by the Schwinger-Dyson renormalize the non-linear coupling. We also find indications of a diffuson-cascade, the appearance of n-loop correction with smaller and smaller exponential suppression.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Spectrally Consistent Scattering, Absorption, and Polarization Properties of Atmospheric Ice Crystals at Wavelengths from 0.2 to 100 um

A data library is developed containing the scattering, absorption, and polarization properties of ice particles in the spectral range from 0.2 to 100 microns. The properties are computed based on a combination of the Amsterdam discrete dipole approximation (ADDA), the T-matrix method, and the improved geometric optics method (IGOM). The electromagnetic edge effect is incorporated into the extinction and absorption efficiencies computed from the IGOM. A full set of single-scattering properties is provided by considering three-dimensional random orientations for 11 ice crystal habits: droxtals, prolate spheroids, oblate spheroids, solid and hollow columns, compact aggregates composed of eight solid columns, hexagonal plates, small spatial aggregates composed of 5 plates, large spatial aggregates composed of 10 plates, and solid and hollow bullet rosettes. The maximum dimension of each habit ranges from 2 to 10,000 microns in 189 discrete sizes. For each ice crystal habit, three surface roughness conditions (i.e., smooth, moderately roughened, and severely roughened) are considered to account for the surface texture of large particles in the IGOM applicable domain. The data library contains the extinction efficiency, single-scattering albedo, asymmetry parameter, six independent nonzero elements of the phase matrix (P11, P12, P22, P33, P43, and P44), particle projected area, and particle volume to provide the basic single-scattering properties for remote sensing applications and radiative transfer simulations involving ice clouds. Furthermore, a comparison of satellite observations and theoretical simulations for the polarization characteristics of ice clouds demonstrates that ice cloud optical models assuming severely roughened ice crystals significantly outperform their counterparts assuming smooth ice crystals.

Optical properties↗

The Potential Impact of Satellite-Retrieved Cloud Parameters on Ground-Level PM2.5 Mass and Composition

Satellite-retrieved aerosol optical properties have been extensively used to estimate ground-level fine particulate matter (PM2.5) concentrations in support of air pollution health effects research and air quality assessment at the urban to global scales. However, a large proportion, approximately 70%, of satellite observations of aerosols are missing as a result of cloud-cover, surface brightness, and snow-cover. The resulting PM2.5 estimates could therefore be biased due to this non-random data missingness. Cloud-cover in particular has the potential to impact ground-level PM2.5 concentrations through complex chemical and physical processes. We developed a series of statistical models using the Multi-Angle Implementation of Atmospheric Correction (MAIAC) aerosol product at 1 km resolution with information from the MODIS cloud product and meteorological information to investigate the extent to which cloud parameters and associated meteorological conditions impact ground-level aerosols at two urban sites in the US: Atlanta and San Francisco. We find that changes in temperature, wind speed, relative humidity, planetary boundary layer height, convective available potential energy, precipitation, cloud effective radius, cloud optical depth, and cloud emissivity are associated with changes in PM2.5 concentration and composition, and the changes differ by overpass time and cloud phase as well as between the San Francisco and Atlanta sites. A case-study at the San Francisco site confirmed that accounting for cloud-cover and associated meteorological conditions could substantially alter the spatial distribution of monthly ground-level PM2.5 concentrations.

MAIAC AO↗

Light scattering by ice clouds in the visible and infrared - A theoretical study.

Computations of the intensity and linear polarization for single scattering by ice clouds have been made based on the assumption that the particles in ice clouds can be approximated by long circular cylinders which are allowed to be polydispersive as well as arbitrarily oriented in space. The results of two models of optically thin ice clouds are presented and compared with those for polydisperse ice spheres. The two models for ice cylinders are assumed to be either uniformly or randomly oriented in a horizontal plane. Four different wavelengths, 0.7, 3, 3.5 and 6.05 microns, are employed in the light scattering computations. It is found that, compared to ice spheres, long ice cylinders scatter more light in the region with scattering angles near 90 deg, at the expense of scattering in both the forward and backward directions. The glory and cloudbows, which occur in light scattered by spherical particles, are either lost (the glory) or largely reduced and distorted (the cloudbows in the case of cylinders). It is probable that for more irregular particles the cloudbows would also disappear. These differences in scattering by spherical and non-spherical scatterers therefore provide useful information for the differentiation between the ice and liquid phase of cloud particles.

Liou, K.-N.↗

Fog Intermittency and Critical Behavior

The intermittency of fog occurrence (the switching between fog and no-fog) is a key stochastic feature that plays a role in its duration and the amount of moisture available. Here, fog intermittency is studied by using the visibility time series collected during the month of July 2022 on Sable Island, Canada. In addition to the visibility, time series of air relative humidity and turbulent kinetic energy, putative variables akin to the formation and breakup conditions of fog, respectively, are also analyzed in the same framework to establish links between fog intermittency and the underlying atmospheric variables. Intermittency in the time series is quantified with their binary telegraph approximations to isolate clustering behavior from amplitude variations. It is shown that relative humidity and turbulent kinetic energy bound many stochastic features of visibility, including its spectral exponent, clustering exponent, and the growth of its block entropy slope. Although not diagnostic, the visibility time series displays features consistent with Pomeau–Manneville Type-III intermittency in its quiescent phase duration PDF scaling (−3/2), power spectrum scaling (−1/2), and signal amplitude PDF scaling (−2). The binary fog time series exhibits properties of self-organized criticality in the relation between its power spectrum scaling and quiescent phase duration distribution.

54 ENVIRONMENTAL SCIENCES↗

Sonic horizons and causality in phase transition dynamics

A system gradually driven through a symmetry-breaking phase transition is subject to the Kibble-Zurek mechanism (KZM). As a consequence of the critical slowing down, its state cannot follow local equilibrium, and its evolution becomes nonadiabatic near the critical point. In the simplest approximation, that stage can be regarded as an “impulse” where the state of the system remains unchanged. It leads to the correct KZM scaling laws. However, such a “freeze-out” might suggest that the coherence length of the nascent order parameter remains unchanged as the critical region is traversed. By contrast, the original causality-based discussion emphasized the role of the sonic horizon: Domains of the broken symmetry phase can expand with a velocity limited by the speed of the relevant sound. This effect was demonstrated in the quantum Ising chain where the dynamical exponent z = 1 and quasiparticles excited by the transition have a fixed speed of sound. To elucidate the role of the sonic horizon, in this paper we study two systems with z > 1 where the speed of sound is no longer fixed, and the fastest excited quasiparticles set the size of the sonic horizon. Their effective speed decays with increasing transition time. In the extreme case, the dynamical exponent z can diverge such as in the Griffiths region of the random Ising chain where the localization of excited quasiparticles freezes the growth of the correlation range when the critical region is traversed. Of particular interest is an example with z < 1 , the long-range extended Ising chain, where there is no upper limit to the velocity of excited quasiparticles with small momenta. Initially, the power-law tail of the correlation function grows adiabatically, but in the nonadiabatic stage it lags behind the adiabatic evolution.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Bidirectional Reflectance of Flat, Optically Thick Particulate Layers: An Efficient Radiative Transfer Solution and Applications to Snow and Soil Surfaces

We describe a simple and highly efficient and accurate radiative transfer technique for computing bidirectional reflectance of a macroscopically flat scattering layer composed of nonabsorbing or weakly absorbing, arbitrarily shaped, randomly oriented and randomly distributed particles. The layer is assumed to be homogeneous and optically semi-infinite, and the bidirectional reflection function (BRF) is found by a simple iterative solution of the Ambartsumian's nonlinear integral equation. As an exact Solution of the radiative transfer equation, the reflection function thus obtained fully obeys the fundamental physical laws of energy conservation and reciprocity. Since this technique bypasses the computation of the internal radiation field, it is by far the fastest numerical approach available and can be used as an ideal input for Monte Carlo procedures calculating BRFs of scattering layers with macroscopically rough surfaces. Although the effects of packing density and coherent backscattering are currently neglected, they can also be incorporated. The FORTRAN implementation of the technique is available on the World Wide Web at http://ww,,v.giss.nasa.gov/-crmim/brf.html and can be applied to a wide range of remote sensing, engineering, and biophysical problems. We also examine the potential effect of ice crystal shape on the bidirectional reflectance of flat snow surfaces and the applicability of the Henyey-Greenstein phase function and the 6-Eddington approximation in calculations for soil surfaces.

Mishchenko, Michael I.↗

Dynamics and control of coherent structures in the turbulent wall layer: An overview

We expand the velocity field in the vicinity of the wall in empirical eigenfunctions obtained from experiment. Truncating our system and using Galerkin projection, we obtain a closed set of non-linear ordinary differential equations with ten degrees of freedom. We find a rich dynamical behavior, including in particular a heteroclinic attracting orbit giving rise to intermittency. The intermittent jump from one attracting point to the other resembles in many respects the bursts observed in experiments. Specifically, the time between jumps and the duration of the jumps, is approximately that observed in a burst; the jump begins with the formation of a narrowed and intensified updraft, like the ejection phase of a burst, and is followed by a gentle, diffuse downdraft like the sweep phase of a burst. The magnitude of the Reynolds stress spike produced during a burst is limited by our truncation. The behavior is quite robust, much of it being due to the symmetries present (Aubry's group has examined dimensions up to 128 with persistence of the global behavior). We have examined eigenvalues and coefficients obtained from experiment, and from exact simulation, which differ in magnitude. Similar behavior is obtained in both cases; in the latter case, the heteroclinic orbits connect limit cycles instead of fixed points, corresponding to cross-stream waving of the streamwise rolls. The bifurcation diagram remains structurally similar, but somewhat distorted. The role of the pressure term is made clear - it triggers the intermittent jumps, which otherwise would occur at longer and longer intervals, as the system trajectory is attracted closer and closer to the heteroclinic cycle. The pressure term results in the jumps occurring at essentially random times, and the magnitude of the signal determines the average timing. Stretching of the wall region shows that the model is consistent with observations of polymer drag reduction. Change of the third order coefficients, corresponding to acceleration or deceleration of the mean flow, changes the heteroclinic cycles from attracting to repelling, increasing or decreasing the stability, in agreement with observations. The existence of fixed points is an artifact introduced by the projection; however, a decoupled model still displays the rich dynamics. Numerous assumptions made in Aubry et al. (1988) can now be proved exactly. Feeding back eigenfuncitons with the proper phase can delay the bursting, (the heteroclinic jump to the other fixed point), decreasing the drag. It is also possible to speed up the bursting, increasing mixing to control separation. Our approach is optimal for short time tracking in control.

Berkooz, Gal↗

Solar protons E greater than 100 Mev incident over Antarctica during January- February 1967.

Commencing at 0825 +a -1 UT on January 28, 1967, a large and prolonged increase in the intensity of penetrating charged particles was observed by balloon-borne instruments floating over Byrd Station, Antarctica. (80°S, 120°W). A peak intensity of approximately 50 protons per cm 2 -sec-steradian with E> 100 MeV occurred at about 1230 UT on the 28th. The event was under observation almost continuously over a period of about 100 hours until the intensity decayed below cosmic-ray background on February 1. The initial decay was rapid but, some 40 hours after onset, went over into a slow exponential decay characterized by a 20 hour time-constant. The decay phase of an additional, though considerably less intense, event was observed on February 3 and 4. Presumably both events had their origins in major disturbances on the far side of the sun since nether event has been definitely linked to any feature which existed on the visible disk within an appropriate time interval. Results pertaining to the time-intensity profile and to the energy spectrum for protons E> 100 MeV are presented for the January 29 event. Comparison of the balloon results with neutron-monitor and satellite measurements and with models of interplanetary diffusion has led to some conclusions regarding the role of small-angle scattering by irregularities and by the random walk of magnetic lines of force relative to the mean interplanetary field within the orbit of earth.

Energy spectrum↗

A phase response curve to single bright light pulses in human subjects

The circadian pacemaker is differentially sensitive to the resetting effects of retinal light exposure, depending upon the circadian phase at which the light exposure occurs. Previously reported human phase response curves (PRCs) to single bright light exposures have employed small sample sizes, and were often based on relatively imprecise estimates of circadian phase and phase resetting. In the present study, 21 healthy, entrained subjects underwent pre- and post-stimulus constant routines (CRs) in dim light (approximately 2-7 lx) with maintained wakefulness in a semi-recumbent posture. The 6.7 h bright light exposure stimulus consisted of alternating 6 min fixed gaze (approximately 10 000 lx) and free gaze (approximately 5000-9000 lx) exposures. Light exposures were scheduled across the circadian cycle in different subjects so as to derive a PRC. Plasma melatonin was used to determine the phase of the onset, offset, and midpoint of the melatonin profiles during the CRs. Phase shifts were calculated as the difference in phase between the pre- and post-stimulus CRs. The resultant PRC of the midpoint of the melatonin rhythm revealed a characteristic type 1 PRC with a significant peak-to-trough amplitude of 5.02 h. Phase delays occurred when the light stimulus was centred prior to the critical phase at the core body temperature minimum, phase advances occurred when the light stimulus was centred after the critical phase, and no phase shift occurred at the critical phase. During the subjective day, no prolonged 'dead zone' of photic insensitivity was apparent. Phase shifts derived using the melatonin onsets showed larger magnitudes than those derived from the melatonin offsets. These data provide a comprehensive characterization of the human PRC under highly controlled laboratory conditions.

Randomized Controlled Trial↗

Simulations of DSB Yields and Radiation-induced Chromosomal Aberrations in Human Cells Based on the Stochastic Track Structure Induced by HZE Particles

The formation of double‐strand breaks (DSBs) and chromosomal aberrations (CAs) is of great importance in radiation research and, specifically, in space applications. We are presenting a new particle track and DNA damage model, in which the particle stochastic track structure is combined with the random walk (RW) structure of chromosomes in a cell nucleus. The motivation for this effort stems from the fact that the model with the RW chromosomes, NASARTI (NASA radiation track image) previously relied on amorphous track structure, while the stochastic track structure model RITRACKS (Relativistic Ion Tracks) was focused on more microscopic targets than the entire genome. We have combined chromosomes simulated by RWs with stochastic track structure, which uses nanoscopic dose calculations performed with the Monte‐Carlo simulation by RITRACKS in a voxelized space. The new simulations produce the number of DSBs as function of dose and particle fluence for high‐energy particles, including iron, carbon and protons, using voxels of 20 nm dimension. The combined model also calculates yields of radiation‐induced CAs and unrejoined chromosome breaks in normal and repair deficient cells. The joined computational model is calibrated using the relative frequencies and distributions of chromosomal aberrations reported in the literature. The model considers fractionated deposition of energy to approximate dose rates of the space flight environment. The joined model also predicts of the yields and sizes of translocations, dicentrics, rings, and more complex‐type aberrations formed in the G0/G1 cell cycle phase during the first cell division after irradiation. We found that the main advantage of the joined model is our ability to simulate small doses: 0.05‐0.5 Gy. At such low doses, the stochastic track structure proved to be indispensable, as the action of individual delta‐rays becomes more important.

Ponomarev, Artem↗

Simulations of DSB Yields and Radiation-induced Chromosomal Aberrations in Human Cells Based on the Stochastic Track Structure iIduced by HZE Particles

The formation of double‐strand breaks (DSBs) and chromosomal aberrations (CAs) is of great importance in radiation research and, specifically, in space applications. We are presenting a new particle track and DNA damage model, in which the particle stochastic track structure is combined with the random walk (RW) structure of chromosomes in a cell nucleus. The motivation for this effort stems from the fact that the model with the RW chromosomes, NASARTI (NASA radiation track image) previously relied on amorphous track structure, while the stochastic track structure model RITRACKS (Relativistic Ion Tracks) was focused on more microscopic targets than the entire genome. We have combined chromosomes simulated by RWs with stochastic track structure, which uses nanoscopic dose calculations performed with the Monte‐Carlo simulation by RITRACKS in a voxelized space. The new simulations produce the number of DSBs as function of dose and particle fluence for high‐energy particles, including iron, carbon and protons, using voxels of 20 nm dimension. The combined model also calculates yields of radiation‐induced CAs and unrejoined chromosome breaks in normal and repair deficient cells. The joined computational model is calibrated using the relative frequencies and distributions of chromosomal aberrations reported in the literature. The model considers fractionated deposition of energy to approximate dose rates of the space flight environment. The joined model also predicts of the yields and sizes of translocations, dicentrics, rings, and more complex‐type aberrations formed in the G0/G1 cell cycle phase during the first cell division after irradiation. We found that the main advantage of the joined model is our ability to simulate small doses: 0.05‐0.5 Gy. At such low doses, the stochastic track structure proved to be indispensable, as the action of individual delta‐rays becomes more important.

Ponomarev, Artem↗