# Culture Medium Geometry: The Dominant Factor Affecting In Vitro RF Exposure Dosimetry

Authors: Alessandra Paffi; Francesca Apollonio; Micaela Liberti; Asher Sheppard; Giorgi Bit-Babik; Quirino Balzano
Journal: International Journal of Antennas and Propagation
 (2015)

---

## Abstract

Biological experiments that expose living cells or tissues to RF energy must have an aqueous medium to provide essential water, ions, nutrients, and growth factors. However, as we show here, the medium inherently functions as a receiving antenna that conveys RF energy to the biological entity in a manner entirely determined by exposure vessel geometry, orientation to the incident RF flux, frequency, and dielectric properties of the medium. We show for two common experimental arrangements that basic antenna theory can predict electromagnetic energy patterns that agree well with those otherwise obtained by computationally intensive methods that require specialized resources.


---

## Body

## 1. Introduction

Over the last half century researchers have exposed cell preparations to RF electromagnetic (EM) energy to explore mechanisms of interaction of living cells with EM fields [1]. Early experiments frequently were conducted by exposing biological preparations to an incident EM wave of given power density. However, the pattern of the RF field exposure of cells within the medium was not investigated for lack of analytical tools, leaving investigators with little knowledge of the distribution of energy absorbed by the cells under test. Unless there is a clear understanding of the physics of EM fields, it is possible to miss the overall character of RF absorption by a biological preparation.The last two decades [2–4] have seen substantial advancements in theoretical and experimental dosimetry. Minimal requirements for exposure systems have been established [5] and are now widely accepted. The detailed distribution of absorbed energy inside an exposure vessel requires computer algorithms that accurately model the fine details of exposure systems and exposed samples. One such numerical method is the finite difference time domain (FDTD) algorithm based on numerical procedures to solve the EM problem within a region of space [6].The FDTD method performs a point-by-point evaluation of the EM fields whereas the analytical approach proceeds from the first principle, Faraday’s law, in describing the interaction between the fields and the medium. Consequently, the analytic models reveal the physical events underlying the computational results.This paper aims to show that the culture medium of a biological preparation not only supports cell metabolism but also conveys the EM energy absorbed by the tissues and cells under test, illustrating the usefulness of the concept that the medium is an antenna whose geometric and dielectric properties entirely determine the biologically relevant exposure. The cells themselves have a passive role as absorbers of EM energy from the local environment. Equivalently, the cells can be considered as a distributed load of a dielectric antenna consisting of the medium in which they are embedded. This paper also presents analytical approaches that are simpler than numerical FDTD computations and sufficiently accurate for fast, preliminary estimation of the SAR in cylindrical vessels containing a biological preparation.


## 2. Exposure of Preparations in Petri Dishes and Test Tubes

Test tubes, flasks, and dishes of various sizes are the three vessels commonly used for the exposure of cells to RF EM energy. We limit our considerations to Petri dishes of 34 mm inner diameter (ID) and test tubes, both of which have cylindrical (rotational) symmetry. Exposure systems with other vessel geometries can be analyzed by a similar methodology. The simplest (canonical) exposure of a cylindrical structure to a propagating wave can be performed with the direction of propagation either oriented along the axis of the cylinder or orthogonal to it. Axial propagation with incidence from the bottom or the top of the vessel is designated ask+ and k-, respectively. Waves with incidence orthogonal to the axis of the cylinder have either the E- or the H-field oriented along the axis. All other exposures are a combination of these four conditions. For stationary waves, there are only three independent geometric exposures: both E- and H-field parallel to the cross section plane and H or E parallel to the axis of the vessel. While all four types of propagating wave exposure have been used in laboratory research [7–10], only the stationary waves with E and H parallel to the cylindrical cross section and the stationary H-field exposure have been employed [11, 12].


## 3. RF EM Field Propagation in Dielectric Cylinders

