Open access peer-reviewed chapter

Fano Resonance in High-Permittivity Objects

By Xianghong Kong, Lina Qiu and Gaobiao Xiao

Submitted: March 28th 2017Reviewed: August 4th 2017Published: November 29th 2017

DOI: 10.5772/intechopen.70518

Downloaded: 1163


In this chapter, Fano resonances in simple structures with high permittivity such as spheres or core-shell particles are analyzed by Mie theory. The Mie scattering coefficients can be decomposed into slow varying backgrounds and narrow resonances, which cause the Fano resonances in scattered field. For structures of arbitrary shapes, temporal coupled-mode theory is applied to explain the Fano resonances found in the scattering cross section. At last, we analyze the periodic structures by using band diagram, and it shows that the Fano resonances can be viewed as the superposition of the Bloch wave and the Mie scattering wave.


  • Fano resonance
  • Mie theory
  • temporal coupled-mode theory
  • photonic crystal
  • sensor

1. Introduction

Fano resonance was first discovered in quantum systems to describe the asymmetrically shaped ionization spectral lines of atoms [1]. The asymmetric profile is caused by the interference between a broad background state and a narrow discrete state. The interference phenomenon also exists in electromagnetic system and was first observed by Wood [2]. With the development of metamaterials, Fano resonances have been observed in many classical oscillator systems, such as nonconcentric ring/disk cavities [3], asymmetric split rings, and dolmen structures [4]. Such Fano systems are caused by symmetry breaking of the geometry and are usually consisted of metal and dielectric. Recently, metamaterials composed of high refractive index materials have attracted researchers’ attentions since they can enhance efficiency significantly [5]. Fano resonances occur in these metamaterials usually have larger quality factor since metal is replaced by lossless high-permittivity dielectric, which makes Fano curve sharper compared with conventional metamaterials.

In this chapter, we will investigate the Fano resonances in high-permittivity objects theoretically and do the simulations to verify the accuracy of theories, which will provide a guidance for the further study and applications. The chapter is organized as follows:

In the second section, we analyze the Fano resonances in high-permittivity spheres, which are the simplest structures that can be analyzed by applying Mie theory easily.

In the third section, we will use the Mie theory to investigate the Fano resonances in core-shell particles. With more degrees of freedom for design, core-shell structures are more suitable for applications such as sensors.

In the fourth section, a new method called temporal coupled-mode theory (TCMT) is used to explain the Fano resonances found in high permittivity arbitrarily shaped objects. Combined with cylindrical wave expansion (2D) or spherical wave expansion (3D), we can use TCMT to model the Fano resonances in scattering by an arbitrary object.

In the fifth section, we do some numerical simulation on periodic array of cylinders and show that Fano resonances can be observed in transmission spectra as a result of interference of leaky guided modes of cylinders with an incident electromagnetic wave.

In the last section, we will draw a conclusion briefly.


2. Fano resonances in high-permittivity spheres

2.1. Mie theory

Mie scattering was first discovered by Mie in 1908 [6]. In spite of the long history, Mie theory still governs the forefront optical devices such as nanoantennas [7] and metamaterials [8]. It describes the scattering of a plane wave by a homogeneous sphere. The solution takes form of an infinite series of spherical multipole partial waves. For different electromagnetic modes, the positions of resonances which can be calculated by Mie theory are different. Resonance arises when the incident wave reaches an eigenmode frequency and excites localized modes in the sphere.

Let us assume the radius of sphere is a. The relative permittivities and permeabilities of sphere (r ≤ a) and embedding medium (r > a) are (ϵ1, μ1) and (ϵ, μ), respectively. The Mie scattering coefficients are [6]:


where x = ka(k = ω/cis the wavenumber of incident wave) is the size parameter of the sphere, and m=ϵ1μ1ϵμis the relative refractive index. jn(x) and hn1xstand for spherical Bessel functions and Hankel functions of the first kind, respectively.

2.2. Decomposition of Mie scattering coefficients

For simplicity, the embedding medium is considered to be vacuum in the following analysis, so we have ϵ = 1, μ = 1. We assume the relative permeability of sphere to be 1 and the relative permittivity ϵ1 to be a purely real number, so the relative refractive index m=ϵ1. Eq. (1), which represents electric scattering coefficients can be rewritten by


