# Stokes-Mueller polarization-based analysis of model SARS-CoV-2 virions

K. U. Spandana, Bhagesh Basavraj Hunakunti, Aymeric Le Gratiet, Ankur

Gogoi, Nirmal Mazumder

► To cite this version:

K. U. Spandana, Bhagesh Basavraj Hunakunti, Aymeric Le Gratiet, Ankur Gogoi, Nirmal Mazumder. Stokes-Mueller polarization-based analysis of model SARS-CoV-2 virions. Lasers in Medical Science, 2023, 38 (1), pp.35. 10.1007/s10103-022-03680-3. hal-04229208

HAL Id: hal-04229208

https://hal.science/hal-04229208v1

Submitted on 5 Oct 2023

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons CC BY 4.0 - Attribution - International License
# Stokes–Mueller polarization-based analysis of model SARS-CoV-2 virions

Spandana K U<fnref n="1" /> · Bhagesh Basavraj Hunakunti<fnref n="1" /> · Aymeric Le Gratiet<fnref n="2" /> · Ankur Gogoi<fnref n="3" /> · Nirmal Mazumder<fnref n="1" />

Received: 6 May 2022 / Accepted: 12 October 2022
© The Author(s) 2022

## Abstract

Understanding the virology of the coronavirus at the structural level has gained utmost importance to overcome the constant and long-term health complications induced by them. In this work, the light scattering properties of SARS-CoV-2 of size 140 nm were simulated by using discrete dipole approximation (DDA) for two incident wavelengths 200 nm and 350 nm, respectively. Three different 3-dimensional (3D) models of SARS-CoV-2 corresponding to 15, 20, and 40 numbers of spike proteins on the viral capsid surface were constructed as target geometries for the DDA calculations. These models were assessed by employing Stokes–Mueller polarimetry to obtain individual polarization properties such as degree of polarization (DOP), degree of linear polarization (DOLP), and degree of circular polarization (DOCP). Irrespective of its spike numbers, all the coronavirus models were found to display higher DOP and DOCP values and negligibly small DOLP values for circularly polarized incident light, indicating the presence of chiral structures. On the other hand, the lack of understanding about the dependence of the Mueller matrix on its microstructural properties was overcome by transforming 16 Mueller elements into sub-matrices with specific structural and physical properties using Lu–Chipman-based Mueller matrix polar decomposition method. The obtained properties such as retardance, diattenuation, and depolarization were used for investigating the composition and microstructural information. The approach presented in this work has the potential to understand the virology of the coronavirus at the structural level and, therefore, will be beneficial in developing effective detection strategies by exploiting their characteristic electromagnetic scattering signatures.

**Keywords** SARS-CoV-2 · Discrete dipole approximation · Mueller matrix · Stokes vector · Polarimetry · Lu-Chipman decomposition

## Introduction

The ongoing corona pandemic crisis is caused by an infection of severe acute respiratory syndrome-related coronavirus-2 (SARS-CoV-2), leading to inflammatory conditions in human lungs. An evolved beta-corona virus similar to human COVID-19 causing SARS-CoV-2 in *Rhinolophus sinicus* (horseshoe bats) [1] was discovered in China soon after the first transmission emergence of SARS-CoV from animals to humans [2]. Epidemiologically, pneumonia incidents in human hosts are a result of varied viral strains, which include adenovirus, influenza virus, Middle East respiratory syndrome virus (MERS-V), parainfluenza virus, respiratory syncytial virus (RSV), SARS-CoV, and enteric

<footnote>☒ Aymeric Le Gratiet
aymericlegratiet@gmail.com</footnote>

<footnote>☒ Ankur Gogoi
ankurgogoi@gmail.com</footnote>

<footnote>☒ Nirmal Mazumder
nirmaluva@gmail.com</footnote>

<footnote n="1">¹ Department of Biophysics, Manipal School of Life Sciences, Manipal Academy of Higher Education, Manipal, Karnataka, India 576104</footnote>

<footnote n="2">² Université de Rennes, CNRS, Institut FOTON - UMR 6082, F-22305 Lannion, France</footnote>

<footnote n="3">³ Department of Physics, Jagannath Barooah College, Jorhat 785001, Assam, India</footnote>
enveloped CoV [3–6]. On 13th January 2020, The World Health Organization (WHO) coined the official terminology 2019-Novel Coronavirus (2019-nCoV), and simultaneously, the disease caused by 2019-nCoV was termed as Coronavirus Disease-2019 (COVID-19) on 11th February 2020 [7]. Due to its similarity with the novel coronavirus and SARS-CoV, the International Committee on Taxonomy of Viruses (ICTV) declared the official nomenclature of the virus as SARS-CoV-2 [8]. SARS-CoV-2 has resulted in unforeseen public health and economic threats worldwide [9]. Globally, as of 25th January 2022, there have been 352,796,704 confirmed cases of COVID-19, including 5,600,434 deaths, reported to WHO [10]. The effect of COVID-19 on humankind has been devastating that the global stock markets experienced their worst crash since 1987, and in the first three months of 2020, the G20 economies fell 3.4% year-on-year [11].

Remarkably, with a size of 60–140 nm and constituting a (+) ssRNA genome of 29,891 bp in size (functions as mRNA directly which encodes 9860 amino acids), SARS-CoV-2 can divide faster and infect rapidly, causing a global health-threatening emergency [12]. Structurally, SARS-CoV-2 constitutes a spike (S) glycoprotein, dimeric HE enzyme, a membrane matrix glycoprotein (M), and RNA [13]. The main point of contact with the host cell is mediated by the S glycoprotein directly and indirectly in the infection cycle. A receptor-binding domain (RBD) is borne on the S glycoprotein of all coronaviruses. The binding of S glycoprotein RBD to its host receptor causes cleavage of this glycoprotein by a host furin-like protease, which releases the spike fusion peptides, thus facilitating entry into the host cell [14]. S glycoprotein has been found to have numerous binding and neutralization epitopes, making it an essential target for vaccine design [15–17]. Notably, even before WHO declared a worldwide pandemic, researchers have been on their toes to develop antiviral drugs against SARS-CoV-2. Nevertheless, coronavirus has been successful in expanding their host ranges, including humans, to their recombination, mutator alleles, and mutational/evolutionary mechanisms [18]. Thus, it has become crucial to understand the virology of the coronavirus at a structural level due to constant as well as long-term health complications from these zoonotic viruses [19].

In this context, polarized light scattering characteristics of particulate matter, including biological entities, carry a lot of information about their physical

and optical properties [20]. Notably, the 4 × 4 Mueller matrix connects the Stokes vectors of the light incident on and scattered by the system of particles under study. The Mueller matrix, which depends on several parameters related to the scattering particle, such as size, shape, and refractive index, can be utilized to extract some of the physical and optical characteristics of the scatterer. In particular, the sample assessment by using Mueller matrix polar decomposition method distinguishes multiply scattered light, revealing individual polarization properties such as retardance, diattenuation, and depolarization [20], which can further be used to extract their compositional and microstructural information. Similarly, approaches based on Stokes–Mueller formalism could be extremely beneficial in understanding the virology of the coronavirus at a structural level. Several studies have been carried out in the past to characterize coronavirus and other virus particles by using light scattering techniques [21–25]. Recently, Petrov calculated the scattering properties of coronavirus particles by using spherical model particles with various spike protein numbers. This work was primarily focused to study the effect of spike proteins on the scattered intensity and degree of linear polarization [26]. On the other hand, the study by Ashraf et al. [27] emphasized the information embedded in the circular polarization state of the scattered light to detect the presence of genome in the virus particle. The group showed that the circular intensity differential scattering is sensitive to the number of turns, handedness, and diameter of the RNA genome. These results greatly improved our understanding of interactions between coronavirus with electromagnetic radiations, whilst using only one or two of the 16 Mueller matrix (MM) elements. Remarkably, it is possible to extract additional optical parameters characteristic to the scattering particle by utilizing other MM elements so that the scatterer of interest (coronavirus in this case) can be properly discriminated from other scattering particles. In this work, Stokes vectors and various polarization parameters of the light scattered from various coronavirus models with different spike numbers are calculated at 200 nm and 350 nm incident wavelengths. Further, Lu–Chipman-based Mueller matrix polar decomposition method is employed in this work to investigate the interaction of coronavirus models with different spike numbers [28] and understand a complete set of physical parameters resulting from the decomposition.

<figure><img src="image_2.png" /></figure>
# Materials and methods

## Stokes–Mueller formalism

