key: cord-0885752-q5az3y1w authors: Shadloo-Jahromi, Alireza; Bavi, Omid; Hossein Heydari, Mohammad; Kharati-Koopaee, Masoud; Avazzadeh, Zakieh title: Dynamics of Respiratory Droplets Carrying SARS-CoV-2 Virus in Closed Atmosphere date: 2020-10-17 journal: Results Phys DOI: 10.1016/j.rinp.2020.103482 sha: ff72d18b1cd732fbe0a823f6cbd55e73206eafa7 doc_id: 885752 cord_uid: q5az3y1w From the epidemiological point of view, the lifetime of cough and sneeze droplets in the ambient atmosphere plays a significant role in the transmission rate of Coronavirus. The lifetime of indoor respiratory droplets, per se, is a function of droplet size, ambient temperature, and humidity. In the attempt to explore the effective factors of droplet lifetime, sufficient knowledge of atomic-scale interactions and dynamics of the droplet with themselves, as well as the airflow molecules in the room space, is necessary. In this study, the vertical traveling of a wide range (100 nm -10 μm) of representative carrier droplets is studied in three ambient temperatures of 258, 298, and 318 K using all-atom molecular dynamics simulation. Our obtained results confirm that by increasing the room temperature, the suspending time of aerosol (suspended droplets carrying virus particles) increases due to the higher dynamics of air and evaporated water molecules in room space. In fact, by increasing the indoor temperature, the collision rate of aerosol and ambient atmosphere molecules increases significantly. Our result shows this higher rate of collision could have a dual effect on the lifetime of aerosol considering the fact of faster deposition of larger (heavier) droplet due to the gravitational force. On one hand, in higher temperatures, the higher collision can split the droplets to smaller ones with a semi-permanent suspension period. On the other hand, the higher dynamics of ambient molecules can lead to meet and coalesce of smaller cough/sneeze droplets making larger (heavier) droplets with faster sediment times. So, the role of indoor humidity to fuel the probability of coalescence phenomenon and lifetime of droplets becomes more determinant in the warmer spaces. The coronavirus discovered firstly in 1966 and isolated from the sick with common cold [1] . In 2019, the spread of new coronavirus (SARS-CoV-2) increased the attention worldwide. On 11 March 2020, the World Health Organization (WHO) notified this virus as a pandemic based on its high spreading rate. SARS-CoV-2 is now known as a new microorganism from the subfamily orthocoronavirinae that causes olfaction failure, respiratory and gastrointestinal infection in hosts. As yet, no effective vaccine has been discovered, thus precautionary measures such as social distancing and home isolation is the best way to minimize the spread of Coronavirus. During the last decade, some epidemiological model for the propagation of viruses has been developed for a verity of virus classes from computer viruses [2, 3] to Coronavirus (2019-nCov) [4] using different mathematical approaches. In case of SARS-CoV-2 virus, specialists believe that this virus spreads mainly from person to person by (i): airborne transmission, (ii): droplets or aerosols transmission, and (iii): surface transmission. Among all of these mechanisms, airborne and droplets transmission is more perilous due to the deprivation of the applicable way to disinfection of the air molecules. When infected people sneeze, cough or, even talk, the droplets or small particles carry the coronavirus into the air and in these modes, the ambient becomes infected by SARS-CoV-2. In 1995 Wells [5] showed that the droplet nuclei (dried-out residual of droplets containing infectious pathogens) are the conveyance for airborne respiratory disease transmission. Therefore, people who are in such an ambient could breathe the droplets into their lungs. This is the major issue faced by governments and the public around the world. Therefore, one of the main concerns related to COVID-19 is the longevity of droplets in the environment. On the other hand, the longevity of droplets in coastal cities or in cold regions could be different from the room condition. Therefore, analyzing the behavior of droplets with different sizes in such conditions could release a new perspective in order to confront COVID-19. Much literature accomplished to study the horizontal spread of droplets. Zhu et al [6] by employing the particle image velocimetry (PIV) method showed that during each coughing, the droplet with a size of 30 μm affect more than 2 m away from the source. Xie et al [7] by using numerical simulation indicated that expelled large droplets (> 60 μm) can travel more than 6 m. Parienta et al [8] employed the mathematical modeling and released that the droplets with a diameter of 16 μm and a coughing velocity of 11.7 m/s, could pass a distance of 7 m. In the mathematical modeling of Bourouiba et al [9] , it has been calculated that droplets with 30 μm diameter, during coughing within a turbulent gas could travel a horizontal distance in range of 2.5 m. Wei and Li [10] with using mathematical modeling released that with mouth opening of 2 cm, large droplets in range of 96 μm could pass a distance up to 1.4 m. Lee et all [11] used optical particle spectrometer to calculate the cough particle concentration of 10 patients. They indicated that transition could spread more than 3 m from the person. In the mathematical study of Chen [12] , a one-dimensional droplet evaporation model was invoked to estimate the droplet lifetime from evaporation in air. He showed that the ambient temperature and ambient relative humidity have the potential effects on virus-carrying drops. Asadi et al. [13] ; Jones and Brosseau [14] indicated that the competing influences of gravity, evaporation, and inertia specify the commendation of the expelled droplets. Although previous studies have released efficacious data in the context of coronavirus and droplets, we believe that studying the longevity and vertical traveling of droplets could be served in precautionary measures which have not been studied clearly hitherto. In this study, molecular dynamics (MD) simulation has been employed to investigate the behavior of water droplets in a room environment. The main finding of this study is to predict the vertical traveling and lifetime of aerosol at different thermodynamics conditions. Different sizes of water droplets were considered to assess the effect of air temperature on the behavior of hovering droplets. We anticipate our study could lead to a better understanding of expelled droplets as a significant factor in the transmission of SARS-CoV-2. In this study, an infinite silicon surface is considered as the floor and some water droplets with wide-range sizes are placed on the above of the surface. The droplets are located in the different coordinates of the surface. Figure 1 shows the schematic configuration of geometry simulated in this study. In figure 1 , the silicon surface is indicated and the droplets which solved in the bulk air box, are put on tho above the silicon surface. The dimension of the simulation box is 43.4 × 43.4 × 120 nm. In this study, the large scale atomic/molecular massively parallel simulator (LAMMPS) [15] was employed to accomplish the atomistic studies. In order to visualize the atomistic output data, we employed OVITO [16] . In MD simulation the classical equations of motion were solved, which for a simple atomic simulation could be written as: Where , and are the mass acceleration of molecule i respectively. In order to determine the , many literature has been reported the governing equation of intermolecular potential. Among all of the potential, Lennard-Jones 6-12 is used extensively as follow: ( ) Where is the distance between two atoms of the system. and are the molecular diameter and well depth of potential function, respectively which known as the Lennard-Jones parameters. The Verlet-velocity algorithm is carried out to integrate the equation of motion. To describe the intermolecular potential for all atoms in this study, such as water molecules, air molecules, and silicon atoms the Lennard-Jones (12-6) potential has been employed. Also, a Lennard-Jones potential with modified parameters of and is implemented for interactions between various Table 1 . In the initial configuration of our system, we consider 22 boxes of water molecules and also air molecules including and molecules with a density of . The air and equilibrate the air molecules at the desired temperature. In the next step, we allow 22 boxes of water molecules to reach the equilibrate state. In this part of the simulation, some small droplets stick together to form larger droplets (see Figure 2 ). As we discussed later, the formation of the large droplets (coalescence phenomenon) in an indoor environment could lead to the faster sediment of droplets. After the equilibration state, the gravitational force is applied to the representative viruscarrying droplets of various sizes. Figure 3 shows the system in different conditions including droplets, air molecules, and evaporated water molecules (room humidity molecules). Figure 3A shows a snapshot of the equilibration state considered as the initial condition ( ) and figures = 0 3B-D indicate the system topology after , for , and respectively. In 780 = 258 298 318 K order to start the nonequilibrium simulation, we exerted the external force of 143E -6 kcal/mol to each water molecule. This force is equivalent to the weight of a drop with a diameter of 10 . µ For a better understanding of temperature effect on vertical traveling of droplets, the biggest droplet and, one of the small droplets are considered as the basis of our simulation. These two droplets are marked with red arrows. As shown in figure 3 , at the big droplet is located at a = 0 higher height, while it reaches to ground sooner than the small one. This issue is more significant in higher values of temperature. For two droplets reach the ground almost = 258 simultaneously while by increasing the air temperature to and the big droplet is more 298 318 interesting to reach the ground. On the other hand, as figure 2 shows, at the low temperature of both the droplets like to hit the surface sooner in comparison to higher temperatures of = 258 and . 298 318 Figure 3 . The configuration of simulated cases in different temperature. In figure 4 , the time-dependent vertical movement of two targeted droplet are plotted for three different temperature of , , and . The results show that throughout the traveling, = 258 298 318 the droplets in higher indoor temperature tend to stay suspended more than those in the lower ambient temperature and this issue becomes more significant by decreasing the droplet size. As shown in Figure 4A , the difference in the traveling distance for small droplets is more considerable than the bigger one (see Figure 4B ). The observed effect of ambient temperature on the vertical traveling distance of droplets reveals that the deposition speed of virus-carrying droplets on the floor increases in higher indoor temperature, especially for smaller droplets (<10 μm). This matter implies that the lifetime of aerosol in the ambient atmosphere of closed space is longer for warmer places. Although Jahangiri et al [21] and Yao et al [22] show that the ambient temperature is no a with the recent study [12] . It has been reported that the lifetime of the droplet does not decrease by increasing the temperature when the ambient relative humidity is set below 37%. The results indicate beyond a critical ambient relative humidity, increasing the ambient temperature can increase the droplet lifetime exponentially [12] . The role of ambient relative humidity can be incorporated to increasing the rate of collision of air/water molecules with droplets as well as the coalescence phenomenon in which, the smaller droplets stick together and form a larger (heavier) droplet with a higher sediment speed. Due to the observation of permanently suspended droplets with 1 µm diameter in our simulation, more studies are invoked to investigate the disinfecting efficiency of large indoor space with fogging machines considering the higher probability of coalescence of droplets and faster deposition rate of larger aerosols. In this study, molecular dynamics simulation has been carried out to analyze the dynamics and lifetime of the virus-carrying water droplets. Our simulation performed for various sizes of droplets and three different ambient temperatures have been considered to assess the vertical traveling of droplets. The major obtained results from the following study are: 1. Although the droplets with a size of less than 10 µm could suspend in an indoor environment for a long time, these small droplets could cohere to each other and form the bigger droplets. This phenomenon could lead to the deposition of expelled small droplets. 2. Droplets could be suspended more time in a warm ambient atmosphere compared to the cold one. Therefore the transmission rate by inhaling the infected aerosols in warm indoor space can be higher than the cold one. 3. Rising the dynamics of ambient air molecules with increasing the temperature or air conditioner, has a more significant effect on small droplets than the bigger one. 4. Increasing the ambient relative humidity can facilitate the coalescence of the smaller droplets and therefore, reduce the suspension time of these droplets. Effects of a 'New' Human Respiratory Virus in Volunteers Modeling the spread of computer virus via Caputo fractional derivative and the beta-derivative A fractional epidemiological model for computer viruses pertaining to a new fractional derivative Modeling the dynamics of novel coronavirus (2019-nCov) with fractional derivative Airborne contagion and air hygiene Jeong-Hoon Yang, Study on transport characteristics of saliva droplets produced by coughing in a calm indoor environment. Building and Environment How far droplets can move in indoor environments-revisiting the Wells evaporation-falling curve. Indoor Air Theoretical analysis of the motion and evaporation of exhaled respiratory droplets of mixed composition Violent expiratory events: on coughing and sneezing Human cough as a two-stage jet and its role in particle transport Quantity, size distribution, and characteristics of cough-generated aerosol produced by patients with an upper respiratory tract infection Effects of ambient temperature and humidity on droplet lifetime -A perspective of exhalation sneeze droplets with COVID-19 virus transmission The coronavirus pandemic and aerosols: Does COVID-19 transmit via expiratory particles? Aerosol Science and Technology Aerosol transmission of infectious disease Fast Parallel Algorithms for Short-Range Molecular Dynamics Visualization and analysis of atomistic simulation data with OVITO-the Open Visualization Tool. Modelling and Simulation in Materials Science and Engineering Intermolecular Forces: Their Origin and Determination The Missing Term In Effectlve Pair Potentlals Friction factor calculation in nanochannels comprising different wall hydrophobicities and superhydrophobic structures : Molecular dynamic simulations. International Communications in Heat and Mass Transfer Vapor-liquid equilibria of mixtures containing alkanes, carbon dioxide, and nitrogen The sensitivity and specificity analyses of ambient temperature and population size on the transmission rate of the novel coronavirus (COVID-19) in different provinces of Iran No association of COVID-19 transmission with temperature or UV radiation in Chinese cities