Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “independent particle 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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
NASA NTRS1991

A model for the frequency fine structure of auroral kilometric radiation

A model for AKR frequency fine structure in which the fine structure is a natural consequence of the emission of electron cyclotron maser radiation in a nonuniform magnetic field is presented. The model is examined by using 1D electromagnetic particle-in-cell simulations. It is found that maser radiation is emitted as individual wave packets combine to form the drifting emission features that make up the fine structure. Rising frequency drifts are associated with wave packets emitted independently of one another. Wave packets that interact with one another produce features with falling frequency drift. The emission features produced have a large spread in both positive and negative frequency drift rates. The average positive drift rate of rapidly drifting features is found to be approximately proportional to the maser growth rate.

Mckean, M. E.↗
DOE OSTI2026

Primordial black hole dark matter: A quantitative parameter sensitivity comparison across formation mechanisms and particle candidates

Primordial black holes (PBHs) in the asteroid-mass window ( 10 17 – 10 22 g ) can account for all of the dark matter without violating any observational constraint, yet are routinely dismissed as fine-tuned. I put that dismissal to the test by applying three complementary sensitivity measures uniformly across a broad landscape: three noninflationary PBH production mechanisms, six classes of inflationary PBH models, and seven particle dark matter benchmarks, all evaluated against the same observable target. Three distinct naturalness universality classes emerge, determined entirely by the analytic structure of the abundance map rather than by the nature of the dark matter candidate. Biased-domain-wall PBHs, in their least model-dependent (free- V b ) form, have the same low sensitivity, Δ = 4.5 , as off-resonance weakly interacting massive particles and freeze-in particles ( Δ = 2 ), a sensitivity that, because it is constant over the entire parameter space of the construction, also coincides trivially with its own Wilson-normalized average within that parameter space (Section Definition and conventions), an equivalence that concerns only the space over which Δ is computed and is not a naturalness statement about the construction as a whole; a further reduction to Δ = 2 is possible only under the additional, independently motivated but not required, assumption that the domain-wall bias is generated by Planck-suppressed operators; early matter-domination PBHs occupy an intermediate tier alongside coannihilating weakly interacting massive particles (WIMPs), unified by a structural identity in which the sensitivity measure equals the logarithm of the ratio of the formation scale to the matter–radiation equality scale; first-order phase transition PBHs, once the more accurate super-exponential collapse probability is used in place of the single-exponential approximation, instead belong to the same highly sensitive tier as resonant WIMP annihilation and single-field inflationary collapse, for a structurally distinct reason; single-field ultraslow-roll inflationary collapse is severely tuned for a distinct reason: a double exponential in which the power spectrum amplitude is itself exponentially sensitive to the inflaton potential coefficients, on top of the exponential collapse sensitivity of the abundance map. My main conclusion is that the claim that PBH dark matter is generically fine-tuned conflates the worst case with a landscape spanning every naturalness tier. The Barbieri-Giudice sensitivity computed here and the Wilson-normalized measure of Iovino and Riotto answer distinct and mutually consistent questions about the same construction, a distinction I clarify within the two-layer decomposition.