Polarization is a fundamental property of electro-magnetic radiation (light), and the interactions of polarized light with biological samples can reveal structural information associated with its pathological condition. Even the slightest variation in structural alignment can induce a significant change in polarization property, which can play a crucial role in the early detection of abnormal morphology [29]. In this regard, Sir George Gabriel represented polarization behavior in terms of observables and described the complete state of polarization of light in terms of four measurable quantities, known as the Stokes parameters. Total optical field intensity is described by the first parameter, and the remaining parameters describe the polarization state. The Stokes parameters are a logical consequence of wave theory and can be arranged in a 4 × 1 column matrix as

$$S = \begin{bmatrix} S_0 \\ S_1 \\ S_2 \\ S_3 \end{bmatrix} = \begin{bmatrix} E_{ox}^2 + E_{oy}^2 \\ E_{ox}^2 - E_{oy}^2 \\ 2E_{ox}E_{oy}\cos\delta \\ 2E_{ox}E_{oy}\sin\delta \end{bmatrix} = \begin{bmatrix} I_0 + I_{90} \\ I_0 - I_{90} \\ I_{45} - I_{135} \\ I_R - I_L \end{bmatrix} \quad (1)$$

where $S_0$ represents the total light intensity, $S_1$ is the difference between 0° and 90° polarization intensities, $S_2$ is the difference between +45° and -45° (135°) polarization intensities, and $S_3$ is a difference between the left and right circular polarization intensities [30]. Several polarization parameters such as degree of polarization (DOP), degree of linear polarization (DOLP), and degree of circular polarization (DOCP) can be constructed from the Stokes parameters as follows,

$$DOP = \frac{\sqrt{S_1^2 + S_2^2 + S_3^2}}{S_0} \quad (2)$$

$$DOLP = \frac{\sqrt{S_1^2 + S_2^2}}{S_0} \quad (3)$$

$$DOCP = \frac{|S_3|}{S_0} \quad (4)$$

In general, fully polarized light gets depolarized when it transmits through a scattering environment [31, 32]. While DOP of a fully polarized light is 1 and is 0 for unpolarized light, DOP for a partially

polarized light is between 0 and 1. Thus, the DOP of the scattered light signifies the amount of incident polarization that persists in the scattered light. Similarly, the DOLP quantitatively describes the amount of linear polarization present in the light beam and essentially illustrates the anisotropy of the distribution of molecular alignment/orientation in the scattering volume. DOCP, on the other hand, is a measure of how effectively chiral molecules interact with the circularly polarized light. The values of both DOLP and DOCP ranges from 0 to 1 [33, 34]. Although the Stokes parameters describe the polarization properties of light, Stokes–Mueller formalism is suitable for detailed analysis of the sample. The transformation of polarization property of light due to its interaction with an optical system can be described by Mueller matrix and characterized by Mueller formalism. When Stokes vector, $S_i$ of incoming light propagates through a scattering media, the output Stokes vector, $S_i'$ can be obtained as a linear combination of the four Stokes parameters of the incident beam in terms of matrix form,

$$\begin{bmatrix} S_0' \\ S_1' \\ S_2' \\ S_3' \end{bmatrix} = \begin{bmatrix} m_{00} m_{01} m_{02} m_{03} \\ m_{10} m_{11} m_{12} m_{13} \\ m_{20} m_{21} m_{22} m_{23} \\ m_{30} m_{31} m_{32} m_{33} \end{bmatrix} \begin{bmatrix} S_0 \\ S_1 \\ S_2 \\ S_3 \end{bmatrix} \quad (5)$$

or

$$S' = M.S \quad (6)$$

where $M = \begin{bmatrix} m_{00} & m_{01} & m_{02} & m_{03} \\ m_{10} & m_{11} & m_{12} & m_{13} \\ m_{20} & m_{21} & m_{22} & m_{23} \\ m_{30} & m_{31} & m_{32} & m_{33} \end{bmatrix}$ is the Mueller matrix.

where $m_{ij}$ are experimentally measurable quantities, $i, j = 0$ to 3.

The Mueller matrix provides a comprehensive explanation of the sample’s optical and structural information related to polarization, thus playing an important role in biomedical research and sample characterization. However, there exists a lack of explicit association of Mueller matrix with microstructural properties. Mueller matrix elements reflect lumped effects due to the simultaneous occurrence of several polarization effects, hence hindering the unique interpretation [35–38]. The Mueller matrix of the analyzed sample encloses its combined polarizing properties. However, the 16 Mueller elements can be transformed into sub-parameters with specific structural and physical properties using Lu–Chipman-based Mueller matrix polar decomposition (MMPD) method. Notably, assessment of biological sample by means of
polar decomposition approach, where decomposition of Mueller matrix into three basis matrices can be used as an effective tool to distinguish multiply scattered light and gain individual polarization properties [39]. Polarization properties such as retardance, diattenuation, and depolarization can be used for investigating composition and microstructural information and could become beneficial in understanding the virology of the corona virus at a structural level.

In this work, Mueller matrix representing SARS-CoV-2 was decomposed into a set of three basis matrices representing diattenuator, retarder, and depolarizer, respectively. These sub-matrices were converted into individual parameters associated to diattenuation, retardance, and depolarization properties. Mueller matrix “M” is decomposed as

$$M = M_{\Delta}M_{R}M_{D}$$

where depolarization matrix, $M_{\Delta}$, describes depolarization effect of the medium, a retardance matrix, $M_{R}$, accounts for the effect of optical activity and linear birefringence, and a diattenuation matrix, $M_{D}$, involves the effect of circular and linear dichroism. The basis matrices were examined further to extract individual polarization properties of medium such as, optical rotation ($\Psi$), linear retardance ($\delta$), and its orientation angle ($\theta$), depolarization coefficient ($\Delta$), and diattenuation ($d$).

The diattenuation matrix, $M_{D}$, is defined using a relation

$$M_D = \begin{bmatrix} 1 & \vec{d}^T \\ \vec{d} & m_D \end{bmatrix}$$

where $m_{D}$ is a 3 × 3 submatrix and the diattenuation vector $\vec{d}$ is defined as

$$\vec{d} = \{1/M(0,0)\} \times [M(0,1)M(0,2)M(0,3)]^T$$

The magnitude of diattenuation was determined using a relation

$$d = \{1/M(0,0)\} \times [(M(0,1))^2 + (M(0,2))^2 + (M(0,3))^2]^{1/2}$$

The elements of the Mueller matrix, $M$, are represented by the $M(i,j)$.

The depolarization matrix, $M_{\Delta}$, is defined as

$$M_{\Delta} = \begin{bmatrix} 1 & \vec{0}^T \\ \vec{P}_{\Delta} & m_{\Delta} \end{bmatrix}$$

Here, $m_{\Delta}$ is the 3 × 3 depolarization submatrix. The parameter $P_{\Delta}$ depends on diattenuation ($d$) and polarization ($P$) [39]. The depolarization coefficients were calculated from the diagonal elements of the depolarization matrix $M_{\Delta}$. The net depolarization coefficient, $\Delta$, is measured using a relation:

$$\Delta = 1 - \{ |\text{tr}(M_{\Delta} - 1)|/3 \}$$

The retardance matrix $M_{R}$ was obtained using the relation

$$M_R = \begin{bmatrix} 1 & \vec{0}^T \\ \vec{0} & m_R \end{bmatrix}$$

The total retardance, $R$ (the coupled effect of circular and linear birefringence), was determined from retardance matrix, $M_{R}$, using the relationship:

<figure><img src="image_3.png" /><figcaption>Fig. 1 Schematic representation of Stokes–Mueller analysis with various SARS-CoV-2 models</figcaption></figure>
<figure><img src="image_4.png" /><figcaption>**Fig. 2 (A)** The blender representation of the target shapes of SARS-COV-2 with different spike numbers used in DDA computations. Model references: PDB ID 6VXX, CHARMM-GUI Archive- COVID-19 Proteins Library. **(B)** 3D rendered dipole representation of the target geometries (SARS-CoV-2) with 15, 20, and 40 spike proteins</figcaption></figure>

$$R = \cos^{-1} \left\{ \frac{\text{tr}(M_R)}{2} - 1 \right\}$$

The total retardance matrix, $M_R$, was denoted as a matrix combination for a circular retarder (optical rotation with magnitude of $\psi$) and linear retarder (having a magnitude of linear retardance = $\delta$). The schematic representation of the Stokes–Mueller analysis is shown in Fig. 1.

## Modeling SARS-CoV-2

In this work, light scattering properties of SARS-CoV-2 were simulated in terms of angle-resolved Mueller matrix elements for two incident wavelengths (200 nm and 350 nm) by using DDSCAT 7.3.0 [40–42]

based on discrete dipole approximation (DDA). Wave-lengths higher than these wavelengths (e.g., 500 nm and 1100 nm) are not reported here since at such larger wavelengths the contributions from the spikes to the overall scattering of a virion became very weak that no noticeable variation in the scattering properties could be observed as compared to a solid sphere of same size and refractive index. Notably, DDA (also referred to as the coupled dipole method) is one of the most widely used methods for calculating electromagnetic scattering properties of arbitrary shaped, inhomogeneous, anisotropic, and optically active particles [40, 43]. It employs a volume-based discretization by filling up the target volume with a finite array of polarizable points or dipoles. The interaction of these dipoles with the incident electric field and each

