Engineering PapersSearch

SEARCH · Engineering Papers

Results for “GCR”

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 109 records · Page 6

The Average Quality Factors by TEPC for Charged Particles

The quality factor used in radiation protection is defined as a function of LET, Q(sub ave)(LET). However, tissue equivalent proportional counters (TEPC) measure the average quality factors as a function of lineal energy (y), Q(sub ave)(Y). A model of the TEPC response for charged particles considers energy deposition as a function of impact parameter from the ion s path to the volume, and describes the escape of energy out of sensitive volume by delta-rays and the entry of delta rays from the high-density wall into the low-density gas-volume. A common goal for operational detectors is to measure the average radiation quality to within accuracy of 25%. Using our TEPC response model and the NASA space radiation transport model we show that this accuracy is obtained by a properly calibrated TEPC. However, when the individual contributions from trapped protons and galactic cosmic rays (GCR) are considered; the average quality factor obtained by TEPC is overestimated for trapped protons and underestimated for GCR by about 30%, i.e., a compensating error. Using TEPC's values for trapped protons for Q(sub ave)(y), we obtained average quality factors in the 2.07-2.32 range. However, Q(sub ave)(LET) ranges from 1.5-1.65 as spacecraft shielding depth increases. The average quality factors for trapped protons on STS-89 demonstrate that the model of the TEPC response is in good agreement with flight TEPC data for Q(sub ave)(y), and thus Q(sub ave)(LET) for trapped protons is overestimated by TEPC. Preliminary comparisons for the complete GCR spectra show that Q(sub ave)(LET) for GCR is approximately 3.2-4.1, while TEPC measures 2.9-3.4 for QQ(sub ave)(y), indicating that QQ(sub ave)(LET) for GCR is underestimated by TEPC.

Kim, Myung-Hee Y.

Evaluating Shielding Effectiveness for Reducing Space Radiation Cancer Risks

We discuss calculations of probability distribution functions (PDF) representing uncertainties in projecting fatal cancer risk from galactic cosmic rays (GCR) and solar particle events (SPE). The PDF s are used in significance tests of the effectiveness of potential radiation shielding approaches. Uncertainties in risk coefficients determined from epidemiology data, dose and dose-rate reduction factors, quality factors, and physics models of radiation environments are considered in models of cancer risk PDF s. Competing mortality risks and functional correlations in radiation quality factor uncertainties are treated in the calculations. We show that the cancer risk uncertainty, defined as the ratio of the 95% confidence level (CL) to the point estimate is about 4-fold for lunar and Mars mission risk projections. For short-stay lunar missions (<180 d), SPE s present the most significant risk, however one that is mitigated effectively by shielding, especially for carbon composites structures with high hydrogen content. In contrast, for long duration lunar (>180 d) or Mars missions, GCR risks may exceed radiation risk limits, with 95% CL s exceeding 10% fatal risk for males and females on a Mars mission. For reducing GCR cancer risks, shielding materials are marginally effective because of the penetrating nature of GCR and secondary radiation produced in tissue by relativistic particles. At the present time, polyethylene or carbon composite shielding can not be shown to significantly reduce risk compared to aluminum shielding based on a significance test that accounts for radiobiology uncertainties in GCR risk projection.

Cucinotta, Francis A.

An Analytical Model for the Prediction of a Micro-Dosimeter Response Function