where xjnx/xhn1xmeans a slow varying background and ixhn1x2/ϵ1xjnϵ1xϵ1xjnϵ1xxhn1xxhn1xmeans a narrow resonance when high-permittivity dielectric sphere is considered [9].

As shown in Figure 1, squared norm of Mie coefficient |a1|2 is plotted when ϵ1 = 1000. It can be described by superposition of narrow resonance and slow varying background.

Figure 1.

Squared norm of Mie coefficient |a1|2 (blue curve), slow varying backgroundxj1xxh11x2(green curve), and narrow resonance (red dot-dash line) for a sphere withϵ1 = 1000.

Similarly, magnetic scattering coefficient ∣bn∣ can be decomposed into two parts:


where jnx/hn1xrepresents a slow varying background and ixhn1x2/ϵ1xjnϵ1xxjnϵ1xxhn1xxhn1xrepresents a narrow resonance (as shown in Figure 2).

Figure 2.

Squared norm of Mie coefficient |b1|2 (blue curve), slow varying backgroundj1xh11x2(green curve), and narrow resonance (red dot-dash line) for a sphere withϵ1 = 1000.

The slow varying backgrounds are the same as scattering coefficients of PEC spheres.

2.3. Rewrite Mie coefficients in the form of Fano function

Normalized Fano function can be expressed as 11+q2q+xx0Γ21+xx0Γ2, where x0, Γ, and q represent resonance position, resonance width, and Fano parameter, respectively. Compared with conventional Lorentz resonance, Fano resonance will exhibit asymmetric line shape and usually has sharper resonant curve.

When the permittivity of sphere is high, we can rewrite the Mie coefficients in the form of Fano function. The resonance position, resonance width, and Fano parameter can be achieved by following Eqs. (10):


In Eqs. (5)(7), Re(x) and Im(x) mean the real and imaginary part of x, yn(x) represent the spherical Neumann functions, sign(x) denotes the sign function. These equations are the rewrite of electric scattering coefficients. Similarly, magnetic scattering coefficients can also be rewritten in the form of Fano function [10].

As shown in Figure 3, the approximate model which can be written in the form of Fano function matches well with the exact Mie scattering coefficient.

Figure 3.

The exact value of Mie scattering coefficients (blue line) and Fano curve predicted by approximate model (red dot-dash line) are shown for (a) electric dipole |a1|2 whenω0 = 3 × 1015 rad/s , a = 64.33 nm , ϵ1 = 1000 and (b) magnetic dipole |b1|2 whenω0 = 3 × 1015 rad/s , a = 148.97 nm , ϵ1 = 1000.

3. Fano resonances in core-shell particles with high-permittivity covers

In most researches [11, 12], Fano resonances observed in coated spheres are derived in the Rayleigh limit. However, the approximation may suffer a loss of precision when frequency gets higher. An exact analysis based on Mie theory is proposed to analyze Fano resonances by coated spheres with high-permittivity covers in a precise way [13].

3.1. Theoretical analysis

Let us assume the inner radius of core-shell particle is a1 and the outer radius is a. The ratio of a1 and acan be denoted as η = a1/a. The relative permittivities and permeabilities of core (0 < r ≤ a1), shell (a1 < r ≤ a), and embedding medium (r > a) are (ϵ1, μ1), (ϵ2, μ2), and (ϵ0, μ0), respectively. The solution of scattering by coated spheres can be described by Mie theory. For simplicity, the embedding medium is considered to be vacuum in the following analysis which means ϵ0 = 1, μ0 = 1. Also, we assume μ1 = μ2 = 1, which means both core and shell are nonmagnetic. When the core-shell particles are covered by high-permittivity dielectric shells, we can decompose the scattering coefficients cnTMand cnTEinto slow varying backgrounds and narrow resonances, which are similar to the high-permittivity spheres. For electric scattering coefficients, we have




x = k0ais the size parameter of outer sphere. snTMrepresents the slow varying background and its expression is given in Eq. (9). As we can see, the background is the same as the electric scattering coefficient of a PEC sphere with radius a. rnTMrepresents the narrow resonance, and it can be calculated formally by subtracting snTMfrom cnTM. The expression for narrow resonance can be found in [13].

Similarly, we can decompose magnetic scattering coefficients into two parts:




As shown in Figure 4, the scattering coefficients can be viewed as the cascade of Fano resonances.

Figure 4.

Squared norm of Mie coefficient (a)c1TM2(blue curve), slow varying backgroundxj1xxh11x2(green curve) and narrow resonances (red dot-dash line) (b)c1TE2(blue curve), slow varying backgroundj1xh11x2(green curve), and narrow resonances (red dot-dash line) for a core-shell particle witha = 100 nm , a1 = 80 nm , ϵ1 = 10 , ϵ2 = 1000.

3.2. Application of sensors

Due to the sharp resonances near the resonance frequencies, Fano resonances have great potential applications in sensing problems [1416]. Although some of them may have high sensitivity, the structures which are designed to produce Fano resonances are usually complicated and cannot be analyzed by formula exactly. Because of the simple structure, core-shell particles have the potential to be a great platform for sensing since they can be fabricated easily. In fact, core-shell particles consisted of metal and dielectric, can exhibit Fano resonances due to the hybridization between the plasmon resonances of the core and shell [17, 18]. However, the loss in metal may flatten the shape of Fano curve, which affects the sensitivity of Fano resonance sensor. Hence, we use lossless high-permittivity dielectric to replace the metal.

For the high-permittivity shell sensors, we can fill the core with unknown materials. The permittivity of the unknown material can be varied continuously such as liquid solvents. By detecting the scattering field over a discrete set of frequencies near Fano resonance position, we can achieve the permittivity we want with high accuracy.

The sensitivity for Fano resonance sensing can be examined by comparing the changes in the scattering coefficients between core-shell structure with high-permittivity shell and homogeneous sphere when the permittivity of material changes. We can define the sensitivity as an analogy to [14]


As shown in Figure 5(a), the difference of sensitivity between core-shell structure and sphere is plotted. As for the core-shell structure, the relative permittivity of core is increased by ∆ϵ1 = 0.1. The maximum value of c1TM2occurs at x = 0.49495 (located by a vertical blue line) when 0.493 ≤ x ≤ 0.496, which is c1TM2=0.1971. In order to make a comparison with Fano resonance sensor, we figure out the scattering coefficients of a sphere with different permittivities as given for the core in core-shell structure. As shown in Figure 5(b), the maximum value c1TM2=4.2259×105occurs at x = 0.496, which shows that the Fano resonance sensor offers a high sensitivity.

Figure 5.

(a)c1TM2as a function ofxfor core-shell structure with different core permittivitiesϵ1 = 1.4 (blue line) andϵ1 = 1.5 (red line), high-permittivity shellϵ2 = 1000 , η = 0.8. (b)c1TM2as a function ofxfor sphere structure withϵ1 = 1.4 (blue line) andϵ1 = 1.5 (red line).

Since Fano resonances of high permittivity core-shell particles mentioned above only exist in scattering coefficients of multipole partial waves, it is difficult to achieve these coefficients separately by measuring the electromagnetic field distribution around the scatterers. In fact, by choosing operating frequency range properly, we can achieve the scattering coefficient of a single partial wave (c1TMfor example) without filtering out from the total electromagnetic fields. When size parameter xis small, most of the scattering coefficients cnTM, cnTEare zero, except for several coefficients with small n. We choose the first resonance frequency of c1TMas the operating frequency. Since Fano resonance is usually sharp and narrow, we find that overlap between different modes can be avoided if the frequency range narrows down. To explain it explicitly, we define the scattering cross section as [19]


The contribution of c1TMto Eq. (13) can be defined as


As shown in Figure 6, when we choose the first Fano resonance position of c1TMas the operating frequency, we find that the scattering cross section of multipole partial waves is the same as the scattering cross section of c1TMfor a narrow frequency range. As shown in Figure 6(a), Fano resonance can also be observed in scattering cross section which can be used to sense the permittivity.

Figure 6.

For the core-shell particle with high-permittivity shellϵ2 = 1000, η = 0.8, (a)log10(Qsca) is plotted as a function ofxandϵ1 when the summation inEq. (13)is truncated ton = 5. (b)log10Q1TM(calculated byEq. (14)) is plotted as a function ofxandϵ1.

The asymmetry parameter ⟨cosΘ⟩ is defined as the average cosine of the scattering angle Θ. For a spherical particle, the asymmetry parameter can be calculated by [20]


