1. Introduction
Viruses are intracellular parasites composed of a nucleic acid surrounded by a protein coat, the capsid. Some viruses contain a lipid envelope, derived from the host, surrounding the capsid. The nucleic acid found in viruses can consist of either RNA or DNA. The Coronaviridae, for example contain a single molecule of RNA consisting of about 30 Kilobases in length where Kb stands for unit of measurement of DNA or RNA length used in genetics [
1,
2,
3,
4,
5,
6]. RNA is composed of nucleotides, each containing a sugar (deoxyribose), a nitrogen containing base (adenine, uracil, guanine, and cytosine), and a phosphate group PO
4. Members of the family Coronoviridae measure 80–160 nm in diameter. The phosphate group it is also present in ambient bacterial DNA and RNA which measure more than 1 μm in diameter, and containing different chemical components in fraction term, thus giving us the ability to distinguish between the two, weighting the different spectrum and flux contributions.
Most of the phosphorus is present in the genetic strand of RNA and in particular in the phosphate group but trace amounts are also in viral proteins that contain the amino acid methionine. Phosphorus, on the other hand, is absent in the EPA pollutants listed in
Table 1.
Therefore, phosphorus, as part of the PO
4 group, has been assumed as a marker of no living matter, potentially a virus. A physical model has been designed in order to detect PO
4 groups amount changes related to different viral clusters. The present study aims to detect, by means of the particle Monte Carlo computer code Fluka 2020 [
7,
8] a cluster element spectrum with concentrations from 1 ppm in air to 0.1%, according to a composition of air in
Table 1 and a simplified chemical form of a virus as a biological matter in air as reported in
Figure 1 [
1] with its own atom modelling.
The proposed methodology is based on a subatomic coupled particle analysis (photon-electron) and the subsequent detection on top of their spectra (i.e., fluxes as a function of PO
4 group contamination ppm). The present study refers to previous works by the same authors [
1,
2] regarding RNA–DNA modeling and their identification through the interactions of particles such as muons and photons.
The purpose of this work is among some preliminary activities in view of experimental work to be performed using a photon beam with energy between 110–150 keV energy range due to a cross sections evaluation and its consequent discrimination ability on top of shielding requirements.
In order to justify the involved physics and the results in the 130 keV beam energy configuration a cross sections analysis for nitrogen and phosphorus content in the PO
4 group has been performed taking into account the EPDL97 Library [
8]. Those two elements have a primary role from a photo atomic contribution point of view in the microscopic and macroscopic cross sections acting in the cluster mixture. The model proposed is based on multiple spherical geometry clusters in grid cells whose spheres (11) have different sizes (as a function of ppm viral contamination). However, those clusters are immersed in a 75% nitrogen composition system which has a weighting factor as far as macroscopic cross section evaluation.
The microscopic cross section represents the effective target area of a single target nucleus for an incident particle and units are given in barn (1. barn = 10
−24 cm
2) or cm
2, while the macroscopic cross-section represents the effective target area of all of the nuclei contained in the volume of the material and units are given in cm
−1. The macroscopic cross section can be obtained from the microscopic according to the following equation where N stands for nuclei density
In a mixture case with different chemical elements, it is necessary to determine the macroscopic cross section for each isotope and therefore sum all the individual macroscopic cross-sections. Moreover, both factors (different atomic densities and different cross-sections) must be considered in the calculation of the macroscopic cross-section of the mixture taken into account.
The Avogadro’s number N
0 = 6.022 × 10
23, is the number of particles that is contained in the amount of substance given by one mole. Thus, if M is the molecular weight, the ratio N
0/M equals to the number of molecules in 1 g of the mixture. The number of molecules per cm
3 in the material of density ρ and the macroscopic cross-section for mixtures are given by the equations
In the photon energy (Ep) interval 0.01 keV < EP < 1 MeV the main contribution to the total photon cross section for the Nitrogen is due to Rayleigh, Compton, and photoelectric effects. Moreover, in that energy interval the total photon cross section has a “higher trend” in terms of photon incident probability (from 0.1 to 5 × 105 cm2/g) due to the combination of the previous three effects with the dominant photoelectric component.
The probability involved is about a photon beam interacting with nuclei in the mixture in different layers such as
where Σ
t stands for total cross section and dl(x) is the number of photons interacting in dx, while −dl(x)/I(x) is its probability to interact in the next layer.
From the previous equation follows the photon probability of interaction in dx such as
At 80 keV the photoelectric and the Rayleigh effect drop leaving the Compton effect as dominant component with a cross section in the interval 0.0001–0.01 barn and decreasing until the pair production becomes the dominant component in the total photon cross section flatting the curve to 0.005 barn.
The photon beam source can be provided from an extraction line of a linear/circular electron accelerator. In this case, a simple bremsstrahlung target interaction with the electron beam produces the requested photons. Another possibility is to obtain the photons of fixed energy from a selected radioactive source by natural radioisotope decay.
2. Materials and Methods
The physical model under analysis and its simulation by Fluka 2020 particles computer code is based on a photon beam source of 130 keV accessible from an extraction line of an electron accelerator. The beam interacts with a cylindrical sample volume—with the axis on x—of air of radius r = 5 cm and height h = 10 cm as
s sample tank (
Figure 2,
Figure 3 and
Figure 4), which is analyzed at x = 10 cm through a double plates’ ionization chamber detector or with a fully depleted pnCCD, a special type of charge coupled devices detector with a CsI (Tl) scintillator.
The sample air chemical composition considered for the present and previous studies is shown in
Table 2 [
1,
9].
Among all possible subatomic secondary particles generated, only photons (coming from first and second fluorescence, bremsstrahlung) and electrons were considered (coming from Compton recoil, photo electric, photon Auger, electron Auger, and knock-on), as other particles are negligible.
The photoelectric effect, which is the predominant reaction in the electron creation, consists of the absorption of the incident photon, with energy E, with emission of several fluorescent photons and the ejection or excitation of an orbital electron of binding energy e < E. Photons of first fluorescence are emitted with energy greater than 1 keV; those ones of second fluorescence are still greater than 1 keV and are caused by residual excitation of the first fluorescence process, leading to a second emission.
The analysis considered both electrons and photons, as said before, without neglecting any secondary photon production by performing a photon/electron coupled calculation and by keeping track of the electron/photon mean free path in the media. All the results proposed concern the photon fluxes and spectra (i.e., fluxes as a function of photon energy) of interest, where also all the possible primary and secondary electrons into the sample volume were taken into account.
The virion particles have dimension 20–250 nm and have been assumed incapsulated in lipidic coat as well as in a water film in order to model a ‘human expectoration process’. The virions are described in 11 cluster configurations through a volumetric cell grid (
Figure 5,
Figure 6,
Figure 7,
Figure 8,
Figure 9 and
Figure 10); each cluster is composed of microspheres with a radius of 5 µm and a volume of 5.24 × 10
−7 mm
3 per incapsulated virion particle, with a mutual distance of 1< d
i < 9 cm among the clusters along all the axes and evaluated in the air sample tank at different concentrations from 1 ppm up to 1000 ppm (
Table 3,
Table 4,
Table 5,
Table 6,
Table 7 and
Table 8); this has been taken into account in the physical modeling with multiple layers radius: ∑ (a
i, b
i, c
i) as a function of the number of particles composing each cluster as index i = 1–1000 ppm (
Figure 6,
Figure 7,
Figure 8,
Figure 9 and
Figure 10). In the geometrical model the PO
4 group analysis outcome has been analyzed by control check volumes/surfaces in order to evaluate energy distributions and particle mean free path (photon and electron) (yellow squares,
Figure 5).
The simulations were performed step by step in different cluster stages: Stage 1, with 0 ppm contamination to investigate the physics involved in the basic case without viral particles; Stage 2, evaluating an escalating contamination grade as maximum stress test of 1 ppm, 10 ppm, 100 ppm, 1000 ppm to determine the subatomic particles’ stopping power and the shielding effects that give the photon/electron fluences and energy spectra in different cluster contaminations with an average relative error of 3% counted on seven independent runs. Each run has been tuned, in term of number histories, based on the deviation of the average relative error which can be detected by error spikes in few energy bins due to lack of sufficient particle histories. The tuning process just mentioned had the capability to converge our statistic and to assure seven independent runs are adequate. Although the 1000 ppm case was also simulated, we did not report it in the main results because a case of non-interest inasmuch as there was interest in low contamination detection lower than 10 ppm.
3. Results
Here we present the results of the analysis showing the photon/electron fluences and energy spectra of the Monte Carlo simulations in the presence of virion cluster contaminations and without it in the detector chamber, located at
x = 10 cm on the top of the sample tank on the
x-axis, as shown in
Figure 2,
Figure 3,
Figure 4,
Figure 5,
Figure 6,
Figure 7,
Figure 8 and
Figure 9.
The study analyzed the photon/electron fluences and their contributions at different virion cluster grades with different energy spectrum peaks due to cross-section considerations and energy spectrum degradation.
Every photon spectrum has been subdivided in three energy bands: low, medium, high in order to highlight the differences in the channel acquisition as a function of the contamination.
As far as the photon spectra as a function of viral contamination in the air sample, the channels:
- -
2 keV, 3 keV, 4 keV, 20 keV, 90 keV are markers as a function of the PO4 group ppm contamination
Concerning electron spectra as a function of viral contamination in the air sample, the channels:
- -
7 keV, 8 keV, 10 keV, 20 keV are markers as a function of the PO4 group ppm contamination
Figure 10 shows the photon/electron fluences as a function of viral contamination coming from the PO
4 group ppm contamination.
4. Discussion
Given the current pandemic as a result of SaRSCoV-2, causing CoVID-19, and the possibility of more viruses crossing species boundaries resulting in yet more viral epidemics, the need for detection of these infectious agents is critical to controlling current and future outbreaks of infectious diseases caused by viruses. Here we have reported that the technology for detecting these infectious agents could exist and could result in quantifiable measurements of ambient viral particles [
2].
Using the physics described in this paper, it could be possible to detect specific viral particles in the environment. This can provide crucial, necessary for the healthcare community to detect viral levels and correlate those to the probability of individuals in a defined area to become infected. Detection of particles in ambient air is crucial in order to lower infection rates in places such as hospitals, retirement homes, and other healthcare facilities. Furthermore, this technology would enable business and government official to make informed decisions about the appropriate measures necessary to take in order to lower infection rates.
As climate change continues and arthropod (vectors) populations migrate, this technology could enable the tracking and identification of the movement of infectious viruses to different geographical regions populated by humans. “According to a report by the U.S. Centers for Disease Control and Prevention. It found vector-borne diseases spread by parasitic insects and arachnids more than tripled in the U.S. over 12 years—from 27,388 cases in 2004 to 96,075 in 2016” [
10].
The photon/electron fluences and spectra can discriminate the amount of virion particles contamination by using its own “particle signature” in terms of photon/electron counts at the detector point combined with the spectrum analysis, as reported. This will enable humans to track the possibility of future viral outbreaks.
The model proposed is based on multiple spherical geometry clusters in grid cells whose spheres have different sizes (as a function of ppm viral contamination). However, those clusters are immersed in a 75% nitrogen composition system which has weighting factor as far as macroscopic cross section contribution.
On
Figure 12,
Figure 13 and
Figure 18 at low contamination 1–100 ppm the photon/electron count profiles show self-shielding and stopping powers behaviors as a function volumetric content as shown in
Figure 20.
As shown in
Figure 18,
Figure 19 and
Figure 20, the electron fluences and spectra associated with the sample of air at different concentrations of the viral particles show a trend in term of electron/cm
2 on different detection energy channels: 7-8-10-20 keV are identifiers as function of PO
4 ppm contamination as much as 120–130 keV. In all the other electron channels superimposed conditions are present, from an energy spectrum point of view, and not enough to discriminate and evaluate their own contribution as function of ppm contamination.
As shown in the cross-section graphs (
Figure 21 and
Figure 22), a 130 keV photon beam allows the identification of phosphorus, as part of the PO
4 group in the cluster configurations.
Analyzing all the significative channels as ‘particle signature’ found in the Fluka simulations it is possible to confirm that: at 2–3 keV channels, nitrogen has a 12 barn cross-section while phosphorus has a cross-section of about 54 barn; therefore, at low energy between 1–10 keV the identification of a virus (that contains phosphorus as a PO4 group) counterbalances a low concentration in ppm as deficit (in proportion) of the first and second fluorescence. Those results are also confirmed by the electron spectrum on the same channels due to photoelectric effect, first and second fluorescence.
By increasing the energy channels to 20 keV, nitrogen has 0.4 barn while phosphorus has a value of 4 barn. As far as 90 keV nitrogen has a value of 0.002 barn while phosphorus 0.03 barn.
To summarize, for primary photon beam energies less than 1 keV the nitrogen total photon cross section has an overall higher probability (in barn) compared to the phosphorus in the same energy range, leading to difficulties in the discrimination process. For energies greater than 1 MeV the steady decreasing of the Compton effect and increasing of pair production make a potential analysis subtle and almost imperceptible. To avoid this problem, it is necessary to proceed with a photon beam between 130 and 150 keV at an experimental level in order to have the ability to discriminate the PO
4 group since the dominant volumetric components would be nitrogen and oxygen. In
Figure 23, a potential new application design and set-up is also shown.
5. Conclusions
This study proposes a new approach to identify low contaminations of viral particles clusters mixed in air using Monte Carlo simulations. The study evaluates the primary/secondary photon/electron energy spectra and fluxes revealed by a potential experimental detector on different energy channels.
Different types of contamination grades can be discriminated using their trends vs. photon/s·cm2—electron/s·cm2 evaluated on different energy intervals as a function of the energy photon beam primary source. Every single contamination is unique in its own spectrum photon/electron signature, and the flux acts as a unique identifier in the detection process so that it can give the ppm amount of virion particles in the air.
Finally, this work has shown that there is still a need for experimental results to confirm the computational analysis to identify RNA in air. Experimental, theoretical, and computational work are all needed for a better understanding of the link between the physical nature of a radiation and its effects for the RNA identification and detection as aerosol in air samples by means of photon/electron interactions.