A rapid analytical procedure for the prediction of a micro-dosimeter response function in low Earth orbit (LEO), correlated with the Space Transportation System (STS, shuttle) Tissue Equivalent Proportional Counter (TEPC) measurements is presented. The analytical model takes into consideration the energy loss straggling and chord length distribution of the detector, and is capable of predicting energy deposition fluctuations in a cylindrical micro-volume of arbitrary aspect ratio (height/diameter) by incoming ions through both direct and indirect (ray) events. At any designated (ray traced) target point within the vehicle, the model accepts the differential flux spectrum of Galactic Cosmic Rays (GCR) and/or trapped protons at LEO as input. On a desktop PC, the response function of TEPC for each ion in the GCR/trapped field is computed at the average rate of 30 seconds/ion. The ionizing radiation environment at LEO is represented by O'Neill fs GCR model (2004), covering charged particles in the 1 less than or equal to Z less than or equal to 28. O'Neill's free space GCR model is coupled with the Langley Research Center (LaRC) angular dependent geomagnetic cutoff model to compute the transmission coefficient in LEO. The trapped proton environment is represented by a LaRC developed time dependent procedure which couples the AP8MIN/AP8MAX, Deep River Neutron Monitor (DRNM) and F10.7 solar radio frequency measurements. The albedo neutron environment is represented by the extrapolation of the Atmospheric Ionizing Radiation (AIR) measurements. The charged particle transport calculations correlated with STS 51 and 114 flights are accomplished by using the most recent version (2005) of the LaRC deterministic High charge (Z) and Energy TRaNsport (HZETRN) code. We present the correlations between the TEPC model predictions (response function) and TEPC measured differential/integral spectra in the lineal energy (y) domain for both GCR and trapped protons, with the conclusion that the model correctly accounts for the increase in flux at low y values where energetic ions are the primary contributor. We further discuss that, even with the incorporation of angular dependency in the cutoffs, comparison of the GCR differential/integral flux between STS 51 and 114 TEPC measured data and current calculations indicates that there still exists an underestimation by the simulations at low to mid range y values. This underestimation is partly related the exclusion of the secondary pion particle production from the current version of HZETRN.

Badavi, Francis F.

Simulation and Comparison of Martian Surface Ionization Radiation

The spectrum of energetic particle radiation and corresponding doses at the surface of Mars is being characterized by the Radiation Assessment Detector (RAD), one of ten science instruments on the Mars Science Laboratory (MSL) Curiosity Rover. The time series of dose rate for the first 300 Sols after landing on Mars on August 6, 2012 is presented here. For the comparison to RAD measurements of dose rate, Martian surface ionization radiation is simulated by utilizing observed space quantities. The GCR primary radiation spectrum is calculated by using the Badhwar-O'Neill 2011 (BO11) galactic cosmic ray (GCR) model, which has been developed by utilizing all balloon and satellite GCR measurements since 1955 and the newer 1997-2012 Advanced Composition Explorer (ACE) measurements. In the BO11 model, solar modulation of the GCR primary radiation spectrum is described in terms of the international smoothed sunspot number and a time delay function. For the transport of the impingent GCR primary radiation through Mars atmosphere, a vertical distribution of atmospheric thickness at each elevation is calculated using the vertical profiles of atmospheric temperature and pressure made by Mars Global Surveyor measurements. At Gale Crater in the southern hemisphere, the seasonal variation of atmospheric thickness is accounted for the daily atmospheric pressure measurements of the MSL Rover Environmental Monitoring Station (REMS) by using low- and high-density models for cool- and warm-season, respectively. The spherically distributed atmospheric distance is traced along the slant path, and the resultant directional shielding by Martian atmosphere is coupled with Curiosity vehicle for dose estimates. We present predictions of dose rate and comparison to the RAD measurements. The simulation agrees to within +/- 20% with the RAD measurements showing clearly the variation of dose rate by heliospheric conditions, and presenting the sensitivity of dose rate by atmospheric pressure, which has been found from the RAD experiments and driven by thermal tides on Martian surface.

Kim, Myung-Hee Y.

Radiation Transport Models in Space: from Supernovae to Cells

Humans embarking on deep space exploration missions will encounter persistent exposure to galactic cosmic rays (GCR) - an energetic and highly complex radiation field that is unlike anything found on Earth. Exposure to such radiation fields is attributed to various adverse health effects, including cancer, cardiovascular disease, and cognitive impairment and is identified by NASA as one of the five main hazards of human spaceflight. It is therefore critical to be able to fully characterize the exposure received by humans behind shielding in space and project consequent health risks. A wide variety of computational models have been developed over the years to help meet this requirement. In this talk, an overview of the GCR environment in deep space is provided. Methods of propagating GCR fields through the shielding that protects humans in space are described along with simulation tools used to assess biological damage at the cellular scale. Finally, the NASA cancer risk model is briefly described, and risk projections are provided for various mission scenarios. Radiation transport models and solution methods pervade many aspects of this talk. For example, the GCR spectrum impinging on spacecraft is determined by solving the Fokker-Planck transport equation to propagate cosmic rays (believed to originate from supernovae) from the edge of the heliosphere to the vicinity of Earth. The Boltzmann transport equation is solved to transport this GCR spectrum through shielding and human tissue. Monte Carlo methods are used to simulate the transport of low energy electrons that dominate biological damage at the cellular scale. Progress and challenges in each of these areas will be highlighted.

