CfAI Physics Thesis 1
CfAI Physics Thesis 1
response in chemically-prepared
Gold-Silver Nanostructures
Aishik Panja
A thesis presented for the degree of
Undergraduate Department
Bengaluru, India
June 2021
Unconventional electrical response in
chemically-prepared Gold-Silver
Nanostructures
Aishik Panja
Abstract
prepared Au-Ag NSs using glancing angle vapour deposition and had
we have tried to reproduce the results of Anshu Pandey et. al. and,
common motif. There are also several samples where the resistance had
huge jumps and drops and it even becomes negative at some points.
meters.
i
Acknowledgements
I would like to thank Prof. Ambarish Ghosh for giving me the opportunity
to work in his group. I am also indebted to my seniors - graduate students
Ms. Priyanka Suri, Mr. Eklayva Vashisth and project workers Ms. Paramita
Modak and Mr. Ankur Jaiswar for their valuable guidance and support.
I am also grateful to the Kishore Vaigyanik Protsahan Yojana (KVPY) gran-
ted by the Government of India for my undergraduate stipend that helped
keep my finances running.
ii
Contents
Declaration iv
3 Optical Characterization 38
3.0.1 Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . 38
3.0.2 Measurements . . . . . . . . . . . . . . . . . . . . . . . . . . 41
4 Appendix A 45
5 Appendix B 48
6 Bibliography 53
iii
Declaration
The work in this thesis is based on research carried out at the Centre for Nanos-
cience and Engineering, Indian Institute of Science, Bengaluru, India. No part of
this thesis has been submitted elsewhere for any other degree or qualification, and
it is the work of the author unless referenced to the contrary in the text.
“The copyright of this thesis rests with the author. No quotation from it should be
published without the author’s prior written consent and information derived from
it should be acknowledged”.
iv
Chapter 1
The plasmonics of Au-Ag NSs
result of absorption and scattering of light from nano-particles that are well within
the dipole limit (i.e., have size much lower than the wavelength of incident light)
later chapter). The data is then fit to a Mie Theoretic model that also incorporates
electromagnetic field. When metal NSs are excited with electromagnetic radiation,
typically visible light in the case of gold and silver, the metal electrons are driven to
resonance and this leads to an enhancement in electric field both inside and outside
external electric field E0 . Next, we solve the Laplace Equation in the system. This
1
1. The plasmonics of Au-Ag NSs
has been has been workedf in Appendix A of the author’s Bachelor of Science thesis
[1] as well as several textbooks [3,4] on plasmonics and optics.
that gives the dipole moment and polarizability. Noting[5] that the electric field
(2cosθr̂ + sinθθ̂)
~ =p
due to a dipole of moment p, at ~r and θ is given by E , we have
4πm r3
(ω) − m
p = 4πm a3 E0 (2a)
(ω) + 2m
The polarizability (p/E0 ) is given by :
(ω) − m
α(ω) = 4πm a3 (2b)
(ω) + 2m
The condition for resonance, also called the Frolich Condition, is given by Re((ω)) =
−2m . This is because the medium is a dielectric (Im(m ) = 0), implying Re((ω)+
2
1. The plasmonics of Au-Ag NSs
cuvette in case of a solution) that is irradiated with light of intensity Ii (ω) only to
be attenuated to It (ω). Assuming that the light does not suffer multiple scattering
sample in the cuvette made of quartz which has near 100% transmission.
Also, It /Ii is the transmittance. The absorbance is given as −log10 (It /Ii ). There-
Here, we fit the above equation to our absorbance data with several parameters
√
- fAu - the volume fraction of gold, N - the number density, nmed (= m ) - the
We use the Maxwell Garnet Method[8,9] to calculate (ω) from fAu and fAg (=
1 − fAu ).
3
1. The plasmonics of Au-Ag NSs
dispersed. The effective dielectric constant ef f is given by solving the following
equation :
ef f − h i − h
. = Σfi (6)
ef f + 2h i + 2h
Thus for a two component system (gold and silver), we have
ef f − Ag Au − Ag
. = fAu (7a)
ef f + 2Ag Au + 2Ag
2fAu (Au − Ag )
Ag
1 +
Au + 2Ag
To be noted is that MGM is valid only at extremely low volume fractions (< 10−5 )
such that the inclusions are spatially separated from the matrix and do not interact
as it can account for small volume fractions, places the matrix and the inclusions
erential directions. However, it has been shown in this author’s Bachelor’s Thesis
that the ef f calculated from the two methods is second order in f , where f is the
derived that all systems - irrespective of the shape and orientation of the inclusions
silver. There are several sources at our perusal - we use the emperical data from
4
1. The plasmonics of Au-Ag NSs
Ciesielski et. al.[13,14] . We have departed from the Johnson and Christy data[15]
that is conventionally used in these cases as it pertains to the bulk metal whereas
the former describes thin films grown on SiO2 subtstrates. The chart is provided in
the appendix. The of gold and silver are plotted below as a function of wavelength.
Figure 1.2: The dielectric constant of silver in the 200 nm (0.2µm) - 900nm
range. From Ciesielski et. al.
Figure 1.3: The dielectric constant of gold from Ciesielski et. al.
5
1. The plasmonics of Au-Ag NSs
position (more details in Chapter 3). from plasmonics, the plasmon resonances
of gold and silver are at 2.3 eV and 3.2 eV respectively; therefore, we expect the
composite to have a peak in both absorption and extinction to lie in this range.
• The extinction has no sharp peak anywhere - only a broad hump-like curve
However, in our experiment, we did not notice any such unusual feature. In fact,
our measured absorption matched with theory in the visible regime to a great
extent.
6
Chapter 2
Preparation and Electrical
Characterization
As discussed in the author’s bachelor’s thesis [1] , samples were thermally deposited
However, in this project, we have prepared our samples using the wet chemical
1. Chemicals Required
2. Procedure
Aqueous solutions of 0.1M CTAB (10 mL), 0.1M NaOH(0.1 ml), 1 M ammonium
bromide(1 mL), 0.1 M potassium iodide (40 µL) and 1 mM silver nitrate(5 mL)
7
2.0.1. Sample Preparation
were mixed. The mixture was continuously stirred at 1000 rpm till slight cloud-
Then 800 µL of 1mM aqueous solution of HAuCl4 was added to the silver nan-
change from golden brown/ maroon to greyish black. It is then centrifuged for 3-4
min at 6500 revolutions per minute (RPM). After centrifugation, the supernatant
Note that the amount of each reagent may be varied at each run; in later runs, we
scale it by 2.
make a saturated solution of Potassium Hydroxide in IPA. The drop casting drill
is the following - Drop cast a small quantity (2µl to 10µl, depending on the desired
size of the spot) of the chloroform solution → add the KOH/IPA solution to it
and wait for 1 minute → rinse thoroughly with Methanol. This process is to be
The drop cast is carried out inside a makeshift plastic container with sliding walls
that is covered in yellow sheet to reduce exposure to ambient light and its top
the later samples are covered with the epoxy AralditeT M to protect the film from
the elements while they are removed from the box to be prepared for electrical
measurements.
Note that we keep on making changes to the reaction conditions such as quantity
Substrates
8
2.0.2. Electrical Measurements
We use quartz substrates that has electrodes/pads made of chromium and gold for
the electrical connections. The pads are fabricated at the Annexe of the National
optical windows are located. Inside this chamber, there is a smaller liquid nitrogen
chamber - one end of that is attached to the lid of the main chamber and the other
end attached to a copper stub that fits into the aforementioned protrusion. A prin-
ted circuit board (PCB) is attached to this stub that has several electrodes which
are connected to the substrate pads via copper wires. These wires are tethered
to the substrate with indium wires. There is a Pt-100 thermometer and a heater
We have a resistance box made of aluminium that has several resistors of resist-
ances between 10 Ohm and 106 Ohm connected between Bayonet Neill–Concelman
(BNC) terminals. We send current and voltage to the PCB, and therefore our
sample, using coaxial cables via a Military connector. The other end of these
(call it the terminal box) that has BNC connectors on the other side.
We have two Stanford Research SR 830 DSP lock-in amplifiers. One of them (call
this Lock-in 1) supplies a 1 V reference signal of frequency 227.7 Hertz across the
10,000 ohm resistor (also called load resistance) on the resistance box and the res-
the terminal box. The I- terminal is grounded using a 50 ohm BNC terminator. For
adapters and BNC cables to Lock-in 1. Similarly, for the four-probe resistance, we
9
2.0.2. Electrical Measurements
connect the V+ and V- terminals to the other lock-in (Lock-in 2). We connect the
the circuit.
The values of voltage displayed on the Channel 1 of the lock-ins are divided by the
The cables from the heater and thermometer also pass through the military con-
The lock-ins and the temperature controller are operated on LabVIEW via a GPIB
interface. The lock-ins are interconnected using an IEEE 488 GPIB which also con-
nects Lock-in 2 with the the temperature controller. The Lock-in 1 is connected to
of the experiment, the phase of both the lock-ins as displayed in Channel 2 are set
We attach the substrate to the stub with strips of carbon tape. The entire stub,
After loading the sample, the refrigerator is pumped using a Pfeiffer HiCube 80
pump for a few hours at the maximum frequency of 1500 Hz. Subsequently, we set
from the top. After the desired temperature is reached, we turn on the heater to
10
2.0.2. Electrical Measurements
a)
b)
Figure 2.1: a) The stub with carbon tapes stuck. Also seen is the PCB and
the copper wires attached to it.
b) The upper part of the cryogenic chamber resting on the lower part. It
contains i) the LN2 chamber (shown wrapped in super-insulation), ii) the
tubes on the top of the lid for electrical connections and an inlet for liquid
nitrogen (white funnel) and iii) the stub attached to the LN2 chamber.
Results Sample 1
This was our very first attempt at making the sample according to the above
protocol. Earlier, we had tried another protocol (Protocol 1[16,17,19,20] which had
yielded our desired Ag-Au solution (as confirmed by spectroscopy) but the drop-
11
2.0.2. Electrical Measurements
batim, except for the amount of HAuCl4 (2.5 ml instead of 800 µl). The drop-cast
was performed on a hot plate that allows the liquids (IPA, methanol and chloro-
form) However, the sample was loaded into the vacuum refrigerator 13 days after it
was prepared. During this time, it was often kept undesiccated and understandably
got oxidised.
Sample 2
Sample 2 was the first sample that was put for electrical measurement right after
it was prepared. During preparation, the solution was heated at 600 C while being
stirred. The amount of the constituents was as mentioned in Protocol 3 except for
HAuCl4 - we added 2 ml of HAuCl4 instead of the prescribed 800 µl. We had also
The substrate has six parallel pads that are labelled 1 through 6 sequentially. The
drop-cast was carried out on a hot plate that was heated to 60-70 0 C to accelerate
the drying process. The dropcast film covers the four pads numbered 2 through
12
2.0.2. Electrical Measurements
a b
c d
Figure 2.2: (a) and (b) are respectively the 2 probe and 4 probe resistance values
as a function of temperature over the entire cycle encompassing the cool-down and
heat-up processes(c) and (d) are the 2-probe and 4-probe values plotted together
in the cool-down and heat-up processes, respectively .
It is evident that none of these processes point to a typical metal that follows a
linear R vs T profile with a positive slope. Moreover, the heat-up and cool-down
profiles are totally different and they do not even converge at room temperature,
implying a permanent change in the film structure over a thermal cycle. In the
cool-down process, the 4-probe and 2-probe resistances follow a similar profile,
sets in.
13
2.0.2. Electrical Measurements
a major portion of the pellets ended up getting stuck to the walls of the tubes and
could not be dissolved in chloroform. Despite this, we had drop-cast the solution
and were able to get electrical continuity between the electrodes numbered 3 and
a)
b)
Figure 2.3: The a) film after it was unloaded from the refrigerator and the b)
2 probe resistance (R34 ) over one thermal cycle (as a function of temperature).
The 2-probe resistance decreases with decreasing temperature till ≈ 225 K, after
which it starts [Link] is quite similar to the Sample 2 profile. During heat-
up, the resistance continuously increases until it reaches 250 K when it flattens out
and eventually starts decreasing. There is a sharp drop in the 300 K - 320 K range
and when cooled back to the room temperature (297 K), it ends up at a much
lower resistance than where it started. Thus, at room temperature, we have three
heat-up (≈ 51.2Ω) and the final resistance (≈ 47.7Ω). In hindsight, these values
may not seem too far off - these are indeed only 3% off each other. However, these
points lie at the extremities of the R vs T graph and are therefore significant.
14
2.0.2. Electrical Measurements
This time, the drop-cast was performed inside a box that had a hairdryer attached
a) b)
Figure 2.4: The a) Sample 4 after it was unloaded from the refrigerator and
the b) The chamber with the hairdryer attached at the top.
a) b)
Figure 2.5: The a)The 2-probe (R25 ) and b)4-probe resistance (R2534 ) of
Sample 4 as a function of temperature over 3 thermal cycles.
The two-probe resistance follows linear (metallic) curves that also show hysteresis.
We confirm this by carrying out 3 thermal cycles. However, the four-probe resist-
15
2.0.2. Electrical Measurements
ance was very similar in nature to the one seen in Sample 2 in that the cool-down
and heat-up processes are too far off with the resistance at room temperature being
treated using a Harrick PlasmaBasic Plasma Cleaner and drop-cast carried inside
Figure 2.6: The Sample 5 thin film. Photograph taken after electrical meas-
urements.
16
2.0.2. Electrical Measurements
a) b)
Figure 2.7: The a)The 2-probe (R25 ) and b)4-probe resistance (R2534 ) of
Sample 5 as a function of temperature over 3 thermal cycles.
The 2-probe resistance profile, though not strictly linear, can be approximated to
that of a metal. The 4-probe has extreme hysteresis - in the cool-down process,
the resistance initially decreases with temperature and starts to increase later. The
A common feature in Samples 4 and 5 that we did not see in the previous samples
is that the 2-probe and 4-probe profiles are totally different in nature. In other
words, the sample is (almost) metal-like in the 2-probe length scale (≈ 1 cm), but
scale across the sample. The substrate had eight concentric pads that, at the cen-
scribed 0.8 ml, while other components were added as per the protocol. We did
17
2.0.2. Electrical Measurements
a) b)
Figure 2.8: The a) The drop cast thin film in Sample 6, photograph taken
after electrical measurements b) The final solution of Sample 6, identifiable
by its golden maroon color.
a) b)
Figure 2.9: The a) The 2-probe and b) 4-probe resistance of Sample 6. The
current was sent through the outer pads and voltage measured over the two
middle pads (marked with indium wires in the previous figure).
Here, both the 2-probe and 4-probe profiles are linear except for temperatures
above 300 K in the heat-up process where the resistance drops sharply. We had
18
2.0.2. Electrical Measurements
strate that has six pads. During the drop-cast, we vented the chamber with Nitro-
a) b)
Figure 2.10: The a) The drop cast thin film in Sample 7, photograph taken
after electrical measurements, the carbon tape and indium wires still stuck
to the substrate b) The final solution of Sample 7.
Clearly, the film is not smooth at all despite the substrate being plasma cleaned.
a) b)
Figure 2.11: The a) The 2-probe and b) 4-probe resistance of Sample 7. For
the 4-probe resistance, we send current through pads numbered 2 and 6 and
measured voltage between 3 and 4.
The two probe was linear apart from a few blips. The four-probe had two linear
19
2.0.2. Electrical Measurements
regions with different slopes that exhibited hypothesis. To check how the differ-
ence in resistance between the cool-down and heat-up processes (due to hysteresis)
depends on temperature, we took the sample through cd3 to 229K (yellow curve)
any given point and start heating up, the two curves start diverging.
7. However, probably due to over-etching caused by IPA and methanol, the film is
not very smooth and the gold pads tend to peel off. Therefore, we only manage to
a)
b)
Figure 2.12: The a) The drop cast thin film in Sample 8, photograph taken
after electrical measurements, the pads looking very worn out b) The final
solution of Sample 8.
20
2.0.2. Electrical Measurements
21
2.0.2. Electrical Measurements
a) b)
Figure 2.14: The a) The drop cast thin film in Sample 9, photograph taken
after electrical measurements, b) The final solution of Sample 9.
b)
a)
Figure 2.15: The a) The 2-probe and b) 4-probe resistance of Sample 9. For
the 4-probe resistance, we send current through pads numbered 2 and 6 and
measured voltage between 3 and 5.
heated the sample up, just to see the deviation (CD2 and HU2 in the figures). We
22
2.0.2. Electrical Measurements
also observe sharp drops and jumps at several temperatures, including at temper-
10:8 (solution:IPA) ratio before centrifuging it. After drop-cast, the film is covered
a)
b)
Figure 2.16: The a) The drop cast thin film in Sample 10, photograph taken
after electrical measurements, the pads looking very worn out b) The final
solution of Sample 10.
a) b)
Figure 2.17: The a) The 2-probe and b) 4-probe resistance of Sample 10. For
the 4-probe resistance, we send current through pads numbered 1 and 5 and
measured voltage between 2 and 3.
23
2.0.2. Electrical Measurements
a chamber that was vented with gaseous nitrogen and the film was covered with
Araldite.
a) b)
Figure 2.18: The a) The drop cast thin film in Sample 11, photograph taken
after electrical measurements, b) The final solution of Sample 11.
a) b)
Figure 2.19: The a) The 2-probe and b) 4-probe resistance of Sample 11.
During the heat-up process, the voltage leads must have got disconnected -
hence the 4-probe was overload.
24
2.0.2. Electrical Measurements
a) b)
Figure 2.20: The a) The drop cast thin film in Sample 12, photograph taken
after electrical measurements, b) The final solution of Sample 12.
a) b)
25
2.0.2. Electrical Measurements
a) b)
a)
b)
Figure 2.23: The a) The drop cast thin film in Sample 14, photograph taken
after electrical measurements, b) The final solution of Sample 14.
26
2.0.2. Electrical Measurements
a) b)
Figure 2.24: The a) The 2-probe and b) 4-probe resistance of Sample 14.
The R vs T curves for Sample are mostly linear in both the 2-probe and 4-probe.
Moreover, the heat-up and cool-down paths are very close together - no sign of
b)
a)
Figure 2.25: The a) The drop cast thin film in Sample 15, photograph taken
after electrical measurements, b) The final solution of Sample 15.
27
2.0.2. Electrical Measurements
a)
b)
Figure 2.26: The a) The 2-probe and b) 4-probe resistance of Sample 15.
Parts of the graph are zoomed in the inset.
In Sample 15, there is a huge spike at temperatures below 150 K that can be
explained as follows - the two phases of the lock-in amplifiers showing the current
and voltage are set to zero. However, at the onset of these sharp spikes, the phase
28
2.0.2. Electrical Measurements
of the voltage lock-in (that measures the 4-probe resistance) departs from zero
and mostly hovered around 80 degrees. As a result the input current (generated
from the reference voltage sent across the 10kΩ resistor) and the output voltage are
always out of phase leading to the resistance taking these unusual values. However,
as seen in the inset, the sample behaves ’normally’ in the 300K-150K - resistance
We are unsure why this happens. Physically this means the input current and
a) b)
Figure 2.27: The a) The drop cast thin film in Sample 16, photograph taken
after electrical measurements, b) The final solution of Sample 16.
29
2.0.2. Electrical Measurements
a ) R (2 probe)
b ) R (4 probe)
Figure 2.28: The a) The 2-probe and b) 4-probe resistance of Sample 16.
Here, the phase difference between the two lock-ins departs from zero at certain
points in the cycle and it is only at those points that the resistance shoots. The
30
2.0.2. Electrical Measurements
a) b)
Figure 2.29: The a) The drop cast thin film in Sample 17, photograph taken
after electrical measurements, b) The final solution of Sample 17.
a ) R (2 probe)
b ) R (4 probe)
Figure 2.30: The a) The 2-probe and b) 4-probe resistance of Sample 17.
31
2.0.2. Electrical Measurements
Sample 17 This time, the phase difference between the signal and reference
exceeds 90 degrees. As a result, the resistance becomes negative (as the time-
Sample 18 Sample 18 was prepared according to the protocol and the re-
agents were scaled by a factor of 2. The drop-cast film was covered with araldite.
a)
b)
Figure 2.31: The a) The drop cast thin film in Sample 18, photograph taken
after electrical measurements, b) The final solution of Sample 18.
32
2.0.2. Electrical Measurements
a ) R (2 probe)
b ) R (4 probe)
Figure 2.32: The a) The 2-probe and b) 4-probe resistance of Sample 18.
As seen in samples 15, 16 and 17, there is a phase mismatch in Sample 18. The
phase difference between the input and signal oscillates between 0 and 180 degrees.
As a result the recorded ’resistance’ is all over the place, even taking negative values
Sample 19 Sample 19
33
2.0.2. Electrical Measurements
a)
b)
Figure 2.33: The a) The drop cast thin film in Sample 19, photograph taken
after electrical measurements, b) The final solution of Sample 19.
a ) R (2 probe)
b ) R (4 probe)
Figure 2.34: The a) The 2-probe and b) 4-probe resistance of Sample 19.
34
2.0.2. Electrical Measurements
Sample 19 does not show the unusual and abrupt phase difference between the ref-
erence and the signal that is accompanied by the resistance taking unconventional
values. This behaves more like Samples 2, 3, 4 etc. where the common theme is an
U or a V shaped R vs T curve.
a) b)
Figure 2.35: The a) The drop cast thin film in Sample 20, photograph taken
after electrical measurements, b) The final solution of Sample 20.
35
2.0.2. Electrical Measurements
a ) R (2 probe)
b ) R (2 probe)
Figure 2.36: The a) The 2-probe resistance of Sample 20. Fig b) is the
zoomed-in section of Fig a) in the 250-300 K range that captures the R vs T
behaviour before the sharp drop. During CD, The resistance decreases and
subsequently increases (a staggered V-shaped curve) before the drop. In the
HU process, the phase difference between the signal and reference stays above
90 degrees hence the sharp vertical lines.
As seen in samples 15-18, the phase difference between the output and reference
signal deviates from zero and takes values greater than 90 degrees (between 170
Samples 15, 16, 17, 18 and 20 are unusual in that the signal and reference fall
out of phase at several points in the experiment; as a result the resistance takes
36
2.0.2. Electrical Measurements
Summary Most samples depart from the conventional metallic nature where
the R vs T is a linear curve with a positive slope. The notable features in the
• During cooldown, in several samples, the resistance first decreases with tem-
perature and later starts to increase. In many cases the upward (increasing)
• For most samples, the R vs T curves for the heat-up and cooldown process
• In some samples, there are sharp drops in resistance during heat-up in the 300
(at the Nanoparticle level) as a result of continuous thermal cycles, the fact
• In Samples 15-20 (with the exception of Sample 19), the phase of the input
becomes negative. These phase changes are abrupt- the sample behaves
37
Chapter 3
Optical Characterization
3.0.1 Introduction
As elaborated in Chapter 1, Anshu Pandey et. al. had seen an unconventional op-
tical response in the samples that had superconducting transitions[18] . The absorb-
ance had been near-zero throughout the spectrum, whereas the extinction spectrum
Both these features are unusual. It was shown in this author’s Bachelors’ Thesis
[1] (pages 60-64) that this corresponds to unnatural values in the Drude Lorentz fit
[3,21,22] .
The first task is to extract the from the optical data. We use the equations (Equa-
tions 3a,b) relating the extinction and absorption coefficients with the (ω), med
and the size and number density of the particles. As the absorbance is close to zero
, (ω) is a real number. Solving for (ω) from the expression of σext and plotting
38
3.0.1. Introduction
Figure 3.1: The effective dielectric constant of the Au-Ag NSs calculated
using the parameters a = 7.5 nm and N = 1015 .
itions into account. As these IB transitions occur at a much higher frequency than
the UV-Vis-NIR regime, they are all clubbed together in the ∞ [Link] relax-
term that is the natural frequency of oscillation of the bound particles. Usually,
this term takes into account the lattice vibrations in non-metals with high-energy
39
3.0.1. Introduction
Figure 3.2: A plot of the effective dielectric constant against energy for the
Au-Ag NSs. The fit data is also plotted. We take the upward branch as that
is what fits with the Drude-Lorentz Model. Code presented in Appendix B.
∞ 9 5 6.7 3.01
ω0 0 0 0.098 0.05,0.08
MgO is modelled a double oscillator which explains the two sets of frequencies in
the table.
Clearly, the values are unphysical, especially the fact that ∞ is negative which
40
3.0.2. Measurements
3.0.2 Measurements
We measure the absorbance of the final solution as the drop-cast film is balckish in
colour and the absorbance would be very high anyway; thereby not providing any
additional information.
We use a simple, purpose-built setup for our measurements. We have a light source
(Ocean Optics HL 2000 or Ocean Optics D 2000) that is connected via an optical
fibre cable to a cuvette holder. Another fibre cable connects the cuvette holder to
We set an appropriate integration time so that the spectrometer does not get sat-
urated. To reduce the noise, we increase the number of scans per measurement.
We then set the baseline with the lamp switched off, thereby eliminating the ef-
fects of ambient light. If we find the data particularly noisy, we do a Fast Fourier
Next, we run a code that takes the data files and fits (sum of least squares fit)
the data to our plasmonic model described in Chapter 1. In most of the fitting
codes, we have considered arrays of size > 100 for each fitting parameter. The
fAu array goes from 0 to 1, the nmed from 1 to 2. The number density range is
chosen rather arbitrarily, using the trial-and-error method, and is usually in the
[1015 , 1018 ] range. The particle size, is usually varied from 10 nm to 20 nm but is
not included in the fitting loop. This is mostly because the particles have a wide
range of size in the solution (also confirmed from TEM images) and it is convenient
to stick with a value and change it before running the code if the situation arises.
41
3.0.2. Measurements
a b
c d
Figure 3.3: (a), (b) , (c) abd (d) are respectively the Abasorbance vs Wavelength
and the least sum of squares fit curves for Samples 3,4,6 and 8. The data was
taken using the Halogen-tungsten light (Ocean Optics Hl 2000) and the Ocean
Optics USB 4000 spectrometer. Both these devices work in the 350-900 nm range.
We have removed the noise using FFT Smoothing on OriginPro. The fit parameters
are given in the text boxes. To be noted is that the particle size is 15 mm in all
these cases unless otherwise mentioned.
Clearly, Samples 3,4 and 6 have a decent agreement with the theoretical data,In all
these samples, the fraction of gold (fAu ) roughly lies between 10 and 15 percent. It
would be seen in the later samples that fAu is in the 10%-20% range. TEM images
of some other samples also confirm that Ag is the major constituent. The ratios
obtained from TEM does not exactly match with the fit parameters but the trend
is well exhibited.
42
3.0.2. Measurements
In Sample 8, the tail does not go to zero, something that cannot be explained by
this theoretical model. We see this behaviour in several other samples as well.
a b
c d
Figure 3.4: (a), (b) , (c) abd (d) are respectively the Absorbance vs Wavelength
for Samples 12,14,15 and 19. The light source was a Halogen-tungsten lamp and
the spectrometer was Ocean Optics Flame.A particle size of 15 mm in all these
cases unless otherwise mentioned.
In other words, lsqs = min(Σi (gi − fi )2 ), where fi ’s are the absorbance values as i
goes over all the data [Link] ’s are the data points in the fitting curve. For every
set of the fitting parameters, we have a unique set of gi ’s. ‘lsqsr’, which is defined
lsqs
as lsqsr = gives the accuracy of fit.
Σi fi2
43
3.0.2. Measurements
Summary
• We analyse 8 samples here and all but two have their lsqsr values less than
amount of gold salt added in the reaction mixture and the peak position. We
attribute this to the fact that the constitution of the nanoparticles may not
44
Chapter 4
Appendix A
(m ). A constant electric field E0 k̂ is switched on. From the Laplace Equation
(∇2 φ = 0) the potential due to this field φ(r, θ) takes the following forms inside
Al rl Pl (cosθ)
P
Φin = l (A1b)
where Pl are Legendre Polynomials of order l(≥ 0) and Al , Bl and Cl are constants.
To be noted is that the ‘−(l + 1)’ is term is absent in Φin as it leads to a singularity
at r = 0.
At large r, E = E0 , implying φ(r → ∞)=-E0 rcosθ. This leads to, from equation
(A1a)
and the fact that P1 (cosθ) = cosθ,
Bl = 0 for l 6= 1
(A2)
B1 = −E0
45
4. Appendix A
Now, from the equality of tangential components of the electric field, i.e.
δφin δφout
δθ r=a = δθ , one has the following :
r=a
P∞ dPl P∞ Cl dPl
l l
l=0 Al a = l=0 Bl a + l+1 which leads to
dθ a dθ
Cl
Al = for l 6= 1
a2l+1
C1
A1 = − E0 (A3)
a3
δφin δφout
From the equality of radial components, i.e. (ω) δr r=a = m δr one
r=a
has :
P∞ P∞ Cl (l + 1)
(ω) l−1 P l−1 −
l=0 Al la = m l=0 Bl la
l
Pl .
al+2
Using
equation (A2), this leads to :
m Cl (l + 1)
A l = − for l 6= 1
(ω) la2l+1
m 2C1
A = − E0 + (A4)
1
3
(ω) a
Equating the corresponding parts of equations (A3) and (A4), one eliminates Al
and
gets the following :
Cl = 0 for l6= 1
(ω) − m
C1 = a3 E0 (A5)
(ω) + 2m
Using (A5) one finds the values of Al :
Al = 0 for l 6= 1
(ω) − m
A1 = −E0 (A6)
(ω) + 2m
Combining Equations (A6) and (A5) with Equation (A1), one has :
(ω) − m E0 cosθ
Φout = −E0 rcosθ + a3 (A7a)
(ω) + 2m r2
3m
Φin = − E0 rcosθ (A7b)
(ω) + 2m
46
4. Appendix A
3m
Clearly, E~in = E0 (Equation (2)) and E~out can be also be written as
(ω) + 2m
(ω) − m a3
E~out = E0 (−cosθr̂ + sinθθ̂) − E0 (2cosθr̂ + sinθθ̂)
(ω) + 2m r3
The first part of the above equation is the external field and the second part is due
to the dipole induced on the metallic particle. The dipole moment is given by
(ω) − m
p = 4πm a3 E0 (A9a)
(ω) + 2m
and the polarizability by
(ω) − m
α(ω) = 4πm a3 (A9b)
(ω) + 2m
When extinction spectroscopic measurements are performed on the sample, one ob-
This is also known as the Frolich’s Condition. The scattering and absorption cross-
sections which together constitute the extinction cross-section are respectively given
k 4 |α|2
by σscat (k) = and σabs (k) = kImag(α) (A11a,b)
6π
where k is the wave vector.
47
Chapter 5
Appendix B
nm = 10ˆ(-9);
r = 15*nm/2;
fclose(samplef);
counts = length(lambdas)
abssnew = [];
%silver
agf = fopen(’[Link]’);
fclose(agf);
ns = p2; ks = p3;
48
5. Appendix B
%gold
auf = fopen(’[Link]’);
fclose(auf);
ng = q2; kg = q3;
count = length(p1)
lambdanm = p1*1000 ;
lambda = p1/106 ;
eg = erg + 1j ∗ eig;
es = ers + 1j ∗ eis;
nm = 10( − 9);
%%%%%%%%%%%%%
sp = 1;
sp = 1;
else
for i = 2:count
if lambdanm(i-1)<s1(1) lambdanm(i)>s1(1)
sp = i;
end
end
end
49
5. Appendix B
for i = sp:count
lambdanmnew(i-sp+1) = lambdanm(i);
end
countnew = length(lambdanmnew)
S = 0;
for i = 1:countnew
for j = 2:counts
abssnew(i) = (s2(j-1)+s2(j))/2;
S = S+abssnew(i)2 ;
end
end
end
length(abssnew)
plot(lambdanmnew,abssnew);hold on
absor = zeros(length(f),length(nmed),length(N),length(abssnew));
diff = zeros(length(f),length(nmed),length(N),length(abssnew));
sumsqdiff=zeros(length(f),length(nmed),length(N));
for j=1:length(f)
50
5. Appendix B
for k = 1:length(nmed)
for m = 1:length(N)
for i = 1:length(abssnew)
(f (j) + 2) ∗ esnew(i)) −
end
end
end
end
for j=1:length(f)
for k = 1:length(nmed)
for m = 1:length(N)
for i = 1:countnew
end
end
end
end
for j=1:length(f)
for k = 1:length(nmed)
for m = 1:length(N)
51
5. Appendix B
end
end
end
end
fmin
nmedmin
Nmin
sumsqdiffmin
sumsqdiffminratio
Absorfinal = [];
for j=1:length(f)
for k = 1:length(nmed)
for m = 1:length(N)
for i = 1:countnew
Absorfinal(i) = absor(jmin,kmin,mmin,i);
end
end
end
end
plot(lambdanmnew,Absorfinal);
52
Chapter 6
Bibliography
6. Beer (1852). "Bestimmung der Absorption des rothen Lichts in farbigen Flüssigkeiten"
[Determination of the absorption of red light in colored liquids]. Annalen der Physik
et umbrae [Photometry, or, On the measure and gradations of light intensity, col-
8, Maxwell Garnett, J. C., "Colours in metal glasses and in metal films," Transac-
9. Sihvola, A., "Two Main Avenues Leading To the Maxwell Garnett Mixing Rule",
53
6. Bibliography
der Mischkörper aus isotropen Substanzen". Annalen der Physik. 416 (7): 636–664.
11. Kreibig, U., Vollmer, M., Optical Properties of Metal Clusters, Springer (1995)
12. Serge, B., "Anisotropic effective medium theories", Journal de Physique I 4(2),
DOI: 10.1051/jp1:1994139
13. Ciesielski, A., Skowronski, L, Trzinski, M., Szoplik, T., Controlling the optical
parameters of self-assembled silver films with wetting layers and annealing, Appl.
14. Ciesielski, A., Skowronski, L, Trzinski, M., Szoplik, T., Evidence of germanium
15. P. B. Johnson and R. W. Christy, Optical Constants of the Noble Metals, Phys.
16. Thapa, D.K., et. al., Coexistence of Diamagnetism and Vanishingly Small
arXiv:1807.08572 (2019)
arXiv:1912.05428
electrical resistance of Ag-Au thin films near room temperature, arXiv preprint
arXiv:1906.00708
21. Ashcroft, N.W., Mermin, N.D., Solid State Physics, Harcourt College Publish-
ing, 1976
54
6. Bibliography
22. Shalaev, V., Cai, W., "Optical Metamaterials- Fundamentals and Applica-
55