1
|
Xu XB, Wang ZH, Xu XN, Fang GY, Gu M. Structural transitions for 2D systems with competing interactions in logarithmic traps. J Chem Phys 2020; 152:054906. [DOI: 10.1063/1.5140816] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/12/2023] Open
Affiliation(s)
- X. B. Xu
- Collaborative Innovation Center of Advanced Microstructures, School of Physics, Nanjing University, Nanjing 210093, People’s Republic of China
| | - Z. H. Wang
- Collaborative Innovation Center of Advanced Microstructures, School of Physics, Nanjing University, Nanjing 210093, People’s Republic of China
| | - X. N. Xu
- Collaborative Innovation Center of Advanced Microstructures, School of Physics, Nanjing University, Nanjing 210093, People’s Republic of China
| | - G. Y. Fang
- Collaborative Innovation Center of Advanced Microstructures, School of Physics, Nanjing University, Nanjing 210093, People’s Republic of China
| | - M. Gu
- Collaborative Innovation Center of Advanced Microstructures, School of Physics, Nanjing University, Nanjing 210093, People’s Republic of China
| |
Collapse
|
2
|
Franzini S, Reatto L, Pini D. Formation of cluster crystals in an ultra-soft potential model on a spherical surface. SOFT MATTER 2018; 14:8724-8739. [PMID: 30357231 DOI: 10.1039/c8sm00997j] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/08/2023]
Abstract
We investigate the formation of cluster crystals with multiply occupied lattice sites on a spherical surface in systems of ultra-soft particles interacting via repulsive, bounded pair potentials. Not all interactions of this kind lead to clustering: we generalize the criterion devised in C. N. Likos et al., Phys. Rev. E, 2001, 63, 031206 to spherical systems in order to distinguish between cluster-forming systems and fluids which display reentrant melting. We use both DFT and Monte Carlo simulations to characterize the behavior of the system, and obtain semi-quantitative agreement between the two. We find that the number of clusters is determined by the ratio between the size σ of the ultra-soft particles and the radius R of the sphere in such a way that each stable configuration spans a certain interval of σ/R. Furthermore, we study the effect of topological frustration on the system due to the sphere curvature by comparing the properties of disclinations, i.e., clusters with fewer than six neighbors, and non-defective clusters. Disclinations are shown to be less stable, contain fewer particles, and be closer to their neighbors than other lattice points: these properties are explained on the basis of geometric and energetic considerations.
Collapse
Affiliation(s)
- Stefano Franzini
- Dipartimento di Fisica "A. Pontremoli", Università di Milano, Via Celoria 16, 20133 Milano, Italy.
| | | | | |
Collapse
|
3
|
Damasceno PF, Glotzer SC, Engel M. Non-close-packed three-dimensional quasicrystals. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2017; 29:234005. [PMID: 28401877 DOI: 10.1088/1361-648x/aa6cc1] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/07/2023]
Abstract
Quasicrystals are frequently encountered in condensed matter. They are important candidates for equilibrium phases from the atomic scale to the nanoscale. Here, we investigate the computational self-assembly of four quasicrystals in a single model system of identical particles interacting with a tunable isotropic pair potential. We reproduce a known icosahedral quasicrystal and report a decagonal quasicrystal, a dodecagonal quasicrystal, and an octagonal quasicrystal. The quasicrystals have low coordination number or occur in systems with mesoscale density variations. We also report a network gel phase.
Collapse
Affiliation(s)
- Pablo F Damasceno
- Applied Physics Program, University of Michigan, Ann Arbor, MI 48109, United States of America. Department of Cellular and Molecular Pharmacology, University of California San Francisco, San Francisco, CA 94158, United States of America
| | | | | |
Collapse
|
4
|
Reichhardt C, Olson Reichhardt CJ. Depinning and nonequilibrium dynamic phases of particle assemblies driven over random and ordered substrates: a review. REPORTS ON PROGRESS IN PHYSICS. PHYSICAL SOCIETY (GREAT BRITAIN) 2017; 80:026501. [PMID: 27997373 DOI: 10.1088/1361-6633/80/2/026501] [Citation(s) in RCA: 89] [Impact Index Per Article: 12.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/06/2023]
Abstract
We review the depinning and nonequilibrium phases of collectively interacting particle systems driven over random or periodic substrates. This type of system is relevant to vortices in type-II superconductors, sliding charge density waves, electron crystals, colloids, stripe and pattern forming systems, and skyrmions, and could also have connections to jamming, glassy behaviors, and active matter. These systems are also ideal for exploring the broader issues of characterizing transient and steady state nonequilibrium flow phases as well as nonequilibrium phase transitions between distinct dynamical phases, analogous to phase transitions between different equilibrium states. We discuss the differences between elastic and plastic depinning on random substrates and the different types of nonequilibrium phases which are associated with specific features in the velocity-force curves, fluctuation spectra, scaling relations, and local or global particle ordering. We describe how these quantities can change depending on the dimension, anisotropy, disorder strength, and the presence of hysteresis. Within the moving phase we discuss how there can be a transition from a liquid-like state to dynamically ordered moving crystal, smectic, or nematic states. Systems with periodic or quasiperiodic substrates can have multiple nonequilibrium second or first order transitions in the moving state between chaotic and coherent phases, and can exhibit hysteresis. We also discuss systems with competing repulsive and attractive interactions, which undergo dynamical transitions into stripes and other complex morphologies when driven over random substrates. Throughout this work we highlight open issues and future directions such as absorbing phase transitions, nonequilibrium work relations, inertia, the role of non-dissipative dynamics such as Magnus effects, and how these results could be extended to the broader issues of plasticity in crystals, amorphous solids, and jamming phenomena.
Collapse
Affiliation(s)
- C Reichhardt
- Theoretical Division, Los Alamos National Laboratory, Los Alamos, NM 87545, USA
| | | |
Collapse
|
5
|
McDermott D, Olson Reichhardt CJ, Reichhardt C. Structural transitions and hysteresis in clump- and stripe-forming systems under dynamic compression. SOFT MATTER 2016; 12:9549-9560. [PMID: 27834430 DOI: 10.1039/c6sm01939k] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/06/2023]
Abstract
Using numerical simulations, we study the dynamical evolution of particles interacting via competing long-range repulsion and short-range attraction in two dimensions. The particles are compressed using a time-dependent quasi-one dimensional trough potential that controls the local density, causing the system to undergo a series of structural phase transitions from a low density clump lattice to stripes, voids, and a high density uniform state. The compression proceeds via slow elastic motion that is interrupted with avalanche-like bursts of activity as the system collapses to progressively higher densities via plastic rearrangements. The plastic events vary in magnitude from small rearrangements of particles, including the formation of quadrupole-like defects, to large-scale vorticity and structural phase transitions. In the dense uniform phase, the system compresses through row reduction transitions mediated by a disorder-order process. We characterize the rearrangement events by measuring changes in the potential energy, the fraction of sixfold coordinated particles, the local density, and the velocity distribution. At high confinements, we find power law scaling of the velocity distribution during row reduction transitions. We observe hysteresis under a reversal of the compression when relatively few plastic rearrangements occur. The decompressing system exhibits distinct phase morphologies, and the phase transitions occur at lower compression forces as the system expands compared to when it is compressed.
Collapse
Affiliation(s)
- Danielle McDermott
- Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA. and Department of Physics, Wabash College, Crawfordsville, Indiana 47933, USA.
| | | | - Charles Reichhardt
- Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA.
| |
Collapse
|
6
|
Spontaneous symmetry breaking in vortex systems with two repulsive lengthscales. Sci Rep 2015; 5:15569. [PMID: 26492969 PMCID: PMC4616066 DOI: 10.1038/srep15569] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/07/2015] [Accepted: 09/24/2015] [Indexed: 11/12/2022] Open
Abstract
Scanning Hall probe microscopy (SHPM) has been used to study vortex structures in thin epitaxial films of the superconductor MgB2. Unusual vortex patterns observed in MgB2 single crystals have previously been attributed to a competition between short-range repulsive and long-range attractive vortex-vortex interactions in this two band superconductor; the type 1.5 superconductivity scenario. Our films have much higher levels of disorder than bulk single crystals and therefore both superconducting condensates are expected to be pushed deep into the type 2 regime with purely repulsive vortex interactions. We observe broken symmetry vortex patterns at low fields in all samples after field-cooling from above Tc. These are consistent with those seen in systems with competing repulsions on disparate length scales, and remarkably similar structures are reproduced in dirty two band Ginzburg-Landau calculations, where the simulation parameters have been defined by experimental observations. This suggests that in our dirty MgB2 films, the symmetry of the vortex structures is broken by the presence of vortex repulsions with two different lengthscales, originating from the two distinct superconducting condensates. This represents an entirely new mechanism for spontaneous symmetry breaking in systems of superconducting vortices, with important implications for pinning phenomena and high current density applications.
Collapse
|
7
|
Edlund E, Lindgren O, Nilsson Jacobi M. A design path for the hierarchical self-assembly of patchy colloidal particles. SOFT MATTER 2015; 11:3913-3919. [PMID: 25869967 DOI: 10.1039/c5sm00596e] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/04/2023]
Abstract
Patchy colloidal particles are promising candidates for building blocks in directed self-assembly. To be successful the surface patterns need to be simple enough to be synthesized, while feature-rich enough to cause the colloidal particles to self-assemble into desired structures. Achieving this is a challenge for traditional synthesis methods. Recently it has been suggested that surface patterns themselves can be made to self-assemble. In this paper we present a design path for the hierarchical targeted self-assembly of patchy colloidal particles based on self-assembling surface patterns. At the level of the surface structure, we use a predictive method utilizing the universality of stripes and spots, coupled with stoichiometric constraints, to cause highly specific and functional patterns to self-assemble on spherical surfaces. We use a minimalistic model of an alkanethiol on gold as a demonstration, showing that even with limited control over the interaction between surface constituents we can obtain patterns that cause the colloidal particles themselves to self-assemble into various complex geometric structures, such as strings, membranes, cubic aggregates and colloidosomes, as well as various crystalline patterns.
Collapse
Affiliation(s)
- E Edlund
- Chalmers University of Technology, Gothenburg.
| | | | | |
Collapse
|
8
|
Ramola K, Damle K, Dhar D. Columnar order and Ashkin-Teller criticality in mixtures of hard squares and dimers. PHYSICAL REVIEW LETTERS 2015; 114:190601. [PMID: 26024157 DOI: 10.1103/physrevlett.114.190601] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/21/2014] [Indexed: 06/04/2023]
Abstract
We show that critical exponents of the transition to columnar order in a mixture of 2×1 dimers and 2×2 hard squares on the square lattice depends on the composition of the mixture in exactly the manner predicted by the theory of Ashkin-Teller criticality, including in the hard-square limit. This result settles the question regarding the nature of the transition in the hard-square lattice gas. It also provides the first example of a polydisperse system whose critical properties depend on composition. Our ideas also lead to some interesting predictions for a class of frustrated quantum magnets that exhibit columnar ordering of the bond energies at low temperature.
Collapse
Affiliation(s)
- Kabir Ramola
- Laboratoire de Physique Théorique et Modèles Statistiques, UMR 8626, Université Paris-Sud 11 and CNRS, Bâtiment 100, Orsay F-91405, France
| | - Kedar Damle
- Tata Institute of Fundamental Research, 1 Homi Bhabha Road, Mumbai, India 400005
| | - Deepak Dhar
- Tata Institute of Fundamental Research, 1 Homi Bhabha Road, Mumbai, India 400005
| |
Collapse
|
9
|
Jin W, Liu Q, Dougherty DB, Cullen WG, Reutt-Robey JE, Weeks J, Robey SW. C60 chain phases on ZnPc/Ag(111) surfaces: Supramolecular organization driven by competing interactions. J Chem Phys 2015; 142:101910. [DOI: 10.1063/1.4906044] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- W. Jin
- Department of Chemistry and Biochemistry, University of Maryland at College Park, College Park, Maryland 20742, USA
| | - Q. Liu
- Institute for Physical Science and Technology, University of Maryland, College Park, Maryland 20742, USA
| | - D. B. Dougherty
- Department of Physics, North Carolina State University, Raleigh, North Carolina 27695, USA
| | - W. G. Cullen
- Department of Physics, University of Maryland at College Park, College Park, Maryland 20742, USA
| | - J. E. Reutt-Robey
- Department of Chemistry and Biochemistry, University of Maryland at College Park, College Park, Maryland 20742, USA
| | - J. Weeks
- Institute for Physical Science and Technology, University of Maryland, College Park, Maryland 20742, USA
| | - S. W. Robey
- National Institute of Standards and Technology, Gaithersburg, Maryland 20878-8372, USA
| |
Collapse
|
10
|
McDermott D, Olson Reichhardt CJ, Reichhardt C. Stripe systems with competing interactions on quasi-one dimensional periodic substrates. SOFT MATTER 2014; 10:6332-6338. [PMID: 25030212 DOI: 10.1039/c4sm01341g] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/03/2023]
Abstract
We numerically examine the two-dimensional ordering of a stripe forming system of particles with competing long-range repulsion and short-range attraction in the presence of a quasi-one-dimensional corrugated substrate. As a function of increasing substrate strength or period we show that a remarkable variety of distinct orderings can be realized, including modulated stripes, prolate clump phases, two dimensional ordered kink structures, crystalline void phases, and smectic phases. Additionally in some cases the stripes align perpendicular to the substrate troughs. Our results suggest that a new route to self assembly for systems with competing interactions can be achieved through the addition of a simple periodic modulated substrate.
Collapse
Affiliation(s)
- Danielle McDermott
- Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA.
| | | | | |
Collapse
|
11
|
Edlund E, Lindgren O, Jacobi MN. Predicting self-assembled patterns on spheres with multicomponent coatings. SOFT MATTER 2014; 10:2955-2960. [PMID: 24695801 DOI: 10.1039/c3sm52827h] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/03/2023]
Abstract
Patchy colloids are promising candidates for building blocks in directed self-assembly, but large scale synthesis of colloids with controlled surface patterns remains challenging. One potential fabrication method is to self-assemble the surface patterns themselves, allowing complex morphologies to organize spontaneously. For this approach to be competitive, prediction and control of the pattern formation process are necessary. However, structure formation in many-body systems is fundamentally hard to understand, and new theoretical methods are needed. Here we present a theory for self-assembling pattern formation in multi-component systems on the surfaces of colloidal particles, formulated as an analytic technique that predicts morphologies directly from the interactions in an effective model. As a demonstration we formulate an isotropic model of alkanethiols on gold, a suggested system for directed self-assembly, and predict its morphologies and transitions as a function of the interaction parameters.
Collapse
Affiliation(s)
- Erik Edlund
- Complex Systems Group, Department of Energy and Environment, Chalmers University of Technology, SE-41296 Göteborg, Sweden.
| | | | | |
Collapse
|
12
|
Bachelard R, Staniscia F. Non-mean-field effects in systems with long-range forces in competition. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2012; 86:051134. [PMID: 23214765 DOI: 10.1103/physreve.86.051134] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/05/2012] [Indexed: 06/01/2023]
Abstract
We investigate the canonical equilibrium of systems with long-range forces in competition. These forces create a modulation in the interaction potential and modulated phases appear at the system scale. The structure of these phases differentiate this system from monotonic potentials, where only the mean-field and disordered phases exist. With increasing temperature, the system switches from one ordered phase to another through a first-order phase transition. Both mean-field and modulated phases may be stable, even at zero temperature, and the long-range nature of the interaction will lead to metastability characterized by extremely long time scales.
Collapse
Affiliation(s)
- R Bachelard
- Instituto de Física de São Carlos, Universidade de São Paulo, 13560-970 São Carlos, São Paulo, Brazil
| | | |
Collapse
|
13
|
Mendoza-Coto A, Stariolo DA. Coarse-grained models of stripe forming systems: phase diagrams, anomalies, and scaling hypothesis. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2012; 86:051130. [PMID: 23214761 DOI: 10.1103/physreve.86.051130] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/21/2012] [Indexed: 06/01/2023]
Abstract
Two coarse-grained models which capture some universal characteristics of stripe forming systems are studied. At high temperatures, the structure factors of both models attain their maxima on a circle in reciprocal space, as a consequence of generic isotropic competing interactions. Although this is known to lead to some universal properties, we show that the phase diagrams have important differences, which are a consequence of the particular k dependence of the fluctuation spectrum in each model. The phase diagrams are computed in a mean field approximation and also after inclusion of small fluctuations, which are shown to modify drastically the mean field behavior. Observables like the modulation length and magnetization profiles are computed for the whole temperature range accessible to both models and some important differences in behavior are observed. A stripe compression modulus is computed, showing an anomalous behavior with temperature as recently reported in related models. Also, a recently proposed scaling hypothesis for modulated systems is tested and found to be valid for both models studied.
Collapse
Affiliation(s)
- Alejandro Mendoza-Coto
- Departamento de Física, Universidade Federal do Rio Grande do Sul, CP 15051, 91501-970 Porto Alegre, Brazil
| | | |
Collapse
|
14
|
Edlund E, Lindgren O, Jacobi MN. Chiral surfaces self-assembling in one-component systems with isotropic interactions. PHYSICAL REVIEW LETTERS 2012; 108:165502. [PMID: 22680733 DOI: 10.1103/physrevlett.108.165502] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/12/2012] [Indexed: 06/01/2023]
Abstract
We show that chiral symmetry can be broken spontaneously in one-component systems with isotropic interactions, i.e., many-particle systems having maximal a priori symmetry. This is achieved by designing isotropic potentials that lead to self-assembly of chiral surfaces. We demonstrate the principle on a simple chiral lattice and on a more complex lattice with chiral supercells. In addition, we show that the complex lattice has interesting melting behavior with multiple morphologically distinct phases that we argue can be qualitatively predicted from the design of the interaction.
Collapse
Affiliation(s)
- E Edlund
- Complex Systems Group, Department of Energy and Environment, Chalmers University of Technology, SE-41296 Göteborg, Sweden
| | | | | |
Collapse
|
15
|
Jin S, Sen A, Sandvik AW. Ashkin-teller criticality and pseudo-first-order behavior in a frustrated Ising model on the square lattice. PHYSICAL REVIEW LETTERS 2012; 108:045702. [PMID: 22400864 DOI: 10.1103/physrevlett.108.045702] [Citation(s) in RCA: 31] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/26/2011] [Indexed: 05/31/2023]
Abstract
We study the challenging thermal phase transition to stripe order in the frustrated square-lattice Ising model with couplings J(1) < 0 (nearest-neighbor, ferromagnetic) and J(2) > 0 (second-neighbor, antiferromagnetic) for g = J(2)/|J(1| > 1/2. Using Monte Carlo simulations and known analytical results, we demonstrate Ashkin-Teller criticality for g ≥ g*; i.e., the critical exponents vary continuously between those of the 4-state Potts model at g = g* and the Ising model for g → ∞. Thus, stripe transitions offer a route to realizing a related class of conformal field theories with conformal charge c = 1 and varying exponents. The transition is first order for g < g* = 0.67 ± 0.01, much lower than previously believed, and exhibits pseudo-first-order behavior for |g* ≤ g </~1.
Collapse
Affiliation(s)
- Songbo Jin
- Department of Physics, Boston University, 590 Commonwealth Avenue, Boston, Massachusetts 02215, USA
| | | | | |
Collapse
|
16
|
Edlund E, Lindgren O, Jacobi MN. Novel self-assembled morphologies from isotropic interactions. PHYSICAL REVIEW LETTERS 2011; 107:085501. [PMID: 21929172 DOI: 10.1103/physrevlett.107.085501] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/20/2011] [Revised: 06/22/2011] [Indexed: 05/31/2023]
Abstract
We present results from particle simulations with isotropic medium range interactions in two dimensions. At low temperature novel types of aggregated structures appear. We show that these structures can be explained by spontaneous symmetry breaking in analytic solutions to an adaptation of the spherical spin model. We predict the critical particle number where the symmetry breaking occurs and show that the resulting phase diagram agrees well with results from particle simulations.
Collapse
Affiliation(s)
- E Edlund
- Complex Systems Group, Department of Energy and Environment, Chalmers University of Technology, SE-41296 Göteborg, Sweden
| | | | | |
Collapse
|
17
|
Edlund E, Lindgren O, Jacobi MN. Designing isotropic interactions for self-assembly of complex lattices. PHYSICAL REVIEW LETTERS 2011; 107:085503. [PMID: 21929174 DOI: 10.1103/physrevlett.107.085503] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/26/2011] [Indexed: 05/31/2023]
Abstract
We present a direct method for solving the inverse problem of designing isotropic potentials that cause self-assembly into target lattices. Each potential is constructed by matching its energy spectrum to the reciprocal representation of the lattice to guarantee that the desired structure is a ground state. We use the method to self-assemble complex lattices not previously achieved with isotropic potentials, such as a snub square tiling and the kagome lattice. The latter is especially interesting because it provides the crucial geometric frustration in several proposed spin liquids.
Collapse
Affiliation(s)
- E Edlund
- Complex Systems Group, Department of Energy and Environment, Chalmers University of Technology, SE-41296 Göteborg, Sweden
| | | | | |
Collapse
|