Tony C. Slaba

Do supernovae of type 1 paly a role in cosmic-ray production?

A model of cosmic-ray origin is suggested which aims to account for some salient features of the composition. Relative to solar abundances, the Galactic cosmic rays (GCR) are deficient in hydrogen and helim (H and He) by an order of magnitude when the two compositions are normalized at iron. Our conjectural model implicates supernovae of Type I (SN-I) as sources of some of the GCR. SN-I occur approximately as often as SN-II, through their genesis is thought to be different. Recent studies of nucleosynthesis in SN-I based on accreting white dwarfs, find that the elements from Si to Fe are produced copiously. On the other hand, SN-I are virtually devoid of hydrogen, and upper limits deduced for He are low. If SN-I contribute significantly to the pool of GCR by injecting energetic particles into the interstellar medium (ISM), then this could explain why the resulting GCR is relatively deficient in H and He. A test of the model is proposed, and difficulties are discussed.

Shapiro, M. M.

Measurement of the cosmic ray iron spectrum above 60 GeV/n

A measurement of the spectral index of Galactic cosmic ray (GCR) iron has been performed using the High Energy Gas Cerenkov Spectrometer. The spectral index of GCR iron, in the energy range 60 to 200 GeV/n, has been determined to be 2.56 + or - 0.11. This result indicates that the source spectrum of GCR iron is similar to that of other primary GCR nuclei, as expected in the standard leaky box model.

Esposito, J. A.

Mars surface radiation exposure for solar maximum conditions and 1989 solar proton events