<table><caption>**Table 1** Relevant parameters for DDA calculation</caption><thead><tr><th rowspan="2">Effective radius (nm)<sup>1</sup></th><th colspan="2">Complex refractive index</th><th rowspan="2">Wave-length (nm)</th><th rowspan="2">No of spikes</th><th rowspan="2">No of dipoles</th><th rowspan="2">Average computation time (in seconds)</th></tr><tr><th>Virion<sup>2</sup></th><th>Nucleic acid<sup>3</sup></th></tr></thead><tbody><tr><td rowspan="8">140</td><td rowspan="4">Real part = 1.48<br>Imaginary part = 0.00</td><td rowspan="4">Real part = 1.68<br>Imaginary part = 0.00</td><td rowspan="2">200</td><td>15</td><td rowspan="4">100 k</td><td>1721</td></tr><tr><td>20</td><td>1655</td></tr><tr><td rowspan="2">350</td><td>40</td><td>1456</td></tr><tr><td>15</td><td>764</td></tr><tr><td rowspan="3">Real part = 1.48<br>Imaginary part = 0.00</td><td rowspan="3">Real part = 1.68<br>Imaginary part = 0.00</td><td>200</td><td>15</td><td>1257</td></tr><tr><td>40</td><td>1280</td></tr></tbody></table>

<footnote n="1">¹This value has been adopted from the papers of Zhu et al. [47], and Petrov [22]</footnote>

<footnote n="2">²This value has been adopted from the RI value of Influenza A Virus published by Wang et al. [48]</footnote>