We focus our attention on exposures utilizing EM fields in the band 900 MHz–5 GHz, the frequencies most used for cellular communications. Our analyses of the absorption of RF energy by cell preparations in Petri dishes and test tubes treat the medium as if it were a simple cylindrical receiving antenna made of a dissipative dielectric material.Cylindrical dielectric antennas have been analyzed, designed, and used for more than a century [13–16]. The distribution of RF energy within dielectric cylinders depends strictly on the boundary conditions at their surface. Consequently, the field distribution differs significantly for systems with propagation along the cylindrical axis from others with propagation perpendicular to the axis [15–17]. RF propagation along the axis of a dielectric cylinder is presented in detail in [16], where it is shown that there are discrete modes of propagation. In the case of a plane wave incident normal to the cylindrical surface with the E- or H-field parallel to the cylindrical axis, the resulting internal fields were analyzed by Stratton [17] and Harrington [18]. All incident directions can be treated analytically if the dielectric cylinder is long enough that end effects can be disregarded. While test tubes can support axial wave propagation if the medium depth is a sufficiently large part of the wavelength, this certainly is not the case for Petri dishes, which contain only a few milliliters of medium with a typical depth of 2–4 mm, a small fraction of a wavelength in the band of our interest. Before presenting the infinite cylinder solutions for medium in a Petri dish or test tube, it is necessary to explore the influence on their validity of end effects and deviations from perfect geometric form.


### 3.1. Meniscus and Antenna Effect in Petri Dishes,E-Polarization

A concave meniscus is present for all aqueous columns that wet their glass or plastic container. Although a meniscus does not break the cylindrical symmetry of the vessels of interest, it can cause substantial end effects if it contains a significant fraction of the liquid volume. We analyze the meniscus effects treating the liquid medium as an antenna. A meniscus does not affect substantially the internal fields of the medium-antenna fork+, k-, or H-polarized waves, but can be a significant factor for E-polarization if the meniscus height is comparable to the liquid depth at the center of the container, as first reported in [4]. The reason is simple: the presence of the meniscus increases substantially the cross sectional area of the medium/antenna that couples with the magnetic flux from the incident H-field as compared to the area of the cross section of the same total liquid volume without a meniscus (Figure 7).AppendixA presents the detailed analysis of the effect of the meniscus on the average specific absorption rate (SAR) at the bottom of a Petri dish exposed to the peak H-field of a E-polarized stationary wave at 1800 MHz. It is the exposure system investigated in [4] using numerical algorithms. The antenna analysis of Appendix A is simple and direct but does not require the computational resources needed for an FDTD analysis as a hand calculator suffices. The results of neglecting the meniscus effect are shown in Figure 1 for a Petri dish with 34 mm ID. The “same volume” and “same height” labels stand for the case of the dish with the same volume of liquid with or without the meniscus and, alternatively, the same liquid height with or without the meniscus. From Figure 1, it is possible to conclude that for Petri dishes of 34 mm ID with 4 mm depth of medium (without counting the meniscus height) the meniscus end effect can be neglected as the SAR error is less than 10%.Figure 1
SAR error at the bottom of a 34 mm ID Petri dish.


### 3.2. Antenna Effects and Guided Waves in Test Tubes,k+-Polarization

