1
|
Eatson JL, Morgan SO, Horozov TS, A. Buzza DM. Programmable 2D materials through shape-controlled capillary forces. Proc Natl Acad Sci U S A 2024; 121:e2401134121. [PMID: 39163335 PMCID: PMC11363311 DOI: 10.1073/pnas.2401134121] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/17/2024] [Accepted: 07/14/2024] [Indexed: 08/22/2024] Open
Abstract
In recent years, self-assembly has emerged as a powerful tool for fabricating functional materials. Since self-assembly is fundamentally determined by the particle interactions in the system, if we can gain full control over these interactions, it would open the door for creating functional materials by design. In this paper, we exploit capillary interactions between colloidal particles at liquid interfaces to create two-dimensional (2D) materials where particle interactions and self-assembly can be fully programmed using particle shape alone. Specifically, we consider colloidal particles which are polygonal plates with homogeneous surface chemistry and undulating edges as this particle geometry gives us precise and independent control over both short-range hard-core repulsions and longer-range capillary interactions. To illustrate the immense potential provided by our system for programming self-assembly, we use minimum energy calculations and Monte Carlo simulations to show that polygonal plates with different in-plane shapes (hexagons, truncated triangles, triangles, squares) and edge undulations of different multipolar order (hexapolar, octopolar, dodecapolar) can be used to create a rich variety of 2D structures, including hexagonal close-packed, honeycomb, Kagome, and quasicrystal lattices. Since the required particle shapes can be readily fabricated experimentally, we can use our colloidal system to control the entire process chain for materials design, from initial design and fabrication of the building blocks, to final assembly of the emergent 2D material.
Collapse
Affiliation(s)
- Jack L. Eatson
- Department of Physics and Astrophysics, George William Gray Centre for Advanced Materials, University of Hull, HullHU6 7RX, United Kingdom
| | - Scott O. Morgan
- Department of Physics and Astrophysics, George William Gray Centre for Advanced Materials, University of Hull, HullHU6 7RX, United Kingdom
| | - Tommy S. Horozov
- Department of Chemistry and Biochemistry, George William Gray Centre for Advanced Materials, University of Hull, HullHU6 7RX, United Kingdom
| | - D. Martin A. Buzza
- Department of Physics and Astrophysics, George William Gray Centre for Advanced Materials, University of Hull, HullHU6 7RX, United Kingdom
| |
Collapse
|
2
|
Eatson JL, Gordon JR, Cegielski P, Giesecke AL, Suckow S, Rao A, Silvestre OF, Liz-Marzán LM, Horozov TS, Buzza DMA. Capillary Assembly of Anisotropic Particles at Cylindrical Fluid-Fluid Interfaces. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2023; 39:6006-6017. [PMID: 37071832 PMCID: PMC10157885 DOI: 10.1021/acs.langmuir.3c00016] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/03/2023]
Abstract
The unique behavior of colloids at liquid interfaces provides exciting opportunities for engineering the assembly of colloidal particles into functional materials. The deformable nature of fluid-fluid interfaces means that we can use the interfacial curvature, in addition to particle properties, to direct self-assembly. To this end, we use a finite element method (Surface Evolver) to study the self-assembly of rod-shaped particles adsorbed at a simple curved fluid-fluid interface formed by a sessile liquid drop with cylindrical geometry. Specifically, we study the self-assembly of single and multiple rods as a function of drop curvature and particle properties such as shape (ellipsoid, cylinder, and spherocylinder), contact angle, aspect ratio, and chemical heterogeneity (homogeneous and triblock patchy). We find that the curved interface allows us to effectively control the orientation of the rods, allowing us to achieve parallel, perpendicular, or novel obliquely orientations with respect to the cylindrical drop. In addition, by tuning particle properties to achieve parallel alignment of the rods, we show that the cylindrical drop geometry favors tip-to-tip assembly of the rods, not just for cylinders, but also for ellipsoids and triblock patchy rods. Finally, for triblock patchy rods with larger contact line undulations, we can achieve strong spatial confinement of the rods transverse to the cylindrical drop due to the capillary repulsion between the contact line undulations of the particle and the pinned contact lines of the sessile drop. Our capillary assembly method allows us to manipulate the configuration of single and multiple rod-like particles and therefore offers a facile strategy for organizing such particles into useful functional materials.
Collapse
Affiliation(s)
- Jack L Eatson
- Department of Physics & Mathematics, University of Hull, Hull HU6 7RX, U.K
| | - Jacob R Gordon
- Department of Chemistry & Biochemistry, University of Hull, Hull HU6 7RX, U.K
| | | | - Anna L Giesecke
- AMO GmbH, Otto-Blumenthal-Str. 25, Aachen 52074, Germany
- University of Duisburg-Essen, Bismarckstr. 81, Duisburg 47057, Germany
| | - Stephan Suckow
- AMO GmbH, Otto-Blumenthal-Str. 25, Aachen 52074, Germany
| | - Anish Rao
- Center for Cooperative Research in Biomaterials (CIC BiomaGUNE), Basque Research and Technology Alliance (BRTA), Paseo de Miramón 182, Donostia-San Sebastián 20014, Spain
| | - Oscar F Silvestre
- Center for Cooperative Research in Biomaterials (CIC BiomaGUNE), Basque Research and Technology Alliance (BRTA), Paseo de Miramón 182, Donostia-San Sebastián 20014, Spain
- Centro de Investigación Biomédica en Red, Bioingeniería, Biomateriales y Nanomedicina (CIBER-BBN), Paseo de Miramón 182, Donostia-San Sebastián 20014, Spain
| | - Luis M Liz-Marzán
- Center for Cooperative Research in Biomaterials (CIC BiomaGUNE), Basque Research and Technology Alliance (BRTA), Paseo de Miramón 182, Donostia-San Sebastián 20014, Spain
| | - Tommy S Horozov
- Department of Chemistry & Biochemistry, University of Hull, Hull HU6 7RX, U.K
| | - D Martin A Buzza
- Department of Physics & Mathematics, University of Hull, Hull HU6 7RX, U.K
| |
Collapse
|
3
|
Morgan SO, Muravitskaya A, Lowe C, Adawi AM, Bouillard JSG, Horozov TS, Stasiuk GJ, Buzza DMA. Using adsorption kinetics to assemble vertically aligned nanorods at liquid interfaces for metamaterial applications. Phys Chem Chem Phys 2022; 24:11000-11013. [PMID: 35467675 DOI: 10.1039/d1cp05484h] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Vertically aligned monolayers of metallic nanorods have a wide range of applications as metamaterials or in surface enhanced Raman spectroscopy. However the fabrication of such structures using current top-down methods or through assembly on solid substrates is either difficult to scale up or have limited possibilities for further modification after assembly. The aim of this paper is to use the adsorption kinetics of cylindrical nanorods at a liquid interface as a novel route for assembling vertically aligned nanorod arrays that overcomes these problems. Specifically, we model the adsorption kinetics of the particle using Langevin dynamics coupled to a finite element model, accurately capturing the deformation of the liquid meniscus and particle friction coefficients during adsorption. We find that the final orientation of the cylindrical nanorod is determined by their initial attack angle when they contact the liquid interface, and that the range of attack angles leading to the end-on state is maximised when nanorods approach the liquid interface from the bulk phase that is more energetically favorable. In the absence of an external field, only a fraction of adsorbing nanorods end up in the end-on state (≲40% even for nanorods approaching from the energetically favourable phase). However, by pre-aligning the metallic nanorods with experimentally achievable electric fields, this fraction can be effectively increased to 100%. Using nanophotonic calculations, we also demonstrate that the resultant vertically aligned structures can be used as epsilon-near-zero and hyperbolic metamaterials. Our kinetic assembly method is applicable to nanorods with a range of diameters, aspect ratios and materials and therefore represents a versatile, low-cost and powerful platform for fabricating vertically aligned nanorods for metamaterial applications.
Collapse
Affiliation(s)
- S O Morgan
- Department of Physics & Mathematics, University of Hull, Hull HU6 7RX, UK.
| | - A Muravitskaya
- Department of Physics & Mathematics, University of Hull, Hull HU6 7RX, UK.
| | - C Lowe
- Department of Physics & Mathematics, University of Hull, Hull HU6 7RX, UK.
| | - A M Adawi
- Department of Physics & Mathematics, University of Hull, Hull HU6 7RX, UK.
| | - J-S G Bouillard
- Department of Physics & Mathematics, University of Hull, Hull HU6 7RX, UK.
| | - T S Horozov
- Department of Chemistry & Biochemistry, University of Hull, Hull HU6 7RX, UK
| | - G J Stasiuk
- Imaging Chemistry & Biology, King's College London, Strand, London WC2R 2LS, UK
| | - D M A Buzza
- Department of Physics & Mathematics, University of Hull, Hull HU6 7RX, UK.
| |
Collapse
|
4
|
Morgan SO, Fox J, Lowe C, Adawi AM, Bouillard JSG, Stasiuk GJ, Horozov TS, Buzza DMA. Adsorption trajectories of nonspherical particles at liquid interfaces. Phys Rev E 2021; 103:042604. [PMID: 34005913 DOI: 10.1103/physreve.103.042604] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/05/2020] [Accepted: 03/22/2021] [Indexed: 11/07/2022]
Abstract
The adsorption of colloidal particles at liquid interfaces is of great importance scientifically and industrially, but the dynamics of the adsorption process is still poorly understood. In this paper we use a Langevin model to study the adsorption dynamics of ellipsoidal colloids at a liquid interface. Interfacial deformations are included by coupling our Langevin dynamics to a finite element model while transient contact line pinning due to nanoscale defects on the particle surface is encoded into our model by renormalizing particle friction coefficients and using dynamic contact angles relevant to the adsorption timescale. Our simple model reproduces the monotonic variation of particle orientation with time that is observed experimentally and is also able to quantitatively model the adsorption dynamics for some experimental ellipsoidal systems but not others. However, even for the latter case, our model accurately captures the adsorption trajectory (i.e., particle orientation versus height) of the particles. Our study clarifies the subtle interplay between capillary, viscous, and contact line forces in determining the wetting dynamics of micron-scale objects, allowing us to design more efficient assembly processes for complex particles at liquid interfaces.
Collapse
Affiliation(s)
- S O Morgan
- Department of Physics & Mathematics, University of Hull, Hull HU6 7RX, United Kingdom
| | - J Fox
- Department of Physics & Mathematics, University of Hull, Hull HU6 7RX, United Kingdom.,School of Physics & Astronomy, University of Leeds, Leeds LS2 9JT, United Kingdom
| | - C Lowe
- Department of Physics & Mathematics, University of Hull, Hull HU6 7RX, United Kingdom
| | - A M Adawi
- Department of Physics & Mathematics, University of Hull, Hull HU6 7RX, United Kingdom
| | - J-S G Bouillard
- Department of Physics & Mathematics, University of Hull, Hull HU6 7RX, United Kingdom
| | - G J Stasiuk
- Imaging Chemistry & Biology, King's College London, Strand, London WC2R 2LS, United Kingdom
| | - T S Horozov
- Department of Chemistry & Biochemistry, University of Hull, Hull HU6 7RX, United Kingdom
| | - D M A Buzza
- Department of Physics & Mathematics, University of Hull, Hull HU6 7RX, United Kingdom
| |
Collapse
|
5
|
Newton B, Mohammed R, Davies GB, Botto L, Buzza DMA. Capillary Interaction and Self-Assembly of Tilted Magnetic Ellipsoidal Particles at Liquid Interfaces. ACS OMEGA 2018; 3:14962-14972. [PMID: 31458162 PMCID: PMC6644019 DOI: 10.1021/acsomega.8b01818] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/30/2018] [Accepted: 10/24/2018] [Indexed: 05/04/2023]
Abstract
Magnetic ellipsoidal particles adsorbed at a liquid interface provide exciting opportunities for creating switchable functional materials, where self-assembly can be switched on and off using an external field [Davies et al., Adv. Mater., 2014, 26, 6715]. In order to gain a deeper understanding of this novel system in the presence of an external field, we study the capillary interaction and self-assembly of tilted ellipsoids using analytical theory and finite element simulations. We derive an analytical expression for the dipolar capillary interaction between tilted ellipsoids in elliptical polar coordinates, which exhibits a 1/r 2 power law dependence in the far field (i.e., large particle separations r) and correctly captures the orientational dependence of the capillary interactions in the near field. Using this dipole potential and finite element simulations, we further analyze the energy landscape of particle clusters consisting of up to eight tilted ellipsoids in contact. For clusters of two particles, we find that the side-to-side configuration is stable, whereas the tip-to-tip configuration is unstable. However, for clusters of more than three particles, we find that circular loops of side-to-side particles become globally stable, whereas linear chains of side-to-side particles become metastable. Furthermore, the energy barrier for the linear-to-loop transition decreases with increasing particle number. Our results explain both thermodynamically and kinetically why tilted ellipsoids assemble side-to-side locally but have a strong tendency to form loops on larger length scales.
Collapse
Affiliation(s)
- Bethany
J. Newton
- Nano3
Group, School of Mathematics & Physical Sciences, University of Hull, Hull HU6 7RX, U.K.
| | - Rizwaan Mohammed
- Nano3
Group, School of Mathematics & Physical Sciences, University of Hull, Hull HU6 7RX, U.K.
- Clare
College, Trinity Lane, Cambridge CB2 1TL, U.K.
| | - Gary B. Davies
- Institute
for Computational Physics, Allmandring 3, 70569 Stuttgart, Germany
| | - Lorenzo Botto
- School
of Engineering and Materials Science, Queen
Mary, University of London, London E1 4NS, U.K.
| | - D. Martin A. Buzza
- Nano3
Group, School of Mathematics & Physical Sciences, University of Hull, Hull HU6 7RX, U.K.
- E-mail: (D.M.A.B.)
| |
Collapse
|
6
|
Xie Q, Davies GB, Harting J. Direct Assembly of Magnetic Janus Particles at a Droplet Interface. ACS NANO 2017; 11:11232-11239. [PMID: 29035521 DOI: 10.1021/acsnano.7b05550] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/07/2023]
Abstract
Self-assembly of nanoparticles at fluid-fluid interfaces is a promising route to fabricate functional materials from the bottom-up. However, directing and controlling particles into highly tunable and predictable structures, while essential, is a challenge. We present a liquid interface assisted approach to fabricate nanoparticle structures with tunable properties. To demonstrate its feasibility, we study magnetic Janus particles adsorbed at the interface of a spherical droplet placed on a substrate. With an external magnetic field turned on, a single particle moves to the location where its position vector relative to the droplet center is parallel to the direction of the applied field. Multiple magnetic Janus particles arrange into reconfigurable hexagonal lattice structures and can be directed to assemble at desirable locations on the droplet interface by simply varying the magnetic field direction. We develop an interface energy model to explain our observations, finding excellent agreement. Finally, we demonstrate that the external magnetic field allows one to tune the particle deposition pattern obtained when the droplet evaporates. Our results have implications for the fabrication of varied nanostructures on substrates for use in nanodevices, organic electronics, or advanced display, printing, and coating applications.
Collapse
Affiliation(s)
- Qingguang Xie
- Department of Applied Physics, Eindhoven University of Technology , P.O. Box 513, 5600MB Eindhoven, The Netherlands
| | - Gary B Davies
- St Paul's Girls' School , Brook Green, Hammersmith, London W6 7BS, United Kingdom
| | - Jens Harting
- Department of Applied Physics, Eindhoven University of Technology , P.O. Box 513, 5600MB Eindhoven, The Netherlands
- Forschungszentrum Jülich, Helmholtz Institute Erlangen-Nürnberg for Renewable Energy (IEK-11) , Fürther Straße 248, 90429 Nürnberg, Germany
| |
Collapse
|
7
|
Dasgupta S, Auth T, Gompper G. Nano- and microparticles at fluid and biological interfaces. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2017; 29:373003. [PMID: 28608781 PMCID: PMC7104866 DOI: 10.1088/1361-648x/aa7933] [Citation(s) in RCA: 21] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/30/2016] [Revised: 04/12/2017] [Accepted: 06/13/2017] [Indexed: 05/05/2023]
Abstract
Systems with interfaces are abundant in both technological applications and biology. While a fluid interface separates two fluids, membranes separate the inside of vesicles from the outside, the interior of biological cells from the environment, and compartmentalize cells into organelles. The physical properties of interfaces are characterized by interface tension, those of membranes are characterized by bending and stretching elasticity. Amphiphilic molecules like surfactants that are added to a system with two immiscible fluids decrease the interface tension and induce a bending rigidity. Lipid bilayer membranes of vesicles can be stretched or compressed by osmotic pressure; in biological cells, also the presence of a cytoskeleton can induce membrane tension. If the thickness of the interface or the membrane is small compared with its lateral extension, both can be described using two-dimensional mathematical surfaces embedded in three-dimensional space. We review recent work on the interaction of particles with interfaces and membranes. This can be micrometer-sized particles at interfaces that stabilise emulsions or form colloidosomes, as well as typically nanometer-sized particles at membranes, such as viruses, parasites, and engineered drug delivery systems. In both cases, we first discuss the interaction of single particles with interfaces and membranes, e.g. particles in external fields, non-spherical particles, and particles at curved interfaces, followed by interface-mediated interaction between two particles, many-particle interactions, interface and membrane curvature-induced phenomena, and applications.
Collapse
Affiliation(s)
- S Dasgupta
- Mechanobiology Institute, National University of Singapore, Singapore 117411, Singapore
- Institut Curie, CNRS, UMR 168, 75005 Paris, France
- Present address: Department of Physics, University of Toronto, Toronto, Ontario M5S1A7, Canada
| | - T Auth
- Theoretical Soft Matter and Biophysics, Institute of Complex Systems and Institute for Advanced Simulation, Forschungszentrum Jülich, 52425 Jülich, Germany
| | - G Gompper
- Theoretical Soft Matter and Biophysics, Institute of Complex Systems and Institute for Advanced Simulation, Forschungszentrum Jülich, 52425 Jülich, Germany
| |
Collapse
|
8
|
Shi W, Lee YH, Ling XY, Li S. Quantitative prediction of the position and orientation for an octahedral nanoparticle at liquid/liquid interfaces. NANOSCALE 2017; 9:11239-11248. [PMID: 28753214 DOI: 10.1039/c7nr02194a] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/07/2023]
Abstract
Shape-controlled polyhedral particles and their assembled structures have important applications in plasmonics and biosensing, but the interfacial configurations that will critically determine their resultant assembled structures are not well-understood. Hence, a reliable theory is desirable to predict the position and orientation of a polyhedron at the vicinity of a liquid/liquid interface. Here we demonstrate that the free energy change theory can quantitatively predict the position and orientation of an isolated octahedral nanoparticle at a liquid/liquid interface, whose vertices and facets can play crucial roles in biosensing. We focus on two limiting orientations of an octahedral nanoparticle, vertex up and facet up. Our proposed theory indicates that the surface wettability (hydrophilic/hydrophobic ratio) of the nanoparticle determines its most stable position and the preferred orientation at a water/oil interface. The surface wettability of an octahedron is adjusted from extremely hydrophobic to extremely hydrophilic by changing the amount of charge on the Ag surface in molecular dynamics (MD) simulations. The MD simulations results are in excellent agreement with our theoretical prediction for an Ag octahedral nanoparticle at a hexane/water interface. Our proposed theory bridges the gap between molecular-level simulations and equilibrium configurations of polyhedral nanoparticles in experiments, where insights from nanoparticle intrinsic wettability details can be used to predict macroscopic superlattice formation experimentally. This work advances our ability to precisely predict the final structures of the polyhedral nanoparticle assemblies at a liquid/liquid interface.
Collapse
Affiliation(s)
- Wenxiong Shi
- Center for Programmable Materials, School of Materials Science and Engineering, Nanyang Technological University, Singapore 639798.
| | | | | | | |
Collapse
|
9
|
Li T, Brandani G, Marenduzzo D, Clegg PS. Colloidal Spherocylinders at an Interface: Flipper Dynamics and Bilayer Formation. PHYSICAL REVIEW LETTERS 2017; 119:018001. [PMID: 28731771 DOI: 10.1103/physrevlett.119.018001] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/10/2016] [Indexed: 06/07/2023]
Abstract
We study the response of a film of colloidal spherocylinders to compression by combining pressure-area isotherm measurements, microscopy, and computer simulations. We find that the behavior of the film depends strongly on the geometry of the particles. For a small aspect ratio, a uniform monolayer forms and then buckles. For a higher aspect ratio, particles flip to orient perpendicular to the interface; we show that flipping occurs in locations where the nematic ordering is low. Our experiments and simulations further demonstrate that the longest particles rearrange to self-assemble a colloidal bilayer, which is stable due to the unique geometry of spherocylinders at an interface.
Collapse
Affiliation(s)
- T Li
- SUPA, School of Physics and Astronomy, University of Edinburgh, Edinburgh EH9 3FD, United Kingdom
| | - G Brandani
- SUPA, School of Physics and Astronomy, University of Edinburgh, Edinburgh EH9 3FD, United Kingdom
| | - D Marenduzzo
- SUPA, School of Physics and Astronomy, University of Edinburgh, Edinburgh EH9 3FD, United Kingdom
| | - P S Clegg
- SUPA, School of Physics and Astronomy, University of Edinburgh, Edinburgh EH9 3FD, United Kingdom
| |
Collapse
|
10
|
Coertjens S, De Dier R, Moldenaers P, Isa L, Vermant J. Adsorption of Ellipsoidal Particles at Liquid-Liquid Interfaces. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2017; 33:2689-2697. [PMID: 28241120 DOI: 10.1021/acs.langmuir.6b03534] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/20/2023]
Abstract
The adsorption of particles at liquid-liquid interfaces is of great scientific and technological importance. In particular, for nonspherical particles, the capillary forces that drive adsorption vary with position and orientation, and complex adsorption pathways have been predicted by simulations. On the basis of the latter, it has been suggested that the timescales of adsorption are determined by a balance between capillary and viscous forces. However, several recent experimental results point out the role of contact line pinning in the adsorption of particles to interfaces and even suggest that the adsorption dynamics and pathways are completely determined by the latter, with the timescales of adsorption being determined solely by particle characteristics. In the present work, the adsorption trajectories of model ellipsoidal particles are investigated experimentally using cryo-SEM and by monitoring the altitudinal orientation angle using high-speed confocal microscopy. By varying the viscosity and the viscosity jump across the interfaces, we specifically interrogate the role of viscous forces.
Collapse
Affiliation(s)
- Stijn Coertjens
- Department of Chemical Engineering, KU Leuven , B-3001 Leuven, Belgium
| | | | - Paula Moldenaers
- Department of Chemical Engineering, KU Leuven , B-3001 Leuven, Belgium
| | | | | |
Collapse
|