THREE-DIMENSIONAL MODEL OF QUANTUM DOTS' SELF-ASSEMBLY UNDER THE ACTION OF LASER RADIATION
V.A. Tkachenko1, A.S. Tsipotan1, A.S. Aleksandrovsky ',2, V. V. Slabko1 1 Siberian Federal University, Krasnoyarsk, Russia, 2 Kirensky Institute of Physics, Russian Academy of Sciences, Krasnoyarsk, Russia
Abstract
This study considered a process of quantum dots' self-assembly into nanostructure arrays with predefined geometry, which proceeds in the external resonant laser field. We considered the simplest case of assembling a stable structure of two particles. The problem was solved numerically using a three-dimensional model of Brownian dynamics. The idea of the method is that the attraction of the dots occurs due to the interaction of resonantly induced dipole moments, with the dots being then captured by the Van der Waals force. Finally, a three-dimensional model was considered; the average nanoparticle aggregation time as a function of the laser radiation wavelength was calculated; the probability of such structures' being formed was estimated for the calculated average aggregation time and for the laser pulse duration used in the experiment. The wavelength of the maximum probability was found to be shifted from the single particle resonance wavelength of 525 nm to the red area of 535 nm, which is in qualitative agreement with the redshift of the resonance wavelength of interacting particles.
Keywords: Nanostructure fabrication; plasmonics; optical tweezers or optical manipulation.
Citation: Tkachenko VA, Tsipotan AS, Aleksandrovsky AS, Slabko VV. Three-dimensional model of quantum dots' self-assembly under the action of laser radiation. Computer Optics 2017; 41(4): 577-580. DOI: 10.18287/2412-6179-2017-41-4-577-580.
Acknowledgements: The work was funded by the Russian Foundation for Basic Research (RFBR) and Krasnoyarsk Krai administration under research project No. 16-42-240410r_a, RFBR research project No. 16-32-00129 and by the Ministry of Education and Science of the Russian Federation (Grant 3.6341.2017/VU).
Introduction
Colloidal quantum dots as well as nanostructures based on them attract considerable interest of researchers due to their unique properties and possibilities of there numerous applications. [1-4]. Hindrance to wide applications of nanostructures composed of quantum dots is the absence of facile route for their production that would enable stable reproduction of the same topology of a nanostructure. The development of self-organization techniques for the formation of quantum dots' structures with pre-defined configuration are under way [5, 6]. The authors of the present communication suggested an universal method of nanostructures' formation under action of external quasi-resonant field due to the self-organization of nanoparticles [7-9]. Topology and content of forming nanostructure depend on the wavelength and the polarization of external field. The first experiments [10] have shown the principal possibility of pairs' formation from an ensemble of chaotic colloidal quantum dots. It was demonstrated that if the wavelength is chosen to correspond the maximum of dipole-dipole attraction of a pair of particles then up to 50 % of initially solitary particles become aggregated into pairs. However, 2D model calculations that were performed earlier do not allow ultimate analysis of the results obtained in the above cited experiment. Present communication reports on the results of 3D modeling of the quantum dots' self-assembly in the external field that allow complementing and generalizing its peculiarities.
Calculation model
As it is known, important role in the process of quantum dots' aggregation is played by the shape of
potential curve [11-12]. Herein, the cumulate potential of pair interaction includes van der Waals forces (responsible for attraction of particles) as well as forces preventing spontaneous aggregation of nanoparticles (electrostatic and elastic repulsion). [13]. Since spontaneous aggregation leads to formation of structures with uncontrolled geometry, it is necessary that repulsion between particles prevailed the attraction in the absence of laser field. On the other hand, the presence of external action, in our case, quasi-resonant laser field allows controllable reducing of the potential barrier height, and since the process is of resonant nature, the geometry and content of the forming structure can be controlled, too.
Preliminary modeling of this process allows estimating the optimal parameters of experiment and revealing the peculiarities of nanoparticles' self-assembly process.
Interaction potential of a pair of nanoparticles
Modeling of the process was done using Brownian dynamics approach. Potential energy of pair interaction included electrostatic repulsion, van der Waals interacion and induced dipole-dipole interaction. Potential of elastic repulsion was not accounted for since it arises only when stabilizing layers of two nanoparticles become overlapped, that is after electrostatic repulsion is overrided, which moment was considered by us as the moment of two-paricle structure formation.
For describing van der Waals energy, let us employ the approximation of Hamaker - de Burgh theory, according to which the interaction energy of two spherical particles with the radius ri made of the same material is expressed [14] as
W
2r.2
2r.2
h2 + 4rh + h2 + 4rh + 2r2
+ ln-
h + 4rh ' h1 + 4rh + 2r2
A
(1)
where Ah is effective Hamaker constant (Ah » 50kT, here T = 300 K, k is Boltzmann constant), h = r--- 2r,, is the in-terparticle spacing, and r- is the distance between particle centers.
General case of electrostatic interaction of two spherical nanoparticles is described on the base of DLVO theory. Every particle in a solution is surrounded by electric double layer (EDL). When spacing between particle surfaces exceeds the thickness of these layers lo (known also
3r„r,
as Debye-Huckel radius or the size of diffuse part of antipolar ions layer), electrostatic interaction is practically absent, while at smaller distance the layers are overlapped.
Potential energy We of Coulomb interaction of overlapping EDLs is:
We = 2nee0r i f ln[1+exp(-hk0)],
(2)
1/ C d1-
3r- - r2
3rr
5 d1x
3r- - r2
1/ Co d2
3r r
5 U2x
0
3r r
__x-JLd
„5 d1y
where ko = 1 / lo is the screening constant; e, eo are the dielectric permeability of the medium and that of vacuum; f o is the potential at the boundary of Helmholtz layer, the thickness of which 5 does not exceed the size of hydrated ions constituting it.
The switch-over from 2D case [14] to 3D one requires essential modification of the equation system for calculation of dipole moments that will have the form:
3r r
d
2y
0
3r r
r
1/ Cod]
0"1y
3ry - r2
1y
3ry - r2
1/ %0d:
2y
0 2y
3r r
__y_Ld
r5 d2y
0
3rr
r5 ^
3r r
1/ %0 d1z
3rz2 - r2
d
3r r
d
3rz2 - r2
1/ C d 2
^0: E
0y
0y
^0; En.
(3)
where %o = |di2|2 / h(Q + /T), W is the detuning of external field frequency from the resonance, h is the Planck constant, and T is the homogeneous line width, |d12|2 electric dipole moment of the first excitonic transition; r, rw are the lengths of the vector connecting the particle, and of its projections; d(i,2xW) are the projections of dipole moments of particles 1 and 2; Eo, £o(w) are the modulus of field strength of external field and its projections.
Expression for dipole-dipole interaction energy of two particles can be presented as:
W =
d1d2r2 -3(d1r)(d2r)
(4)
from which the formula for dipole-dipole interaction force is obtained
Fdd =
9 (d1r )(d2r )-3 d1d2
(5)
Accounting all forces of pair interaction allows correct description of Brownian motion in the field of laser radiation and modeling of the self-assembly of nanoparticles into pre-defined structures.
Methods
Within the frame of above-described model, the interaction of two semiconductor nanoparticles suspended in water at limited distance under action of external laser radiation was considered.
The problem was investigated numerically by simple time-stepping technique. At the beginning of every pass the first particle was placed into the center of bounding cubic box while second particle was placed into a random point inside the box. The dimension of the box was chosen from the requirement that at the maximum initial distance between particles Lmax = 21.7 nm the energy of their attraction was equal to the thermal motion energy kT [14]; consequently, the cube edge is defined to be a = 2Lmax / V3 » 25 nm.
At every calculation step the dipole moments, cumulate forces acting on the particles, the velocities of particles and their renewed spatial coordinates were calculated. The viscous friction force, van der Waals force (1), electrostatic repulsion of surface potentials (2), and dipole-dipole interaction force (4) were accounted for. Besides, at every step the particle acquired additional random momentum in accordance with Brownian dynamics.
Calculations terminated either at the step when in-terparticle distance had become smaller than the sum of radii or when one of the particles leaved the bounding box. First case corresponds to particles' collision and pair formation while the second case means the fail of aggregation. Final output values calculated at every pass are the successfulness of aggregation and the required time that corresponds to the number of steps. Based on 1oo repetitions of modeling pass, the cumulate probability of
+
6
r
r
r
d
d
0
0
5
r
r
0
0
r
r
r
0
0
r
r
r
0
0
r
0
r
r
2
r
6
r
aggregation, the probability of aggregation per time (10 ns) and average aggregation time were determined. Every pass differed by random initial position of the second particle and by random action of Brownian forces.
Calculations were done for several values of wavelength in the vicinity of resonance, that for the particles under consideration was equal to 1r = 525 nm. Radius of particles including the absorbed layer was r=6 nm, the density r = 5.81 g/cm3, electric dipole moment of the first excitonic transition |d^|2 = 2-10-46 J-m3. These particles' parameters corresponded to those of CdTe quantum dots used in the experimental study of the authors [10]. Water viscosity was taken to be h = 8.9-10-5 P, dielectric constant e = 81; laser field strength in all calculations was taken to be E0 = 6-106 V/m.
Results and Discussion
Fig. 1 presents the dependence of aggregation probability on the laser wavelength. The maximum of the probability is close to 100% at the wavelength l = 535.5 nm, somewhat redshifted with respect to the resonant value for isolated particle 1r = 525 nm. This redshift evidently is due to the corresponding redshift of the resonance of closely spaced interacting pair of particles. [7-9]. Probability of aggregation steadily diminishes with the increase of the wavelength above 545 nm and approaches zero at
525 530 535 540 545 550 X(nm)
Fig. 1. Dependence of aggregation probability on the laser wavelength
Average aggregation time of a pair experiences minimum at the wavelengths close to the maximum aggregation probability and grows with the detuning off these wavelengths as plotted in Fig. 2. The aggregation wavelength is maximum this region and approaches 100 %. t(hs)
525 530 535 540 545 550 X(nm)
Fig. 2. Dependence of average aggregation time for a pair of particles on the laser wavelength
Average aggregation time is of order of microseconds for all studied wavelengths. Therefore, high aggregation probability can be obtained using cw laser or pulsed laser
with the pulse duration of order of several microsecond. However, use of cw laser with the intensity necessary for aggregation of nanoparticles (I = 106 W/cm2) will lead to strong heating of both nanoparticles and local environment, and negative influence on the obtained structures must be expected. Pulsed lasers with microsecond pulse duration with the possibility of tuning in the broad wavelength range necessary for efficient assembly do not exist. Due to this fact, the experiments in [10] were conducted using tunable optical parametric oscillator with the required power and pulse duration 10 nm. Therefore, it is worthy to estimate the probability of pair formation during a single 10 nm pulse (Fig. 3).
Pl0ns(%)
8
525 530 535 540 545 550 X(nm)
Fig. 3. Wavelength dependence of aggregation probability during a single pulse with the 10 ns duration
Maximum probability is at l = 525 nm and near to 8 %. Therefore, using repeated pulsed action allows achieving formation of acceptable number of structures, which was observed in the experiment. [10] (after 3000 pulses approximately 40 % of particles were assembled into pairs).
Conclusion
3D modeling of nanoparticles' self-assembly under the action of quasi-resonant laser field resulted in estimation of probability of this process and of characteristic time required. Maximum probability and shortest time required for self-assembly is achieved at the external field wavelength 535.5 nm, which is redshifted with respect to the resonance of isolated nanoparticle, in accordance with earlier 2D calculations and 3D experiments. Average assembly time is as high as several microseconds. Since 10 nm laser was used in the experiment, the probability of aggregation of a pair of nanoparticles during a single pulse was calculated, being near to 8 %.
References
[1] Yang Y, Zheng Y, Cao W, Titov A, Hyvonen J, Manders JR, Xue J, Holloway PH, Qian L. High-efficiency light-emitting devices based on quantum dots with tailored nanostructures. Nature Photonics 2014; 9: 259-266. DOI: 10.1038/nphoton.2015.36.
[2] Volkov Y. Quantum dots in nanomedicine: Recent trends, advances and unresolved issues. Biochemical and Biophysical Research Communications 2015; 468(3): 419-427. DOI: 10.1016/j.bbrc.2015.07.039.
[3] Zhang J, Cheng F, Li J, Zhu J, Lu Y. Fluorescent nano-probes for sensing and imaging of metal ions: Recent advances and future perspectives. Nanotoday 2016; 11(3): 309-329. DOI: 10.1016/j.nantod.2016.05.010.
[4] Frecker T, Bailey D, Arzeta-Ferrer X, McBride J, Rosenthal SJ. Review - Quantum dots and their application in
lighting, displays, and biology. ESC Journal of Solid State Science and Technology 2o16; 5(1): R3o19-R3o31. DOI: 1o.1149/2.oo316o1jss.
[5] Tikhomirov G, Hoogland S, Lee PE, Fischer A, Sargent EH, Kelley SO. DNA-based programming of quantum dot valency, self-assembly and luminescence. Nature Nano-technology 2o11; 6; 485-49o. DOI: 1o.1o38/nna-no.2o11.1oo.
[6] Wang T, LaMontagne D, Lynch J, Zhuangc J, Cao YC. Colloidal superparticles from nanoparticle assembly. Chem Soc Rev 2o13; 42(7): 28o4-2823. DOI: 1o.1o39/c2cs35318k.
[7] Slabko VV, Khachatryan GG, Aleksandrovsky AS. Self-organized aggregation of small metal particles controlled by an external light field. JETP Letters 2oo6; 84(6): 3oo-3o4. DOI: 1o.1134/Soo21364oo618oo56.
[8] Slabko VV, Tsipotan AS, Aleksandrovsky AS. Resonant light-controlled self-assembly of ordered nanostructures. Photonics and Nanostructures - Fundamentals and Applications 2o12; 1o(4): 636-643. DOI: 1o.1o16/j.photo-nics.2o12.o6.oo2.
[9] Slabko VV, Tsipotan AS, Aleksandrovsky AS. Self-organised aggregation of a pair of particles with different resonant frequencies and electric dipole moments of transitions, controlled by an external quasi-resonant field. Quan-
tum Electronics 2o13; 43(5): 458-462. DOI: 1o.1 o7o/QE2o 13vo43no5ABEHo 15o21.
[10] Tsipotan AS, Gerasimova MA, Slabko VV, Aleksandrov-sky AS. Laser-induced wavelength-controlled self-assembly of colloidal quasi-resonant quantum dots. Optics Express 2o16; 24(1o): 11145-1115o. DOI: 1o.1364/OE.24.o11145.
[11] Wang C-L, Fang M, Han J-S, Zhang H, Cui Y-P, Yang B. Two opposite effects of alcohols in the precipitation of aqueous nanocrystals. The Journal of Physical Chemistry C 2oo9; 113(45): 19445-19451. DOI: 1o.1o21/jp9o54o7y.
[12] Xu J, Xu S, Lv C, Wang C, Cui Y. Excellent dynamic stability under saturated salt solution for aqueous quantum dots capped by multi-branched ligands. Materials Research Express 2o16; 3(9): o959o3. DOI: 1o.1o88/2o53-1591/3/9/o959o3.
[13] Zhang H, Liu Y, Zhang J, Wang C, Li M, Yang B. Influence of interparticle electrostatic repulsion in the initial stage of aqueous semiconductor nanocrystal growth. The Journal of Physical Chemistry C 2oo8; 112(6): 1885-1889. DOI: 1o.1o21/jp71o11on.
[14] Slabko VV, Tsipotan AS, Aleksandrovsky AS, Slyusareva EA. Dynamics of self-organized aggregation of resonant nanoparticles in a laser field. Applied Physics B 2o14; 117(1): 271-278. DOI: 1o.1oo7/soo34o-o14-5831-o.
Authors' information
Viktor Aleksandrovich Tkachenko (b. 1989) graduated from Siberian Federal University (SFU) in 2o14, majoring in Physics. Currently he is a PhD student in SFU. Research interests are nonlinear optics, self assembling of nanostructures, programming. E-mail: [email protected] .
Aleksey Sergeevich Tsipotan (b. 1988), PhD (Siberian Federal University, 2o15). Currently he works as scientific worker at laboratory for nonlinear optics and photonics Siberian Federal University. Research interests are nanophoton-ics, nanotechnology. E-mail: [email protected].
Vitaliy Vasilyevich Slabko (b. 1949), Doctor, Professor (Siberian Federal University), Currently he works as leading scientific researcher, professor in Siberian Federal University. Research interests are nonlinear optics, nanophoton-ics. E-mail: [email protected] .
Code of State Categories Scientific and Technical Information (in Russian - GRNTI): 29.31.27. Received June 6, 2017. The final version - August 20, 2017.