Profumo, Stefano [University of California, Santa ↗
NASA NTRS1975

The effect of suprathermal protons on the physical conditions in Seyfert galaxy nuclei II

The radiative transfer of Ly-alpha, Ly-beta, and H-alpha in a hydrogen gas containing dust was considered with application to the nuclear gas in Seyfert galaxies. By neglecting the direct escape of line radiation and by averaging over the gas, the radiation transfer in space is suppressed and the frequency transfer only is considered. The dust degrades the line radiation via frequency-independent absorption, converting the energy to infrared luminosity. The source functions in the lines were solved, using appropriate approximations, in order to determine under what conditions the narrow component of the Balmer line radiation from the gas can be self-absorbed and degraded without similar degradation of the broad component, which originates from the suprathermals themselves. The results are used to find self-consistent values for the temperature and ionization of the gas for various amounts of dust and various concentrations of suprathermal particles.

Stoner, R.↗
DOE OSTI2020

First implementation of gyrokinetic exact linearized Landau collision operator and comparison with models

Gyrokinetic simulations are fundamental to understanding and predicting turbulent transport in magnetically confined fusion plasmas. Previous simulations have used model collision operators with approximate field-particle terms of unknown accuracy and/or have neglected collisional finite Larmor radius (FLR) effects. We have implemented the linearized Fokker–Planck collision operator with exact field-particle terms and full FLR effects in a gyrokinetic code (GENE). The new operator, referred to as “exact” in this paper, allows the accuracy of model collision operators to be assessed. The conservative Landau form is implemented because its symmetry underlies the conservation laws and the H-theorem, and enables numerical methods to preserve this conservation, independent of resolution. The implementation utilizes the finite-volume method recently employed to discretize the Sugama collision model in GENE, allowing direct comparison between the two operators. Results show that the Sugama model appears accurate for the growth rates of trapped electron modes (TEMs) driven only by density gradients, but appreciably underestimates the growth rates as the collisionality and electron temperature gradient increase. The TEM turbulent fluxes near the nonlinear threshold using the exact operator are similar to the Sugama model for the n e = d ln T e /d ln n e = 0 case, but substantially larger than the Sugama model for the n e = 1 case. The FLR effects reduce the growth rates increasingly with wavenumber deepening a “valley” at the intermediate binormal wavenumber as the unstable mode extends from the TEM regime to the electron temperature gradient instability regime. Application to the Hinton–Rosenbluth problem shows that zonal flows decay faster as the radial wavenumber increases and the exact operator yields weaker decay rates.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Non-asymptotic analysis of ensemble Kalman updates: effective dimension and localization

Many modern algorithms for inverse problems and data assimilation rely on ensemble Kalman updates to blend prior predictions with observed data. Ensemble Kalman methods often perform well with a small ensemble size, which is essential in applications where generating each particle is costly. This paper develops a non-asymptotic analysis of ensemble Kalman updates, which rigorously explains why a small ensemble size suffices if the prior covariance has moderate effective dimension due to fast spectrum decay or approximate sparsity. Here, we present our theory in a unified framework, comparing everal implementations of ensemble Kalman updates that use perturbed observations, square root filtering and localization. As part of our analysis, we develop new dimension-free covariance estimation bounds for approximately sparse matrices that may be of independent interest.

Mathematics↗
DOE OSTI2020

Average nonlinear dynamics of particles in gravitational pulses: Effective Hamiltonian, secular acceleration, and gravitational susceptibility

Particles interacting with a prescribed quasimonochromatic gravitational wave (GW) exhibit secular (average) nonlinear dynamics that can be described by Hamilton’s equations. We derive the Hamiltonian of this “ponderomotive” dynamics to the second order in the GW amplitude for a general background metric. For the special case of vacuum GWs, we show that our Hamiltonian is equivalent to that of a free particle in an effective metric, which we calculate explicitly. We also show that already a linear GW pulse displaces a particle from its unperturbed trajectory by a finite distance that is independent of the GW phase and proportional to the integral of the pulse intensity.This effect is independent from the nonlinear memory effects that has been known. We calculate the particle displacement analytically and show that our result is in agreement with numerical simulations. We also show how the Hamiltonian of the nonlinear averaged dynamics naturally leads to the concept of the linear gravitational susceptibility of a particle gas with an arbitrary phase-space distribution. We calculate this susceptibility explicitly to apply it, in a follow-up paper, toward studying self-consistent GWs in inhomogeneous media within the geometrical-optics approximation

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1974

Quantum mechanical streamlines. I - Square potential barrier

Exact numerical calculations are made for scattering of quantum mechanical particles hitting a square two-dimensional potential barrier (an exact analog of the Goos-Haenchen optical experiments). Quantum mechanical streamlines are plotted and found to be smooth and continuous, to have continuous first derivatives even through the classical forbidden region, and to form quantized vortices around each of the nodal points. A comparison is made between the present numerical calculations and the stationary wave approximation, and good agreement is found between both the Goos-Haenchen shifts and the reflection coefficients. The time-independent Schroedinger equation for real wavefunctions is reduced to solving a nonlinear first-order partial differential equation, leading to a generalization of the Prager-Hirschfelder perturbation scheme. Implications of the hydrodynamical formulation of quantum mechanics are discussed, and cases are cited where quantum and classical mechanical motions are identical.

Hirschfelder, J. O.↗
NASA NTRS1977

Test particle propagation in magnetostatic turbulence. 3: The approach to equilibrium

The asymptotic behavior, for large time, of the quasi-linear diabatic solutions and their local approximations is considered. A time averaging procedure is introduced which yields the averages of these solutions over time intervals which contain only large time values. A discussion of the quasi-linear diabatic solutions which is limited to those solutions that are bounded from below as functions of time is given. It is shown that as the upper limit of the time averaging interval is allowed to approach infinity the time averaged quasi-linear diabatic solutions must approach isotropy (mu-independence). The first derivative with respect to mu of these solutions is also considered. This discussion is limited to first derivatives which are bounded functions of time. It is shown that as the upper limit of the time averaging interval is allowed to approach infinity, the time averaged first derivative must approach zero everywhere in mu except at mu = 0 where it must approach a large value which is calculated. The impact of this large derivative on the quasi-linear expansion scheme is discussed. An H-theorem for the first local approximation to the quasi-linear diabatic solutions is constructed. Without time averaging, the H-theorem is used to determine sufficient conditions for the first local approximate solutions to asymptote, with increasing time, to exactly the same final state which the time averaged quasi-linear diabatic solutions must approach as discussed above.

Klimas, A. J.↗
DOE OSTI2022

Modified gravitational wave propagation with higher modes and its degeneracies with lensing

Low-energy alternatives to General Relativity (GR) generically modify the phase of gravitational waves (GWs) during their propagation. As detector sensitivities increase, it becomes key to understand how these modifications affect the GW higher modes and to disentangle possible degeneracies with astrophysical phenomena. We apply a general formalism — the WKB approach — for solving analytically wave propagation in the spatial domain with a modified dispersion relation (MDR). We compare this WKB approach to applying a stationary phase approximation (SPA) in the temporal domain with time delays associated to the group or particle velocity. To this end, we extend the SPA to generic signals with higher modes, keeping careful track of reference phases and arrival times. We find that the WKB approach coincides with the SPA using the group velocity, in agreement with the principles of wave propagation. We then explore the degeneracies between a GW propagation with an MDR and a strongly-lensed GW in GR, since the latter can introduce a frequency-independent phase shift which is not degenerate with source parameters in the presence of higher modes. We find that for a particular MDR there is an exact degeneracy for wave propagation, unlike with the SPA for particle propagation. For the other cases, we search for the values of the MDR parameters that minimize the χ 2 and conclude that strongly-lensed GR GWs could be misinterpreted as GWs in modified gravity. As a result, future MDR constraints with higher mode GWs should include the possibility of frequency-independent phase shifts, allowing for the identification of modified gravity and strong lensing distortions at the same time.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1999

Cloud Thickness from Diffusion of Lidar Pulses in Clouds

Measurements of the distribution of reflected light from a laser beam incident on an aqueous suspension of particles or "cloud" with known thickness and particle size distribution are reported. The distribution is referred to as the "cloud radiative Green's function", G. In the diffusion domain, G is sensitive to cloud thickness, allowing that important quantity to be retrieved. The goal of the laboratory simulation is to provide preliminary estimates of sensitivity of G to cloud thickness,for use in the optimal design of an offbeam Lidar instrument for remote sensing of cloud thickness (THOR, Thickness from Offbeam Returns). These clouds of polystyrene microspheres suspended in water are analogous to real clouds of water droplets suspended in air. The microsphere size distribution is roughly lognormal, from 0.5 microns to 25 microns, similar to real clouds. Density of suspended spheres is adjusted so mean-free-path of visible photons is about 10 cm, approximately 1000 times smaller than in real clouds. The light source is a ND:YAG laser at 530 nm. Detectors are flux and photon-counting Photomultiplier Tube (PMTS), with a glass probe for precise positioning. A Labview 5 VI controls positioning, and data acquisition, via an NI Motion Control board connected to a stepper motor driving an Edmund linear slider, and a 16-channel 16-bit NI-DAQ board. The stepper motor is accurate to 10 microns, and step size is selectable from the VI software. Far from the incident beam, the rate of exponential increase as the direction of the incident beam is approached scales as expected from diffusion theory, linearly with the cloud thickness, and inversely as the square root of the reduced optical thickness, and is independent of particle size. Near the beam the signal begins to increase faster than exponential, due to single and low-order scattering near the backward direction, and here the distribution depends on particle size. Results are being used to verify 3D Monte Carlo radiative transfer simulations, used to estimate signal-to-noise ratios for remotely sensed off beam returns, for both homogeneous and inhomogeneous clouds. Signal-to-noise estimates show that unfiltered observations are straight forward at night, while narrow band pass filters are being studied for day.

Cahalan, Robert F.↗
DOE OSTI2023

Methane, ethane, and propane production in Greenland ice core samples and a first isotopic characterization of excess methane

Air trapped in polar ice provides unique records of the past atmospheric composition ranging from key greenhouse gases such as methane (CH 4 ) to short-lived trace gases like ethane (C 2 H 6 ) and propane (C 3 H 8 ). Recently, the comparison of CH 4 records obtained using different extraction methods revealed disagreements in the CH 4 concentration for the last glacial in Greenland ice. Elevated methane levels were detected in dust-rich ice core sections measured discretely, pointing to a process sensitive to the melt extraction technique. To shed light on the underlying mechanism, we performed targeted experiments and analyzed samples for methane and the short-chain alkanes ethane and propane covering the time interval from 12 to 42 kyr. Here, we report our findings of these elevated alkane concentrations, which scale linearly with the amount of mineral dust within the ice samples. The alkane production happens during the melt extraction step of the classic wet-extraction technique and reaches 14 to 91 ppb of CH 4 excess in dusty ice samples. We document for the first time a co-production of excess methane, ethane, and propane, with the observed concentrations for ethane and propane exceeding their past atmospheric background at least by a factor of 10. Independent of the produced amounts, excess alkanes were produced in a fixed molar ratio of approximately 14 : 2 : 1, indicating a shared origin. The measured carbon isotopic signature of excess methane is (-47.0 ± 2.9) ‰ and its deuterium isotopic signature is (-326 ± 57) ‰. With the co-production ratios of excess alkanes and the isotopic composition of excess methane we established a fingerprint that allows us to constrain potential formation processes. This fingerprint is not in line with a microbial origin. Moreover, an adsorption–desorption process of thermogenic gas on dust particles transported to Greenland does not appear very likely. Instead, the alkane pattern appears to be indicative of abiotic decomposition of organic matter as found in soils and plant leaves.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1977

Techniques of Global analysis applied to gravitation theories: A cosmological black hole?

An elementary model of freely falling observers and emitters within a black hole's radius is examined to determine the redshift spectrum reaching a typical observer. The model is independent of scale, the fundamental unit being the radius (mass) of the black hole. The observers/emitters all follow the same kinds of trajectories: radially inward and starting from rest at spatial infinity. The test-particle role is assumed throughout; i.e., the observers/emitters do not themselves contribute to the gravitational field of the system. By means of redshift formulas and luminosity distance to the emitters, a picture of actual redshifts and blueshifts, with their intensities, emerges for an observer within the black hole's radius. No luminosity distances greater than approximately one-half the radius are considered in this particular study; nevertheless, redshifts and blueshifts up to approximately 0.6 are seen in portions of the observer's celestial sphere. An exotic application can be made, as a curiosity, to a black hole the size of the universe, resulting in a particular anisotropic "cosmology."

Debney, G.↗
DOE OSTI2024

A Vibrational Energy Harvesting Sensor Based on Linear and Rotational Electromechanical Effects

In this investigation, a magnetically coupled double-spring design is presented for harvesting low-level non-stationary random vibrational energy. The sensor relies on multimodal coupling between the translation and rotation of a two-spring magnet and coil system to widen the harvesting bandwidth. Energy methods are used to develop a model to characterize the electromechanical response of the system, the solution of which is obtained using stochastic techniques based on a particle swarm algorithm. This approach provides an efficient method to estimate system parameters that otherwise are difficult or impossible to determine with independent measurements. The experimental results demonstrate agreement with the theoretical predictions over a limited bandwidth. The sensor can effectively harvest non-stationary vibration energy down to 10 -4 g within a limited bandwidth of 130–150 Hz. The sensor prototype has an operational volume of 2.6 cm 3 with a calculated power density of 0.2 W/cm 3 . The sensor’s small size results in a coupling efficiency of approximately 6% across the tested bandwidth.

42 ENGINEERING↗
NASA NTRS1995

Variation of fractional electron density fluctuations inside 40 R(sub 0) observed by Ulysses ranging measurements

The first measurements of fractional electron density fluctuations delta-n(sub e)/n(sub e), where delta-n(sub e) is rms electron density fluctuation and n(sub e) is the mean electron density, have been carried out inside 40 R(sub 0) using 1991 Ulysses dual-frequency S- and X-band (13 and 3.6 cm) ranging (time delay) measurements. In the frequency band of approximately 6 x 10(exp -5) - 8 x 10(exp -4) Hz (periods of 20 min to 5 hr), delta-n(sub e)/n(sub e) varies from a high near 20% in the slow wind close to the neutral line to a low of 1% in the fast wind far from the neutral line. For spatial wavenumber K approximately = 1.4 x 10(exp -6)/km (period of 5 hr at 250 km/s), delta-n(sub e)/n(sub e) is essentially independent of heliocentric distance over 0.03-1.0 AU in the slow wind; it is a factor of 30 lower in the fast wind than in the slow wind inside 0.1 AU, but exhibits dramatic growth with heliocentric distance inside 0.3 AU. This latter result reinforces current views of the evolution of MHD turbulence and the association of Alfven waves with high speed streams based on in situ fields and particles measurements beyond 0.3 AU. That regions of enhanced density fluctuations near or above the neutral line coincide with regions of enhanced density confirms previous conclusions that they are the interplanetary manifestation of the heliospheric current sheet and extensions of coronal streamers. While the regions of enhanced density fluctuations lie within those of enhanced density, they have boundaries that are distinctly more abrupt, suggesting the separation of plasma of different nature and origin.

Woo, Richard↗
NASA NTRS2019

New Horizons in Cosmology with Spectral Distortions of the Cosmic Microwave Background

Following the pioneering observations with COBE in the early 1990s, studies of the cosmic mi- crowave background (CMB) have primarily focused on temperature and polarization anisotropies. CMB spectral distortions – tiny departures of the CMB energy spectrum from that of a perfect blackbody – provide a second, independent probe of fundamental physics, with a reach deep into the primordial Universe. The theoretical foundation of spectral distortions has seen major advances in recent years, highlighting the immense potential of this emerging field. Spectral distortions probe a fundamental property of the Universe – its thermal history – thereby providing additional insight into processes within the cosmological standard model(I) (CSM) as well as new physics beyond. Spectral distortions are an important tool for understanding inflation and the nature of dark matter. They shed new light on the physics of recombination and reionization, both prominent stages in the evolution of our Universe, and furnish critical information on baryonic feedback processes, in addition to probing primordial correlation functions at scales inaccessible to other tracers. In principle the range of signals is vast: many orders of magnitude of discovery space can be explored by detailed observations of the CMB energy spectrum. Several CSM signals are predicted and provide clear experimental targets that are observable with present-day technology. Confirmation of these signals would extend the reach of the CSM by orders of magnitude in physical scale as the Universe evolves from the initial stages to its present form. Their absence would pose a huge theoretical challenge, immediately pointing to new physics. Here, we advocate for a dedicated effort to measure CMB spectral distortions at the largest angular scales (greater than approximately 1°) within the ESA Voyage 2050 Program. We argue that an L-class mission with a pathfinder would allow a precise measurement of all the expected CSM distortions. With an M-class mission, the primordial distortions (created at z >~ 10(exp 3)) would still be detected at modest significance, while the late-time distortions will continue to be measured to high accuracy. Building on the heritage of COBE/FIRAS, a spectrometer that consists of multiple, cooled (approximately equal to 0.1 K), absolutely-calibrated Fourier Transform Spectrometers (FTS) with wide frequency coverage (ν approximately equal to 10 GHz to a few x THz) and all-sky spectral sensitivity at the level of 0.1 0.5 Jy/sr would be the starting point for the M-class option. A scaled and further optimized version of this concept is being envisioned as the L-class option. Such measurements can only be done from space and would deliver hundreds of absolutely-calibrated maps of the Universe at large scales, opening numerous science opportunities for cosmology and astrophysics. This will provide independent probes of inflation, dark matter and particle physics, recombination and the energy output of our Universe from at late times, turning the long-standing spectral distortion limits of COBE/FIRAS into clear detections.

Chluba, J.↗
NASA NTRS2003

Combined Radiometer-Radar Microphysical Profile Estimations with Emphasis on High Frequency Brightness Temperature Observations

Information about the vertical microphysical cloud structure is useful in many modeling and predictive practices. Radiometers and radars are used to observe hydrometeor properties. This paper describes an iterative retrieval algorithm that combines the use of airborne active and wideband (10 to 340 GHz) passive observations to estimate the vertical content and particle size distributions of liquid and frozen hydrometeors. The physically-based retrieval algorithm relies on the high frequencies (greater than 89 GHz) to provide details on the frozen hydrometeors. Neglecting the high frequencies yielded acceptable estimates of the liquid profiles, but the ice profiles were poorly retrieved. Airborne radar and radiometer observations from the third Convection and Moisture EXperiment (CAMEX-3) were used in the retrieval algorithm as constraints. Nadir profiles were estimated for a minute each of flight time (approximately 12.5 km along track) from an anvil, convection, and quasi- stratiform rain. The complex structure of the frozen hydrometeors required the most iterations for convergence for the anvil cloud type. The wideband observations were found to more than double the estimated frozen hydrometeor content as compared to retrievals using only 90-GHz and below. The convective and quasi-stratiform quickly reached convergence (minimized difference between observations and calculations using the estimated profiles). A qualitative validation using coincident in situ CAMEX-3 observations shows that the retrieved particle size distributions are well corroborated with independent measurements.

Gail Skofronick-Jackson↗
DOE OSTI2021

Scale-fixed predictions for γ + η c production in electron-positron collisions at NNLO in perturbative QCD

In the paper, we present QCD predictions for γ + η c production at an electron-positron collider up to next-to-next-to-leading order (NNLO) accuracy without renormalization scale ambiguities. The NNLO total cross-section for e + + e – → γ + η c using the conventional scale-setting approach has large renormalization scale ambiguities, usually estimated by choosing the renormalization scale to be the e + e – center-of-mass collision energy $\sqrt{s}$. The Principle of Maximum Conformality (PMC) provides a systematic way to eliminate such renormalization scale ambiguities by summing the nonconformal β contributions into the QCD coupling α s ( Q 2 ). The renormalization group equation then sets the value of α s for the process. The PMC renormalization scale reflects the virtuality of the underlying process, and the resulting predictions satisfy all of the requirements of renormalization group invariance, including renormalization scheme invariance. After applying the PMC, we obtain a renormalization scale-and-scheme independent prediction, σ | NNLO,PMC ≃ 41.18 fb for $\sqrt{s}$=10.6 GeV. The resulting pQCD series matches the series for conformal theory and thus has no divergent renormalon contributions. The large K factor which contributes to this process reinforces the importance of uncalculated NNNLO and higher-order terms. Using the PMC scale-and-scheme independent conformal series and the Padé approximation approach, we predict σ | NNNLO,PMC+Pade ≃ 18.99 fb, which is consistent with the recent BELLE measurement ${\sigma}^{\mathrm{obs}}={16.58}_{-9.93}^{+10.51}$ fb at $\sqrt{s}$ ≃ 10.6 GeV. This procedure also provides a first estimate of the NNNLO contribution.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS2001

Monte Carlo predictions of DNA fragment-size distributions for large sizes after HZE particle irradiation

DSBs (double-strand breaks) produced by densely ionizing space radiation are not located randomly in the genome: recent data indicate DSB clustering along chromosomes. DSB clustering at large scales, from >100 Mbp down to approximately 2 kbp, is modeled using a Monte-Carlo algorithm. A random-walk model of chromatin is combined with a track model, that predicts the radial distribution of energy from an ion, and the RLC (randomly-located-clusters) formalism, in software called DNAbreak. This model generalizes the random-breakage model, whose broken-stick fragment-size distribution is applicable to low-LET radiation. DSB induction due to track interaction with the DNA volume depends on the radiation quality parameter Q. This dose-independent parameter depends only weakly on LET. Multi-track, high-dose effects depend on the cluster intensity parameter lambda, proportional to fluence as defined by the RLC formalism. After lambda is determined by a numerical experiment, the model reduces to one adjustable parameter Q. The best numerical fits to the experimental data, determining Q, are obtained. The knowledge of lambda and Q allows us to give biophysically based extrapolations of high-dose DNA fragment-size data to low doses or to high LETs.

NASA Center JSC↗