key: cord-325471-fc152wm3
authors: Zhang, Xiaqi; Wang, Fang; Weng, Ming; Cao, Meng
title: Calculations of energy deposition and ionization in the 2019 novel coronavirus by electron beam irradiation
date: 2020-09-04
journal: Radiat Phys Chem Oxf Engl 1993
DOI: 10.1016/j.radphyschem.2020.109169
sha: 
doc_id: 325471
cord_uid: fc152wm3

Using Monte Carlo methods, this study investigates energy deposition of energetic electrons and ionization in the 2019 novel coronavirus by electron irradiation, which are important characteristic quantities related with biological damage formation. The inelastic scattering of low-energy electrons (<10 keV) was calculated by dielectric theory. The optical energy-loss functions of viral proteins and RNA were derived from an empirical method in the energy-loss range <40 eV and the calculation of optical parameters of the biomolecules. The densities and distributions of energy deposition and ionization were calculated from the stopping power and inelastic cross-sections in the electron-cascade simulation. Electrons with primary energies of approximately 1–3 keV produced significant energy deposition and ionization in the target coronavirus. More energetic electrons were less effective due to the larger electron range and fewer scattering events in the coronavirus.

The 2019 novel coronavirus (2019-noCoV) was responsible for the pandemic of 2020 (Sohrabi et al., 2020; WHO, 2020) . To avoid exposure to COVID-19, effective methods that kill or inactivate the virus on contaminated media or in specific environments are demanded. Electron beam irradiation is widely used for killing microorganisms, bacteria, and viruses (Berovic et al., 2002; Benson, 2002; Feng et al., 2020) , because its penetration depth and dose can be controlled to comply with environmental protection requirements. Moreover, under controlled irradiation, the virus can be partially or totally inactivated, providing useful information on its properties (Lidzey et al., 1995; Smolko and Lombardo, 2005) .

Electron irradiation damage on biological macromolecules is mainly characterized by inactivity and the irreversibly chain broken (Tan and Liu, 2014) . The severity of these effects is highly correlated with the energy deposition and ionization in the molecules, which occur by inelastic interactions between the electrons and molecules (Boudaïffa et al., 2000; Kempner, 2001; Kempner, 2011; Tan et al., 2018; Zhang and Tan, 2010) . To properly understand electron irradiation effects on the novel coronavirus, the energy deposition and ionization must be calculated, because experimental evaluation with currently available methods is very difficult.

High-energy electrons (above tens of keV) are usually applied for antiseptic purposes. When primary electrons interact with a target molecule by inelastic scattering, they lose energy and excite the shell electrons as secondary electrons (SEs).

Low-energy electrons, including decelerated primary electrons and SEs with energies below 50 eV, significantly contribute to the large number of SEs during the electron-cascade process, exerting a considerable ionization effect (Emfietzoglou et al., 2012; Nikjoo et al., 1997; Zhang and Tan, 2010) . Therefore, the interaction of low-energy electrons with the main proteins and RNA macromolecules of the novel coronavirus was the key issue of our calculation. The stopping power was suitable for describing the energy-deposition probability in the target, and the ionization events were determined as the inelastic cross-sections of collisions with the shell electrons of J o u r n a l P r e -p r o o f an atom. The stopping power and ionization events of low-energy electrons are usually calculated by dielectric theory, which requires the energy-loss function (ELF). The ELF can be derived from the optical ELF (OELF) through a proper dispersion relation.