We analyzed the exposure of biological preparations in test tubes to a 10 W/m2 plane wave at 837 MHz using CST Microwave Studio [19] and FDTDLab, an independently developed computational suite [20]. Both algorithms yield similar results. The meshing size used to represent the liquid in the computations was 0.3 mm. Figure 2 shows the SAR distribution in the medium. The depth of the biological liquid inside the test tube is L=62.7 mm and the radius is a=8.45 mm. The dielectric characteristics of the medium (permittivity εr and conductivity σ) are given in the figure legends. The liquid column behaves as a loop antenna in the plane orthogonal to the magnetic field (Figure 2, left side). Applying Faraday’s law to the outermost closed path that can be drawn in the medium, one finds EAV, the average value of the E-field, and SARAV, the average value of the SAR along the periphery of the cross section normal to the H-field. They are(1)EAV=ωμ0H2a×L-2-π/2a222a+L+a×2π-4=6.21V/m,SARAV=70mW/kg.The peak SAR can be approximated by applying Lenz’s law, according to which minimization of magnetic flux concatenated with the antenna requires that current induction is greater in the two long sides of the loop. The RF currents can be represented roughly by a sinusoid that peaks in the middle of the side. This procedure yields an estimated peak SAR value of π2/4 × 70 = 173 mW/kg, which is about 0.5 dB higher than the peak SAR value from FDTD methods shown in Figure 2. The detailed distribution of the SAR in the medium mass requires the use of numerical methods or a rigorous application of Lenz’s law, a lengthy procedure involving variational calculus. In contrast, the very simple antenna model just described gives an acceptable approximation of the peak value of the exposure. The E-field has negligible coupling with the liquid column, only end effects, as expected given the very small dimension of its diameter (≈0.05 λ) with respect to wavelength.Figure 2
SAR patterns in a liquid medium.The approach just described has limitations: the rectangular loop perimeter must be less than a half wavelength and two sides of the loop must be much longer than the other two in order to ignore the currents induced in the short sides.The physical details of EM energy absorption by a cylindrical liquid mass can be analyzed by varying its height (or depth) from a very small to a substantial fraction of the wavelength. For this purpose, the absorption patterns for various depths of a dielectric liquid exposed to a plane wave at 1.8 GHz were computed. The same computational tools discussed above were used and the incident EM fields again propagated in the direction of the cylindrical axis. The patterns in the planes orthogonal to theH- and E-field are shown in Figures 3 and 4, respectively. One feature is clearly observed comparing Figures 2 and 3: the depth of penetration of the RF energy is much greater at the higher frequency, in agreement with the results reported in [2]. This fact shows that there are modal excitation and wave propagation inside the cylindrical mass. The phenomenon is clearly visible in Figure 4 for the higher liquid columns where end effects do not mask the attenuated propagation of EM fields at the core of the liquid. The theory of the modes in a dielectric rod is found in [16], where the following equation is given for the fundamental TM mode: (2)-εrcJ1qaqaJ0qa+K1papaK0pa=0.In (2), q and p are the radial propagation constants inside and outside the liquid, εrc is its complex relative dielectric constant, the J’s are Bessel functions of the first kind and the K’s are modified Bessel functions of the second kind of orders 0 and 1.Figure 3
SAR patterns in the plane orthogonal to theH-field for various liquid depths with incident 1800 MHz plane wave.Figure 4
SAR patterns in the plane orthogonal to theE-field for various liquid depths with incident 1800 MHz plane wave.The axial propagation constant (β) is given by β2=εrc2π/λ2-p2, with λ = free space wavelength.The following relation also holds:pa2+qa2=εrc-12πa/λ2. Equation (2) can be solved using MATLAB [20] and the axial attenuation constant of the fundamental TM mode is found to be 2π/λ×1.51 per/wavelength [15]. This value matches the axial attenuation of the SAR in the rightmost panel of Figure 4. Computations of the propagation constant of the dielectric cylinder modes can be carried out with readily available mathematical software. If the axial attenuation constants of the modes are not very large, then the simple low frequency approximation of a loop in a constant H-field is not valid because the cylindrical mass behaves like a dielectric waveguide, a device often used as an antenna [21].


### 3.3. SAR in Petri Dishes,H-Polarization

