WEBVTT

00:00:05.550 --> 00:00:07.950
HA DANG: Today, I am
going to use Mathematica

00:00:07.950 --> 00:00:11.640
to visualize a nanoparticle
polymer network

00:00:11.640 --> 00:00:14.450
and to investigate how
the polymer length affects

00:00:14.450 --> 00:00:18.360
of probability that polymer
would cross the nanoparticles

00:00:18.360 --> 00:00:20.400
and how the polymer
length affects

00:00:20.400 --> 00:00:23.000
the nanoparticle polymer
network formation.

00:00:25.590 --> 00:00:29.400
Polymers and macromolecules
composed of hundred to thousand

00:00:29.400 --> 00:00:30.540
subunits.

00:00:30.540 --> 00:00:33.660
They can adopt many
different configuration,

00:00:33.660 --> 00:00:37.170
depending on the identity
of their subunits.

00:00:37.170 --> 00:00:39.135
Common everyday
examples of polymers

00:00:39.135 --> 00:00:44.400
are PTC, polystyrene,
Teflon, and so on.

00:00:44.400 --> 00:00:46.650
When particles and
polymers of are

00:00:46.650 --> 00:00:50.680
favorable interaction are
mixed, it is not always apparent

00:00:50.680 --> 00:00:53.410
whether they will form
a cross-linked network.

00:00:53.410 --> 00:00:56.100
If the polymer
configuration makes the end

00:00:56.100 --> 00:00:58.290
to end distance of
the polymer shorter

00:00:58.290 --> 00:01:00.630
than the distance
between two particles,

00:01:00.630 --> 00:01:03.870
no polymer-particle
network will form.

00:01:03.870 --> 00:01:07.230
This project will focus on
visualizing the interaction

00:01:07.230 --> 00:01:10.980
between particles and polymers
given polymer's length,

00:01:10.980 --> 00:01:12.130
or the number of subunits.

00:01:14.960 --> 00:01:19.060
In stage 1, we make polymer
by use the random walk method.

00:01:24.110 --> 00:01:27.470
We consider a particle
as many lattice sites.

00:01:27.470 --> 00:01:32.170
Polymer with start a
random lattice site

00:01:32.170 --> 00:01:36.490
This initiation part just
ensure that the polymer will

00:01:36.490 --> 00:01:38.830
truly random walk at
a certain distance

00:01:38.830 --> 00:01:43.110
away from the surface
of the particle.

00:01:43.110 --> 00:01:46.650
In the random walk
simulaton, the step size

00:01:46.650 --> 00:01:50.100
is the length of the
subunit of the polymer,

00:01:50.100 --> 00:01:56.590
which is about 10 times smaller
than the size of the particles.

00:01:56.590 --> 00:01:59.100
The step angle tell
you the direction

00:01:59.100 --> 00:02:01.920
that the poem will walk.

00:02:01.920 --> 00:02:05.460
The step tell you how
much the polymer will walk

00:02:05.460 --> 00:02:08.650
in the x and the y direction.

00:02:08.650 --> 00:02:12.930
The next position of the polymer
will be the current position

00:02:12.930 --> 00:02:16.350
of the polymer plus step.

00:02:16.350 --> 00:02:20.340
Each time the new polymer
position is computed,

00:02:20.340 --> 00:02:26.030
it is stored inside on
list which we call path.

00:02:26.030 --> 00:02:28.010
Once we have the
polymer path, we

00:02:28.010 --> 00:02:31.970
can visualize it inside the
particle matrix using Graphic.

00:02:38.020 --> 00:02:41.890
In stage 2, we want to
generate many polymers

00:02:41.890 --> 00:02:44.500
and visualize how they
distribute and interact

00:02:44.500 --> 00:02:48.620
with particles inside a
lattice of many particles.

00:02:48.620 --> 00:02:52.270
And in order to do so,
we combine all the steps

00:02:52.270 --> 00:02:56.590
in state 1 into one
function, polymerRandomwalk,

00:02:56.590 --> 00:02:59.020
which takes in
the polymer length

00:02:59.020 --> 00:03:02.890
and the lattice of particles
and returns a polymer path.