<footnote n="3">³This value has been adopted from Inagaki et al. [49]</footnote>
other results in the total secondary scattered radiation. In essence, DDSCAT is the numerical implementation of DDA in the form of a freely available open-source Fortran-90 software package developed by Bruce T. Draine and Piotr J. Flatau [41, 42]. In addition to a variety of pre-defined standard and regular target geometries (e.g., cylinders, ellipsoids, and hexagonal prisms), the DDSCAT can also import DDA compatible user defined arbitrary target geometries for the calculation of their light scattering properties. Here, three different shapes of SARS-CoV corresponding to 15, 20, and 40 numbers of spike proteins on the viral capsid surface were constructed as target geometries for the calculations as shown in Fig. 2A. It was reported that about 40 spikes decorate the virion in a SARS-CoV 2 [44]. Open-source modeling tools Blender (V2.90.1) and Autodesk meshmixer (V3.5.474) were used to construct a SARS-CoV-2 3-dimensional (3D) model. Protein data bank (PDB) structure of the SARS-CoV-2 spike glycoprotein (closed state) [45] along with SARS (2002) Spike S-2P with receptor binding domain (RBD) down in the inactive conformation and SARS-CoV-2 Open Spike (S-Protein) was used as a reference to construct the spike protein. Subsequently, these target geometries were converted into DDSCAT compatible dipole arrays by using the open-source target generation tool DDSCAT Convert (https://nanohub.org/resources/ddaconvert), of which the 3D rendered representations were realized by using ParaView software (https://www.paraview.org/), as shown in Fig. 2B. The converted dipole arrays were subsequently used as the shape files in DDSCAT to calculate the light scattering properties.

## Simulation of light scattering properties of SARS-CoV 2 using discrete dipole approximation (DDA)

The target geometries were represented in the form of dipole arrays consisting of ~100 k dipoles for all the three SARS-CoV-2 models with 15, 20, and 40 spikes. These dipole numbers satisfied the so-called $|m|kd$ criterion required for the accurate calculation of Mueller matrix elements [40, 42]. The size of SARS-CoV-2 was considered to be 140 nm, and the refractive indices of the virion and the nucleic acid were considered to be 1.48 and 1.68, respectively. The relevant parameters of DDA calculation are given in Table 1. Notably, the calculated light scattering properties were averaged over the angles $β$ (from 0° to 360°), $Θ$ (from 0° to 180°), and $Φ$ (from 0° to 360°) to account for the random orientation of the target geometries

with respect to the scattering plane. However, orientation averaging significantly increased the computation time. Therefore, computations were restricted to 27 target orientations ($β = 0°, 120°, 240°; Θ = 0°, 60°, 120°; \text{and } Φ = 0°, 120°, 240°$) due to the limitations in computational resources and the enormous time required for light scattering calculations. The graphical representations of the calculated Mueller matrix elements are shown in Supplementary Fig. 1. The results are also compared with Mie calculations for a solid sphere of radius 140 nm and refractive index 1.48, obtained by using the code TUMiescat.c [46]. Remarkably, significant variations in the scattering properties were observed for all the three target geometries as compared to that of the solid sphere, which may be attributed to their structural (surface roughness and presence of spikes in virus geometries) and compositional (refractive indices for virus body and nucleic acid) differences.

## Results and discussion

The Mueller matrices of the SARS-CoV-2 models with 15, 20, and 40 spikes were illuminated with light of 0°, 45°, and 90°; subsequently, RCP polarization and output polarization states were determined. The simulation was performed using light of two different wavelengths corresponding to 200 nm and 350 nm. Using MATLAB Software, Stokes vector parameters, $S_0$, $S_1$, $S_2$, and $S_3$ were reconstructed for both input and output states. The output Stokes vectors and polarization parameters at 0°, 45°, and 90° scattering angles, calculated for these three SARS-CoV-2 models, for various input polarization states are shown in Table 2. For these scattering angles, all the three target geometries were found to exhibit greater values of Stokes polarization parameters (DOP, DOLP, and DOCP) at an incident wavelength of 350 nm than that of 200 nm. Further, for 0° input polarization, the SARS-CoV-2 models display higher DOP and DOLP values and negligibly small DOCP values irrespective of its spike numbers. Remarkably, the models follow reverse trend by exhibiting greater DOP, DOCP values and reduced DOLP values for RCP input. Further, the graph of DOP, DOLP, and DOCP as a function of spike numbers for two different wavelengths, 200 nm and 350 nm, is shown in Fig. 3. The results clearly depict the gradual decrease in DOP and DOLP values with an increase in scattering angle for the coronavirus model with 15 spikes in the case of 0° input polarization, for both 200 nm and 350 nm wavelengths.
<figure><img src="image_5.png" /><figcaption>Mueller Matrix and Lu-Chipman decomposition matrices for spike models (15, 20 and 40 spikes) at various scattering angles (0°, 45° and 90°), at two different wavelengths (200 nm and 350 nm)</figcaption></figure>
15 Spikes, 350nm

M

$$ \begin{bmatrix} 28.856 & -0.00004 & 0 & -0.0009 \\ -0.0004 & 28.86 & 0.00019 & 0 \\ 0 & -0.00019 & 28.86 & 0.00005 \\ -0.00092 & 0 & -0.00005 & 28.86 \end{bmatrix} $$

Polar decomposition

<figure><img src="image_6.png" /></figure>

$$ \begin{bmatrix} 1 & 0 & 0 & 0 \\ 0 & 1.0001 & 0 & 0 \\ 0 & 0 & 1.0001 & 0 \\ 0 & 0 & 0 & 1.0001 \end{bmatrix} $$

$$ \begin{bmatrix} 1 & 0 & 0 & 0 \\ 0 & 1 & 0 & 0 \\ 0 & 0 & 1 & 0 \\ 0 & 0 & 0 & 0 \end{bmatrix} $$

0° scattering angle

M

$$ \begin{bmatrix} 8.4282 & -1.1575 & 0.00099 & -0.00105 \\ -1.158 & 8.428 & 0.00114 & 0.00489 \\ 0.0011 & -0.00196 & 8.275 & 1.102 \\ -0.0005 & -0.0045 & -1.102 & 8.275 \end{bmatrix} $$

Polar decomposition

<figure><img src="image_7.png" /></figure>

$$ \begin{bmatrix} 1 & 0 & 0 & 0 \\ -0.0001 & 1 & 0 & 0 \\ 0 & 0 & 1 & 0 \\ 0 & 0 & 0 & 1 \end{bmatrix} $$

$$ \begin{bmatrix} 1 & 0 & 0 & 0 \\ 0 & 1 & 0.0001 & 0.0006 \\ 0 & -0.0002 & 0.9912 & 0.132 \\ 0 & -0.0005 & -0.132 & 0.9912 \end{bmatrix} $$

$$ \begin{bmatrix} 1 & -0.1373 & 0.0001 & -0.0001 \\ -0.1373 & 1 & 0 & 0 \\ 0.0001 & 0 & 0.9905 & 0 \\ -0.0001 & 0 & 0 & 0.9905 \end{bmatrix} $$

45° scattering angle

M

$$ \begin{bmatrix} 0.61439 & 0.29996 & 0.00236 & 0.00067 \\ 0.2999 & 0.6141 & 0.00137 & 0.00077 \\ 0.0018 & -0.00046 & 0.5189 & 0.1173 \\ 0.00007 & -0.00054 & -0.1173 & 0.5189 \end{bmatrix} $$

Polar decomposition

<figure><img src="image_8.png" /></figure>

$$ \begin{bmatrix} 1 & 0 & 0 & 0 \\ 0.0002 & 0.9995 & -0.0001 & -0.0002 \\ -0.0002 & -0.0001 & 0.9922 & 0 \\ 0.0005 & -0.0002 & 0 & 0.9922 \end{bmatrix} $$

$$ \begin{bmatrix} 1 & 0 & 0 & 0 \\ 0 & 1 & 0.0014 & 0.0013 \\ 0 & -0.0017 & 0.9754 & 0.2205 \\ 0 & -0.001 & -0.2205 & 0.9754 \end{bmatrix} $$

$$ \begin{bmatrix} 1 & 0.4882 & 0.0038 & 0.0011 \\ 0.4882 & 1 & 0.001 & 0.0003 \\ 0.0038 & 0.001 & 0.8727 & 0 \\ 0.0011 & 0.0003 & 0 & 0.8727 \end{bmatrix} $$

90° scattering angle
<figure><img src="image_9.png" /></figure>
20 Spikes, 350 nm

M

$$ \begin{bmatrix} 2.86E+01 & 5.73E-05 & -6.18E-08 & -2E-04 \\ 5.73E-05 & 2.86E+01 & 6.12E-05 & -1.79E-08 \\ -5.65E-08 & -5.66E-04 & 2.86E+01 & -6.12E-05 \\ -2.01E-04 & 1.39E-08 & 6.12E-05 & 2.86E+01 \end{bmatrix} $$

Polar decomposition

<figure><img src="image_10.png" /></figure>

0° scattering angle

M

$$ \begin{bmatrix} 8.27 & -1.18 & -1.49E-03 & 1.72E-03 \\ -1.18 & 8.27 & -1.12 & -1.04E-02 \\ -2.02E-03 & 5.07E-03 & 8.11 & 1.12 \\ 4.94E-04 & 9.62E-03 & -1.12 & 8.11 \end{bmatrix} $$

Polar decomposition

<figure><img src="image_11.png" /></figure>

45° scattering angle

M

$$ \begin{bmatrix} 6.07E-01 & 2.93E-01 & -3.61E-03 & 1.95E-03 \\ 2.93E-01 & 6.07E-01 & -1.35E-01 & 2.47E-03 \\ -2.35E-03 & 9.03E-04 & 5.1E-01 & 1.35E-01 \\ 1.29E-03 & -2.53E-03 & -1.35E-01 & 5.10E-01 \end{bmatrix} $$

Polar decomposition

<figure><img src="image_12.png" /></figure>

90° scattering angle
<figure><img src="image_13.png" /></figure>
<figure><img src="image_14.png" /></figure>

Table 3 shows the Lu–Chipman decomposition matrices, such as depolarization ($M_Δ$), retardance ($M_R$), and diattenuation ($M_D$) of the SARS-CoV-2 models at 0°, 45°, and 90° scattering angles for both 200 nm and 350 nm incident wavelengths. It can be observed that all the spike models yielded negative depolarization (Δ) values, prominently at 0°, 45°, and 90° scattering angles for both incident wavelengths. At 0° and 45° scattering angles, all the models were found to exhibit greater diattenuation (d) value for 200 nm incident wavelength, whereas, at 90° scattering angle, a higher diattenuation value was observed at 350 nm for all the spike numbers. Further, the diattenuation value was observed to be the highest (≥ 1) at 90° scattering angle at 350 nm. For the 15-spike model, a negative optical rotation (ψ) was observed for both 45° and 90° scattering angles at 350 nm incident wavelength; however, in case of the models with 20 and 40 spikes, the negative optical rotation was observed at 200 nm incident wavelength. The Lu–Chipman decomposition parameters such as diattenuation (d), depolarization (Δ), retardance (R), optical rotation (ψ), and linear retardance (δ) for spike models (15, 20, and 40 spikes) at various scattering angles (0°, 45°, and 90°), at two different wavelengths (200 nm and 350 nm) are given in Table 3. The graph of Lu–Chipman decomposition parameters such as diattenuation (d), linear retardance (δ), and optical

rotation (ψ) as a function of spike number at 200 nm and 350 nm for 0°, 45°, and 90° scattering angle is shown in Fig. 4.

As the scattering angle increases, decomposition parameter values such as diattenuation and linear retardance are found to exhibit a gradual increase in their values, irrespective of the number of spikes in the SARS-CoV-2 model. At the scattering angle of 0°, all the models show variation in diattenuation values for 200 nm incident wavelength, whereas diattenuation was found to be 0 for 350 nm incident wavelength. Further, at 45° scattering angle, the sample exhibits diattenuation and retardance values for both 200 nm and 350 nm but does not exhibit any optical rotation for 350 nm. For 90° scattering angle, the sample shows variation in every decomposition parameter irrespective of the spike model and incident wavelength.

## Conclusion

During the past two years, humankind has witnessed unprecedented threat from SARS-CoV-2 that caused millions of deaths and adversely modulated lives and livelihoods throughout the world. The spike proteins of SARS-CoV-2 is the key of the virus that enables it to bind and invade human cells, and the virus infectivity increases significantly with a larger number of spike proteins. In this work, we have carried out light scattering calculations on SARS-CoV-2 models with 15, 20, and 40 numbers of spike proteins on the viral capsid surface to investigate how these numbers affect the light scattering properties by using DDSCAT 7.3.0 software package based on discrete dipole approximation. Notably, significant variations in the scattering properties were observed for all the three target geometries compared to that of a solid sphere of equal size, suggesting that the number of spikes is an essential parameter to understand the interaction of light with such virus particles. Further, the SARS-CoV-2 3D model assessed through Stokes polarimetry and Lu–Chipman decomposition approach revealed the individual polarization properties such as DOP, DOLP, and DOCP. Remarkably, the coronavirus models display higher DOP and DOLP values and vanishingly small DOCP values for 0° input polarization irrespective of their spike numbers. This indicates the presence of chiral/helical structures in the SARS-CoV-2. However, the models follow reverse trend for RCP input by exhibiting greater DOP and DOCP values but negligible DOLP
<table><caption>**Table 2** Stokes vectors and various polarization parameters</caption><thead><tr><th>S spikes</th><th>Wavelength</th><th>Scatting angle (θ)</th><th>Input</th><th>S<sub>in</sub></th><th>S<sub>out</sub></th><th>DOP</th><th>DOLP</th><th>DOCP</th></tr></thead><tbody><tr><td rowspan="15">15</td><td rowspan="15">200 nm</td><td rowspan="4">0°</td><td>0°</td><td>[1 1 0 0]<sup>T</sup></td><td>[1 1 0.0001 0]<sup>T</sup></td><td>1</td><td>1</td><td>0</td></tr><tr><td>45°</td><td>[1 0 1 0]<sup>T</sup></td><td>[1 -0.0001 1 0]<sup>T</sup></td><td>1</td><td>1</td><td>0</td></tr><tr><td>90°</td><td>[1 -1 0 0]<sup>T</sup></td><td>[1 -1 0 1]<sup>T</sup></td><td>1</td><td>1</td><td>0</td></tr><tr><td>RCP</td><td>[1 0 0 1]<sup>T</sup></td><td>[1 0 0 1]<sup>T</sup></td><td>1</td><td>0</td><td>1</td></tr><tr><td rowspan="4">45°</td><td>0°</td><td>[1 1 0 0]<sup>T</sup></td><td>[1.3387 1.3371 -0.0024 -0.0049]<sup>T</sup></td><td>0.9988</td><td>0.9988</td><td>3.70E-03</td></tr><tr><td>45°</td><td>[1 0 1 0]<sup>T</sup></td><td>[0.9998 0.3409 0.9275 0.1348]<sup>T</sup></td><td>1.1356</td><td>1.1276</td><td>0.1348</td></tr><tr><td>90°</td><td>[1 -1 0 0]<sup>T</sup></td><td>[0.6613 -0.66 0.0009 0.0008]<sup>T</sup></td><td>0.9988</td><td>0.998</td><td>0.0012</td></tr><tr><td>RCP</td><td>[1 0 0 1]<sup>T</sup></td><td>[0.9989 0.3404 -0.1375 0.9256]<sup>T</sup></td><td>0.9968</td><td>0.3675</td><td>9.27E-01</td></tr><tr><td rowspan="4">90°</td><td>0°</td><td>[1 1 0 0]<sup>T</sup></td><td>[1.4158 1.41 0.0215 -0.0052]<sup>T</sup></td><td>0.996</td><td>0.996</td><td>3.7E-03</td></tr><tr><td>45°</td><td>[1 0 1 0]<sup>T</sup></td><td>[1.0153 0.4169 0.8551 0.2878]<sup>T</sup></td><td>0.9789</td><td>0.937</td><td>0.2835</td></tr><tr><td>90°</td><td>[1 -100]<sup>T</sup></td><td>[0.5842 -0.579 0.0099 -0.0044]<sup>T</sup></td><td>0.9913</td><td>0.9912</td><td>0.0075</td></tr><tr><td>RCP</td><td>[1 0 0 1]<sup>T</sup></td><td>[0.9928 0.4145 -0.2763 0.8359]<sup>T</sup></td><td>0.9801</td><td>0.5018</td><td>0.842</td></tr><tr><td rowspan="15">350 nm</td><td rowspan="15">200 nm</td><td rowspan="4">0°</td><td>0°</td><td>[1 1 0 0]<sup>T</sup></td><td>[1 1.0001 0 -0.0002]<sup>T</sup></td><td>1.0001</td><td>1.0001</td><td>2E-04</td></tr><tr><td>0°</td><td>[1 1 0 0]<sup>T</sup></td><td>[1 0 1.0001 0]<sup>T</sup></td><td>1.0001</td><td>1.0001</td><td>0</td></tr><tr><td>45°</td><td>[1 0 1 0]<sup>T</sup></td><td>[1 -1.0001 0 0]<sup>T</sup></td><td>1.0001</td><td>1.0001</td><td>0</td></tr><tr><td>RCP</td><td>[1 0 0 1]<sup>T</sup></td><td>[1 0 0 1.0001]</td><td>1.0001</td><td>0</td><td>1.0001</td></tr><tr><td rowspan="4">45°</td><td>0°</td><td>[1 1 0 0]<sup>T</sup></td><td>[0.8627 0.8626 -0.0001 -0.0006]<sup>T</sup></td><td>0.9999</td><td>0.9999</td><td>6.95E-04</td></tr><tr><td>45°</td><td>[1 0 1 0]<sup>T</sup></td><td>[1.0001 -0.1373 0.982 -0.1308]<sup>T</sup></td><td>1</td><td>0.9915</td><td>1.31E-01</td></tr><tr><td>90°</td><td>[1 -100]<sup>T</sup></td><td>[1.1373 -1.1374 0.0004 0.0005]<sup>T</sup></td><td>1.0001</td><td>1.0001</td><td>4.40E-04</td></tr><tr><td>RCP</td><td>[1 0 0 1]<sup>T</sup></td><td>[0.9999 -0.1368 0.1309 0.9818]<sup>T</sup></td><td>1</td><td>0.1894</td><td>0.9819</td></tr><tr><td rowspan="4">90°</td><td>0°</td><td>[1 1 0 0]<sup>T</sup></td><td>[1.4882 1.4877 0.0022 -0.0008]<sup>T</sup></td><td>0.9997</td><td>0.9997</td><td>5E-04</td></tr><tr><td>45°</td><td>[1 0 1 0]<sup>T</sup></td><td>[1.0038 0.4904 0.8475 -0.1908]<sup>T</sup></td><td>0.9938</td><td>0.9755</td><td>0.1901</td></tr><tr><td>90°</td><td>[1 -1 0 0]<sup>T</sup></td><td>[0.5118 -0.5114 0.0037 0.001]<sup>T</sup></td><td>0.9992</td><td>0.9992</td><td>2E-03</td></tr><tr><td>RCP</td><td>[1 0 0 1]<sup>T</sup></td><td>[1.0011 0.4894 0.1939 0.8447]<sup>T</sup></td><td>0.9942</td><td>0.5258</td><td>0.8438</td></tr><tr><td rowspan="15">20</td><td rowspan="15">200 nm</td><td rowspan="4">0°</td><td>0°</td><td>[1 1 0 0]<sup>T</sup></td><td>[1 0.9999 0 -0.0002]<sup>T</sup></td><td>0.9999</td><td>0.9999</td><td>2E-04</td></tr><tr><td>45°</td><td>[1 0 1 0]<sup>T</sup></td><td>[1 0 0.9999 -0.0002]<sup>T</sup></td><td>0.9999</td><td>0.9999</td><td>2E-04</td></tr><tr><td>90°</td><td>[1 -1 0 0]<sup>T</sup></td><td>[1 -0.9999 0 -0.0002]<sup>T</sup></td><td>0.9999</td><td>0.9999</td><td>2E-04</td></tr><tr><td>RCP</td><td>[1 0 0 1]<sup>T</sup></td><td>[0.9998 0 0 0.9997]<sup>T</sup></td><td>0.9999</td><td>0</td><td>1</td></tr><tr><td rowspan="4">45°</td><td>0°</td><td>[1 1 0 0]<sup>T</sup></td><td>[1.3093 1.3079 -0.0019 -0.0045]<sup>T</sup></td><td>0.9989</td><td>0.9989</td><td>3.4E-03</td></tr><tr><td>45°</td><td>[1 0 1 0]<sup>T</sup></td><td>[0.9956 0.3956 0.941 0.0835]<sup>T</sup></td><td>1.0287</td><td>1.0253</td><td>0.0839</td></tr><tr><td>90°</td><td>[1 -1 0 0]<sup>T</sup></td><td>[0.6907 -0.6896 -0.0057 -0.0012]<sup>T</sup></td><td>0.9984</td><td>0.9984</td><td>0.0017</td></tr><tr><td>RCP</td><td>[1 0 0 1]<sup>T</sup></td><td>[0.998 0.31 -0.0898 0.9412]<sup>T</sup></td><td>0.997</td><td>0.3234</td><td>9.43E-01</td></tr><tr><td rowspan="4">90°</td><td>0°</td><td>[1 1 0 0]<sup>T</sup></td><td>[1.4318 1.4239 0.0216 -0.022]<sup>T</sup></td><td>0.9947</td><td>0.9946</td><td>1.54E-02</td></tr><tr><td>45°</td><td>[1 0 1 0]<sup>T</sup></td><td>[1.0058 0.5947 0.8723 0.1518]<sup>T</sup></td><td>1.0604</td><td>1.0496</td><td>0.1509</td></tr><tr><td>90°</td><td>[1 -1 0 0]<sup>T</sup></td><td>[0.5682 -0.5614 -0.0001 -0.0011]<sup>T</sup></td><td>0.988</td><td>0.988</td><td>0.0019</td></tr><tr><td>RCP</td><td>[1 0 0 1]<sup>T</sup></td><td>[0.9917 0.4421 -0.1519 0.8523]<sup>T</sup></td><td>0.9802</td><td>0.4714</td><td>8.59E-01</td></tr><tr><td rowspan="15">350 nm</td><td rowspan="15">200 nm</td><td rowspan="4">0°</td><td>0°</td><td>[1 1 0 0]<sup>T</sup></td><td>[1 1 0 0]<sup>T</sup></td><td>1</td><td>1</td><td>0</td></tr><tr><td>45°</td><td>[1 0 1 0]<sup>T</sup></td><td>[1 0 1 0]<sup>T</sup></td><td>1</td><td>1</td><td>0</td></tr><tr><td>90°</td><td>[1 -1 0 0]<sup>T</sup></td><td>[1 -1 0 0]<sup>T</sup></td><td>1</td><td>1</td><td>0</td></tr><tr><td>RCP</td><td>[1 0 0 1]<sup>T</sup></td><td>[1 0 0 1]<sup>T</sup></td><td>1</td><td>0</td><td>1</td></tr><tr><td rowspan="4">45°</td><td>0°</td><td>[1 1 0 0]<sup>T</sup></td><td>[0.8573 0.8573 0.0004 0.0012]<sup>T</sup></td><td>1</td><td>1</td><td>1.4E-03</td></tr><tr><td>45°</td><td>[1 0 1 0]<sup>T</sup></td><td>[0.9998 -0.2781 0.9804 -0.1354]<sup>T</sup></td><td>1.0282</td><td>1.0193</td><td>1.35E-01</td></tr><tr><td>90°</td><td>[1 -1 0 0]<sup>T</sup></td><td>[1.1427 -1.1427 -0.0009 -0.0011]<sup>T</sup></td><td>1</td><td>1</td><td>9.62E-04</td></tr><tr><td>RCP</td><td>[1 0 0 1]<sup>T</sup></td><td>[0.9999 0 0 -0.9986]<sup>T</sup></td><td>1.0002</td><td>0.1974</td><td>0.9805</td></tr><tr><td rowspan="4">90°</td><td>0°</td><td>[1 1 0 0]<sup>T</sup></td><td>[1.4827 1.4827 -0.00024 -0.002]<sup>T</sup></td><td>1</td><td>1</td><td>1.4E-03</td></tr><tr><td>45°</td><td>[1 0 1 0]<sup>T</sup></td><td>[0.9941 0.2603 0.8303 -0.2203]<sup>T</sup></td><td>0.9085</td><td>0.8811</td><td>0.2216</td></tr><tr><td>90°</td><td>[1 -1 0 0]<sup>T</sup></td><td>[0.5173 -0.5173 -0.0054 0.0063]<sup>T</sup></td><td>1.0001</td><td>1.0001</td><td>1.22E-02</td></tr><tr><td>RCP</td><td>[1 0 0 1]<sup>T</sup></td><td>[1.0032 0.4868 0.2185 0.8423]<sup>T</sup></td><td>0.9939</td><td>0.5319</td><td>0.8396</td></tr></tbody></table>
<table><caption>**Table 2 (continued)**</caption><thead><tr><th>Spikes</th><th>Wavelength</th><th>Scatteing angle (θ)</th><th>Input</th><th>S<sub>in</sub></th><th>S<sub>out</sub></th><th>DOP</th><th>DOLP</th><th>DOCP</th></tr></thead><tbody><tr><td rowspan="12">40</td><td rowspan="12">200 nm</td><td rowspan="4">0°</td><td>0°</td><td>[1 1 0 0]<sup>T</sup></td><td>[1 1.0001 – 0.0001 0.0001]<sup>T</sup></td><td>1.0001</td><td>1.0001</td><td>1E–04</td></tr><tr><td>45°</td><td>[1 0 1 0]<sup>T</sup></td><td>[1 0 1.0001 0.0001]<sup>T</sup></td><td>1.0001</td><td>1.0001</td><td>1E–04</td></tr><tr><td>90°</td><td>[1 – 1 0 0]<sup>T</sup></td><td>[1 – 1.0001 0.0001 0.0001]<sup>T</sup></td><td>1.0001</td><td>1.0001</td><td>1E–04</td></tr><tr><td>RCP</td><td>[1 0 0 1]<sup>T</sup></td><td>[1 0 0 0.9999]<sup>T</sup></td><td>0.9999</td><td>0</td><td>0.9999</td></tr><tr><td rowspan="4">45°</td><td>0°</td><td>[1 1 0 0]<sup>T</sup></td><td>[1.3285 1.3281 0.0005 – 0.0019]<sup>T</sup></td><td>0.9997</td><td>0.9997</td><td>1.4E–03</td></tr><tr><td>45°</td><td>[1 0 1 0]<sup>T</sup></td><td>[0.9996 0.4445 0.9364 0.1163]<sup>T</sup></td><td>1.0435</td><td>1.037</td><td>1.16E–01</td></tr><tr><td>90°</td><td>[1 – 1 0 0]<sup>T</sup></td><td>[0.6715 – 0.6712 0.0009 0.0023]<sup>T</sup></td><td>0.9996</td><td>0.9996</td><td>3.40E–03</td></tr><tr><td>RCP</td><td>[1 0 0 1]<sup>T</sup></td><td>[1.0021 0.3313 – 0.1163 0.9368]<sup>T</sup></td><td>0.9983</td><td>0.3504</td><td>0.9348</td></tr><tr><td rowspan="4">90°</td><td>0°</td><td>[1 1 0 0]<sup>T</sup></td><td>[1.3873 1.3832 – 0.0167 0.0108]<sup>T</sup></td><td>0.9971</td><td>0.9971</td><td>7.8E–03</td></tr><tr><td>45°</td><td>[1 0 1 0]<sup>T</sup></td><td>[0.989 0.6566 0.8632 0.2732]<sup>T</sup></td><td>1.1309</td><td>1.0966</td><td>0.2762</td></tr><tr><td>90°</td><td>[1 – 1 0 0]<sup>T</sup></td><td>[0.6127 – 0.6112 – 0.0058 – 0.0057]<sup>T</sup></td><td>0.9976</td><td>0.9976</td><td>0.0093</td></tr><tr><td>RCP</td><td>[1 0 0 1]<sup>T</sup></td><td>[1.0014 0.3792 – 0.2817 0.8764]<sup>T</sup></td><td>0.9942</td><td>0.4717</td><td>0.9942</td></tr><tr><td rowspan="12">350 nm</td><td rowspan="4">0°</td><td>0°</td><td>[1 1 0 0]<sup>T</sup></td><td>[1 1.0001 0 0]<sup>T</sup></td><td>1.0001</td><td>1.0001</td><td>0</td></tr><tr><td>45°</td><td>[1 0 1 0]<sup>T</sup></td><td>[1 0 1.0001 0]<sup>T</sup></td><td>1.0001</td><td>1.0001</td><td>0</td></tr><tr><td>90°</td><td>[1 – 1 0 0]<sup>T</sup></td><td>[1 – 1.0001 0 0]<sup>T</sup></td><td>1.0001</td><td>1.0001</td><td>0</td></tr><tr><td>RCP</td><td>[1 0 0 1]<sup>T</sup></td><td>[1 0 0 – 1.0001]<sup>T</sup></td><td>1.0001</td><td>0</td><td>1.0001</td></tr><tr><td rowspan="4">45°</td><td>0°</td><td>[1 1 0 0]<sup>T</sup></td><td>[0.8586 0.8587 0.0001 – 0.0002]<sup>T</sup></td><td>1.0001</td><td>1.0001</td><td>2.33E–04</td></tr><tr><td>45°</td><td>[1 0 1 0]<sup>T</sup></td><td>[1.0002 – 0.2746 0.9811 – 0.1333]<sup>T</sup></td><td>1.0273</td><td>1.0186</td><td>0.1333</td></tr><tr><td>90°</td><td>[1 – 1 0 0]<sup>T</sup></td><td>[1.1414 – 1.1413 0.0002 – 0.0002]<sup>T</sup></td><td>0.9999</td><td>0.9999</td><td>1.75E–04</td></tr><tr><td>RCP</td><td>[1 0 0 1]<sup>T</sup></td><td>[1 – 0.1411 0.1334 0.9809]<sup>T</sup></td><td>0.9999</td><td>0.1942</td><td>0.9809</td></tr><tr><td rowspan="4">90°</td><td>0°</td><td>[1 1 0 0]<sup>T</sup></td><td>[1.4813 1.4811 – 0.0008 0.0009]<sup>T</sup></td><td>0.9999</td><td>0.9999</td><td>6.07E–04</td></tr><tr><td>45°</td><td>[1 0 1 0]<sup>T</sup></td><td>[1 0.3045 0.8572 – 0.1761]<sup>T</sup></td><td>0.9266</td><td>0.9097</td><td>0.1761</td></tr><tr><td>90°</td><td>[1 – 1 0 0]<sup>T</sup></td><td>[0.5187 – 0.5814 0.0005 0.0006]<sup>T</sup></td><td>0.9994</td><td>0.9994</td><td>1.20E–03</td></tr><tr><td>RCP</td><td>[1 0 0 1]<sup>T</sup></td><td>[1.0007 0.4814 0.1767 0.8579]<sup>T</sup></td><td>0.9988</td><td>0.5124</td><td>0.8573</td></tr></tbody></table>

<figure><img src="image_15.png" /></figure>
<figure><img src="image_16.png" /><figcaption>**Fig. 3** Graphs of Stokes vector polarization parameters such as DOP, DOLP, and DOCP as a function of spike numbers with 0° and RCP excitation polarization states at two incident wavelengths 200 nm and 350 nm for (a) 0°, (b) 45°, and (c) 90° scattering angle</figcaption></figure>
<table><caption>**Table 3** Lu-Chipman decomposition parameters</caption><thead><tr><th rowspan="2">Spikes</th><th rowspan="2">Wavelength</th><th rowspan="2">Scattering angle (θ)</th><th rowspan="2">Diattenuation (d)</th><th rowspan="2">Depolarization (Δ)</th><th colspan="2">Retardance (R)</th><th rowspan="2">Optical rotation (ψ)</th><th colspan="2">Linear retardance (δ)</th></tr><tr><th>Degree</th><th>Radian</th><th>Degree</th><th>Radian</th></tr></thead><tbody><tr><td rowspan="9">15</td><td rowspan="3">200 nm</td><td>0°</td><td>1.37E-05</td><td>-3.33E-09</td><td>0.0057</td><td>9.94E-05</td><td>0.0057</td><td>0.008i</td><td>0.00013</td></tr><tr><td>45°</td><td>0.3387</td><td>2.60E-03</td><td>8.3857</td><td>0.1463</td><td>0.0086</td><td>8.3891</td><td>0.1464</td></tr><tr><td>90°</td><td>0.4161</td><td>0.0167</td><td>19.1904</td><td>0.3349</td><td>0.1503</td><td>19.1947</td><td>0.335</td></tr><tr><td rowspan="3">350 nm</td><td>0°</td><td>3.12E-05</td><td>-1E-04</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td></tr><tr><td>45°</td><td>0.1373</td><td>3.27E-05</td><td>7.5856</td><td>0.132</td><td>-0.0086</td><td>7.6067</td><td>0.132</td></tr><tr><td>90°</td><td>0.4882</td><td>0.0053</td><td>12.7397</td><td>0.222</td><td>-0.0899</td><td>12.7344</td><td>0.222</td></tr><tr><td rowspan="9">20</td><td rowspan="3">200 nm</td><td>0°</td><td>2.28E-04</td><td>1E-04</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td></tr><tr><td>45°</td><td>0.3094</td><td>1.6E-03</td><td>5.8088</td><td>0.101</td><td>-2.55</td><td>5.2</td><td>0.09</td></tr><tr><td>90°</td><td>0.4319</td><td>0.018</td><td>11.8583</td><td>0.206</td><td>-4.9427</td><td>10.7741</td><td>0.186</td></tr><tr><td rowspan="3">350 nm</td><td>0°</td><td>7.27E-06</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td></tr><tr><td>45°</td><td>0.1427</td><td>-1.70E-03</td><td>8.7671</td><td>0.152</td><td>0.0144</td><td>9.5986</td><td>0.16</td></tr><tr><td>90°</td><td>0.4827</td><td>-6.8E-04</td><td>16.5463</td><td>0.288</td><td>7.3006</td><td>14.8628</td><td>0.244</td></tr><tr><td rowspan="9">40</td><td rowspan="3">200 nm</td><td>0°</td><td>7.87E-05</td><td>-3.33E-05</td><td>0.0029</td><td>5.06E-05</td><td>0</td><td>0</td><td>0</td></tr><tr><td>45°</td><td>0.3285</td><td>-5.77E-04</td><td>7.892</td><td>0.137</td><td>-3.4992</td><td>7.018</td><td>0.124</td></tr><tr><td>90°</td><td>0.3875</td><td>-0.0013</td><td>19.1391</td><td>0.334</td><td>-8.5635</td><td>17.137</td><td>0.299</td></tr><tr><td rowspan="3">350 nm</td><td>0°</td><td>7.24E-06</td><td>-1E-04</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td></tr><tr><td>45°</td><td>0.1414</td><td>-1.5E-03</td><td>8.6389</td><td>0.15</td><td>0.0462</td><td>9.46</td><td>0.165</td></tr><tr><td>90°</td><td>0.4813</td><td>-2.30E-03</td><td>12.9879</td><td>0.226</td><td>5.7437</td><td>11.6391</td><td>0.203</td></tr></tbody></table>

NERO/2181), Assam Science Technology and Environment Council (Grant No.: ASTEC/S&T/1614/8/2018-19/1159), and Science and Engineering Research Board (SERB), Department of Science and Technology (DST) (Grant No. CRG/2019/004868), Government of India, India.