In our analysis, the medium in the Petri dish is treated as a fluid column of infinite extent in the axial direction. Clearly this is an approximate model of the real situation where the liquid is only a few mm deep. Nonetheless, as we will show below, the model provides useful results at a sufficiently high frequency.The magnetic field of a plane wave polarized in the direction of the cylindrical axis (z) can be represented in terms of cylindrical functions [17]:(3)Hzi=H0e-jkx=H0∑n=-∞+∞j-nJnkρejnϕ,where k is the free space propagation constant k=2π/λ, ρ the radial coordinate of a cylindrical reference system with z as the axis, and Jn denotes a Bessel function of first kind and order n.The electric field of the plane wave is found from the equation(4)E→=1-jωε0∇×z0H0e-jkxwhichyields:Eyi=kHoωε0e-jkx=kHoωε0∑n=-∞+∞j-nJnkρejnϕ,Eϕi=kHoωε0cosϕ∑n=-∞+∞j-nejnϕJnkρhh=jkHoωε0∑n=-∞+∞j-nejnϕJn′kρ,Eρi=kHoωε0sinϕ∑n=-∞+∞j-nJnkρejnϕhh=kHoωε0∑n=-∞+∞j-nejnϕnJnkρkρ.In the equation above Jn′ denotes the derivative with respect to the argument of the Bessel function of first kind and order n. The magnetic field Hzs scattered by an infinite dielectric cylinder can be represented by the infinite series [17](5)Hzs=∑n=-∞+∞j-nbnHn2kρejnϕ.In (5), Hn(2)kρ is a cylindrical Bessel function of the third kind of order n and bn is the space harmonic magnetic field amplitude.The electric field components associated with the scattered magnetic field are given by(6)Eϕs=1jωε0ejnϕ∂∂ρ∑n=-∞+∞j-nbnHn2kρ,Eρs=1-jωε0∂ρ∂ϕ∑n=-∞+∞j-nbnHn2kρejnϕ.The fields inside the cylinder can be represented as a superposition of TE (with respect to the axial direction) electromagnetic fields [16]. Specifically(7)Hzcρ,ϕ=∑n=-∞+∞anψn=∑n=-∞+∞anJnkερejnϕ,Eρcρ,ϕ=-μωkε2ρ∑n=-∞+∞annJnkερejnϕ=-Z0εrc2π/λ1ρ∑n=-∞+∞annJnkερejnϕ,Eϕcρ,ϕ=-jμωkε2∑n=-∞+∞an∂∂ρψn=-jZ0εrc∑n=-∞+∞anJn′kερejnϕ.In the equations above kε is the propagation constant of the cylindrical material, given by kε=(2π/λ)εr-jσ/ωε0=(2π/λ)εrc, where σ is conductivity and εr is relative dielectric constant.The details for the evaluation of the coefficientsan are given in Appendix B, where the following equations for the electric field inside the dielectric cylinder are established:(8)Eρcρ,ϕ=∑n=-∞+∞Eρncρ,ϕ;Eϕcρ,ϕ=∑n=-∞+∞Eϕncρ,ϕ.The distribution of the power loss density (PLD) is the plot of the SAR multiplied by the mass density of the medium. The PLD in the medium is evaluated by(9)PLDρ,ϕ=σ2∑n=-∞∞Eρncρ,ϕ2+∑n=-∞∞Eϕncρ,ϕ2.Given the ID of the Petri dish (34 mm), it is worth noting that at 1 GHz ka=0.366, so relatively few harmonics have an amplitude of relevance (n≈3-4). At 5 GHz, ka=1.83, so a large number of harmonics (n≈20–30) enter into the computation of the PLD. Note that the PLD is symmetric with respect to the ϕ coordinate, as expected.Now, let us compare the results from (9) with the computations using a FDTD approach. Figure 5 shows the PLD at 0.2 mm from the bottom of a 34 mm ID Petri dish containing 4 mL of medium at 3, 4, and 5 GHz. A plane wave with electric field amplitude of 100 V/m is incident from the x direction. The PLD was computed using the commercial software CST Microwave Studio 2010 considering a 34 mm ID Petri dish of Perspex with 0.5 mm wall thickness (i.e., a 35 mm outside dish diameter) filled with 4 mL of Dulbecco’s modified Eagle’s medium (DMEM). The model takes into account the presence of the meniscus. The analysis domain is a cubic box, 40 cm on the side, filled with air. The electric properties of the culture medium are εr=42 and σ=1.6S/m and those of the sample holder are εr=2.6 and σ=0S/m. These values were used for all frequencies and both computational methods.Figure 5
PLD versus frequency at 0.2 mm from the bottom of a 35 mm Petri dish (FDTD).
(a)
f
=
3
G
H
z
(b)
f
=
4
G
H
z
(c)
f
=
5
G
H
zFor each frequency, the module was solved with radiation boundary conditions and mesh of 20 lines per wavelength, resulting in more than3×106 mesh cells.Figure6 plots the PLD patterns in an infinite (in the axial direction) cylindrical column of liquid having the same dielectric characteristics and diameter as the medium in the Petri dish of Figure 5, using the analytic method of (9). The plane wave is H-polarized in the direction of the cylindrical axis and is incident from the same direction as in Figure 5.Figure 6
PLD versus frequency in an infinite cylindrical liquid column exposed to a plane wave (analytic method).
(a)
f
=
3
G
H
z
(b)
f
=
4
G
H
z
(c)
f
=
5
G
H
zFigure 7
Faraday paths in the meniscus forming loops along the bottom of the dish.Comparing the plots in Figures5 and 6, it is clear that the infinite cylindrical model gives results that are more similar to the FDTD computations with increasing frequency, as expected when a larger number of harmonics enter the calculation. The analytical results present sharper contours of the cylindrical stationary waves boundaries because of the absence of voxel meshing and of smoothing that occurs because of field scattering from the boundaries of a liquid of short height. Notably, despite the absence of meshing and scattering, the infinite model correctly predicts the radii of curvature of the stationary cylindrical waves and the progress from single to double focusing of the incident energy. As both calculations show, at 5 GHz, a 34 mm ID Petri dish with 4 mL medium acts as a dual focus lens delivering a highly nonuniform exposure to a cell preparation at the bottom of the vessel. An almost identical pattern would be obtained by a dissipative dielectric antenna of the same diameter and infinite axial length. This last example clearly illustrates that the medium acts as a receiving antenna that delivers the incident EM energy to the biological cells at the bottom of a Petri dish in a pattern depending only on the geometry, frequency, and the dielectric characteristics of the liquid. The living cells under test are a passive load of the receiving antenna, albeit being a load of very small mass.