The Langley heavy-ion/nucleon transport code, HZETRN, and the high-energy nucleon transport code, BRYNTRN, are used to predict the propagation of galactic cosmic rays (GCR's) and solar flare protons through the carbon dioxide atmosphere of Mars. Particle fluences and the resulting doses are estimated on the surface of Mars for GCR's during solar maximum conditions and the Aug., Sep., and Oct. 1989 solar proton events. These results extend previously calculated surface estimates for GCR's at solar minimum conditions and the Feb. 1956, Nov. 1960, and Aug. 1972 solar proton events. Surface doses are estimated with both a low-density and a high-density carbon dioxide model of the atmosphere for altitudes of 0, 4, 8, and 12 km above the surface. A solar modulation function is incorporated to estimate the GCR dose variation between solar minimum and maximum conditions over the 11-year solar cycle. By using current Mars mission scenarios, doses to the skin, eye, and blood-forming organs are predicted for short- and long-duration stay times on the Martian surface throughout the solar cycle.

Simonsen, Lisa C.

Galactic cosmic ray transport methods and radiation quality issues

An overview of galactic cosmic ray (GCR) interaction and transport methods, as implemented in the Langley Research Center GCR transport code, is presented. Representative results for solar minimum, exo-magnetospheric GCR dose equivalents in water are presented on a component by component basis for various thicknesses of aluminum shielding. The impact of proposed changes to the currently used quality factors on exposure estimates and shielding requirements are quantified. Using the cellular track model of Katz, estimates of relative biological effectiveness (RBE) for the mixed GCR radiation fields are also made.

Townsend, L. W.

Numerical simulation of cosmogenic nuclide production in lunar rocks

The production rates of cosmogenic nuclides depend on the primary cosmic-ray particles, the irradiated-body's bulk composition, size, and shape, and the sample's composition and shielding depth. Although much work has been done on some of these dependencies, more detailed studies still need to be done on others. This work describes the influence of irradiation geometry on nuclide production in lunar rocks. In most cases, computer simulations of cosmogenic nuclide production were restricted to spherical objects irradiated with a 4 pi isotropic flux (meteoroids) or in lunar core samples irradiated by a 2 pi flux incident on semi-infinite layers or cylinders of huge sizes. Many lunar samples are rocks found on top of the lunar surface. For these rocks, neither of the above-mentioned models correspond to the real conditions. We present results of our simulations of cosmogenic nuclide production in models simulating the irradiation of rocks sitting on top of the lunar surface. The Galactic Cosmic Rays (GCR) production profiles in lunar rocks were calculated using the Los Alamos 3-D Monte Carlo LAHET Code System (LCS). The irradiated object was modeled as the union of a sphere with the radius of the Moon and a small hemisphere with radii varying from 10 to 100 g/sq cm simulating the lunar rock. These calculations for the production of cosmogenic nuclides in lunar rocks by GCR particle show that there are important differences between the results obtained by commonly used geometric irradiation models and the lunar-rock models presented. The steeper GCR production profiles for a rock could help to explain the poor agreement for Be-10 in rock 68815, where slab models give GCR profiles flatter than the observed profiles.

Reedy, R. C.

Mars Surface Ionizing Radiation Environment: Need for Validation

Protection against the hazards from exposure to ionizing radiation remains an unresolved issue in the Human Exploration and Development of Space (HEDS) enterprise [1]. The major uncertainty is the lack of data on biological response to galactic cosmic ray (GCR) exposures but even a full understanding of the physical interaction of GCR with shielding and body tissues is not yet available and has a potentially large impact on mission costs. "The general opinion is that the initial flights should be short-stay missions performed as fast as possible (so-called 'Sprint' missions) to minimize crew exposure to the zero-g and space radiation environment, to ease requirements on system reliability, and to enhance the probability of mission success." The short-stay missions tend to have long transit times and may not be the best option due to the relatively long exposure to zero-g and ionizing radiation. On the other hand the short-transit missions tend to have long stays on the surface requiring an adequate knowledge of the surface radiation environment to estimate risks and to design shield configurations. Our knowledge of the surface environment is theoretically based and suffers from an incomplete understanding of the physical interactions of GCR with the Martian atmosphere, Martian surface, and intervening shield materials. An important component of Mars surface robotic exploration is the opportunity to test our understanding of the Mars surface environment. The Mars surface environment is generated by the interaction of Galactic Cosmic Rays (GCR) and Solar Particle Events (SPEs) with the Mars atmosphere and Mars surface materials. In these interactions, multiple charged ions are reduced in size and secondary particles are generated, including neutrons. Upon impact with the Martian surface, the character of the interactions changes as a result of the differing nuclear constituents of the surface materials. Among the surface environment are many neutrons diffusing from the Martian surface and especially prominent are energetic neutrons with energies up to a few hundred MeV. Testing of these computational results is first supported by ongoing experiments at the Brookhaven National Laboratory but equally important is the validation to the extent possible by measurements on the Martian surface. Such measurements are limited by power and weight requirements of the specific mission and simplified instrumentation by necessity lacks the full discernment of particle type and spectra as is possible with laboratory experimental equipment. Yet, the surface measurements are precise and a necessary requisite to validate our understanding of the surface environment. At the very minimum the surface measurements need to provide some spectral information on the neutron environment. Of absolute necessity is the precise knowledge of the detector response functions for absolute comparisons between the computational model of the surface environment and the detector measurements on the surface.

Wilson, J. W.

Radiation Transport Properties of Potential In Situ-Developed Regolith-Epoxy Materials for Martian Habitats

Mission crews in space outside the Earth s magnetic field will be exposed to high energy heavy charged particles in the galactic cosmic radiation (GCR). These highly ionizing particles will be a source of radiation risk to crews on extended missions to the Moon and Mars, and the biological effects of and countermeasures to the GCR have to be investigated as part of the planning of exploration-class missions. While it is impractical to shield spacecraft and planetary habitats against the entire GCR spectrum, biological and physical studies indicate that relatively modest amounts of shielding are effective at reducing the radiation dose. However, nuclear fragmentation in the shielding materials produces highly penetrating secondary particles, which complicates the problem: in some cases, some shielding is worse than none at all. Therefore the radiation transport properties of potential shielding materials need to be carefully investigated. One intriguing option for a Mars mission is the use of material from the Martian surface, in combination with chemicals carried from Earth and/or fabricated from elements found in the Martian atmosphere, to construct crew habitats. We have measured the transmission properties of epoxy-Martian regolith composites with respect to heavy charged particles characteristic of the GCR ions which bombard the Martian surface. The composites were prepared at NASA Langley Research Center using simulated Martian regolith, in the process also evaluating fabrication methods which could lead to technologies for in situ fabrication on Mars. Initial evaluation of the radiation shielding properties is made using radiation transport models developed at NASA-LaRC, and the results of these calculations are used to select the composites with the most favorable radiation transmission properties. These candidates are then evaluated at particle accelerators which produce beams of heavy charged particles representative in energy and charge of the radiation at the surface of Mars. The ultimate objective is to develop the models into a design tool for use by mission planners, flight surgeons and radiation health specialists.

Miller, Jack

Application of real-time radiation dosimetry using a new silicon LET sensor

A new type of real-time radiation monitoring device, RRMD-III, consisting of three double-sided silicon strip detectors (DSSDs), has been developed and tested on-board the Space Shuttle mission STS-84. The test succeeded in measuring the linear energy transfer (LET) distribution over the range of 0.2 keV/micrometer to 600 keV/micrometer for 178 h. The Shuttle cruised at an altitude of 300 to 400 km and an inclination angle of 51.6 degrees for 221.3 h, which is equivalent to the International Space Station orbit. The LET distribution obtained for particles was investigated by separating it into galactic cosmic ray (GCR) particles and trapped particles in the South Atlantic Anomaly (SAA) region. The result shows that the contribution in dose-equivalent due to GCR particles is almost equal to that from trapped particles. The total absorbed dose rate during the mission was 0.611 mGy/day; the effective quality factor, 1.64; and the dose equivalent rate, 0.998 mSv/day. The average absorbed dose rates are 0.158 mGy/min for GCR particles and 3.67 mGy/min for trapped particles. The effective quality factors are 2.48 for GCR particles and 1.19 for trapped particles. The absorbed doses obtained by the RRMD-III and a conventional method using TLD (Mg(2)SiO(4)), which was placed around the RRMD-III were compared. It was found that the TLDs showed a lower efficiency, just 58% of absorbed dose registered by the RRMD-III.

NASA Discipline Radiation Health

An Algorithm for the Transport of Anisotropic Neutrons

One major obstacle to human space exploration is the possible limitations imposed by the adverse effect of long-term exposure to the space environment. Even before human spaceflight began, the potentially brief exposure of astronauts to the very intense random solar particle events (SPE) were of great concern. A new challenge appears in deep space exploration from exposure to the low-intensity heavy-ion flux of the galactic cosmic rays (GCR) since the missions are of long duration and the accumulated GCR exposures can be high. Because cancer induction rates increase behind low to rather large thicknesses of aluminum shielding, according to available biological data on mammalian exposures to GCR like ions, the shield requirements for a Mars mission are prohibitively expensive in terms of mission launch costs. Therefore, a critical issue in the Human Exploration and Development of Space enterprise is cost effective mitigation of risk associated with ionizing radiation exposure. In order to estimate astronaut risk to GCR exposure and associated cancer risks and health hazards, it is necessary to do shield material studies. To determine an optimum radiation shield material it is necessary to understand nuclear interaction processes such as fragmentation and secondary particle production which is a function of energy dependent cross sections. This requires knowledge of material transmission characteristics either through laboratory testing or improved theoretical modeling. Here ion beam transport theory is of importance in that testing of materials in the laboratory environment generated by particle accelerators is a necessary step in materials development and evaluation for space use. The approximations used in solving the Boltzmann transport equation for the space setting are often not sufficient for laboratory work and those issues are a major emphasis of the present work.

Tweed, J.

Managing Lunar and Mars Mission Radiation Risks: Cancer Risks, Uncertainties, and Shielding Effectiveness - Part 1

This document addresses calculations of probability distribution functions (PDFs) representing uncertainties in projecting fatal cancer risk from galactic cosmic rays (GCR) and solar particle events (SPEs). PDFs are used to test the effectiveness of potential radiation shielding approaches. Monte-Carlo techniques are used to propagate uncertainties in risk coefficients determined from epidemiology data, dose and dose-rate reduction factors, quality factors, and physics models of radiation environments. Competing mortality risks and functional correlations in radiation quality factor uncertainties are treated in the calculations. The cancer risk uncertainty is about four-fold for lunar and Mars mission risk projections. For short-stay lunar missins (<180 d), SPEs present the most significant risk, but one effectively mitigated by shielding. For long-duration (>180 d) lunar or Mars missions, GCR risks may exceed radiation risk limits. While shielding materials are marginally effective in reducing GCR cancer risks because of the penetrating nature of GCR and secondary radiation produced in tissue by relativisitc particles, polyethylene or carbon composite shielding cannot be shown to significantly reduce risk compared to aluminum shielding. Therefore, improving our knowledge of space radiobiology to narrow uncertainties that lead to wide PDFs is the best approach to ensure radiation protection goals are met for space exploration.

Cucinotta, Francis A.

The Projection of Space Radiation Environments with a Solar Cycle Statistical Model

A solar cycle statistical model has been developed to project sunspot numbers which represent the variations in the space radiation environment. The resultant projection of sunspot numbers in near future were coupled to space-related quantities of interest in radiation protection, such as the galactic cosmic radiation (GCR) deceleration potential (f) and the mean occurrence frequency of solar particle event (SPE). Future GCR fluxes have been derived from a predictive model, in which GCR temporal dependence represented by f was derived from GCR flux and ground-based Climax neutron monitor rate measurements over the last four decades. Results showed that the point dose equivalent inside a typical spacecraft in interplanetary radiation fields was influenced by solar modulation up to a factor of three. One important characteristic of sporadic SPEs is their mean frequency of occurrence, which is dependent on solar activity. Projections of future mean frequency of SPE occurrence were estimated from a power law function of sunspot number. Furthermore, the cumulative probabilities of SPE during short-period missions were defined with the continuous database of proton fluences of SPE. The analytic representation of energy spectra of SPE was constructed by the Weibull distribution for different event sizes. The representative exposure level at each event size was estimated for the guideline of protection systems for astronauts during future space exploration missions.

Kim, Myung-Hee

Observations of the Li, Be, and B Isotopes and Constraints on Cosmic-ray Propagation

The abundance of Li, Be, and B isotopes in galactic cosmic rays (GCR) between E=50-200 MeV/nucleon has been observed by the Cosmic Ray Isotope Spectrometer (CRIS) on NASA's ACE mission since 1997 with high statistical accuracy. Precise observations of Li, Be, B can be used to constrain GCR propagation models. We find that a diffusive reacceleration model with parameters that best match CRIS results (e.g. B/C, Li/C, etc) are also consistent with other GCR observations. A approx. 15-20% overproduction of Li and Be in the model predictions is attributed to uncertainties in the production cross-section data. The latter becomes a significant limitation to the study of rare GCR species that are generated predominantly via spallation.

deNolfo, G. A.

Comparisons of Integrated Radiation Transport Models with Microdosimetry Data in Spaceflight

Astronauts are exposed to galactic cosmic rays (GCR), trapped protons, and possible solar particle events (SPE) during spaceflight. For such complicated mixtures of radiation types and kinetic energies, tissue equivalent proportional counters (TEPC's) represent a simple time-dependent approach for radiation monitoring. Of interest in radiation protection is the average quality factor of a radiation field defined as a function of linear energy transfer, LET, Q(sub ave)(LET). However TEPC's measure the average quality factors as a function of lineal energy (y), Q(sub ave)(y) defined as the average energy deposition in a volume divided by the average chord length of the volume. Lineal energy, y deviates from LET due to energy straggling, delta-ray escape or entry, and nuclear fragments produced in the detector. Using integrated space radiation models that includes the transport code HZETRN/BRYNTRN, the quantum nuclear interaction model, QMSFRG, and results from Monte-Carlo track simulations of TEPC's response to ions, we consider comparisons of model calculations to TEPC results from NASA missions in low Earth orbit and make predictions for lunar and Mars missions. Good agreement between the model and measured spectra from past NASA missions is found. A finding of this work is that TEPC's values for trapped or solar protons of Q(sub ave)(y) range from 1.9-2.5, overestimating Q(sub ave)(LET), which ranges from 1.4-1.6 with both quantities increasing with shielding depth due to nuclear secondaries Comparisons for the complete GCR spectra show that Q(sub ave)(LET) for GCR is approximately 3.5-4.5, while TEPC's measure 2.9-3.4 for Q(sub ave)(y) with the GCR values decreasing with depth as heavy ions are absorbed in shielding material. Our results support the use of TEPC's for space radiation environmental monitoring when computational analysis is used for proper data interpretation.

Cucinotta, Francis A.