The asymmetry parameter is positive if the particle scatters more light toward the forward direction while it is negative if more light is scattered toward the backscattering direction.

As shown in Figure 7, the width between maximum and minimum for a fixed core permittivity ϵ1 narrows down compared with scattering cross section which is shown in Figure 6. With the increase of size parameter x, the asymmetry parameter reaches its maximum and decreases sharply to its minimum.

Figure 7.

The asymmetry parameters for high-permittivity shell particles withϵ2 = 1000 , η = 0.8.

To check the average scattering direction changes from front to back, we use numerical simulation software COMSOL 5.0 to simulate the scattering of a plane wave by a core-shell structure. The incident wave travels in the +z-direction and the electric field is oriented in the x-direction. Draw a horizontal line at ϵ1 = 1.5 in Figure 7, and we can find the average cosine of the scattering angle ⟨cosΘ⟩ has a lineshape of Fano resonance as a function of size parameter x. As shown in Figure 8(a), the high-permittivity shell structure scatters more light to the forward direction at x = 0.498672. When xincreases, the asymmetry parameter decreases sharply from positive value to negative value. The minimum value is achieved at x = 0.499276 and the scattering wave is concentrated in the backward direction as shown in Figure 8(c). Among the maximum value and the minimum value of asymmetry parameter, we find ⟨cosΘ⟩ ≈ 0 at x = 0.498980, which means the scattering is symmetric with respect to the plane z = 0.

Figure 8.

Scattering pattern of a high-permittivity shell particle withϵ2 = 1000 , ϵ1 = 1.5 , η = 0.8 , a = 100 nm when (a)x = 0.498672 (b)x = 0.498980 (c)x = 0.499276.

Hence, Fano resonances in core-shell particles can be used to detect the slight changes of core permittivity since they are sensitive in both magnitude and direction.

4. Fano resonances in arbitrary objects with high-permittivity dielectric

When the structure gets more complicated, the Mie theory is no longer valid for the solution of scattering field. We have to use the temporal coupled-mode theory (TCMT) to replace Mie theory when investigating the Fano resonances in arbitrary objects with high-permittivity dielectric.

4.1. Temporal coupled-mode theory

The temporal coupled-mode theory provides a useful general framework to study the interaction of a resonance with external waves. It has been well developed when dealing with particles that have cylindrical or spherical shapes [21]. In [22], TCMT has been generalized to analysis the scattering of arbitrary shape structures. The temporal coupled-mode equations can be expressed as [23]


In Eq. (16), |A|2 corresponds to the energy inside the resonator. s+ and s represent incoming waves and outgoing waves, respectively. They couple directly by the resonant mode Ais coupled with the outgoing waves s through dand is excited by the incoming waves s+ through κT. ω0 is the resonance frequency and 1τis the external leakage rate.

There exists some constrains between B, d,and κ,which are imposed by energy conservation and time-reversal invariance [22]. The constrain conditions are


For a 2D arbitrary object, we can expand scattering field into cylindrical waves


s+ and s in Eq. (16) can be viewed as coefficients of input wave and outgoing wave on the basis of cylindrical waves.

The incident plane wave can also be expanded into cylindrical waves


where θ0 is the incident angle. Combined with cylindrical wave expansion, we can use TCMT to describe the Fano resonances in arbitrary objects with high-permittivity dielectric.

4.2. Numerical simulation

The method we determine the coefficients in Eq. (16) is similar to the method described in [22]. Firstly, we use eigenmode analysis in COMSOL 5.0 to figure out the resonance frequency ω0 and the external leakage rate 1τ. Secondly, being different from the method in [22] where they set B = I, we calculate the Bthrough the simulation results of the arbitrary object covered by PEC illuminated by the plane wave. Since the slow varying background of high-permittivity sphere is the PEC sphere as mentioned above, it is intuitive to assume the slow varying background of arbitrary object which is described by Bis the same as the object covered by PEC. Thirdly, combined with the field distribution of eigenmode simulation and background scattering matrix B, we can figure out the resonant radiation coefficients d. At last, κcan be solved through constrain conditions in Eq. (17).

Once the coefficients in Eq. (16) are determined, we can use the TCMT to predict the scattering fields by different incident frequencies and incident angles.