## 3.1. Meniscus and Antenna Effect in Petri Dishes,E-Polarization

A concave meniscus is present for all aqueous columns that wet their glass or plastic container. Although a meniscus does not break the cylindrical symmetry of the vessels of interest, it can cause substantial end effects if it contains a significant fraction of the liquid volume. We analyze the meniscus effects treating the liquid medium as an antenna. A meniscus does not affect substantially the internal fields of the medium-antenna fork+, k-, or H-polarized waves, but can be a significant factor for E-polarization if the meniscus height is comparable to the liquid depth at the center of the container, as first reported in [4]. The reason is simple: the presence of the meniscus increases substantially the cross sectional area of the medium/antenna that couples with the magnetic flux from the incident H-field as compared to the area of the cross section of the same total liquid volume without a meniscus (Figure 7).AppendixA presents the detailed analysis of the effect of the meniscus on the average specific absorption rate (SAR) at the bottom of a Petri dish exposed to the peak H-field of a E-polarized stationary wave at 1800 MHz. It is the exposure system investigated in [4] using numerical algorithms. The antenna analysis of Appendix A is simple and direct but does not require the computational resources needed for an FDTD analysis as a hand calculator suffices. The results of neglecting the meniscus effect are shown in Figure 1 for a Petri dish with 34 mm ID. The “same volume” and “same height” labels stand for the case of the dish with the same volume of liquid with or without the meniscus and, alternatively, the same liquid height with or without the meniscus. From Figure 1, it is possible to conclude that for Petri dishes of 34 mm ID with 4 mm depth of medium (without counting the meniscus height) the meniscus end effect can be neglected as the SAR error is less than 10%.Figure 1
SAR error at the bottom of a 34 mm ID Petri dish.


## 3.2. Antenna Effects and Guided Waves in Test Tubes,k+-Polarization