**Data availability** Data will be available on request to the corresponding author.

## Declarations

**Conflict of interest** The author declares no conflict of interest.

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

values. This observation is critical as it indicates the capability of light scattering techniques to distinguish viral RNA, a chiral genetic structure. On the other hand, Lu-Chipman-based Mueller matrix polar decomposition method was employed to obtain various polarization properties of samples. Stokes-Mueller polarimetry-based characterization can provide insightful information of virions in a remote sensing architecture.

**Supplementary Information** The online version contains supplementary material available at https://doi.org/10.1007/s10103-022-03680-3.

**Funding** Open access funding provided by Manipal Academy of Higher Education, Manipal. The work is funded by the Global Innovation and Technology Alliance (GITA), Department of Science and Technology (DST) (Project Number- GITA/DST/TWN/P-95/2021), Indian Council of Medical Research (ICMR) (Project Number-ITR/Ad-hoc/43/2020-21, ID No. 2020-3286), University Grants Commission (UGC) (Grant No. F.5-376/2014-15/MRP/

<figure><img src="image_17.png" /></figure>
<figure><img src="image_18.png" /><figcaption>**Fig. 4** Graph of Lu–Chipman decomposition parameters such as diattenuation (d), linear retardance (δ), and optical rotation (ψ), as a function of spikes number at different wavelengths 200 nm and 350 nm for (a) 0°, (b) 45°, and (c) 90° scattering angle</figcaption></figure>

## References

1. Lau SKP, Woo PCY, Li KSM, Huang Y, Tsoi HW, Wong BHL, Wong SSY, Leung SY, Chan KH, Yuen KY (2005) Severe acute respiratory syndrome coronavirus-like virus in Chinese horse-shoe bats. Proc Natl Acad Sci USA 102:14040–14045. https://doi.org/10.1073/pnas.0506735102

2. Wang M, Yan M, Xu H, Liang W, Kan B, Zheng B, Chen H, Zheng H, Xu Y, Zhang E, Wang H, Ye J, Li G, Li M, Cui Z, Liu YF, Guo RT, Liu XN, Zhan LH, Xu J (2005) SARS-CoV infection in a restaurant from palm civet. Emerg Infect Dis 1:1860–1865. https://doi.org/10.3201/eid1112.041293

3. Al-Tawfiq JA, Memish ZA (2014) Middle East respiratory syndrome-coronavirus (MERS-CoV) infection. Emerg Infect Dis 18:1–190. https://doi.org/10.1016/B978-0-12-416975-3.00014-5

4. Cucinotta D, Vanelli M (2020) WHO declares COVID-19 a Pandemic. Acta Bio-medica: Atenei Parmensis 91(1):157–160. https://doi.org/10.23750/abm.v91i1.9397

5. Cunha CB, Opal SM (2014) Middle East respiratory syndrome (MERS): a new zoonotic viral pneumonia. Virulence 5(6):650–654. https://doi.org/10.4161/viru.32077

6. White J (2014) Middle Eastern respiratory syndrome coronavirus (MERS-CoV). Clin Microbiol Newsl 36(15):115–122. https://doi.org/10.1016/j.clinmicnews.2014.07.002

7. Kim CH (2020) Sars-cov-2 evolutionary adaptation toward host entry and recognition of receptor o-acetyl sialylation in virus-host interaction. Int J Mol Sci 21:1–34. https://doi.org/10.3390/ijms21124549

8. Gorbalenya AE, Baker SC, Baric RS, de Groot RJ, Drosten C, Gulyaeva AA, Haagmans BL, Lauber C, Leontovich AM, Neuman BW, Penzar D, Perlman S, Poon LLM, Samborskiy DV, Sidorov IA, Sola I, Ziebuhr J (2020) The species severe acute respiratory syndrome-related coronavirus: classifying 2019-nCoV and naming it SARS-CoV-2. Nat Microbiol 5:536–544. https://doi.org/10.1038/s41564-020-0695-z

9. Mittal A, Manjunath K, Ranjan RK, Kaushik S, Kumar S, Verma V (2020) COVID-19 pandemic: insights into structure, function, and hACE2 receptor recognition by SARS-CoV-2. PLoS Pathog 16:e1008762
10. WHO (2020) Coronavirus Disease (COVID-19) Dashboard [Internet]. [cited 2022 Jan 26]. Available from: https://covid19.who.int/

11. "G20 GDP growth - rirst quarter of 2020, OECD" [Internet]. [cited 2020 Nov 29]. Available from: https://www.oecd.org/economy/g20-gdp-growth-first-quarter-2020-oecd.htm

12. Zhou P, Yang XL, Wang XG, Hu B, Zhang L, Zhang W, Si HR, Zhu Y, Li B, Huang CL, Chen HD (2020) A pneumonia outbreak associated with a new coronavirus of probable bat origin. Nature 579:270–273. https://doi.org/10.1038/s41586-020-2012-7

13. Tu YF, Chien CS, Yarmishyn AA, Lin YY, Luo YH, Lin YT, Lai WY, Yang DM, Chou SJ, Yang YP, Wang ML (2020) A review of SARS-CoV-2 and the ongoing clinical trials. Int J Mol Sci 21:2657. https://doi.org/10.3390/ijms21072657

14. Simmons WK, Burrows K, Avery JA, Kerr KL, Taylor A, Bodurka J, Potter W, Teague TK, Drevets WC (2020) Appetite changes reveal depression subgroups with distinct endocrine, metabolic, and immune states. Mol Psychiatry 25:1457–1468. https://doi.org/10.1038/s41380-018-0093-6

15. Qiu L, Jiao R, Zhang A, Chen X, Ning Q, Fang F, Zeng F, Tian N, Zhang Y, Huang Y, Sun Z (2020) A case of critically ill infant of coronavirus disease 2019 with persistent reduction of T lymphocytes. J Pediatr Infect Dis 39:e87-90. https://doi.org/10.1097/INF.00000000000002720

16. Reguera J, Santiago C, Mudgal G, Ordoño D, Enjuanes L, Casas-novas JM (2012) Structural bases of coronavirus attachment to host aminopeptidase N and its inhibition by neutralizing antibodies. PLoS Pathog 8:e1002859. https://doi.org/10.1371/journal.ppat.1002859

17. Prabakaran P, Gan J, Feng Y, Zhu Z, Choudhry V, Xiao X, Ji X, Dimitrov DS (2006) Structure of severe acute respiratory syndrome coronavirus receptor-binding domain complexed with neutralizing antibody. J Biol Chem 281:15829–15836. https://doi.org/10.1074/jbc.M600697200

18. Peck KM, Burch CL, Heise MT, Baric RS (2015) Coronavirus host range expansion and Middle East respiratory syndrome coronavirus emergence: biochemical mechanisms and evolutionary perspectives. Annu Rev Virol 2:95–117. https://doi.org/10.1146/annurev-virology-100114-055029

19. Mittal R, Ni R, Seo JH (2020) The flow physics of COVID-19. J Fluid Mech 894:F2. https://doi.org/10.1017/jfm.2020.330

20. Ghosh N, Wood MF, Li SH, Weisel RD, Wilson BC, Li RK, Vitkin IA (2009) Mueller matrix decomposition for polarized light assessment of biological tissues. J Biophotonics 2:145–156

21. Valagiannopoulos C, Sihvola A (2021) Maximal interaction of electromagnetic radiation with corona virions. Phys Rev B 103:014114. https://doi.org/10.1103/PhysRevB.103.014114

22. Faez S, Lahini Y, Weidlich S, Garmann RF, Wondraczek K, Zeisberger M, Schmidt MA, Orrit M, Manoharan VN (2015) Fast label-free tracking of single viruses and weakly scattering nanoparticles in a nanofluidic optical fiber. ACS Nano 9:12349–12357. https://doi.org/10.1021/acsnano.5b05646

23. Driskell JD, Jones CA, Tompkins SM, Tripp RA (2011) One-step assay for detecting influenza virus using dynamic light scattering and gold nanoparticles. Analyst 136:3083–3090. https://doi.org/10.1039/c1an15303j

24. Makra I, Terejánszky P, Gyurcsányi RE (2015) A method based on light scattering to estimate the concentration of virus particles without the need for virus particle standards. MethodsX 2:91–99. https://doi.org/10.1016/j.mex.2015.02.003

25. Guo B, Ni B, Jin X, Zhang H, Zhao H, Hou L, Marsh JH, Dong L, Li S, Xiong J, Liu X (2021) Detection of virus particles by

scattering field using the multi perspective polarization modulation imaging method. JOSA B 38:3592–3600. https://doi.org/10.1364/JOSAB.436357

26. Petrov D (2020) Photopolarimetrical properties of coronavirus model particles: spike proteins number influence. J Quant Spectrosc Radiat Transf 248:107005. https://doi.org/10.1016/j.jqsrt.2020.107005

27. Ashraf MW, Ranjan R, Diaspro A (2021) Circular intensity differential scattering of light to characterize the coronavirus particles. JOSA B 38:1702–1709. https://doi.org/10.1364/JOSAB.422646

28. Lu SY, Chipman RA (1996) Interpretation of Mueller matrices based on polar decomposition. JOSA A 13:1106–1113. https://doi.org/10.1364/JOSAA.13.001106

29. Goldstein DH (2003) Polarized light, revised and expanded. CRC press

30. Bickel WS, Bailey WM (1985) Stokes vectors, Mueller matrices, and polarized scattered light. Am J Phys 53:468–478. https://doi.org/10.1119/1.14202

31. Cronin TW, Marshall J (2011) Patterns and properties of polarized light in air and water. Philos Trans R Soc Lond Ser B Biol Sci 366(1565):619–626. https://doi.org/10.1098/rstb.2010.0201

32. Li D, Guo K, Sun Y, Bi X, Gao J, Guo Z (2021) Depolarization characteristics of different reflective interfaces indicated by indices of polarimetric purity (IPPs). Sensors (Basel, Switzerland) 21(4):1221. https://doi.org/10.3390/s21041221

33. Mazumder N, Qiu J, Foreman MR, Romero CM, Hu CW, Tsai HR et al (2012) Polarization-resolved second harmonic generation microscopy with a four-channel Stokes-polarimeter. Opt Express 20:14090–14099. https://doi.org/10.1364/OE.20.014090

34. (2002) Scattering, absorption, and emission of light by small particles. In: Michael I. Mishchenko, Larry D. Travis, Andrew A. Lacis (eds) Cambridge University press

35. Dong Y, He H, Sheng W, Wu J, Ma H (2017) A quantitative and non-contact technique to characterise microstructural variations of skin tissues during photo-damaging process based on Mueller matrix polarimetry. Sci Rep 7(1):14702. https://doi.org/10.1038/s41598-017-14804-z

36. Dong Y, Qi J, He H, He C, Liu S, Wu J, Elson DS, Ma H (2017) Quantitatively characterizing the microstructural features of breast ductal carcinoma tissues in different progression stages by Mueller matrix microscope. Biomed Opt Express 8(8):3643–3655. https://doi.org/10.1364/BOE.8.003643

37. He C, He H, Chang J, Dong Y, Liu S, Zeng N, He Y, Ma H (2015) Characterizing microstructures of cancerous tissues using multispectral transformed Mueller matrix polarization parameters. Biomed Opt Express 6(8):2934–2945. https://doi.org/10.1364/BOE.6.002934

38. Li P, Dong Y, Wan J, He H, Aziz T, Ma H (2021) Polaromics: deriving polarization parameters from a Mueller matrix for quantitative characterization of biomedical specimen. J Phys D Appl Phys 55(3):034002. https://doi.org/10.1088/1361-6463/ac292f

39. Ghosh N, Wood MFG, Vitkin IA (2008) Mueller matrix decomposition for extraction of individual polarization parameters from complex turbid media exhibiting multiple scattering, optical activity, and linear birefringence. J Biomed Opt 13:0440361–0440364. https://doi.org/10.1117/1.2960934

40. Draine BT, Piotr JF (1994) Discrete-dipole approximation for scattering calculations. JOSA A 11:1491–1499. https://doi.org/10.1364/JOSAA.11.001491

41. Draine BT, Flatau PJ (2000) DDSCAT: The discrete dipole approximation for scattering and absorption of light by irregular particles. Ascl pp.ascl-0008
42. Draine, BT, Flatau PJ (2013) User guide to the discrete dipole approximation code DDSCAT 7.3. http://arxiv.org/abs/1305.6497

43. Yurkin MA, Hoekstra AG (2007) The discrete dipole approximation: an overview and recent developments. J Quant Spectrosc Radiat Transf 106:558–589. https://doi.org/10.1016/j.jqsrt.2007.01.034

44. Turoňová B, Sikora M, Schürmann C, Hagen WJ, Welsch S, Blanc FE, von Bülow S, Gecht M, Bagola K, Hörner C, van Zandbergen G (2020) In situ structural analysis of SARS-CoV-2 spike reveals flexibility mediated by three hinges. Science 370(6513):203–208

45. Shajahan A, Archer-Hartmann S, Supekar NT, Gleinich AS, Heiss C, Azadi P (2020) Comprehensive characterization of N- and O- glycosylation of SARS-CoV-2 human receptor angiotensin converting enzyme 2. Glycobiology 31:410–424. https://doi.org/10.1093/glycob/cwaa101

46. Gogoi A, Choudhury A, Ahmed GA (2010) Mie scattering computation of spherical particles with very large size parameters using an improved program with variable speed and accuracy. J Mod Opt 57:2192–2202. https://doi.org/10.1080/09500340.2010.533206

47. Zhu N, Zhang D, Wang W, Li X, Yang B, Song J et al (2020) Novel coronavirus from patients with pneumonia in China, 2019. N Engl J Med 382:727–733. https://doi.org/10.1056/NEJMo a2001017

48. Wang S, Shan X, Patel U, Huang X, Lu J, Li J, Tao N (2010) Label-free imaging, detection, and mass measurement of single viruses by surface plasmon resonance. Proc Natl Acad Sci U S A 107:16028–16032. https://doi.org/10.1073/pnas.10052 64107

49. Inagaki T, Hamm RN, Arakawa ET, Painter LR (1974) Optical and dielectric properties of DNA in the extreme ultraviolet. J Chem Phys 61:4246–4250. https://doi.org/10.1063/1.16817 24

**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.