00:03:06.450 --> 00:03:11.490
We will generate 100 random
random walk polymers of 1,000

00:03:11.490 --> 00:03:14.070
subunit by using Table.

00:03:14.070 --> 00:03:17.385
We visualize the polymers and
particle matrix using Graphic.

00:03:25.510 --> 00:03:30.370
In stage 3, we want to calculate
the probability of crossover

00:03:30.370 --> 00:03:33.310
between nanoparticle
in polymers by counting

00:03:33.310 --> 00:03:37.390
how many polymers that wrap
around more than two particles.

00:03:37.390 --> 00:03:41.260
In a nutshell, we arrive a
function and crosslinkSuccess,

00:03:41.260 --> 00:03:43.870
which test whether each
subunit of the polymer

00:03:43.870 --> 00:03:47.620
would touch a neighboring
particle by using RegionMember.

00:03:58.670 --> 00:04:02.060
If there is a least a subunit
that successfully cross

00:04:02.060 --> 00:04:06.050
a particle, the polymer
connects two or more particles

00:04:06.050 --> 00:04:09.730
and the function return 1.

00:04:09.730 --> 00:04:12.370
If the polymer fails
to cross a particle,

00:04:12.370 --> 00:04:13.900
the function will return 0.

00:04:17.320 --> 00:04:21.010
So we run this function
over many polymers

00:04:21.010 --> 00:04:23.050
by using the PR function.

00:04:23.050 --> 00:04:25.420
And this is a
resource that we got.

00:04:25.420 --> 00:04:30.520
And we really want to know the
ratio of the cross-link polymer

00:04:30.520 --> 00:04:33.460
to a total number of polymers.

00:04:33.460 --> 00:04:35.470
And to do that, we
just do a simple math

00:04:35.470 --> 00:04:40.045
by dividing the
number of successfully

00:04:40.045 --> 00:04:43.345
cross-linked polymer to the
total number of polymers.

00:04:46.940 --> 00:04:48.650
Below is a function
that calculate

00:04:48.650 --> 00:04:51.950
the ratio of polymer that link
two or more particles together.

00:05:01.570 --> 00:05:04.565
Next, we calculate the ratio
of cross linking for polymer

00:05:04.565 --> 00:05:07.530
of various lengths--

00:05:07.530 --> 00:05:15.410
100, 200, 500, 800,
1,000 and 1,500.

00:05:15.410 --> 00:05:17.555
Since the operation
take a while,

00:05:17.555 --> 00:05:18.880
I'll show just the results.

00:05:25.900 --> 00:05:29.110
The most interesting thing that
we can obtain for this data

00:05:29.110 --> 00:05:31.960
is actually plotting
the ratio of cross-link

00:05:31.960 --> 00:05:33.370
versus the polymer length.

00:05:36.400 --> 00:05:39.130
What we find out is that
as the polymer length

00:05:39.130 --> 00:05:42.820
increases, the probability
to cross-link increases.

00:05:42.820 --> 00:05:46.150
This makes sense because
as a polymer gets longer,

00:05:46.150 --> 00:05:49.750
it can walk a larger distance
and touch a particle.

00:05:49.750 --> 00:05:53.020
What is even more
interesting is that the graph

00:05:53.020 --> 00:05:55.300
has a shape of a lock function.

00:05:55.300 --> 00:05:59.680
We go ahead and
use Fit to generate

00:05:59.680 --> 00:06:02.320
least square fit equation.

00:06:02.320 --> 00:06:05.990
And then we plot
it with the data.

00:06:05.990 --> 00:06:10.000
As you can see, the
square fit function

00:06:10.000 --> 00:06:15.470
give you a pretty
good fit of the data.

00:06:15.470 --> 00:06:18.130
What this fit tell
you is that first,

00:06:18.130 --> 00:06:21.380
there's a minimum
number of subunits

00:06:21.380 --> 00:06:24.220
that cross-link will occur.

00:06:24.220 --> 00:06:27.520
In our case, the
minimum polymer length

00:06:27.520 --> 00:06:31.510
will have to be between
100 and 200 subunits

00:06:31.510 --> 00:06:34.930
to be able to provide you
some type of cross-linking