We analyzed the exposure of biological preparations in test tubes to a 10 W/m2 plane wave at 837 MHz using CST Microwave Studio [19] and FDTDLab, an independently developed computational suite [20]. Both algorithms yield similar results. The meshing size used to represent the liquid in the computations was 0.3 mm. Figure 2 shows the SAR distribution in the medium. The depth of the biological liquid inside the test tube is L=62.7 mm and the radius is a=8.45 mm. The dielectric characteristics of the medium (permittivity εr and conductivity σ) are given in the figure legends. The liquid column behaves as a loop antenna in the plane orthogonal to the magnetic field (Figure 2, left side). Applying Faraday’s law to the outermost closed path that can be drawn in the medium, one finds EAV, the average value of the E-field, and SARAV, the average value of the SAR along the periphery of the cross section normal to the H-field. They are(1)EAV=ωμ0H2a×L-2-π/2a222a+L+a×2π-4=6.21V/m,SARAV=70mW/kg.The peak SAR can be approximated by applying Lenz’s law, according to which minimization of magnetic flux concatenated with the antenna requires that current induction is greater in the two long sides of the loop. The RF currents can be represented roughly by a sinusoid that peaks in the middle of the side. This procedure yields an estimated peak SAR value of π2/4 × 70 = 173 mW/kg, which is about 0.5 dB higher than the peak SAR value from FDTD methods shown in Figure 2. The detailed distribution of the SAR in the medium mass requires the use of numerical methods or a rigorous application of Lenz’s law, a lengthy procedure involving variational calculus. In contrast, the very simple antenna model just described gives an acceptable approximation of the peak value of the exposure. The E-field has negligible coupling with the liquid column, only end effects, as expected given the very small dimension of its diameter (≈0.05 λ) with respect to wavelength.Figure 2
SAR patterns in a liquid medium.The approach just described has limitations: the rectangular loop perimeter must be less than a half wavelength and two sides of the loop must be much longer than the other two in order to ignore the currents induced in the short sides.The physical details of EM energy absorption by a cylindrical liquid mass can be analyzed by varying its height (or depth) from a very small to a substantial fraction of the wavelength. For this purpose, the absorption patterns for various depths of a dielectric liquid exposed to a plane wave at 1.8 GHz were computed. The same computational tools discussed above were used and the incident EM fields again propagated in the direction of the cylindrical axis. The patterns in the planes orthogonal to theH- and E-field are shown in Figures 3 and 4, respectively. One feature is clearly observed comparing Figures 2 and 3: the depth of penetration of the RF energy is much greater at the higher frequency, in agreement with the results reported in [2]. This fact shows that there are modal excitation and wave propagation inside the cylindrical mass. The phenomenon is clearly visible in Figure 4 for the higher liquid columns where end effects do not mask the attenuated propagation of EM fields at the core of the liquid. The theory of the modes in a dielectric rod is found in [16], where the following equation is given for the fundamental TM mode: (2)-εrcJ1qaqaJ0qa+K1papaK0pa=0.In (2), q and p are the radial propagation constants inside and outside the liquid, εrc is its complex relative dielectric constant, the J’s are Bessel functions of the first kind and the K’s are modified Bessel functions of the second kind of orders 0 and 1.Figure 3
SAR patterns in the plane orthogonal to theH-field for various liquid depths with incident 1800 MHz plane wave.Figure 4
SAR patterns in the plane orthogonal to theE-field for various liquid depths with incident 1800 MHz plane wave.The axial propagation constant (β) is given by β2=εrc2π/λ2-p2, with λ = free space wavelength.The following relation also holds:pa2+qa2=εrc-12πa/λ2. Equation (2) can be solved using MATLAB [20] and the axial attenuation constant of the fundamental TM mode is found to be 2π/λ×1.51 per/wavelength [15]. This value matches the axial attenuation of the SAR in the rightmost panel of Figure 4. Computations of the propagation constant of the dielectric cylinder modes can be carried out with readily available mathematical software. If the axial attenuation constants of the modes are not very large, then the simple low frequency approximation of a loop in a constant H-field is not valid because the cylindrical mass behaves like a dielectric waveguide, a device often used as an antenna [21].


## 3.3. SAR in Petri Dishes,H-Polarization