As shown in Figure 9, the relative permittivity of rounded-corner triangle is 600. The structure has a resonance frequency of ω0 = 0.12172ωpand the leakage rate is 1τ=1.1771×104, which can be figured out by COMSOL. We use Matlab to set the temporal coupled-mode model as shown in Eq. (16). By comparing with the simulation results of COMSOL at different incident frequencies (near the resonance frequency) and incident angles, we can prove the validity of TCMT.

Figure 9.

(a) Rounded-corner triangle withr=0.15λp,a1=0.3λp,φ=π6,ϵ=600. (b) The real part ofHzfor the eigenmode analysis at the frequencyω0 = 0.12172ωp.

When a TM wave impinges on the scatterer, the scattering cross section can be defined as


where Psctis the rate at which energy is scattered across the circle far away from the scatterer and I0=12μ0ϵ0H02is the intensity of the incident plane wave. For TCMT, the scattering cross section can be calculated by [22]




As shown in Figure 10(b, c), for different incident angles, scattering cross sections predicted by TCMT match well with the results simulated by COMSOL.

Figure 10.

(a) Scattering cross section predicted by TCMT as a function of incident frequencyωand incident angleθ0. Comparison between TCMT and COMSOL simulation results for different incident anglesθ0 = 0 (b) andθ0=π2(c).

As shown in Figure 11, the green line represents the assumption in [22] that the background scattering matrix Bcan be set to Iwhile the blue line represents the assumption that Bis achieved through the simulation results of the arbitrary object covered by PEC. As we can see, when the incident frequency equals to the resonance frequency, both the green line and blue line (TCMT models with different parameters) can match the simulation results. However, when the incident frequency deviates from the resonance frequency, our TCMT model shows a better accuracy compared with the TCMT model in [22], which indicates that the background scattering matrix Bcannot be set to Ieasily when the permittivity of object is high.

Figure 11.

The far-field amplitude of the scattering field with different incident frequencies and angles of (a) ω = ω0 , θ0 = 0, (b)ω=ω0,θ0=π2, (c) ω = ω0 + 2Γ , θ0 = 0, (d)ω=ω0+2Γ,θ0=π2.

5. Fano resonances in periodic structures

Fano resonances have been widely observed in the various periodic structures [24, 25]. The theory of Fano resonance in periodic structures is well developed. In [26], temporal coupled-mode theory is applied to analysis the transmission spectra of photonic crystal slab. According to the TCMT, the Fano resonances existed in transmission spectra are the result of the coupling of leaky mode to the external waves. Recently, the experimental discovery of Fano resonances involving interference between Mie scattering and Bragg scattering is studied in [27]. By comparing the disordered system with the periodic structure, they conclude the sharp resonances in periodic structure are caused by the Bloch waves. In order to study the interference between Mie scattering and Bragg scattering theoretically, the inverse dispersion method is proposed to calculate the photonic band diagram and distinguish unambiguously between Bragg and Mie gaps in the spectra [28]. The method reduces Maxwell’s equations to a problem with the eigenvalue kwhile ω is considered to be a real parameter. It is not so intuitive since conventional approach will reduce the Maxwell’s equations to standard eigenproblem for the frequency [29]. In [30], the author shows that the frequencies of observed Fano resonances existed in a linear array of dielectric cylinders coincide with the position of narrow frequency bands found in the spectra of corresponding two-dimensional photonic crystals. Inspired by [28, 30], we figure out the eigenfrequency of the photonic crystal slab and compare with the band diagram of two-dimensional photonic crystal. We are surprised to find that the occurrence of Fano resonances in photonic crystal slab can be predicted by the band diagram of photonic crystal.

5.1. Transmission spectra of the photonic Crystal slab

The structure of photonic crystal slab is shown in Figure 12. We assume the dielectric cylinders are parallel to the zaxis. When TE waves with different angles incident on the slab, we can calculate the transmission coefficients and plot them in Figure 13.

Figure 12.

Photonic crystal slab with radius of cylindersr = 0.4a(ais the period of the slab), ϵ = 12, μ = 1.

Figure 13.

Transmission coefficient of plane wave incident on the photonic crystal slab. The incident angles are (a) ϕ = 0 and (b)ϕ=π100.

