The profile of the microfluidic chip in a typical on-chip BSI experiment is shown in
Figure 1. For a schematic showing the full experimental implementation of BSI, the reader is directed to the work of Kammer et al. [
10] and Mulkerns et al. [
19]. In cross-section, the chip consists of three distinct pieces: a base, a cap, and the enclosed fluid. The base, consisting of an infinitely thick flat plane with a semicircular groove of radius
r, has a refractive index of
, and is typically fashioned from glass or polydimethylsiloxane. The channel is formed by sealing the base with a rectangular cap of some thickness
t and will contain the fluid of interest with refractive index
. In the majority of the literature, the same material of refractive index
is used for the cap as for the base; this convention will be echoed here. The chip itself is surrounded by a medium of refractive index
, typically assumed to be air. The light incident on the channel is considered to be a plane wave of wavelength
at some angle
to the surface normal, causing the irradiance at any point on the surface to be constant. As the channel is symmetric about the origin (denoted by a black dot in
Figure 1), the problem can be simplified by considering one side of the channel in the case of normally incident light, but this symmetry is broken under oblique illumination. The general case of non-perpendicular incident light is considered here; however, special attention is paid to the special case of
which is typically used in BSI.
As ray tracing is a deterministic technique, the propagation of a light ray through the system can be fully calculated by geometrical arguments. It should be noted that effects such as diffraction can not be simulated using ray tracing; however, in the case analysed here, the length scales probed are much larger than the wavelength of light used, and as such, geometrical optics is appropriate. First, the angle of each ray at a given interface is elucidated, which can be used to determine the path lengths as well as the relative amplitude and phase of each ray. In addition, the bounds for which a given ray is actually defined can be determined by considering these angles. From this, the final fringe pattern at a given position is calculated by interfering the appropriate rays.
2.1. Angular Bounds
Consider a ray that is incident on the surface at some angle
that intersects the cap at some point
x away from the origin. In this manuscript, left is defined to be negative
x and right positive, with
at the origin of the radius of the channel, as shown in
Figure 1. The incident angle of light is confined to be between
. In the special case that
and
, the reflected ray will travel back along the same path it entered through. However, if
and
(for example, see
Figure 1A), then the ray will intercept the channel circumference at
and some angle
which is described by
However, in general (i.e., when
),
x and
are not the same but are linked by
where
,
t is the thickness of the cap, and the other parameters as defined in
Figure 1. As
is dependent on
, Equation (
2) can be fully solved in terms of
. This form, though less intuitively understandable, is given in the
Supporting Information.
Depending upon the magnitude of
and the incident angle of the light
, the ray will go on to intersect the radius of the channel
i times, with
i henceforth referred to as the intersection number of a given ray. In this case, at very small
x values, the intersection number will be 1, but as
, the number of intersections will approach
∞. In the case of
, the ray will be reflected from the channel surface at an angle
, however in the case that
, the ray will be reflected at an angle
as shown in
Figure 1B. Here
can be found by repeated usage of Snell’s law.
The relationship between the intersection number and the incident position of the light ray can be elucidated by considering the value of
as shown in
Figure 1. Taking the example of the ray in
Figure 1B, as the incident position of the ray moves towards the origin, the intersection point
will move closer to the cap until
and
ceases to exist. For some intersection number
i,
can be determined by summing the angles around the origin:
As
for a given
i, rearranging leads to
which sets the bounds on the angles (and therefore the
x values given a known
value) that cause a ray with
i intersections. Equation (
5) defines a continuous function, but the number of intersections of the ray with the channel can only be a natural number. Therefore, it should be amended to
where
denotes the floor function. Equation (
6) can be seen to give the correct values of
at
(i.e.,
) and
as
(i.e.,
) at normal incidence. The bounds for the first 6 intersection numbers for the case of normal incidence (
) are given in
Table 1.
By considering the right-angled triangle formed by the final intersection point and the origin, the angle of any ray that strikes the underside of the cap will always do so at an angle of
, where
By application of Snell’s law, the angle of the ray exiting the cap of the channel will be
As
is typically higher than
, some of these rays will strike the cap-air boundary at angles greater than the critical angle
(for
,
). By combining Equations (
7) and (
8), for a given intersection number
i the minimum and maximum allowed
values are given implicitly by
By substituting Equation (
4) into Equation (
9) and rearranging, in the special case that
, the range of
values that lead to totally internally reflected rays are given by
Thus, for
,
values between 0.431 and 0.542 (correct to 3 decimal places) will not contribute to the final fringe pattern seen.
2.2. Path Lengths
As the angles and intersection number of a given ray are now defined, the optical path length of the ray between each interface can be evaluated [
20]. By calculating the optical path length of a given ray, the phase difference between the rays can be determined when they are interfered.
For clarity, each incident ray can be split into three distinct sub-rays labelled
a,
b, and
c, as shown in
Figure 2. Rays
a and
b are the rays that reflect from the top and the bottom of the cap respectively. Ray
c is defined to be any ray that intersects with the channel radius, no matter how many times.
Taking first the general case of a
c-type ray with intersection number
i, the optical path length can be split into 6 separate parts, as shown in
Figure 2. The path length of the light from the top surface to the first intersection with the channel radius is given by
and
where all the symbols have the same definitions as previously described. From the first intersection point
to the final intersection point
, the path length is simply the sum of the base lengths of the isosceles triangles between consecutive points. Therefore, the subsequent path length contribution can be written as
In the case that
,
as no paths between intersection points are created. The next part of the path length from the final intersection of the ray with the channel to the cap is given by
In the special case that
,
and so
describes the optical path length of the ray between the only intersection point
and the cap. Lastly, the path inside the cap and from the cap to the imaging plane can be written as
and
respectively, where
d is the distance from the cap-air interface to the imaging plane and
is the final exit angle of the ray from the chip given by
The final path length of ray
c is given by
For rays
a and
b, the formulae are simpler. The optical path lengths for these rays are given by
which for a fixed input angle
are independent of
x. No path lengths between the source and the cap are considered for any ray types as this is the same for all rays due to only rays of a single angle being considered. However, this constraint could be relaxed if needed, such as to facilitate simulation of a focused beam.
2.3. Intensity Calculations
To calculate the fringe pattern seen, all that is left to calculate is the relative amplitude of each ray that arises for a given incident angle and position of the illumination source.
A ray of fully
s- or
p-polarised light incident on a surface at some angle
A and reflected/transmitted at some angle
B will have amplitude reflection and transmission coefficients defined by the Fresnel equations. In terms of
A and
B these are
for
s-polarised light, whereas for
p-polarised light they are
where the polarisation is not specifically stated,
R and
T will be used to denote the polarisation dependent versions
and
. In general,
B can be determined using Snell’s law. It is pertinent to note that
R and
T are here defined in terms of amplitude, not power, to allow inclusion of the phase changes of the ray upon reflection/transmission. In addition, the shorthand of, for example,
is used to denote the reflection coefficient of a ray of angle A travelling in a medium with refractive index
incident on a medium of refractive index
.
It is possible to determine the relative amplitude of the rays from every
value by following the path taken through the system and applying the correct reflection or transmission coefficient. In this paper, the incident intensity of the light is always
, though it would be possible to substitute it for any other factor or indeed one that is a function of
(e.g., a Gaussian beam). For any ray, the first two interfaces will give rise to two reflected components from the interfaces encountered (
a and
b), as well as a transmitted ray (
c) that will continue on, as seen in
Figure 2. In general, a
b-type ray will be approximately
lower in intensity than an
a ray for the materials typically used in on-chip BSI. The rays that continue on will strike the circumference of the channel at an angle
a number of times equal to
i. This leads to an extra factor of
if it is assumed that the material of the cap is the same as the channel material. Lastly, the ray must transmit back through the cap to be detected, leading to the final factors of
. Therefore, for a single
value with a corresponding intersection number
i, a maximum of three different rays will be created with relative amplitudes of
where
. It can immediately be seen that for a given
value,
and
are constant in magnitude, but
varies with
x.
The density of incident rays is equal on the surface of the channel; however, the density of rays reaching a certain location on the viewing plane is not constant. Therefore, equally spaced rays in the imaging plane would not lead to evenly spaced rays in the incident plane, potentially under- or over-sampling rays of certain intersection numbers. To correct for this, a multiplicative factor of
f must be introduced to each
value [
15,
21]. Consider a piece of cap extending from
x to
in width that has incident upon it
N rays of intensity
each, which may be a function of position (here only light that eventually makes it to the detector is considered). As these rays are transmitted and reflected, the extremal rays will be detected at
and
, where
g is some function defining the relation between the input and output
x-coordinates. It is assumed here that the
segment is small enough to not contain any discontinuities and that g is not multivalued within the range
to
(i.e., the intersection number is constant). In that case, the number of rays between
and
must also be
N. If
is sufficiently small,
can be written as
and so assuming conservation of energy density holds
which is clearly incorrect. This shows that a factor of
is required to balance the equation. As the number of rays and area cannot be changed, as a reduction of one would reduce the other, the amplitude of a given ray therefore must be multiplied by a factor
to conserve energy density, where the absolute value is used to ensure intensities cannot be negative.
By considering only the change in the
x position of the light as it propagates through the chip, the full analytical equation for the output position on the detector plane
in terms of the input position
x for a type
c ray is found to be
where
is defined to be the
x position change from first to final intersection point within the channel, which is dependent upon the intersection number of the ray. This term is given by
To determine the final intensity
of the fringe pattern at a particular horizontal position
on the image plane at some distance
d, the rays must be interfered pairwise [
14,
16]. The total intensity at some position
is given by [
17]
where
is the wavelength of the light and
l is the total number of rays that can intersect the detector plane at
. In Equation (
32),
and
are defined to be the
x positions of the incident light that give rise to an output position
for the
j and
k rays. The
and
terms here should each include the factor of
f as set out in Equation (
29) for their respective incident positions if appropriate.
By considering a detection plane some distance
d away from the surface of the chip, evenly spaced points at which to evaluate the intensity can be chosen to simulate the pixels of a camera. From this, the incident position of rays that give rise to components that would all intersect the detector plane at some position
can be determined using Equation (
30). These rays can then be interfered pairwise using Equation (
32) and the interference pattern determined. If not otherwise stated, the values used in the following simulations are
,
,
,
,
,
,
, and
, where
is the width of the detector used here.