In our analysis, the medium in the Petri dish is treated as a fluid column of infinite extent in the axial direction. Clearly this is an approximate model of the real situation where the liquid is only a few mm deep. Nonetheless, as we will show below, the model provides useful results at a sufficiently high frequency.The magnetic field of a plane wave polarized in the direction of the cylindrical axis (z) can be represented in terms of cylindrical functions [17]:(3)Hzi=H0e-jkx=H0∑n=-∞+∞j-nJnkρejnϕ,where k is the free space propagation constant k=2π/λ, ρ the radial coordinate of a cylindrical reference system with z as the axis, and Jn denotes a Bessel function of first kind and order n.The electric field of the plane wave is found from the equation(4)E→=1-jωε0∇×z0H0e-jkxwhichyields:Eyi=kHoωε0e-jkx=kHoωε0∑n=-∞+∞j-nJnkρejnϕ,Eϕi=kHoωε0cosϕ∑n=-∞+∞j-nejnϕJnkρhh=jkHoωε0∑n=-∞+∞j-nejnϕJn′kρ,Eρi=kHoωε0sinϕ∑n=-∞+∞j-nJnkρejnϕhh=kHoωε0∑n=-∞+∞j-nejnϕnJnkρkρ.In the equation above Jn′ denotes the derivative with respect to the argument of the Bessel function of first kind and order n. The magnetic field Hzs scattered by an infinite dielectric cylinder can be represented by the infinite series [17](5)Hzs=∑n=-∞+∞j-nbnHn2kρejnϕ.In (5), Hn(2)kρ is a cylindrical Bessel function of the third kind of order n and bn is the space harmonic magnetic field amplitude.The electric field components associated with the scattered magnetic field are given by(6)Eϕs=1jωε0ejnϕ∂∂ρ∑n=-∞+∞j-nbnHn2kρ,Eρs=1-jωε0∂ρ∂ϕ∑n=-∞+∞j-nbnHn2kρejnϕ.The fields inside the cylinder can be represented as a superposition of TE (with respect to the axial direction) electromagnetic fields [16]. Specifically(7)Hzcρ,ϕ=∑n=-∞+∞anψn=∑n=-∞+∞anJnkερejnϕ,Eρcρ,ϕ=-μωkε2ρ∑n=-∞+∞annJnkερejnϕ=-Z0εrc2π/λ1ρ∑n=-∞+∞annJnkερejnϕ,Eϕcρ,ϕ=-jμωkε2∑n=-∞+∞an∂∂ρψn=-jZ0εrc∑n=-∞+∞anJn′kερejnϕ.In the equations above kε is the propagation constant of the cylindrical material, given by kε=(2π/λ)εr-jσ/ωε0=(2π/λ)εrc, where σ is conductivity and εr is relative dielectric constant.The details for the evaluation of the coefficientsan are given in Appendix B, where the following equations for the electric field inside the dielectric cylinder are established:(8)Eρcρ,ϕ=∑n=-∞+∞Eρncρ,ϕ;Eϕcρ,ϕ=∑n=-∞+∞Eϕncρ,ϕ.The distribution of the power loss density (PLD) is the plot of the SAR multiplied by the mass density of the medium. The PLD in the medium is evaluated by(9)PLDρ,ϕ=σ2∑n=-∞∞Eρncρ,ϕ2+∑n=-∞∞Eϕncρ,ϕ2.Given the ID of the Petri dish (34 mm), it is worth noting that at 1 GHz ka=0.366, so relatively few harmonics have an amplitude of relevance (n≈3-4). At 5 GHz, ka=1.83, so a large number of harmonics (n≈20–30) enter into the computation of the PLD. Note that the PLD is symmetric with respect to the ϕ coordinate, as expected.Now, let us compare the results from (9) with the computations using a FDTD approach. Figure 5 shows the PLD at 0.2 mm from the bottom of a 34 mm ID Petri dish containing 4 mL of medium at 3, 4, and 5 GHz. A plane wave with electric field amplitude of 100 V/m is incident from the x direction. The PLD was computed using the commercial software CST Microwave Studio 2010 considering a 34 mm ID Petri dish of Perspex with 0.5 mm wall thickness (i.e., a 35 mm outside dish diameter) filled with 4 mL of Dulbecco’s modified Eagle’s medium (DMEM). The model takes into account the presence of the meniscus. The analysis domain is a cubic box, 40 cm on the side, filled with air. The electric properties of the culture medium are εr=42 and σ=1.6S/m and those of the sample holder are εr=2.6 and σ=0S/m. These values were used for all frequencies and both computational methods.Figure 5
PLD versus frequency at 0.2 mm from the bottom of a 35 mm Petri dish (FDTD).
(a)
f
=
3
G
H
z
(b)
f
=
4
G
H
z
(c)
f
=
5
G
H
zFor each frequency, the module was solved with radiation boundary conditions and mesh of 20 lines per wavelength, resulting in more than3×106 mesh cells.Figure6 plots the PLD patterns in an infinite (in the axial direction) cylindrical column of liquid having the same dielectric characteristics and diameter as the medium in the Petri dish of Figure 5, using the analytic method of (9). The plane wave is H-polarized in the direction of the cylindrical axis and is incident from the same direction as in Figure 5.Figure 6
PLD versus frequency in an infinite cylindrical liquid column exposed to a plane wave (analytic method).
(a)
f
=
3
G
H
z
(b)
f
=
4
G
H
z
(c)
f
=
5
G
H
zFigure 7
Faraday paths in the meniscus forming loops along the bottom of the dish.Comparing the plots in Figures5 and 6, it is clear that the infinite cylindrical model gives results that are more similar to the FDTD computations with increasing frequency, as expected when a larger number of harmonics enter the calculation. The analytical results present sharper contours of the cylindrical stationary waves boundaries because of the absence of voxel meshing and of smoothing that occurs because of field scattering from the boundaries of a liquid of short height. Notably, despite the absence of meshing and scattering, the infinite model correctly predicts the radii of curvature of the stationary cylindrical waves and the progress from single to double focusing of the incident energy. As both calculations show, at 5 GHz, a 34 mm ID Petri dish with 4 mL medium acts as a dual focus lens delivering a highly nonuniform exposure to a cell preparation at the bottom of the vessel. An almost identical pattern would be obtained by a dissipative dielectric antenna of the same diameter and infinite axial length. This last example clearly illustrates that the medium acts as a receiving antenna that delivers the incident EM energy to the biological cells at the bottom of a Petri dish in a pattern depending only on the geometry, frequency, and the dielectric characteristics of the liquid. The living cells under test are a passive load of the receiving antenna, albeit being a load of very small mass.