00:06:34.930 --> 00:06:37.900
between polymer
and nanoparticles.

00:06:37.900 --> 00:06:41.950
Second, the graph plateau
at a certain polymer length.

00:06:41.950 --> 00:06:44.830
Any polymers that
are longer than those

00:06:44.830 --> 00:06:47.410
would not increase the
probability of cross-link.

00:06:51.740 --> 00:06:54.710
Although the analysis of
probability for polymer

00:06:54.710 --> 00:06:57.230
to cross-link it's
useful, it does not

00:06:57.230 --> 00:07:00.150
tell any information about
the network of polymers

00:07:00.150 --> 00:07:03.580
and particle as a whole.

00:07:03.580 --> 00:07:07.860
So in stage 4, I will
use cluster analysis

00:07:07.860 --> 00:07:11.190
to analyze network formation
given polymer length

00:07:11.190 --> 00:07:14.390
and particle distribution.

00:07:14.390 --> 00:07:20.570
In a nutshell, the function
singlePolymerCrossLinkTest will

00:07:20.570 --> 00:07:25.610
return all the particles'
lattice site which the input

00:07:25.610 --> 00:07:27.530
polymer can connect.

00:07:30.120 --> 00:07:38.580
Then we use a built-in function,
UndirectedEdge and Graph

00:07:38.580 --> 00:07:41.990
to visualize this connection.

00:07:41.990 --> 00:07:45.140
We want to do this simulation
for all the polymers

00:07:45.140 --> 00:07:49.160
and particle to visualize the
connectedness with the polymer

00:07:49.160 --> 00:07:53.025
nanoparticle network then use
a cluster analysis to analyze

00:07:53.025 --> 00:07:55.395
the network formed
by mixing polymers

00:07:55.395 --> 00:07:57.005
and nanoparticles together.

00:08:02.920 --> 00:08:09.390
In here, we write
a little function

00:08:09.390 --> 00:08:14.640
that basically run
the single polymer

00:08:14.640 --> 00:08:20.640
cross-link test over many
polymer of the same length.

00:08:20.640 --> 00:08:27.930
What I want to show you
is this GraphCommunityPlot

00:08:27.930 --> 00:08:30.840
and FindGraphCommunities.

00:08:30.840 --> 00:08:32.730
In this graph, you
see that there's

00:08:32.730 --> 00:08:37.650
many clusters which has a
very dense cross-linking

00:08:37.650 --> 00:08:40.650
between polymer
and nanoparticles.

00:08:40.650 --> 00:08:45.720
However, most of the clusters
do not connect it together.

00:08:45.720 --> 00:08:52.620
This suggests that a network
of polymer and nanoparticle

00:08:52.620 --> 00:08:54.690
may not be formed,
even though you

00:08:54.690 --> 00:08:57.840
have many cluster of
densely cross-linking

00:08:57.840 --> 00:08:59.565
between nanoparticles
and polymers.

00:09:02.620 --> 00:09:05.620
What I want to show
next is a comparison

00:09:05.620 --> 00:09:10.270
between a cluster network
of 100 subunit polymer

00:09:10.270 --> 00:09:13.810
and a cluster network of
1,000 subunit polymers.

00:09:13.810 --> 00:09:17.630
As you can see, as a
polymer length increase,

00:09:17.630 --> 00:09:25.100
there is a higher chance that
cluster can connect together.

00:09:25.100 --> 00:09:27.770
So definitely the
polymer length can

00:09:27.770 --> 00:09:30.980
affect the connectedness
of the system

00:09:30.980 --> 00:09:35.120
and give a higher
chance that polymer

00:09:35.120 --> 00:09:38.630
can connect nanoparticle
and form a cohesive network.

00:09:44.410 --> 00:09:48.040
In a realistic system of
polymers and nanoparticles

00:09:48.040 --> 00:09:50.710
which have favorable
interaction, one

00:09:50.710 --> 00:09:53.110
a polymer hits a
nanoparticle, it

00:09:53.110 --> 00:09:55.750
will stick to that nanoparticle.

00:09:55.750 --> 00:09:57.850
In the system that
we used above,

00:09:57.850 --> 00:10:01.360
polymers are allowed to walk
in such a way that does not