However, as important input data for the calculation of the OELF, the optical parameters of macromolecules of proteins, RNA are not available. Tan et al. developed an empirical method that estimates the OELF of organic biomaterials at energy losses below 40 eV (Tan et al., 2004a (Tan et al., , 2004b . The same approach is expected to be applicable to the novel coronavirus.

In this work, we calculated the distributions of energy deposition and ionization in the novel coronavirus irradiated with electrons having primary energies below 10 keV. 

For the calculations, we used a simplified physical model of the novel coronavirus with a multilayer spherical shell configuration. The model is based on published information of the novel coronavirus and similar known coronaviruses (CDC, 2020; Lu, 2020) . The novel coronavirus is comparatively large (with an approximate diameter of 120 nm), and its main proteins are spiked glycoproteins (S) on the outer layer, envelope proteins (E) in the next layer, a water layer, the nucleocapsid (N), and RNA in the centre (see Fig. 1 ). The diameters of each shell are marked in Fig. 1 . The components of each shell were assumed as the typical amino acid sequences or nucleobase sequences of the main proteins within the shell, extracted from the sequence resources published in the National Center for Biotechnology Information (NCBI, 2020). The statistical results J o u r n a l P r e -p r o o f yielded the average molecular formulas of the three types of proteins as follows: 

The energetic electrons interact with the novel coronavirus through elastic and inelastic scattering. Elastic scattering describes the collision between electrons and the nuclei of the target atoms, which changes the direction of the electron trajectories. The elastic-scattering calculation of electrons with energy <10 keV was usually based on the Mott cross-section. In this calculation, the atomic differential elastic-scattering cross-sections were taken from the NIST electron elastic-scattering cross-section database (Jablonski et al., 2010) . Following the additivity rule (Chiari et al., 2013; Jiang et al., 1997) , the differential cross-section of each kind of molecule in the novel coronavirus was calculated by summing the differential cross-sections σ el of the individual atoms in the target molecule as follows:

where n i is the atomic number of the i-th atom in the molecular formula, and Ω is the solid angle determined by the scattering angle θ.

The inelastic cross-sections of low-energy electrons are usually calculated by dielectric theory. This calculation requires the ELF, given by the negative of the reciprocal of the imaginary part of the dielectric function (Im − 1 , ∆ ⁄ ). Penn (1987) and Ashley (1988) obtained the ELF by extending the dielectric function in the optical limit, where the momentum transfer is = 0 . By fitting the OELF Im − 1 0, ∆ ⁄ to the experimental optical data and introducing a suitable dispersion relation, the differential inelastic cross-section σ inel of the incident electrons colliding with the shell electrons can be calculated as (Ashley, 1988) ,△

with

In these expressions, 5 ( is the Bohr radius, 6 is the number of molecules per unit volume in the target, is the electron energy, and Δ is the energy transfer.

To calculate the energy loss along the electron trajectories, we require the stopping power of the electrons, which represents the average energy loss of an electron in the unit distance. The stopping power is given by 

However, the experimental OELF data of most complex compounds are lacking in the literature. For the constituent molecules in the novel coronavirus, which certainly lack available optical data, the OELFs at energy losses below 40 eV were estimated by an empirical approach for bioorganic compounds developed by Tan et al. (2004a; 2004b) .

The OELF was expressed by the following single Drude function:

where the parameters G and H are linearly related to the average atomic number I ̅ of the molecule, and 5 can be determined by inputting the mass density of the target, the molecular weight and the total number of electrons per molecule to the f-sum rule. For J o u r n a l P r e -p r o o f energy losses larger than 50 eV, the atomic photoabsorption data of Henke et al. (1993) are available, so the OELF can be calculated from the refraction index and the extinction coefficient (Ding and Shimizu, 1996) . For energy losses between 40 eV and 50 eV, the OELF was calculated by spline interpolation.

If the electron energy decreases by 2-3 times the bandgap energy during the electron cascade, the main energy-loss mechanism should be electron-phonon interactions, and the cross-section can be treated by Fröhlich's theory (Fröhlich, 1954) . Low-energy electrons are also affected by the polarization field induced by their own motions through the target. This interaction is known as the polaronic effect. In this study, the cross-section of the polaronic effects was calculated by a semi-empirical approach proposed by Ganachaud and Mokrani (1995) .

As shown in Fig. 1 , uniformly distributed electrons with primary energies below 10 keV were irradiated over the 120-nm diameter of the novel coronavirus. The number of primary electrons was 10 > . The densities of protein, RNA and water were set to 1.3, 1.35, and 1.0 g cm 4 ⁄ , respectively.

The electron trajectories were traced by two Monte Carlo strategies. The energy deposition was determined by a simple and fast Monte Carlo method based on the continuous slowing down approximation (CSDA). The CSDA assumes that electrons continuously lose energy along the electron trajectories, and change directions by elastic scattering. In the ionization-characterizing calculation, each energy loss and direction change along the electron trajectories was calculated by the energy straggling strategy. If the amount of energy lost in a single inelastic scattering event exceeded the binding energy of the atomic shell, the site became ionized, producing an inner secondary electron that was traced by the same procedure as the primary electrons (Caleman et al., 2009 ). This simulated electron cascade considered the events of elastic scattering, electron-electron inelastic scatterings, and electron-phonon and electronpolaron interactions (Dapor, 2014) .

To describe the inelastic interactions between the electrons and molecules, the OEFLs of the constituent molecules in the novel coronavirus were calculated by using Eq. (6) and atomic photoabsorption data (Henke et al., 1993) , as mentioned above. The OEFL of water was adopted from Heller et al. (1974) . As shown in Fig. 2 , the OEFL curves of the five components peaked at energies between 20 and 25 eV. The envelope yielded the highest OEFL peak, and RNA has the maximum peak position. From the OEFLs, the probability of inelastic scattering events and the energy-loss distributions were determined. and RNA. The OEFL of water was adopted from Heller et al. (1974) . 

The Table 1 . Like the deposition energy, the ionization was maximized in the envelope.

Both the energy deposition and ionization densities appeared over a larger range of primary energies (~400 eV-3 keV) in the envelope than in the other viral components (see Figs. 3 and 4) . 

J o u r n a l P r e -p r o o f Figure 5 shows the distribution characteristics of the deposited energy and ionization densities through a vertical section of the novel coronavirus. To simplify the presentation, the sparse structure of the spike glycoproteins was neglected, and the deposited energy and ionization were assumed to be uniformly distributed through the shell volume. The primary energy was varied as 500 eV, 1.5 keV, and 5 keV. At 500 eV, the energy deposition and ionization were concentrated in the upper shell (the envelope) because the penetration range of the primary electrons was small.

Approximately 60% of the primary electrons lost all their energy and were deposited in the target. At a primary energy of 1.5 keV, the energy deposition and ionization covered Most of the electrons with energies below hundreds of eV will finally loss their energy and be deposited inside the target. As the primary energy decreases, the total amounts and densities of the energy deposition decrease accordingly. The estimated deposition depth of primary electrons with energies above 2 keV in the thick materials of the envelope, nucleocapsid, and RNA is 180 nm (Fitting, 1974) , larger than the coronavirus diameter. As fast electrons have lower probability of inelastic interactions with molecules than slow electrons, their energy deposition and ionization in the target will decrease, as shown in Figs. 3 and 4 . Therefore, the energy deposition and ionization densities are maximized over an optimal range of primary energies (1-3 J o u r n a l P r e -p r o o f keV). The presented results demonstrate that the characteristics of energy deposition and ionization depend on the energy loss and electron range in the coronavirus, which are mainly determined by the primary energy of the electrons. Furthermore, the electron energy deposition is correlated with the temperature rise imparted by the electron irradiation. High temperature is expected to be an important deactivating factor of the novel coronavirus.

Although electrons with lower energies (1-3 keV) effectively deposit their energies and ionize the viral macromolecules, irradiating the novel coronavirus by low-energy electrons is difficult to achieve in practice. Theoretical calculations may provide fundamental physical information on the electron irradiation of viruses.

In this paper, we performed calculations of energy deposition and ionization on the novel 2019 coronavirus by electron irradiation, using Monte Carlo methods. The 

Interaction of low-energy electrons with condensed matter: stopping powers and inelastic mean free paths from optical data

Use of radiation in biomaterials science

The kinetics of radiation damage to the protein luciferase and recovery of enzyme activity after irradiation

Resonant formation of DNA strand breaks by low-energy (3 to 20 eV) electrons

Radiation damage in biological material: Electronic properties and electron impact ionization in urea

Total, elastic, and inelastic cross sections for positron and electron collisions with tetrahydrofuran

Transport of energetic electrons in solids: Computer simulation with applications to materials analysis and characterization

A Monte Carlo modeling of electron interaction with solids including cascade secondary electron production

Inelastic scattering of low-energy electrons in liquid water computed from optical-data models of the Bethe surface

Electron beam irradiation on novel coronavirus (COVID-19): A Monte-Carlo simulation

Transmission, energy distribution, and SE excitation of fast electrons in thin solid films

Electrons in lattice fields

Theoretical study of the secondary electron emission of insulating targets

Collective oscillation in liquid water

X-Ray Interactions: Photoabsorption, Scattering, Transmission and Reflection at E = 50-30,000 eV, Z= 1-92. At. Data Nucl

NIST Electron Elastic-Scattering Cross-Section Database, National Institute of Standards and Technology

Geometric shielding effects of electron scattering from polyatomic molecules

Effects of high-energy electrons and gamma rays directly on protein molecules

Direct Effects effects of ionizing radiation on macromolecules

A critical analysis of the use of radiation inactivation to measure the mass of protein

Genomic characterisation and epidemiology of 2019 novel coronavirus: implications for virus origins and receptor binding

Computational modelling of low-energy electron-induced DNA damage by early physical and chemical events

Electron mean-free-path calculations using a model dielectric function

Virus inactivation studies using ion beams, electron and gamma irradiation

World Health Organization declares global emergency: A review of the 2019 novel coronavirus (COVID-19)

Simple and universal model for electron-impact ionization of complex biomolecules

Monte Carlo calculations of energy deposition distributions of electrons below 20 keV in protein

Cross sections of electron inelastic interactions in DNA

Electron stopping power and mean free path in organic compounds over the energy range of 20-10,000 eV

WHO World Health Organization, 2020. WHO statement on cases of COVID-19

A new calculation on spectrum of direct DNA damage induced by low-energy electrons