## 4. Conclusion

This paper has shown that the details of the geometry of the medium inside a Petri dish or a test tube are extremely important to the determination of the SAR or the PLD in the vessel. Based on these facts, it should be clear that the medium acts as a dielectric receiving antenna for the incident electromagnetic energy to which the living cells are exposed. At frequencies below 2 GHz, depending on polarization and wavelength, the liquids in Petri dishes and test tubes behave as a loop antenna or a dielectric waveguide, which in fact is just another type of antenna. The PLD patterns in a shallow liquid (few mm in depth) contained in a relatively wide dish can be modeled with good approximation as an infinitely long dissipative cylindrical antenna with a diameter equal to the ID of the vessel, provided that the liquid depth is more thanλ/30. Viewing the medium as the antenna of the biological preparations has important ramifications for the results of in vitro experiments using different vessels. The same incident EM fields may result in very different exposure patterns to the test cells or tissues depending on the dielectric characteristics of the medium, its geometry (shape and size), and the material of the vessel walls (which affect the meniscus shape). Obviously, different spatial power absorption may cause different biological outcomes for in vitro experiments.Finally, the analytic methods presented here may be valuable to research biologists for evaluating the behavior of the medium as an antenna for fast, preliminary selection of the best vessels for the intended research, using widely available computational tools, if FDTD software and computational power are not accessible. In future efforts, the present analytical model (an infinitely long cylinder of dissipative dielectric) for the medium of a biological preparation exposed in a Petri dish will be refined to a liquid disk of finite height. Such a model could be used as a canonical model for testing FDTD computations. The antenna models and methods presented are not intended to replace the use of numerical algorithms for the accurate evaluation of the EM exposure of biopreparations, but they can be used for a quick, preliminary evaluation of an exposure system and as validation of FDTD computations.



---

*Source: 438962-2015-05-19.xml*