00:10:01.360 --> 00:10:03.190
take into account
the interaction

00:10:03.190 --> 00:10:05.950
between nanoparticles
and polymers.

00:10:05.950 --> 00:10:08.380
In state 5, we'll
attempt to take

00:10:08.380 --> 00:10:13.410
into account that interaction
by using this algorithm step.

00:10:13.410 --> 00:10:18.340
First, we make a lattice
of nanoparticle sites.

00:10:18.340 --> 00:10:21.570
One lattice is
randomly selected.

00:10:21.570 --> 00:10:25.240
And a random walker will spawn
from the site and can walk.

00:10:25.240 --> 00:10:27.740
If polymer walker
hits a charge site,

00:10:27.740 --> 00:10:31.960
including the original one, then
the polymer stick to that side.

00:10:35.150 --> 00:10:38.000
The modifiedPolymerRandomWalk
function

00:10:38.000 --> 00:10:42.600
is very similar to the random
walk polymer function above,

00:10:42.600 --> 00:10:49.790
except for its use RegionMember
to test whether or not

00:10:49.790 --> 00:10:52.260
a subunit touched a particle.

00:10:52.260 --> 00:10:55.790
If the subunit
touched a particle,

00:10:55.790 --> 00:10:58.670
the function will stop and
return the current path.

00:11:01.390 --> 00:11:08.960
Here I generate hundreds
polymer of 1,000 subunit

00:11:08.960 --> 00:11:14.720
by using this modified
polymer random walk function

00:11:14.720 --> 00:11:17.330
and using Graphic
to visualize it.

00:11:17.330 --> 00:11:20.780
As you can look
into this picture,

00:11:20.780 --> 00:11:24.830
the probability that polymers
will cross-link nanoparticles

00:11:24.830 --> 00:11:32.930
decreases because you restricted
it to stopping once the polymer

00:11:32.930 --> 00:11:37.040
hitting a nanoparticle.

00:11:37.040 --> 00:11:40.010
Although this approach takes
into account the interaction

00:11:40.010 --> 00:11:42.380
between nanoparticles
and polymers,

00:11:42.380 --> 00:11:44.960
it is still not a
realistic representation

00:11:44.960 --> 00:11:46.430
of a real system.

00:11:46.430 --> 00:11:48.775
In this approach, we
let a polymer subunit

00:11:48.775 --> 00:11:53.320
spawn on the particle surface
and grow into a polymer.

00:11:53.320 --> 00:11:57.650
In a real system, polymers are
thrown into a disperse solution

00:11:57.650 --> 00:11:58.880
of particles.

00:11:58.880 --> 00:12:02.630
The confirmation of polymers
in a dynamic process in which

00:12:02.630 --> 00:12:04.910
the entire polymer
twist, turns so

00:12:04.910 --> 00:12:08.660
that the two ends of a polymer
find nanoparticle surfaces.

00:12:08.660 --> 00:12:10.760
This project has not
taken into account

00:12:10.760 --> 00:12:15.770
the dynamic confirmation of
polymers in a dilute solution.

00:12:15.770 --> 00:12:18.140
For future research,
I would like

00:12:18.140 --> 00:12:22.160
to improve this project by
taking into account polymer

00:12:22.160 --> 00:12:25.220
confirmation when doing
polymer random walk

00:12:25.220 --> 00:12:27.660
and finding how
polymer concentration

00:12:27.660 --> 00:12:31.670
and nanoparticle concentration
affect the nanoparticle polymer

00:12:31.670 --> 00:12:34.460
network by doing
cluster analysis.

00:12:34.460 --> 00:12:37.370
I would like to thank Dr.
Keane for your guidance

00:12:37.370 --> 00:12:38.610
on this project.

00:12:38.610 --> 00:12:45.080
I would like to thank Professor
Carter for teaching 3016, which

00:12:45.080 --> 00:12:50.030
I found very helpful
for the other project

00:12:50.030 --> 00:12:52.250
I'm working on my
UROP right now.

00:12:52.250 --> 00:12:54.840
And I would like to tell
you all the TA's in 3.016

00:12:54.840 --> 00:12:58.389
say for helping me
during office hour.