As shown in Figure 13, Fano resonances with narrow resonance width can be observed. The permittivity of photonic crystal slab does not need to be as high as single cylinder in order to achieve same quality factor and such materials may be easily found in nature.

Let us assume the Fano resonances in Figure 13 satisfy the Fano function 11+q2q+ωω0Γ21+ωω0Γ2. Firstly, we use eigenmode analysis in COMSOL to figure out the eigenfrequency of photonic crystal slab. The real and imaginary part of eigenfrequency represent the resonance frequency ω0 and resonance width Γ respectively. Secondly, we use the fitting method in Matlab to get the optimal Fano parameter qin Fano function. Thirdly, with given ω0, Γ, and q, we can plot the Fano function with respect to frequency ω. As shown in Figure 14, the Fano curve matches well with the transmission coefficient simulated by COMSOL. The horizontal ordinate is chosen as a/λfor convenience, which is proportional to frequency ω.

Figure 14.

The Fano curve and the simulation result of photonic crystal slab are plotted.

5.2. Band diagram of photonic Crystal

The photonic crystal slab is periodic in only one direction while two-dimensional photonic crystal is periodic in two directions. For a photonic crystal as shown in Figure 15, the band diagram for ky=0,0kxπais plotted in Figure 16(a). The eigenfrequencies of the photonic crystal are real while the eigenfrequencies of the photonic crystal slab are complex due to the existence of radiation loss. Hence, only the real parts of eigenfrequencies are plotted as shown in Figure 16(b). By comparing the resonance frequencies shown in Figure 13 and eigenfrequencies in Figure 16(b), we can conclude that the occurrence of Fano resonances in transmission spectra of photonic crystal slab can be predicted by the real parts of the eigenfrequencies of the system. In addition, for the Fano resonances, which are observed in Figure 13(b) but cannot be observed in Figure 13(a), they all have the eigenfrequencies with Q→∞. Hence, the resonance widths tend to zero and the resonances cannot be observed.

Figure 15.

Photonic crystal with radius of cylindersr = 0.4a(R1 = R2 = a), ϵ = 12, μ = 1.

Figure 16.

(a) Band diagram of two-dimensional photonic crystal as shown inFigure 15whenky=0,0kxπa. (b) Real parts of eigenfrequencies of photonic crystal slab as shown inFigure 12withQ > 70 are plotted.

As shown in Figure 16, the Fano resonances of the transmission spectra coincide with the band diagram of the two-dimensional photonic crystal, which further explains that Fano resonances in periodic structures can be viewed as the superposition of the Bloch wave, which provides the narrow resonances and the Mie scattering wave which provides the slow varying background.


6. Conclusion

In this chapter, we have presented various structures with high permittivity, which have Fano resonances, such as spheres, core-shell particles, arbitrary shape objects, and periodic structures. For each structure, different theoretical methods together with numerical analysis have been presented. Compared with conventional Fano resonances observed in structures consisted of metal and dielectric, high-permittivity structures can enhance the quality factor significantly, which may open up new opportunities for applications such as sensors, switches, and permittivity measuring technique.

© 2017 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution 3.0 License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

How to cite and reference

Link to this chapter Copy to clipboard

Cite this chapter Copy to clipboard

Xianghong Kong, Lina Qiu and Gaobiao Xiao (November 29th 2017). Fano Resonance in High-Permittivity Objects, Resonance, Jan Awrejcewicz, IntechOpen, DOI: 10.5772/intechopen.70518. Available from:

chapter statistics

1163total chapter downloads

More statistics for editors and authors

Login to your personal dashboard for more detailed statistics on your publications.

Access personal reporting

Related Content

This Book

Next chapter

Laser-Induced Fano Resonance in Condensed Matter Physics

By Ken-ichi Hino, Yohei Watanabe, Nobuya Maeshima and Muneaki Hase

Related Book

First chapter

Numerical Solution of Many-Body Wave Scattering Problem for Small Particles and Creating Materials with Desired Refraction Coefficient

By M. I. Andriychuk and A. G. Ramm

We are IntechOpen, the world's leading publisher of Open Access books. Built by scientists, for scientists. Our readership spans scientists, professors, researchers, librarians, and students, as well as business professionals. We share our knowledge and peer-reveiwed research papers with libraries, scientific and engineering societies, and also work with corporate R&D departments and government entities.

More About Us