Analysis of pattern formation in systems with competing range interactions
Identifieur interne : 001713 ( Istex/Corpus ); précédent : 001712; suivant : 001714Analysis of pattern formation in systems with competing range interactions
Auteurs : H J Zhao ; V R Misko ; F M PeetersSource :
- New Journal of Physics [ 1367-2630 ] ; 2012-06.
Abstract
We analyzed pattern formation and identified various morphologies in a system of particles interacting through a non-monotonic potential with a competing range interaction characterized by a repulsive core (rrc), using molecular-dynamics simulations. Depending on parameters, the interaction potential models the inter-particle interaction in various physical systems ranging from atoms, molecules and colloids to vortices in low type-II superconductors and in recently discovered type-1.5 superconductors. We constructed a morphology diagram in the plane critical radius rcdensity n and proposed a new approach to characterizing the different types of patterns. Namely, we elaborated a set of quantitative criteria in order to identify the different pattern types, using the radial distribution function (RDF), the local density function and the occupation factor.
Url:
DOI: 10.1088/1367-2630/14/6/063032
Links to Exploration step
ISTEX:A36699F86F9CD28199495900E74448C936502F68Le document en format XML
<record><TEI wicri:istexFullTextTei="biblStruct"><teiHeader><fileDesc><titleStmt><title>Analysis of pattern formation in systems with competing range interactions</title><author><name sortKey="Zhao, H J" sort="Zhao, H J" uniqKey="Zhao H" first="H J" last="Zhao">H J Zhao</name><affiliation><mods:affiliation>Department of Physics, University of Antwerpen, Groenenborgerlaan 171, B-2020 Antwerpen, Belgium</mods:affiliation></affiliation><affiliation><mods:affiliation>E-mail: haijun.zhao@ua.ac.be</mods:affiliation></affiliation></author><author><name sortKey="Misko, V R" sort="Misko, V R" uniqKey="Misko V" first="V R" last="Misko">V R Misko</name><affiliation><mods:affiliation>Department of Physics, University of Antwerpen, Groenenborgerlaan 171, B-2020 Antwerpen, Belgium</mods:affiliation></affiliation><affiliation><mods:affiliation>E-mail: vyacheslav.misko@ua.ac.be</mods:affiliation></affiliation></author><author><name sortKey="Peeters, F M" sort="Peeters, F M" uniqKey="Peeters F" first="F M" last="Peeters">F M Peeters</name><affiliation><mods:affiliation>Department of Physics, University of Antwerpen, Groenenborgerlaan 171, B-2020 Antwerpen, Belgium</mods:affiliation></affiliation><affiliation><mods:affiliation>E-mail: francois.peeters@ua.ac.be</mods:affiliation></affiliation></author></titleStmt><publicationStmt><idno type="wicri:source">ISTEX</idno><idno type="RBID">ISTEX:A36699F86F9CD28199495900E74448C936502F68</idno><date when="2012" year="2012">2012</date><idno type="doi">10.1088/1367-2630/14/6/063032</idno><idno type="url">https://api.istex.fr/document/A36699F86F9CD28199495900E74448C936502F68/fulltext/pdf</idno><idno type="wicri:Area/Istex/Corpus">001713</idno></publicationStmt><sourceDesc><biblStruct><analytic><title level="a">Analysis of pattern formation in systems with competing range interactions</title><author><name sortKey="Zhao, H J" sort="Zhao, H J" uniqKey="Zhao H" first="H J" last="Zhao">H J Zhao</name><affiliation><mods:affiliation>Department of Physics, University of Antwerpen, Groenenborgerlaan 171, B-2020 Antwerpen, Belgium</mods:affiliation></affiliation><affiliation><mods:affiliation>E-mail: haijun.zhao@ua.ac.be</mods:affiliation></affiliation></author><author><name sortKey="Misko, V R" sort="Misko, V R" uniqKey="Misko V" first="V R" last="Misko">V R Misko</name><affiliation><mods:affiliation>Department of Physics, University of Antwerpen, Groenenborgerlaan 171, B-2020 Antwerpen, Belgium</mods:affiliation></affiliation><affiliation><mods:affiliation>E-mail: vyacheslav.misko@ua.ac.be</mods:affiliation></affiliation></author><author><name sortKey="Peeters, F M" sort="Peeters, F M" uniqKey="Peeters F" first="F M" last="Peeters">F M Peeters</name><affiliation><mods:affiliation>Department of Physics, University of Antwerpen, Groenenborgerlaan 171, B-2020 Antwerpen, Belgium</mods:affiliation></affiliation><affiliation><mods:affiliation>E-mail: francois.peeters@ua.ac.be</mods:affiliation></affiliation></author></analytic><monogr></monogr><series><title level="j">New Journal of Physics</title><idno type="eISSN">1367-2630</idno><imprint><publisher>IOP Publishing</publisher><date type="published" when="2012-06">2012-06</date><biblScope unit="volume">14</biblScope><biblScope unit="issue">6</biblScope></imprint></series><idno type="istex">A36699F86F9CD28199495900E74448C936502F68</idno><idno type="DOI">10.1088/1367-2630/14/6/063032</idno><idno type="href">http://stacks.iop.org/NJP/14/063032</idno><idno type="ArticleID">nj427523</idno></biblStruct></sourceDesc></fileDesc><profileDesc><textClass></textClass><langUsage><language ident="en">en</language></langUsage></profileDesc></teiHeader><front><div type="abstract">We analyzed pattern formation and identified various morphologies in a system of particles interacting through a non-monotonic potential with a competing range interaction characterized by a repulsive core (rrc), using molecular-dynamics simulations. Depending on parameters, the interaction potential models the inter-particle interaction in various physical systems ranging from atoms, molecules and colloids to vortices in low type-II superconductors and in recently discovered type-1.5 superconductors. We constructed a morphology diagram in the plane critical radius rcdensity n and proposed a new approach to characterizing the different types of patterns. Namely, we elaborated a set of quantitative criteria in order to identify the different pattern types, using the radial distribution function (RDF), the local density function and the occupation factor.</div> </front></TEI><istex><corpusName>iop</corpusName><author><json:item><name>H J Zhao</name><affiliations><json:string>Department of Physics, University of Antwerpen, Groenenborgerlaan 171, B-2020 Antwerpen, Belgium</json:string><json:string>E-mail: haijun.zhao@ua.ac.be</json:string></affiliations></json:item><json:item><name>V R Misko</name><affiliations><json:string>Department of Physics, University of Antwerpen, Groenenborgerlaan 171, B-2020 Antwerpen, Belgium</json:string><json:string>E-mail: vyacheslav.misko@ua.ac.be</json:string></affiliations></json:item><json:item><name>F M Peeters</name><affiliations><json:string>Department of Physics, University of Antwerpen, Groenenborgerlaan 171, B-2020 Antwerpen, Belgium</json:string><json:string>E-mail: francois.peeters@ua.ac.be</json:string></affiliations></json:item></author><subject><json:item><lang><json:string>eng</json:string></lang><value>Paper</value></json:item></subject><articleId><json:string>nj427523</json:string></articleId><language><json:string>eng</json:string></language><originalGenre><json:string>research-article</json:string></originalGenre><abstract>We analyzed pattern formation and identified various morphologies in a system of particles interacting through a non-monotonic potential with a competing range interaction characterized by a repulsive core (r>rc) and an attractive tail (r>rc), using molecular-dynamics simulations. Depending on parameters, the interaction potential models the inter-particle interaction in various physical systems ranging from atoms, molecules and colloids to vortices in low type-II superconductors and in recently discovered type-1.5 superconductors. We constructed a morphology diagram in the plane critical radius rcdensity n and proposed a new approach to characterizing the different types of patterns. Namely, we elaborated a set of quantitative criteria in order to identify the different pattern types, using the radial distribution function (RDF), the local density function and the occupation factor.</abstract><qualityIndicators><score>6.976</score><pdfVersion>1.4</pdfVersion><pdfPageSize>595.276 x 841.89 pts (A4)</pdfPageSize><refBibsNative>true</refBibsNative><keywordCount>1</keywordCount><abstractCharCount>896</abstractCharCount><pdfWordCount>8045</pdfWordCount><pdfCharCount>42583</pdfCharCount><pdfPageCount>21</pdfPageCount><abstractWordCount>123</abstractWordCount></qualityIndicators><title>Analysis of pattern formation in systems with competing range interactions</title><genre><json:string>research-article</json:string></genre><host><volume>14</volume><publisherId><json:string>nj</json:string></publisherId><issue>6</issue><genre><json:string>journal</json:string></genre><language><json:string>unknown</json:string></language><eissn><json:string>1367-2630</json:string></eissn><title>New Journal of Physics</title></host><publicationDate>2012</publicationDate><copyrightDate>2012</copyrightDate><doi><json:string>10.1088/1367-2630/14/6/063032</json:string></doi><id>A36699F86F9CD28199495900E74448C936502F68</id><score>0.21893549</score><fulltext><json:item><original>true</original><mimetype>application/pdf</mimetype><extension>pdf</extension><uri>https://api.istex.fr/document/A36699F86F9CD28199495900E74448C936502F68/fulltext/pdf</uri></json:item><json:item><original>false</original><mimetype>application/zip</mimetype><extension>zip</extension><uri>https://api.istex.fr/document/A36699F86F9CD28199495900E74448C936502F68/fulltext/zip</uri></json:item><istex:fulltextTEI uri="https://api.istex.fr/document/A36699F86F9CD28199495900E74448C936502F68/fulltext/tei"><teiHeader><fileDesc><titleStmt><title level="a">Analysis of pattern formation in systems with competing range interactions</title></titleStmt><publicationStmt><authority>ISTEX</authority><publisher>IOP Publishing</publisher><availability><p>IOP Publishing and Deutsche Physikalische Gesellschaft</p></availability><date>2012-06-25</date></publicationStmt><sourceDesc><biblStruct type="inbook"><analytic><title level="a">Analysis of pattern formation in systems with competing range interactions</title><author xml:id="author-1"><persName><forename type="first">H J</forename><surname>Zhao</surname></persName><email>haijun.zhao@ua.ac.be</email><affiliation>Department of Physics, University of Antwerpen, Groenenborgerlaan 171, B-2020 Antwerpen, Belgium</affiliation></author><author xml:id="author-2"><persName><forename type="first">V R</forename><surname>Misko</surname></persName><email>vyacheslav.misko@ua.ac.be</email><affiliation>Department of Physics, University of Antwerpen, Groenenborgerlaan 171, B-2020 Antwerpen, Belgium</affiliation></author><author xml:id="author-3"><persName><forename type="first">F M</forename><surname>Peeters</surname></persName><email>francois.peeters@ua.ac.be</email><affiliation>Department of Physics, University of Antwerpen, Groenenborgerlaan 171, B-2020 Antwerpen, Belgium</affiliation></author></analytic><monogr><title level="j">New Journal of Physics</title><idno type="eISSN">1367-2630</idno><imprint><publisher>IOP Publishing</publisher><date type="published" when="2012-06"></date><biblScope unit="volume">14</biblScope><biblScope unit="issue">6</biblScope></imprint></monogr><idno type="istex">A36699F86F9CD28199495900E74448C936502F68</idno><idno type="DOI">10.1088/1367-2630/14/6/063032</idno><idno type="href">http://stacks.iop.org/NJP/14/063032</idno><idno type="ArticleID">nj427523</idno></biblStruct></sourceDesc></fileDesc><profileDesc><creation><date>2012-06-25</date></creation><langUsage><language ident="en">en</language></langUsage><abstract><p>We analyzed pattern formation and identified various morphologies in a system of particles interacting through a non-monotonic potential with a competing range interaction characterized by a repulsive core (rrc), using molecular-dynamics simulations. Depending on parameters, the interaction potential models the inter-particle interaction in various physical systems ranging from atoms, molecules and colloids to vortices in low type-II superconductors and in recently discovered type-1.5 superconductors. We constructed a morphology diagram in the plane critical radius rcdensity n and proposed a new approach to characterizing the different types of patterns. Namely, we elaborated a set of quantitative criteria in order to identify the different pattern types, using the radial distribution function (RDF), the local density function and the occupation factor.</p> </abstract><textClass><keywords scheme="keyword"><list><head>article-type</head><item><term>Paper</term></item></list></keywords></textClass></profileDesc><revisionDesc><change when="2012-06-25">Created</change><change when="2012-06">Published</change></revisionDesc></teiHeader></istex:fulltextTEI><json:item><original>false</original><mimetype>text/plain</mimetype><extension>txt</extension><uri>https://api.istex.fr/document/A36699F86F9CD28199495900E74448C936502F68/fulltext/txt</uri></json:item></fulltext><metadata><istex:metadataXml wicri:clean="corpus iop not found" wicri:toSee="no header"><istex:xmlDeclaration>version="1.0" encoding="ISO-8859-1" </istex:xmlDeclaration><istex:docType PUBLIC="-//NLM//DTD Journal Publishing DTD v3.0 20080202//EN" URI="http://ej.iop.org/dtd/nlm-3.0/journalpublishing3.dtd" name="istex:docType"></istex:docType><istex:document><article article-type="research-article"><front><journal-meta><journal-id journal-id-type="publisher-id">nj</journal-id><journal-id journal-id-type="coden">NJOPFM</journal-id><journal-title-group><journal-title>New Journal of Physics</journal-title><abbrev-journal-title abbrev-type="IOP">NJP</abbrev-journal-title><abbrev-journal-title abbrev-type="publisher">New J. Phys.</abbrev-journal-title></journal-title-group><issn pub-type="epub">1367-2630</issn><publisher><publisher-name>IOP Publishing</publisher-name></publisher></journal-meta><article-meta><article-id pub-id-type="publisher-id">nj427523</article-id><article-id pub-id-type="doi">10.1088/1367-2630/14/6/063032</article-id><article-id pub-id-type="manuscript">427523</article-id><article-categories><subj-group subj-group-type="article-type"><subject>Paper</subject></subj-group></article-categories><title-group><article-title>Analysis of pattern formation in systems with competing range interactions</article-title><alt-title alt-title-type="ascii">Analysis of pattern formation in systems with competing range interactions</alt-title><alt-title alt-title-type="short">Analysis of pattern formation in systems with competing range interactions</alt-title><alt-title alt-title-type="short-ascii">Analysis of pattern formation in systems with competing range interactions</alt-title></title-group><contrib-group content-type="all"><contrib contrib-type="author"><name><surname>Zhao</surname><given-names>H J</given-names></name><xref ref-type="aff" rid="nj427523af1"></xref><xref ref-type="aff" rid="nj427523em1"></xref></contrib><contrib contrib-type="author"><name><surname>Misko</surname><given-names>V R</given-names></name><xref ref-type="aff" rid="nj427523af1"></xref><xref ref-type="aff" rid="nj427523em2"></xref></contrib><contrib contrib-type="author"><name><surname>Peeters</surname><given-names>F M</given-names></name><xref ref-type="aff" rid="nj427523af1"></xref><xref ref-type="aff" rid="nj427523em3"></xref></contrib><aff id="nj427523af1">Department of Physics, <institution>University of Antwerpen</institution>, Groenenborgerlaan 171, B-2020 Antwerpen, <country>Belgium</country></aff><ext-link ext-link-type="email" id="nj427523em1">haijun.zhao@ua.ac.be</ext-link><ext-link ext-link-type="email" id="nj427523em2">vyacheslav.misko@ua.ac.be</ext-link><ext-link ext-link-type="email" id="nj427523em3">francois.peeters@ua.ac.be</ext-link><author-comment content-type="short-author-list"><p>H J Zhao <italic>et al</italic></p></author-comment></contrib-group><pub-date pub-type="ppub"><month>6</month><year>2012</year></pub-date><pub-date pub-type="epub"><day>25</day><month>6</month><year>2012</year></pub-date><volume>14</volume><issue>6</issue><elocation-id content-type="artnum">063032</elocation-id><history><date date-type="received"><day>23</day><month>3</month><year>2012</year></date></history><permissions><copyright-statement>© IOP Publishing and Deutsche Physikalische Gesellschaft</copyright-statement><copyright-year>2012</copyright-year></permissions><self-uri xlink:href="http://stacks.iop.org/NJP/14/063032"></self-uri><abstract><title>Abstract</title><p>We analyzed pattern formation and identified various morphologies in a system of particles interacting through a non-monotonic potential with a competing range interaction characterized by a repulsive core (<italic>r</italic> < <italic>r</italic><sub>c</sub>) and an attractive tail (<italic>r</italic> > <italic>r</italic><sub>c</sub>), using molecular-dynamics simulations. Depending on parameters, the interaction potential models the inter-particle interaction in various physical systems ranging from atoms, molecules and colloids to vortices in low <italic>κ</italic> type-II superconductors and in recently discovered ‘type-1.5’ superconductors. We constructed a ‘morphology diagram’ in the plane ‘critical radius <italic>r</italic><sub>c</sub>–density <italic>n</italic>’ and proposed a new approach to characterizing the different types of patterns. Namely, we elaborated a set of quantitative criteria in order to identify the different pattern types, using the radial distribution function (RDF), the local density function and the occupation factor.</p></abstract><counts><page-count count="21"></page-count></counts><custom-meta-group><custom-meta><meta-name>ccc</meta-name><meta-value>1367-2630/12/063032+21$33.00</meta-value></custom-meta><custom-meta><meta-name>table-of-contents</meta-name><meta-value>yes</meta-value></custom-meta></custom-meta-group></article-meta></front><body><sec id="nj427523s1"><label>1.</label><title>Introduction</title><p>In a variety of systems pattern formation is governed by competing interactions [<xref ref-type="bibr" rid="nj427523bib01">1</xref>]. Examples of such systems are: Langmuir monolayers [<xref ref-type="bibr" rid="nj427523bib02">2</xref>, <xref ref-type="bibr" rid="nj427523bib03">3</xref>], colloids and gels [<xref ref-type="bibr" rid="nj427523bib04">4</xref>–<xref ref-type="bibr" rid="nj427523bib08">8</xref>], ferrofluids [<xref ref-type="bibr" rid="nj427523bib09">9</xref>–<xref ref-type="bibr" rid="nj427523bib12">12</xref>], magnetic garnet thin films, type-I superconductors [<xref ref-type="bibr" rid="nj427523bib13">13</xref>], the pasta phase in neutron stars [<xref ref-type="bibr" rid="nj427523bib14">14</xref>], etc. In general, there is a strong correlation between pattern formation and inter-particle interaction. Thus, attraction favors aggregation while repulsion favors low local densities. This competition between repulsive and attractive interactions leads to very rich morphologies and phases, such as stripes, clusters, bubbles, etc [<xref ref-type="bibr" rid="nj427523bib15">15</xref>]. Those different patterns were observed in many diverse systems. In ferrofluid systems, rich phases due to the competition between dipolar forces and short-range forces opposing density variations were found experimentally [<xref ref-type="bibr" rid="nj427523bib09">9</xref>] and theoretically [<xref ref-type="bibr" rid="nj427523bib10">10</xref>]. In neutron stars, the competition between the short-range nuclear attraction and the long-range Coulomb repulsion leads to complex pasta phases [<xref ref-type="bibr" rid="nj427523bib14">14</xref>].</p><p>Pattern formation was extensively studied in colloidal systems where the colloid–colloid interaction is characterized by the competition of the hardcore excluded volume interaction, on the one hand, and the polarization of the particles, on the other. The interaction potential can be further controlled by introducing other contributions to the interaction. As a result, rich configurations, such as clusters, repulsive or attractive glassy states and gels, were found numerically [<xref ref-type="bibr" rid="nj427523bib16">16</xref>–<xref ref-type="bibr" rid="nj427523bib21">21</xref>], analytically [<xref ref-type="bibr" rid="nj427523bib22">22</xref>, <xref ref-type="bibr" rid="nj427523bib23">23</xref>] and experimentally [<xref ref-type="bibr" rid="nj427523bib05">5</xref>, <xref ref-type="bibr" rid="nj427523bib07">7</xref>] in colloidal systems with short-range attractive interaction. The properties of isolated clusters formed by the short-range attractive and long-range repulsive interaction were studied recently [<xref ref-type="bibr" rid="nj427523bib24">24</xref>]. By controlling the interaction, the growth of the cluster was shown to change from nearly spherical to one-dimensional (1D) patterns. The 1D growth of the clusters facilitated the collective packing into columnar or lamellar phases [<xref ref-type="bibr" rid="nj427523bib24">24</xref>]. The columnar and lamellar phases in three dimensions were analyzed using molecular dynamics (MD) simulations [<xref ref-type="bibr" rid="nj427523bib25">25</xref>].</p><p>In superconductors, the vortex–vortex interaction is usually considered to be either repulsive (in type-II superconductors, where the Ginzburg–Landau parameter <italic>κ</italic>, i.e. the ratio of the magnetic field penetration depth <italic>λ</italic> to the coherence length <italic>ξ</italic>, is <inline-formula><tex-math></tex-math><inline-graphic xlink:href="nj427523ieqn1.gif"></inline-graphic></inline-formula>) or attractive (in type-I superconductors, where <inline-formula><tex-math></tex-math><inline-graphic xlink:href="nj427523ieqn2.gif"></inline-graphic></inline-formula> and vortices are unstable), while vortices do not interact with each other at the so-called ‘dual point’ when <inline-formula><tex-math></tex-math><inline-graphic xlink:href="nj427523ieqn3.gif"></inline-graphic></inline-formula> [<xref ref-type="bibr" rid="nj427523bib26">26</xref>]. However, a somewhat deeper analysis of the inter-vortex interaction in type-II superconductors near the dual point revealed an attractive tail [<xref ref-type="bibr" rid="nj427523bib27">27</xref>, <xref ref-type="bibr" rid="nj427523bib28">28</xref>]. This repulsive–attractive inter-vortex interaction was used for the interpretation of unusual patterns in the intermediate state in low-<italic>κ</italic> superconductors (e.g. Nb): islands of the Meissner phase surrounded by the vortex phase or vice versa, i.e. vortex clusters surrounded by the Meissner phase [<xref ref-type="bibr" rid="nj427523bib27">27</xref>, <xref ref-type="bibr" rid="nj427523bib28">28</xref>].</p><p>The recent discovery of ‘type-1.5’ superconductors [<xref ref-type="bibr" rid="nj427523bib29">29</xref>] induced a new wave of interest in systems with non-monotonic interactions (see, e.g., [<xref ref-type="bibr" rid="nj427523bib30">30</xref>, <xref ref-type="bibr" rid="nj427523bib31">31</xref>]), due to the fact that the observed vortex patterns in the new superconductors revealed a clear signature of repulsive–attractive inter-vortex interaction. In particular, some typical characteristics of the observed vortex patterns (e.g. vortex lattice with voids and the nearest-neighbor distribution) were explained using a simple model involving a non-monotonic inter-vortex interaction based on a more general approach to multi-order-parameter condensates [<xref ref-type="bibr" rid="nj427523bib32">32</xref>].</p><p>In this paper, we consider a model competing range interaction potential which is repulsive for the short range and attractive for the longer range. In principle, this type of interaction potential could be used as a model for various systems with non-monotonic inter-particle interaction, e.g. atoms or molecules (i.e. the Lennard–Jones potential), or vortices in two-band superconductors, depending on specific parameters of the potential. We investigate the pattern formation in 2D systems for various interaction potential profiles; for example, we distinguish ‘soft-core’ and ‘hard-core’ interactions and analyze the transitions between different types of patterns. Based on this analysis, we construct a ‘morphology diagram’ for various interaction parameters and particle densities. We propose a new approach to characterizing the different morphologies: instead of qualitative characterization of the patterns (e.g. clusters or labyrinths), we introduce a number of quantitative criteria to distinguish these. In particular, the obtained patterns are analyzed in terms of the radial distribution function (RDF) and additional quantities characterizing, e.g., the local density of particles in clusters.</p><p>The paper is organized as follows. In section <xref ref-type="sec" rid="nj427523s2">2</xref>, we describe the model. The pattern formation for different interaction parameters is discussed in section <xref ref-type="sec" rid="nj427523s3">3</xref>. In section <xref ref-type="sec" rid="nj427523s4">4</xref>, we analyze different types of patterns using the RDF and discuss criteria for the identification of different patterns. The conclusions are given in section <xref ref-type="sec" rid="nj427523s5">5</xref>.</p></sec><sec id="nj427523s2"><label>2.</label><title>The model</title><p>The model inter-particle interaction potential is chosen in the following form: <disp-formula id="nj427523eqn1"><label>1</label><tex-math></tex-math><graphic xlink:href="nj427523eqn1.gif"></graphic></disp-formula> Here, <italic>K</italic><sub>0</sub> is the zeroth-order modified Bessel function, <italic>r</italic><sub><italic>ij</italic></sub> = |<bold>r</bold><sub><italic>i</italic></sub> − <bold>r</bold><sub><italic>j</italic></sub>| is the inter-vortex distance, <italic>V</italic><sub>0</sub> and <italic>λ</italic> are the units of energy and length, respectively. These units can be specified, depending on the system. For example, in the case of a type-II superconductor, the appropriate units of length <italic>λ</italic> and energy <italic>V</italic><sub>0</sub>, correspondingly, are the magnetic field penetration depth <italic>λ</italic> and <italic>V</italic><sub>0</sub> = Φ<sup>2</sup><sub>0</sub>/8<italic>π</italic><sup>2</sup><italic>λ</italic><sup>2</sup>, where Φ<sub>0</sub> = <italic>hc</italic>/2<italic>e</italic> (see, e.g., [<xref ref-type="bibr" rid="nj427523bib33">33</xref>]).</p><p>In the dimensionless form, the interaction potential (<xref ref-type="disp-formula" rid="nj427523eqn1">1</xref>) reads as <disp-formula id="nj427523eqn2"><label>2</label><tex-math></tex-math><graphic xlink:href="nj427523eqn2.gif"></graphic></disp-formula> where the dimensionless length is defined as <italic>r</italic>'<sub><italic>ij</italic></sub> = <italic>r</italic><sub><italic>ij</italic></sub>/<italic>λ</italic>. Further on, we will omit the primes and use the dimensionless form of the potential (<xref ref-type="disp-formula" rid="nj427523eqn2">2</xref>). The interaction force is then given by
<disp-formula id="nj427523eqn3"><label>3</label><tex-math></tex-math><graphic xlink:href="nj427523eqn3.gif"></graphic></disp-formula> where <italic>K</italic><sub>1</sub> is the first-order modified Bessel function, <inline-formula><tex-math></tex-math><inline-graphic xlink:href="nj427523ieqn4.gif"></inline-graphic></inline-formula>, and <italic>a</italic> and <italic>b</italic> are two positive coefficients.</p><p>The first term of equation (<xref ref-type="disp-formula" rid="nj427523eqn3">3</xref>) is repulsive, while the second term describes an attractive interaction force. Indeed, for <inline-formula><tex-math></tex-math><inline-graphic xlink:href="nj427523ieqn5.gif"></inline-graphic></inline-formula>. The interaction force (<xref ref-type="disp-formula" rid="nj427523eqn3">3</xref>) has a repulsive (attractive) tail when <italic>b</italic> < 1 (<italic>b</italic> > 1), while for <italic>r</italic> → 0,<italic>K</italic><sub>1</sub>(<italic>x</italic>) → 1/<italic>x</italic> and thus <italic>F</italic><sub><italic>ij</italic></sub> → (<italic>a</italic>/<italic>b</italic> − 1)/<italic>r</italic>. Therefore, for the short range the interaction force (<xref ref-type="disp-formula" rid="nj427523eqn3">3</xref>) is repulsive (attractive) when <italic>a</italic> > <italic>b</italic> (<italic>a</italic> < <italic>b</italic>), and we only consider the case <italic>a</italic> > <italic>b</italic> since an attractive interaction for short distances would result in a collapse of our system of point particles. When <italic>a</italic> > <italic>b</italic> and <italic>b</italic> < 1, the interaction is always repulsive, and particles form a Wigner crystal structure. The most interesting case is realized when <italic>a</italic> > <italic>b</italic> and <italic>b</italic> > 1, and the interaction has a <italic>repulsive core</italic> and <italic>attractive tail</italic> <xref ref-type="fn" rid="nj427523fn1"><sup>1</sup></xref><fn id="nj427523fn1"><label>1</label><p>Note that we do not characterize the interaction force as ‘short-range repulsive and long-range attractive’ since both the contributions are rather ‘short-range’, due to the Bessel functions. The fact that the interaction is short-range allows us to use a cut-off for long distances.</p></fn>. In this case, there exists a critical distance <italic>r</italic><sub>c</sub>, where the inter-particle interaction energy (<xref ref-type="disp-formula" rid="nj427523eqn1">1</xref>) reaches a minimum (and the interaction force (<xref ref-type="disp-formula" rid="nj427523eqn3">3</xref>) changes sign). By setting the force equal to zero, the coefficient <italic>a</italic> is given by
<disp-formula id="nj427523eqn4"><label>4</label><tex-math></tex-math><graphic xlink:href="nj427523eqn4.gif"></graphic></disp-formula> The pattern formation is determined by the coefficients <italic>b</italic>, <italic>r</italic><sub>c</sub> and the particle density <italic>n</italic>.</p><p>We study pattern formation in a system of interacting point-like particles using the Langevin equations. The dynamics of a single particle <italic>i</italic> obeys the following overdamped equation of motion [<xref ref-type="bibr" rid="nj427523bib16">16</xref>–<xref ref-type="bibr" rid="nj427523bib18">18</xref>, <xref ref-type="bibr" rid="nj427523bib33">33</xref>]:
<disp-formula id="nj427523eqn5"><label>5</label><tex-math></tex-math><graphic xlink:href="nj427523eqn5.gif"></graphic></disp-formula> Here, <italic>η</italic> is the viscosity, which is set to unity, <bold>F</bold><sub><italic>ij</italic></sub> is the inter-particle interaction force defined by equation (<xref ref-type="disp-formula" rid="nj427523eqn3">3</xref>) and <bold>F</bold><sup>p</sup><sub><italic>i</italic></sub>, given by
<disp-formula id="nj427523eqn6"><label>6</label><tex-math></tex-math><graphic xlink:href="nj427523eqn6.gif"></graphic></disp-formula> is the pinning force (although often neglected, in the general case we introduced random pinning sites, see section <xref ref-type="sec" rid="nj427523s3">3</xref>). Here, <italic>N</italic><sub>p</sub> is the number of pinning sites, <italic>f</italic><sub>p</sub> is the maximum pinning force of each potential well, <italic>r</italic><sub>p</sub> is the range of the pinning potential, Θ is the Heaviside step function, and <inline-formula><tex-math></tex-math><inline-graphic xlink:href="nj427523ieqn6.gif"></inline-graphic></inline-formula>. The thermal stochastic term <bold>F</bold><sup>T</sup><sub><italic>i</italic></sub> in equation (<xref ref-type="disp-formula" rid="nj427523eqn5">5</xref>) obeys the following conditions: <disp-formula id="nj427523eqn7"><label>7</label><tex-math></tex-math><graphic xlink:href="nj427523eqn7.gif"></graphic></disp-formula> and <disp-formula id="nj427523eqn8"><label>8</label><tex-math></tex-math><graphic xlink:href="nj427523eqn8.gif"></graphic></disp-formula></p><p>We consider a 2D square simulation region <italic>L</italic><sub><italic>x</italic></sub> × <italic>L</italic><sub><italic>y</italic></sub> in the <italic>xy</italic>-plane and apply periodic boundary conditions in the <italic>x</italic>- and <italic>y</italic>-directions. For the interaction force given by equation (<xref ref-type="disp-formula" rid="nj427523eqn3">3</xref>), which exponentially decays for large distances, we use a cut-off for <italic>r</italic> > 8. Note that in this way the interaction in our system is finite-range. The length of the square cell <italic>L</italic><sub><italic>x</italic></sub> = <italic>L</italic><sub><italic>y</italic></sub> = <italic>L</italic> has been varied from 60 to 180 to examine the finite-size effects, and we set <italic>L</italic> = 120 to optimize the calculations, without appreciable influence on the results. To obtain stable particle patterns, we performed simulated annealing simulations of interacting particles. For this purpose, particles were initially randomly distributed inside the simulation region at some suitable non-zero temperature. Then the temperature was gradually reduced to zero, and the simulation was continued until the total force acting on any single particle became much smaller than typical forces in the system.</p></sec><sec id="nj427523s3"><label>3.</label><title>Pattern formation</title><p>In this section, we study pattern formation and identify different morphologies depending on parameters of the interaction <italic>r</italic><sub>c</sub>, <italic>b</italic> and the particle density <italic>n</italic>. In figure <xref ref-type="fig" rid="nj427523fig1">1</xref>, we illustrate the change of the interaction force profile due to the increase of <italic>r</italic><sub>c</sub> and <italic>b</italic>. The interaction force is very sensitive to <italic>r</italic><sub>c</sub>, since its increase directly leads to an increase in the repulsive part in equation (<xref ref-type="disp-formula" rid="nj427523eqn3">3</xref>) (see the main panel of figure <xref ref-type="fig" rid="nj427523fig1">1</xref>). On the other hand, an increase of <italic>b</italic> leads to an increase of both the repulsive and attractive parts (see the inset of figure <xref ref-type="fig" rid="nj427523fig1">1</xref>). However, the repulsive interaction increases much faster for <italic>r</italic> < <italic>r</italic><sub>c</sub> than the attractive interaction for <italic>r</italic> > <italic>r</italic><sub>c</sub>.</p><fig id="nj427523fig1" position="float"><label>Figure 1.</label><caption id="nj427523fc1"><p>The profile of the model inter-particle interaction force, given by equation (<xref ref-type="disp-formula" rid="nj427523eqn2">2</xref>), versus the distance. The main panel and the inset panel show the change in the force profile due to an increase of <italic>r</italic><sub>c</sub> and <italic>b</italic>, respectively.</p></caption><graphic id="nj427523f1_eps" content-type="print" xlink:href="nj427523f1_pr.eps"></graphic><graphic id="nj427523f1_online" content-type="online" xlink:href="nj427523f1_online.jpg"></graphic></fig><p>Thus, increasing <italic>b</italic> facilitates stabilization of the minimum inter-particle distance at <italic>r</italic><sub>c</sub>. In other words, an increase of <italic>b</italic> results in hardening of the core in the interaction force (<xref ref-type="disp-formula" rid="nj427523eqn3">3</xref>), i.e. the interaction changes continuously from a ‘soft-core’ to a ‘hard-core’ regime.</p><sec id="nj427523s3-1"><label>3.1.</label><title>Soft-core interaction</title><p>In this subsection, we analyze pattern formation in the absence of additional pinning. The influence of random pinning on the pattern formation will be discussed in the next subsection.</p><p>To study pattern formation in the soft-core regime, we set the coefficient <italic>b</italic> = 1.1. We define the critical density as <inline-formula><tex-math></tex-math><inline-graphic xlink:href="nj427523ieqn7.gif"></inline-graphic></inline-formula>, which is the density of an ideal (hexagonal) Wigner crystal with the lattice constant <italic>a</italic> = <italic>r</italic><sub>c</sub>. The density is defined as <italic>n</italic> = <italic>N</italic>/<italic>S</italic>, where <italic>S</italic> = 120 × 120 is the area of the simulation region and <italic>N</italic> is the number of particles (further on, in our analysis of various patterns, we will refer either to the number of particles <italic>N</italic> in the simulation cell or to the density which is <italic>n</italic> = <italic>N</italic>/14 400). For the case of <italic>n</italic> < <italic>n</italic><sub>c</sub> which is considered here, the Wigner crystal is not stable due to the attractive interaction. In figure <xref ref-type="fig" rid="nj427523fig2">2</xref>, we present patterns formed by <italic>N</italic> = 1500 particles when <italic>r</italic><sub>c</sub> increases from 1.9 (a) to 2.9 (f). Note that the condition <italic>n</italic> < <italic>n</italic><sub>c</sub> is always fulfilled for all the values of <italic>r</italic><sub>c</sub> in this range. We found that for <italic>r</italic><sub>c</sub> < 2.1 particles form clusters similar to the formation of the ‘clump’ phase found in colloids and 2D electrons [<xref ref-type="bibr" rid="nj427523bib01">1</xref>, <xref ref-type="bibr" rid="nj427523bib16">16</xref>–<xref ref-type="bibr" rid="nj427523bib18">18</xref>] (see figures <xref ref-type="fig" rid="nj427523fig2">2</xref>(a) and (b)). The main difference to the patterns found in [<xref ref-type="bibr" rid="nj427523bib16">16</xref>–<xref ref-type="bibr" rid="nj427523bib18">18</xref>] is that the relatively softer core in our case is compressed due to the attractive interaction, and the clusters acquire a circular shape. In addition, the interaction between the clusters (decaying exponentially for long distances within the interaction range) becomes negligible for inter-cluster distance of the order of a few to several <italic>r</italic><sub>0</sub>. Therefore, the clusters can be considered as non-interacting for low densities (although they still do not approach each other), contrary to the situation described in [<xref ref-type="bibr" rid="nj427523bib16">16</xref>, <xref ref-type="bibr" rid="nj427523bib18">18</xref>] where a super-lattice is formed due to the long-range cluster–cluster repulsion. As <italic>r</italic><sub>c</sub> increases, the clusters expand. In particular, for <italic>r</italic><sub>c</sub> > 2.1, the clusters start elongating, which is an indication of the instability with respect to the transition to stripes. For 2.1 < <italic>r</italic><sub>c</sub> < 2.3, a mixed state with both stripes and clusters is observed (see figure <xref ref-type="fig" rid="nj427523fig2">2</xref>(c)). Further increase of <italic>r</italic><sub>c</sub> gradually destroys the clusters and leads to the formation of labyrinths (see figures <xref ref-type="fig" rid="nj427523fig2">2</xref>(d)–(f)).</p><fig id="nj427523fig2" position="float"><label>Figure 2.</label><caption id="nj427523fc2"><p>Patterns for <italic>N</italic> = 1500 particles in a unit cell <italic>L</italic> × <italic>L</italic> with <italic>L</italic> = 120 for varying <italic>r</italic><sub>c</sub>: <italic>r</italic><sub>c</sub> = 1.9 (a), 2.1 (b), 2.3 (c), 2.5 (d), 2.7 (e) and 2.9 (f).</p></caption><graphic id="nj427523f2_eps" content-type="print" xlink:href="nj427523f2_pr.eps"></graphic><graphic id="nj427523f2_online" content-type="online" xlink:href="nj427523f2_online.jpg"></graphic></fig><p>In order to reveal the influence of the density on the pattern formation, we gradually increase the number of particles from 1500 to 10 500 in our simulation cell. First, in figure <xref ref-type="fig" rid="nj427523fig3">3</xref>, we present patterns formed by <italic>N</italic> = 5500 particles for varying <italic>r</italic><sub>c</sub>. Note that <italic>n</italic> > <italic>n</italic><sub>c</sub> for <italic>N</italic> = 5500. As compared to the lower density case shown in figure <xref ref-type="fig" rid="nj427523fig2">2</xref>, in the case of clusters, the additional particles lead to the expansion of the clusters (see figures <xref ref-type="fig" rid="nj427523fig3">3</xref>(a) and (b)). A mixture of clusters and stripes is formed when the additional particles form bridges connecting the clusters (see figure <xref ref-type="fig" rid="nj427523fig3">3</xref>(c)), which will be discussed in detail below. For labyrinths (shown in figures <xref ref-type="fig" rid="nj427523fig2">2</xref>(d)–(f) for <italic>N</italic> = 1500), the additional particles fill the empty regions (voids), resulting in the formation of the hexagonal lattice with varying local density (figure <xref ref-type="fig" rid="nj427523fig3">3</xref>(d)) and, finally, a regular hexagonal lattice (see figures <xref ref-type="fig" rid="nj427523fig3">3</xref>(e) and (f)). Note that the lattice with varying local density (figure <xref ref-type="fig" rid="nj427523fig3">3</xref>(d)) is not stable in systems with pure repulsive interaction such as vortices in type-II superconductors.</p><fig id="nj427523fig3" position="float"><label>Figure 3.</label><caption id="nj427523fc3"><p>Patterns for <italic>N</italic> = 5500 particles in a unit cell <italic>L</italic> × <italic>L</italic> with <italic>L</italic> = 120 for varying <italic>r</italic><sub>c</sub>: <italic>r</italic><sub>c</sub> = 1.9 (a), 2.1 (b), 2.3 (c), 2.5 (d), 2.7 (e) and 2.9 (f).</p></caption><graphic id="nj427523f3_eps" content-type="print" xlink:href="nj427523f3_pr.eps"></graphic><graphic id="nj427523f3_online" content-type="online" xlink:href="nj427523f3_online.jpg"></graphic></fig><p>We analyzed in detail the intermediate regime (i.e. corresponding to the transition from clusters to stripes, see figure <xref ref-type="fig" rid="nj427523fig2">2</xref>(c)) where <italic>r</italic><sub>c</sub> ≈ 2.3. The resulting patterns for varying density are shown in figure <xref ref-type="fig" rid="nj427523fig4">4</xref>. Figure <xref ref-type="fig" rid="nj427523fig4">4</xref>(a) displays a configuration at low density with <italic>N</italic> = 1500 when many individual clusters are formed. With increasing the density (see figure <xref ref-type="fig" rid="nj427523fig4">4</xref>(b)), the clusters connect to each other and form stripes. A further increase of the number of particles (figure <xref ref-type="fig" rid="nj427523fig4">4</xref>(c)) results in the formation of a complex mixture of interconnected stripes with voids. A very interesting and counter-intuitive evolution is observed when the number of particles increases from <italic>N</italic> = 5500 to <italic>N</italic> = 6500 (see figure <xref ref-type="fig" rid="nj427523fig4">4</xref>(d)): in contrast to the gradual transition from a void-rich configuration to a lattice-rich configuration when <italic>N</italic> changes from <italic>N</italic> = 3500 (b) to <italic>N</italic> = 5500 (c), in the case of <italic>N</italic> = 6500 we observe ‘re-entrant’ behavior, i.e. the void-rich configuration starts recovering, which is compensated for by an increase of the local density in the stripes. However, further increasing the density results in the expansion of the stripes to the empty regions, which is accompanied by a decrease of the local density in the stripes, as shown in figure <xref ref-type="fig" rid="nj427523fig4">4</xref>(e). The distribution of particles becomes more uniform, with only a few small voids. Finally, for <italic>N</italic> = 9500 (figure <xref ref-type="fig" rid="nj427523fig4">4</xref>(f)), we obtain a deformed hexagonal lattice characterized by a varying local density with only one small void.</p><fig id="nj427523fig4" position="float"><label>Figure 4.</label><caption id="nj427523fc4"><p>Patterns for <italic>r</italic><sub>c</sub> = 2.3 and varying number of particles <italic>N</italic> (in the computational unit cell <italic>L</italic> × <italic>L</italic> with <italic>L</italic> = 120): <italic>N</italic> = 1500 (a), 3500 (b), 5500 (c), 6500 (d), 7500 (e) and 9500 (f).</p></caption><graphic id="nj427523f4_eps" content-type="print" xlink:href="nj427523f4_pr.eps"></graphic><graphic id="nj427523f4_online" content-type="online" xlink:href="nj427523f4_online.jpg"></graphic></fig><p>Our calculations show that the obtained patterns are very sensitive to variations in <italic>r</italic><sub>c</sub>. Thus, if <italic>r</italic><sub>c</sub> slightly decreases (e.g. <italic>r</italic><sub>c</sub> = 2.25), the number of clusters greatly increases as compared to the case <italic>r</italic><sub>c</sub> = 2.3, for the same density of particles. For <italic>r</italic><sub>c</sub> = 2.25, we also observe the transition from a void-rich configuration to a lattice-rich configuration. However, since the decrease of <italic>r</italic><sub>c</sub> increases the attractive component of the inter-particle interaction this occurs at a much higher density. For even smaller values of <italic>r</italic><sub>c</sub>, i.e. <italic>r</italic><sub>c</sub> < 2.1, we do not observe a lattice pattern, even for an extremely large number of particles (up to <italic>N</italic> = 20 000).</p></sec><sec id="nj427523s3-2"><label>3.2.</label><title>Analysis of the stability of the patterns</title><p>To further study the stability of the patterns, here we analyze the influence of the cut-off procedure and the annealing time used in our simulations. Alternatively, we consider the effect of additional weak pinning and modify the inter-particle interaction by introducing a weak repulsive tail.</p><p>First, we consider the case when clusters are formed (i.e. <italic>r</italic><sub>c</sub> < 2.3). We choose <italic>r</italic><sub>c</sub> = 1.9 and increase the cut-off distance from 8 to 16<italic>λ</italic>. We also increase the annealing time from typically (10 000–20 000)Δ<italic>t</italic> to 200 000Δ<italic>t</italic>, where Δ<italic>t</italic> = 0.04 is the time step of the MD simulation. In figure <xref ref-type="fig" rid="nj427523fig5">5</xref>, we plot the evolution of the energy with time, and we show some representative snapshots in the insets. We start from a random configuration of particles (see figure <xref ref-type="fig" rid="nj427523fig5">5</xref>(b)). During the annealing, the energy first drops rapidly due to the clusterization of the particles. After a relatively short time, small clusters start to merge and form larger clusters (see figure <xref ref-type="fig" rid="nj427523fig5">5</xref>(c)). The average inter-cluster distance increases when the total number of clusters decreases. The merging becomes more difficult and considerably slows down when the clusters are well separated. As a result, short-living metastable states are formed (see figure <xref ref-type="fig" rid="nj427523fig5">5</xref>(d)). However, during the merging of the clusters, the lifetime of the metastable states increases very rapidly, and gradually the short-living metastable states evolve into <italic>long-living metastable states</italic> when the clusters are well separated (see figure <xref ref-type="fig" rid="nj427523fig5">5</xref>(e)). Comparing with the relatively faster annealing case (i.e. 20 000Δ<italic>t</italic>), we see that the slower annealing process results in a decrease of the number of clusters inside the simulation cell, which increases the lifetime of the metastable states. However, even for the shorter annealing times (i.e. 20 000Δ<italic>t</italic>), we arrive at long-living metastable configurations which can have either slightly higher energy or even the same energy as those states obtained during the long-time annealing (note that the energy of the configuration does not change over the time interval from ≈15 000Δ<italic>t</italic> to 50 000Δ<italic>t</italic>). It is worth noting that the <italic>morphology</italic> is similar to that found in the ‘fast’ annealing process.</p><fig id="nj427523fig5" position="float"><label>Figure 5.</label><caption id="nj427523fc5"><p>The evolution of the energy <italic>E</italic> versus time <italic>t</italic> (a); snapshots corresponding to various annealing times (b–e). Two merging clusters are shown in red (gray) in (d).</p></caption><graphic id="nj427523f5_eps" content-type="print" xlink:href="nj427523f5_pr.eps"></graphic><graphic id="nj427523f5_online" content-type="online" xlink:href="nj427523f5_online.jpg"></graphic></fig><p>Thermodynamically stable states can be obtained in the presence of weak random pinning. To verify this, we introduce weak random pinning sites in the system, with the maximum pinning force <italic>f</italic><sub>p</sub> = 0.2 and radius <italic>r</italic><sub>p</sub> = 0.3 and the density of random pinning sites being the same as the density of the particles. We found that such a weak pinning produces no appreciable impact on stripes, labyrinths or lattices. In the case of clusters, the bonding energy due to the inter-particle interaction within every individual cluster is strong, and the weak pinning is not sufficient to destroy the structure of the clusters. However, the inter-cluster interaction is weak, therefore, those clusters can be easily stabilized by weak pinning. Thus, comparing to the patterns obtained without pinning, we find that slightly smaller clusters are formed in the presence of pinning, and the merging of clusters is prevented even for two relatively close but small neighbor clusters (see figure <xref ref-type="fig" rid="nj427523fig6">6</xref>). Note that in the absence of pinning, clusters have nearly circular shape while their shape is less regular in the case of weak pinning. In the latter case, some of the particles were pinned individually, i.e. outside the clusters, see figure <xref ref-type="fig" rid="nj427523fig6">6</xref>(b).</p><fig id="nj427523fig6" position="float"><label>Figure 6.</label><caption id="nj427523fc6"><p>Calculated patterns for <italic>r</italic><sub>c</sub> = 2.1,<italic>N</italic> = 4500 and the annealing time 50 000Δ<italic>t</italic>: without pinning (a) and in the presence of weak random pinning with the maximum pinning force <italic>f</italic><sub>p</sub> = 0.2 and radius <italic>r</italic><sub>p</sub> = 0.3 (b).</p></caption><graphic id="nj427523f6_eps" content-type="print" xlink:href="nj427523f6_pr.eps"></graphic><graphic id="nj427523f6_online" content-type="online" xlink:href="nj427523f6_online.jpg"></graphic></fig><p>An alternative way to stabilize clusters (and approach the thermodynamically stable configurations) consists in using a modified force with an additional weak <italic>repulsive tail</italic>, e.g. taken in the form: <bold>F</bold>'<sub><italic>ij</italic></sub> = <bold>F</bold><sub><italic>ij</italic></sub> cosh(<italic>r</italic><sub>e</sub> − <italic>r</italic>), where <italic>r</italic><sub>e</sub> is a parameter chosen as <italic>r</italic><sub>e</sub> = 10, and the cut-off of the interaction is set at 16<italic>λ</italic>. Then the resulting modified force coincides with that given by equation (<xref ref-type="disp-formula" rid="nj427523eqn3">3</xref>) for <italic>r</italic> < 8 but the weak attractive tail smoothly evolves to repulsive when <italic>r</italic> > 10. The modified interaction force leads to very similar patterns as the original one (given by equation (<xref ref-type="disp-formula" rid="nj427523eqn3">3</xref>)) because the weak repulsive tail is not efficient for the clusters to form super-structures (see, e.g., [<xref ref-type="bibr" rid="nj427523bib16">16</xref>, <xref ref-type="bibr" rid="nj427523bib18">18</xref>, <xref ref-type="bibr" rid="nj427523bib21">21</xref>]), but is sufficient for stabilizing the configurations of particles. When we increase the inter-cluster <italic>repulsive</italic> interaction by increasing the particle density, the resulting patterns are similar to the previously calculated long-living metastable states. The reason for this behavior is that the radius of a cluster is insensitive to the particle numbers within the cluster. As a result, the inter-cluster repulsion increases very slowly as the density increases.</p></sec><sec id="nj427523s3-3"><label>3.3.</label><title>Morphology diagram</title><p>Based on the above analysis of the patterns, we constructed a ‘morphology diagram’ (as we explained above, the calculated patterns are long-living metastable states, i.e. in the absence of additional pinning or the long-range repulsive tail—rather than real thermodynamic phases; thus the diagram presents various morphologies) in the plane of <italic>r</italic><sub>c</sub> and the number density <italic>n</italic> (see figure <xref ref-type="fig" rid="nj427523fig7">7</xref>). For extremely low density (<italic>n</italic> < 0.03), particles form small clusters that are well separated, and the patterns are rather insensitive to variations in <italic>r</italic><sub>c</sub>. However, morphologies become richer when density <italic>n</italic> increases (i.e. <italic>n</italic> > 0.1). Low values of <italic>r</italic><sub>c</sub> (i.e. <italic>r</italic><sub>c</sub> < 2.1) still favor cluster formation for a broad range of densities, although the particle density in clusters increases for large <italic>n</italic>. With increasing <italic>r</italic><sub>c</sub>, clusters become unstable with respect to the formation of the ‘bridges’ between separate clusters, which is a precursor of the formation of stripes. Stripes are formed in a rather narrow range of <italic>r</italic><sub>c</sub> when <italic>r</italic><sub>c</sub> > 2.1 (see figure <xref ref-type="fig" rid="nj427523fig7">7</xref>). The stripe patterns can be divided into two sub-groups, I and II, i.e. void-rich (see figure <xref ref-type="fig" rid="nj427523fig7">7</xref>(i)) and void-poor (lattice-rich) stripes (figure <xref ref-type="fig" rid="nj427523fig7">7</xref>(j)). For larger <italic>r</italic><sub>c</sub> and <italic>n</italic> < <italic>n</italic><sub>c</sub>, particles form labyrinth structures. However, when increasing the density, additional particles fill the empty regions and finally they form a deformed hexagonal lattice with varying density. Depending on <italic>r</italic><sub>c</sub> and <italic>n</italic> < <italic>n</italic><sub>c</sub>, deformed lattice is characterized by a varying local density or by the appearance of voids. Correspondingly, we distinguish two regions (1 and 2 in figure <xref ref-type="fig" rid="nj427523fig7">7</xref>). Note that in the vicinity of the morphology boundaries, patterns are always mixtures of the two morphologies (e.g. clusters and stripes) except the transition from the stripes to the deformed hexagonal lattice. Near this morphology boundary, particles form either stripes with high local density or deformed hexagonal lattices with lower local density. The probability of finding the lattice configuration greatly increases when <italic>r</italic><sub>c</sub> or <italic>n</italic> increases. Finally, for very high density, when the average inter-particle distance becomes smaller than <italic>r</italic><sub>c</sub> the repulsive interaction prevails, resulting in the formation of an almost perfect hexagonal lattice.</p><fig id="nj427523fig7" position="float"><label>Figure 7.</label><caption id="nj427523fc7"><p>The ‘morphology diagram’ in the plane ‘critical radius <italic>r</italic><sub>c</sub>–density <italic>n</italic>’ (n) and representative patterns (a–m). For <italic>r</italic><sub>c</sub> < 2.1 and extremely low density (<italic>n</italic> < 0.03), particles form small, well-separated clusters. For larger <italic>r</italic><sub>c</sub>, the size of the clusters increases (see the change from (d) to (e) and from (b) to (c)). Increase of the density results in increasing size of the clusters (see the change from (d) to (b) and from (e) to (c)). Further increase of <italic>r</italic><sub>c</sub> leads to elongation of the clusters and the formation of ‘bridges’ between them (a). Thus, the configurations change gradually from clusters to stripes (i). At larger <italic>r</italic><sub>c</sub> or <italic>n</italic>, stripes interconnect and form patterns with voids (j). Thus the stripes are divided into two sub-groups, I and II. When <italic>r</italic><sub>c</sub> becomes even larger, labyrinth-like configurations are formed for <italic>n</italic> < <italic>n</italic><sub>c</sub> (see (g) and (f)). With further increase of <italic>r</italic><sub>c</sub> or <italic>n</italic>, deformed hexagonal lattices with (k) or without voids (l) are formed. The deformation of a hexagonal lattice is reduced for even larger values of <italic>r</italic><sub>c</sub> and <italic>n</italic> (m).</p></caption><graphic id="nj427523f7_eps" content-type="print" xlink:href="nj427523f7_pr.eps"></graphic><graphic id="nj427523f7_online" content-type="online" xlink:href="nj427523f7_online.jpg"></graphic></fig><p>The stability of clusters for varying parameter of the inter-particle interaction, <italic>r</italic><sub>c</sub>, can be examined by applying an additional weak <italic>confinement potential</italic> and calculating the energy corresponding to the varying number of particles in the system. For this purpose, we use a parabolic potential, similar to <italic>F</italic><sup><italic>i</italic></sup><sub>p</sub> in equation (<xref ref-type="disp-formula" rid="nj427523eqn6">6</xref>), to confine all the particles and gradually decrease <italic>f</italic><sub>p</sub> to zero during the annealing. This allows us to obtain one large cluster for various <italic>r</italic><sub>c</sub>. Then we study the stability of these induced clusters, by analyzing the average energy per particle <italic>U</italic> for various <italic>r</italic><sub>c</sub> and <italic>N</italic><sub>c</sub>, where <italic>N</italic><sub>c</sub> is the total number of particles composing the cluster. We found that for <italic>r</italic><sub>c</sub> < 2.3, the clusters are stable: <italic>U</italic> is approximately proportional to <italic>N</italic><sup>2</sup><sub>c</sub> and decreases while <italic>N</italic><sub>c</sub> increases (see figure <xref ref-type="fig" rid="nj427523fig8">8</xref>(a)). The size of the cluster grows with the number of constituent particles: it goes to infinity for <italic>N</italic><sub>c</sub> → ∞, which corresponds to <italic>phase separation</italic>. Note that this result could be understood from a comparison of the model potential (equation (<xref ref-type="disp-formula" rid="nj427523eqn2">2</xref>)) with the Lennard-Jones potential: for small <italic>r</italic><sub>c</sub>, both the potentials become very similar, in particular, for <italic>r</italic><sub>c</sub> ≈ 1. On the other hand, it is well known that the Lennard-Jones potential leads to phase separation.</p><fig id="nj427523fig8" position="float"><label>Figure 8.</label><caption id="nj427523fc8"><p>The potential energy per particle <italic>U</italic> versus the total number of particles within the cluster, <italic>N</italic><sub>c</sub>, for <italic>r</italic><sub>c</sub> = 1.8 (a), 2.3 (b) and 2.5 (c).</p></caption><graphic id="nj427523f8_eps" content-type="print" xlink:href="nj427523f8_pr.eps"></graphic><graphic id="nj427523f8_online" content-type="online" xlink:href="nj427523f8_online.jpg"></graphic></fig><p>For <italic>r</italic><sub>c</sub> = 2.3, which is a typical value for stripe formation, there is a saddle point in the function <italic>U</italic> versus <italic>N</italic><sub>c</sub>, which shows that the clusters have an <italic>optimal</italic> size, e.g. <italic>N</italic><sub>c</sub> ≈ 280 (see figure <xref ref-type="fig" rid="nj427523fig8">8</xref>(b)). Finally, for <italic>r</italic><sub>c</sub> > 2.3,<italic>U</italic> changes sign when <italic>N</italic><sub>c</sub> > 7, resulting in a positive energy, indicating that the clusters in this case are overheated by the parabolic potential and therefore are not stable. Labyrinths are more favorable in this case.</p><p>The patterns presented in figure <xref ref-type="fig" rid="nj427523fig7">7</xref> are found in many real systems. For example, clusters are found in such systems as colloids [<xref ref-type="bibr" rid="nj427523bib04">4</xref>, <xref ref-type="bibr" rid="nj427523bib06">6</xref>, <xref ref-type="bibr" rid="nj427523bib19">19</xref>, <xref ref-type="bibr" rid="nj427523bib20">20</xref>] and 2D electrons [<xref ref-type="bibr" rid="nj427523bib16">16</xref>–<xref ref-type="bibr" rid="nj427523bib18">18</xref>]. The obtained stripe patters are very similar to those in Langmuir monolayers [<xref ref-type="bibr" rid="nj427523bib03">3</xref>].</p><p>Furthermore, several of the calculated morphologies were found in superconductors, e.g. the stripe phase in the mixed state of type-I superconductors, clusters of the Meissner phase or vortex clusters in the intermediate state in low-<italic>κ</italic> type-II superconductors [<xref ref-type="bibr" rid="nj427523bib27">27</xref>, <xref ref-type="bibr" rid="nj427523bib28">28</xref>], the labyrinth phase (i.e. vortex lattice with voids) or vortex clusters in type-1.5 superconductors [<xref ref-type="bibr" rid="nj427523bib29">29</xref>]. In spite of the variety of physical systems and length scales ranging from nano- and micro-objects to cosmic objects, the common feature of all these systems is a competing attractive–repulsive inter-particle interaction, which allows analyzing them within the same approach.</p></sec><sec id="nj427523s3-4"><label>3.4.</label><title>Hard-core interaction</title><p>Let us now analyze the influence of the <italic>b</italic> parameter on the pattern formation. As mentioned above (see figure <xref ref-type="fig" rid="nj427523fig1">1</xref>), an increase in parameter <italic>b</italic> changes the potential from soft-core to hard-core. Hardening the repulsive core in the interaction potential generally leads to a decreasing compressibility of the inner parts of the patterns where particles are closely packed. Clearly, in this case the patterns change as compared to the soft-core regime.</p><p>As shown in figure <xref ref-type="fig" rid="nj427523fig9">9</xref>, for <italic>b</italic> = 4 and low density (<italic>N</italic> = 1500 per simulation region), all the patterns are clusters of different shapes. Thus for <italic>r</italic><sub>c</sub> = 1.15 and <italic>r</italic><sub>c</sub> = 1.3, the clusters are of a circular shape similar to that found in the soft-core regime, although for <italic>r</italic><sub>c</sub> = 1.3 some clusters, composed of a small number of particles, have a polygonal shape. With increasing <italic>r</italic><sub>c</sub>, polygonal-shaped clusters become more favorable, which allows us to identify them as ‘short stripes’. For even larger <italic>r</italic><sub>c</sub>, the clusters become much larger. They represent separate islands of the hexagonal lattice.</p><fig id="nj427523fig9" position="float"><label>Figure 9.</label><caption id="nj427523fc9"><p>Patterns for <italic>N</italic> = 1500 particles in a unit cell <italic>L</italic> × <italic>L</italic> with <italic>L</italic> = 120,
<italic>b</italic> = 4 and for varying <italic>r</italic><sub>c</sub>: <italic>r</italic><sub>c</sub> = 1.15 (a), 1.3 (b), 1.45 (c), 1.6 (d), 1.75 (e) and 1.9 (f).</p></caption><graphic id="nj427523f9_eps" content-type="print" xlink:href="nj427523f9_pr.eps"></graphic><graphic id="nj427523f9_online" content-type="online" xlink:href="nj427523f9_online.jpg"></graphic></fig><p>For higher density (<italic>N</italic> = 6500), the variety of patterns is much richer. Although for <italic>r</italic><sub>c</sub> = 1.15 circular-shaped clusters are still observed (see figure <xref ref-type="fig" rid="nj427523fig10">10</xref>(a)), which grow in size for <italic>r</italic><sub>c</sub> = 1.3 and slightly deviate in shape from circular (figure <xref ref-type="fig" rid="nj427523fig10">10</xref>(b)), for larger <italic>r</italic><sub>c</sub> = 1.45, the deviations of the cluster shape from circular become more pronounced (see figure <xref ref-type="fig" rid="nj427523fig10">10</xref>(c)). With further increasing <italic>r</italic><sub>c</sub>, i.e. <italic>r</italic><sub>c</sub> = 1.6, stripes start to form (see figure <xref ref-type="fig" rid="nj427523fig10">10</xref>(d)) followed by lattices with voids for <italic>r</italic><sub>c</sub> = 1.75 and <italic>r</italic><sub>c</sub> = 1.9. In the hard-core regime, deformations of lattices occur via the appearance of voids instead of varying local density (as happens in the soft-core regime).</p><fig id="nj427523fig10" position="float"><label>Figure 10.</label><caption id="nj427523fc10"><p>Patterns for <italic>N</italic> = 6500 particles in a unit cell <italic>L</italic> × <italic>L</italic> with <italic>L</italic> = 120,<italic>b</italic> = 4 and for varying <italic>r</italic><sub>c</sub>: <italic>r</italic><sub>c</sub> = 1.15 (a), 1.3 (b), 1.45 (c), 1.6 (d), 1.75 (e) and 1.9 (f).</p></caption><graphic id="nj427523f10_eps" content-type="print" xlink:href="nj427523f10_pr.eps"></graphic><graphic id="nj427523f10_online" content-type="online" xlink:href="nj427523f10_online.jpg"></graphic></fig><p>The obtained patterns presented in figures <xref ref-type="fig" rid="nj427523fig9">9</xref> and <xref ref-type="fig" rid="nj427523fig10">10</xref> resemble the nanoscale domain patterns formed by Cu and Pb atoms on a surface of Cu [<xref ref-type="bibr" rid="nj427523bib34">34</xref>, <xref ref-type="bibr" rid="nj427523bib35">35</xref>]. In particular, the lattice islands (figures <xref ref-type="fig" rid="nj427523fig9">9</xref>(d)–(f)), the stripes (figure <xref ref-type="fig" rid="nj427523fig10">10</xref>(d)) and void-rich configurations (figure <xref ref-type="fig" rid="nj427523fig10">10</xref>(e)) are similar to the patterns evolving from islands to stripes and ‘inverted’ islands formed by Cu and Pb atoms arranged in a single atomic layer on a Cu(111) surface [<xref ref-type="bibr" rid="nj427523bib34">34</xref>, <xref ref-type="bibr" rid="nj427523bib35">35</xref>], with increasing atomic coverage.</p></sec></sec><sec id="nj427523s4"><label>4.</label><title>Analysis of the patterns</title><sec id="nj427523s4-1"><label>4.1.</label><title>Radial distribution function</title><p>Although the different patterns and morphologies studied above are qualitatively well distinguished, it is highly desirable to build a set of solid criteria that would allow us to identify them in a quantitative manner. For this purpose, we analyze here different morphologies by means of the RDF. The RDF, <italic>g</italic>(<italic>r</italic>), describes the variation of the atomic (particle) density as a function of the distance from one particular atom (particle). If we define an average density as <italic>n</italic> = <italic>N</italic>/<italic>V</italic> (where <italic>V</italic> is the volume or surface area in the 2D case), then the local density at distance <italic>r</italic> is <italic>ng</italic>(<italic>r</italic>). The knowledge of RDF is important since one can measure <italic>g</italic>(<italic>r</italic>) experimentally using digital video microscopy [<xref ref-type="bibr" rid="nj427523bib36">36</xref>]. Moreover, macroscopic thermodynamic quantities can be calculated using <italic>g</italic>(<italic>r</italic>) in thermodynamics [<xref ref-type="bibr" rid="nj427523bib37">37</xref>].</p><p>In our calculations, we define the RDF <italic>g</italic><sub><italic>i</italic></sub>(<italic>r</italic>) as follows:
<disp-formula id="nj427523eqn9"><label>9</label><tex-math></tex-math><graphic xlink:href="nj427523eqn9.gif"></graphic></disp-formula> Here, the lower index indicates that the RDF centers at the position of the <italic>i</italic>th particle; △<italic>N</italic> is the number of particles whose distance to the <italic>i</italic>th particles is between <italic>r</italic> and <italic>r</italic> + △<italic>r</italic>. The average <italic>g</italic>(<italic>r</italic>) is given by <disp-formula id="nj427523eqn10"><label>10</label><tex-math></tex-math><graphic xlink:href="nj427523eqn10.gif"></graphic></disp-formula></p><p>In figure <xref ref-type="fig" rid="nj427523fig11">11</xref>(a), we plot the function <italic>g</italic>(<italic>r</italic>) calculated for the low-density cluster configuration shown in figure <xref ref-type="fig" rid="nj427523fig2">2</xref>(b). The function <italic>g</italic>(<italic>r</italic>) has two well-pronounced peaks. The first peak corresponds to the average distance to the first coordination sphere (nearest neighbors), while the second one is located at a distance approximately twice the distance to the first peak, which shows short-range periodicity. For <italic>r</italic> larger than the size of the cluster and smaller than the inter-cluster distance, <italic>g</italic>(<italic>r</italic>) < 1, since only a few particles are located inside this range. The decreasing tail in the RDF shows strong aggregation of the particles forming clusters. The minimum of <italic>g</italic>(<italic>r</italic>) is very close to zero, indicating that most of the clusters have circular symmetry and are well separated. The position of the minimum of <italic>g</italic>(<italic>r</italic>) (marked by the gray arrow in figure <xref ref-type="fig" rid="nj427523fig11">11</xref>(a)) gives an estimate of the average diameter of the clusters.</p><fig id="nj427523fig11" position="float"><label>Figure 11.</label><caption id="nj427523fc11"><p>The average RDF of the patterns formed at low density (<italic>N</italic> = 1500). Panels (a), (b) and (c) correspond to the patterns shown in figures <xref ref-type="fig" rid="nj427523fig2">2</xref>(b), (d) and (f), respectively. The peaks marked by dark arrows show the short-range periodicity. The gray arrows show the minimum of the RDF.</p></caption><graphic id="nj427523f11_eps" content-type="print" xlink:href="nj427523f11_pr.eps"></graphic><graphic id="nj427523f11_online" content-type="online" xlink:href="nj427523f11_online.jpg"></graphic></fig><p>For stripes (see figure <xref ref-type="fig" rid="nj427523fig11">11</xref>(b) for the pattern shown in figure <xref ref-type="fig" rid="nj427523fig2">2</xref>(d)), there are also two peaks indicating the short-range periodicity, similar to clusters. However, unlike the case of clusters, the minimum of <italic>g</italic>(<italic>r</italic>) is non-zero since the stripes are generally continuous. For labyrinths (see figure <xref ref-type="fig" rid="nj427523fig11">11</xref>(c) for the pattern shown in figure <xref ref-type="fig" rid="nj427523fig2">2</xref>(f)), only short-range periodic ordering exists, and <italic>g</italic>(<italic>r</italic>) becomes uniform for larger <italic>r</italic>.</p><p>In figure <xref ref-type="fig" rid="nj427523fig12">12</xref>, we plot <italic>g</italic>(<italic>r</italic>) for the patterns formed at high density (<italic>N</italic> = 5500). For clusters, increase in the density leads to an increasing variation of the cluster size. Note, however, that the sizes of individual clusters are rather insensitive to moderate variations of the number of particles in the clusters, due to the strong compressibility of the core. As a result, the function <italic>g</italic>(<italic>r</italic>) for the high-density clusters misses the short-range periodicity and the second peak disappears (see figure <xref ref-type="fig" rid="nj427523fig12">12</xref>(a)). Strikingly, this peak still exists for the void-rich patterns (see figure <xref ref-type="fig" rid="nj427523fig12">12</xref>(b)). The decrease of the tail of <italic>g</italic>(<italic>r</italic>) becomes very slow, which shows a minor aggregation of the particles.</p><fig id="nj427523fig12" position="float"><label>Figure 12.</label><caption id="nj427523fc12"><p>The average RDF of the patterns formed at high density (<italic>N</italic> = 5500). Panels (a), (b) and (c) correspond to the patterns shown in figures <xref ref-type="fig" rid="nj427523fig3">3</xref>(b), (d) and (f), respectively. The arrows have the same meaning as in figure <xref ref-type="fig" rid="nj427523fig11">11</xref>.</p></caption><graphic id="nj427523f12_eps" content-type="print" xlink:href="nj427523f12_pr.eps"></graphic><graphic id="nj427523f12_online" content-type="online" xlink:href="nj427523f12_online.jpg"></graphic></fig><p>For lattices (see figure <xref ref-type="fig" rid="nj427523fig12">12</xref>(b)), the position of the first peak is <inline-formula><tex-math></tex-math><inline-graphic xlink:href="nj427523ieqn8.gif"></inline-graphic></inline-formula>, where <italic>a</italic> is the distance between two neighboring particles in the ideal hexagonal lattice with the density <italic>n</italic>. The positions of the second and third peaks are at <inline-formula><tex-math></tex-math><inline-graphic xlink:href="nj427523ieqn9.gif"></inline-graphic></inline-formula> and at <italic>r</italic><sub>2</sub> ≈ 2<italic>a</italic>, respectively, corresponding to those in a hexagonal lattice. The fourth, fifth and even sixth peaks appear, which show that this structure is much more ordered. However, due to the variation of the local density, these peaks are inhomogeneously broadened: they are actually a combination of many neighboring peaks. These peaks become more pronounced for larger <italic>r</italic><sub>c</sub> or <italic>n</italic>, which implies that the lattices become more regular.</p><p>It is interesting to compare the results for the function <italic>g</italic>(<italic>r</italic>) for hard-core interaction with those for soft-core interaction. Although the clusters shown in figures <xref ref-type="fig" rid="nj427523fig10">10</xref>(a) and (b) have shapes very similar to those of the clusters shown in figures <xref ref-type="fig" rid="nj427523fig3">3</xref>(a) and (b), the analysis using the RDF shows that most of the clusters shown in figure <xref ref-type="fig" rid="nj427523fig10">10</xref>(b) have hexagonal ordered cores. The peaks in <italic>g</italic>(<italic>r</italic>) appear to be much better separated than, e.g., for the deformed hexagonal lattice formed in the soft-core case (see figure <xref ref-type="fig" rid="nj427523fig12">12</xref>(c)). For the stripe and void-rich patterns, the RDF shows a much clearer hexagonal lattice ordering.</p><p>Therefore, the analysis using the RDF allowed us to establish quantitative criteria for the different patterns and reveal the differences in the structure of the patterns (which look similar, e.g. clusters) in the case of soft- and hard-core inter-particle interactions (figure <xref ref-type="fig" rid="nj427523fig13">13</xref>).</p><fig id="nj427523fig13" position="float"><label>Figure 13.</label><caption id="nj427523fc13"><p>The average RDF of the patterns formed in the hard-core case. Panels (a), (b), (c), (d) and (e) correspond to patterns (a), (b), (c), (d) and (e) in figure <xref ref-type="fig" rid="nj427523fig10">10</xref>, respectively.</p></caption><graphic id="nj427523f13_eps" content-type="print" xlink:href="nj427523f13_pr.eps"></graphic><graphic id="nj427523f13_online" content-type="online" xlink:href="nj427523f13_online.jpg"></graphic></fig></sec><sec id="nj427523s4-2"><label>4.2.</label><title>Local density</title><p>Using the RDF, let us define another useful quantity, the ‘order parameter’, to characterize various types of patterns. By definition, the local density is <disp-formula id="nj427523eqn11"><label>11</label><tex-math></tex-math><graphic xlink:href="nj427523eqn11.gif"></graphic></disp-formula> Here, <italic>N</italic><sub><italic>ξ</italic></sub> ≈ <italic>π</italic><italic>ξ</italic><sup>2</sup><italic>n</italic> is the average number of particles within the circle centered at the <italic>i</italic>th particle with radius <italic>R</italic> = <italic>ξ</italic>. In an ideal hexagonal lattice, one particle has six nearest neighbors, therefore <italic>N</italic><sub><italic>ξ</italic></sub> = 7 and <inline-formula><tex-math></tex-math><inline-graphic xlink:href="nj427523ieqn10.gif"></inline-graphic></inline-formula>. Thus, for any configuration characterized by a small local density fluctuation, the average local density <italic>I</italic> = ⟨<italic>I</italic><sub><italic>i</italic></sub>⟩ = 6. From the definition given by equation (<xref ref-type="disp-formula" rid="nj427523eqn11">11</xref>), we can see that the presence of a large fraction of empty regions can result in a considerable increase of <italic>I</italic> since only a ‘shell’ of those regions of thickness <italic>ξ</italic> is considered. In figure <xref ref-type="fig" rid="nj427523fig14">14</xref>(a), we plot the function <italic>I</italic> for different <italic>r</italic><sub>c</sub> and <italic>N</italic> for the soft-core interaction with <italic>b</italic> = 1.1. We see that the function <italic>I</italic> can be used to characterize the different types of patterns. Thus we found that, for clusters, <italic>I</italic> is always large (<italic>I</italic> > 20) since aggregation is strong. For stripes, aggregation is weaker, and thus <italic>I</italic> is smaller (10 < <italic>I</italic> < 17). For labyrinths, the regions free of particles are relatively small. Therefore, <italic>I</italic> ranges between 6 and 9. Finally, for a hexagonal lattice, <italic>I</italic> is always 6. Therefore, the function <italic>I</italic> serves as a measure of the aggregation of particles and allows us to effectively distinguish the different types of patterns.</p><fig id="nj427523fig14" position="float"><label>Figure 14.</label><caption id="nj427523fc14"><p>The average local density <italic>I</italic> versus the particle density <italic>n</italic> for <italic>b</italic> = 1.1 (soft-core) (a). Curves A, B, C, D, E and F correspond to <italic>r</italic><sub>c</sub> = 1.9,2.1,2.3,2.5,2.7 and 2.9, respectively. For clusters, <italic>I</italic> > 20, for stripes, 10 < <italic>I</italic> < 20, for labyrinths, 6 < <italic>I</italic> < 10;<italic>I</italic> = 6 for lattice and deformed lattice configurations. The local density <italic>I</italic> versus the particle density <italic>n</italic> for <italic>b</italic> = 4 (hard-core) (b). Curves A, B, C, D, E and F correspond to <italic>r</italic><sub>c</sub> = 1.15,1.3,1.45,1.6,1.75 and 1.9, respectively.</p></caption><graphic id="nj427523f14_eps" content-type="print" xlink:href="nj427523f14_pr.eps"></graphic><graphic id="nj427523f14_online" content-type="online" xlink:href="nj427523f14_online.jpg"></graphic></fig><p>The function <italic>I</italic> also provides a tool for analyzing the stability of patterns with increasing density. This is demonstrated in figure <xref ref-type="fig" rid="nj427523fig14">14</xref>(b), where we plot <italic>I</italic> for the hard-core case when <italic>b</italic> = 4. For <italic>r</italic><sub>c</sub> = 1.15,<italic>I</italic> > 20 for all the values of the density, and thus the clusters are stable. However, the situation changes for <italic>r</italic><sub>c</sub> = 1.3: while clusters are formed for low densities up to <italic>N</italic> ≈ 7500, for larger density <italic>I</italic> decreases below the value <italic>I</italic> = 20, which means that clusters elongate and interconnect, giving rise to the onset of stripe formation. For <italic>r</italic><sub>c</sub> = 1.45, the particles start to form stripes at even lower density (<italic>N</italic> = 3500). For higher density, the additional particles fill in the empty regions and finally they form a regular lattice. For <italic>r</italic><sub>c</sub> > 1.6, the island-like clusters formed at low densities are very unstable with respect to an increase in density. The additional particles rapidly occupy the empty regions, and the patterns change from clusters to stripes and from stripes to lattice with increasing density.</p><p>Therefore, the phenomenological description of pattern evolution with increasing density we revealed above in section <xref ref-type="sec" rid="nj427523s3-1">3.1</xref> has been verified in a quantitative manner using the solid basis of the local density function approach.</p></sec><sec id="nj427523s4-3"><label>4.3.</label><title>Occupation factor</title><p>We demonstrated that the (average) RDF <italic>g</italic>(<italic>r</italic>) (equation (<xref ref-type="disp-formula" rid="nj427523eqn10">10</xref>)) and the local density function <italic>I</italic> (equation (<xref ref-type="disp-formula" rid="nj427523eqn11">11</xref>)) allowed us to unambiguously characterize various types of patterns. At the same time, it is still hard to distinguish, using the above tools, between labyrinths (or a lattice with voids) and lattices, especially in the soft-core regime when the corresponding RDF do not differ much from each other, and both types of patterns are characterized by small values of <italic>I</italic>. However, these can be easily distinguished by employing another additional criterion, namely, the occupation factor which characterizes the fraction of the space occupied by particles to the total space.</p><p>Let us assume that all the patterns are formed by islands of hexagonal lattice with the average distance between two nearest-neighbor particles <italic>r</italic> = <italic>r</italic><sub>1</sub>, where <italic>r</italic><sub>1</sub> is the position of the first peak of the RDF. Then the ratio of the area occupied by the particles and the whole simulation area is <italic>A</italic> = (<italic>r</italic><sub>1</sub>/<italic>a</italic>)<sup>2</sup>. As shown in figure <xref ref-type="fig" rid="nj427523fig15">15</xref>(a), this ratio can be used to distinguish between labyrinths and lattices, since for lattices <italic>A</italic> ≈ 1.</p><fig id="nj427523fig15" position="float"><label>Figure 15.</label><caption id="nj427523fc15"><p>The occupation factor <italic>A</italic> = (<italic>r</italic><sub>1</sub>/<italic>a</italic>)<sup>2</sup> for <italic>b</italic> = 1.1 (soft-core) (a). Curves A, B, C, D, E and F correspond to <italic>r</italic><sub>c</sub> = 1.9,2.1,2.3,2.5,2.7 and 2.9, respectively. The occupation factor <italic>A</italic> for <italic>b</italic> = 4 (hard-core) (b). Curves A, B, C, D, E and F correspond to <italic>r</italic><sub>c</sub> = 1.15,1.3,1.45,1.6,1.75 and 1.9, respectively.</p></caption><graphic id="nj427523f15_eps" content-type="print" xlink:href="nj427523f15_pr.eps"></graphic><graphic id="nj427523f15_online" content-type="online" xlink:href="nj427523f15_online.jpg"></graphic></fig><p>We also found that circular clusters in the case of soft-core interaction are stable with respect to an increase in the density. In this case, the additional particles cannot efficiently increase the occupation factor due to the increase in the local density of the core. However, for stripes and labyrinths, the occupation factor <italic>A</italic> increases with the density (see figure <xref ref-type="fig" rid="nj427523fig15">15</xref>(a)). These patterns are not stable for high density. Thus the large jump in <italic>A</italic>(<italic>N</italic>) at <italic>N</italic> = 6500 in curve C in figure <xref ref-type="fig" rid="nj427523fig15">15</xref>(a) shows the transition from stripes to lattice.</p><p>As one can expect, in the case of hard-core interaction, the occupation factor increases nearly linearly versus the density (see figure <xref ref-type="fig" rid="nj427523fig15">15</xref>(b)). Therefore, the polygon shaped and island-like clusters are not stable: with increasing density, they evolve to the lattice configuration (the plateau in <italic>A</italic>(<italic>N</italic>) in figure <xref ref-type="fig" rid="nj427523fig15">15</xref>(b)).</p><p>In addition, let us introduce another useful quantity which characterizes the degree of ‘perfection’ of a lattice. Let us define the particles with inter-particle distance shorter than <italic>ξ</italic> as neighboring particles. Then <disp-formula id="nj427523eqn12"><label>12</label><tex-math></tex-math><graphic xlink:href="nj427523eqn12.gif"></graphic></disp-formula> is the average displacement of the inter-distance between two neighbor particles (measured in units of <italic>a</italic>), which is independent of the density. The function <italic>ϵ</italic> is non-zero for a deformed hexagonal lattice and zero for the ideal hexagonal lattice. Thus, <italic>ϵ</italic> quantitatively measures the degree of perfection of a lattice. Note that <italic>ϵ</italic> is only used for lattices as an auxiliary tool to distinguish ordered lattices from less ordered ones (see table <xref ref-type="table" rid="nj427523t1">1</xref>).</p><table-wrap id="nj427523t1" position="float"><label>Table 1.</label><caption id="nj427523tc1"><p>The set of criteria used to quantitatively identify various morphologies in terms of the RDF <italic>g</italic>(<italic>r</italic>), the local density function <italic>I</italic>, the occupation factor <italic>A</italic> and the parameter <italic>ϵ</italic> (for details, see the text).</p></caption><table frame="hsides"><colgroup><col align="left"></col><col align="center"></col><col align="center"></col><col align="center"></col><col align="center"></col></colgroup><thead><tr><th>Patterns</th><th><italic>g</italic>(<italic>r</italic>)</th><th><italic>I</italic></th><th><italic>A</italic></th><th><italic>ϵ</italic></th></tr></thead><tbody><tr><td>Clusters</td><td>Peak at <italic>r</italic><sub>1</sub>(<italic>r</italic><sub>2</sub>),</td><td>>20</td><td></td><td></td></tr><tr><td></td><td><inline-formula><tex-math></tex-math><inline-graphic xlink:href="nj427523ieqn11.gif"></inline-graphic></inline-formula></td><td></td><td></td><td></td></tr><tr><td>Stripes</td><td>Peak at <italic>r</italic><sub>1</sub>(<italic>r</italic><sub>2</sub>),</td><td>10–20</td><td></td><td></td></tr><tr><td></td><td><inline-formula><tex-math></tex-math><inline-graphic xlink:href="nj427523ieqn12.gif"></inline-graphic></inline-formula></td><td></td><td></td><td></td></tr><tr><td>Labyrinths</td><td>Peak at <italic>r</italic><sub>1</sub>(<italic>r</italic><sub>2</sub>),</td><td>6–10</td><td><1</td><td></td></tr><tr><td></td><td><italic>g</italic>(<italic>r</italic>)≈ const, <italic>r</italic>><italic>r</italic><sub>2</sub></td><td></td><td></td><td></td></tr><tr><td>Lattice</td><td>Several peaks: <italic>r</italic><sub>1</sub>,<italic>r</italic><sub>2</sub>, …</td><td>≈6</td><td>≈1</td><td>0</td></tr><tr><td>Deformed</td><td>Several peaks: <italic>r</italic><sub>1</sub>,<italic>r</italic><sub>2</sub>, …</td><td>≈6</td><td>≈1</td><td>>0</td></tr><tr><td>lattice</td><td></td><td></td><td></td><td></td></tr></tbody></table></table-wrap></sec></sec><sec id="nj427523s5"><label>5.</label><title>Conclusions</title><p>Using MD simulations, we analyzed the pattern formation and identified different morphologies in a system of particles interacting via a non-monotonic competing range potential, with a repulsive short-range part and an attractive long-range part. The form of the interacting potential is rather general: it describes, depending on the specific parameters, the inter-particle interaction in a variety of physical systems ranging from, e.g., atoms and molecules (Lennard-Jones potential) to colloids and vortices in superconductors. In particular, it can be used as a model of inter-vortex interaction in low <italic>κ</italic> type-II superconductors and in the recently discovered so-called ‘type-1.5’ superconductors.</p><p>Since the non-monotonic range potential involves interaction frustration, the obtained patterns are in general long-living metastable states rather than equilibrium configurations. We demonstrated that thermodynamically stable patterns can be obtained if we either use a modified potential with an additional repulsive tail or introduce weak pinning.</p><p>The obtained morphologies were summarized in a ‘morphology diagram’ in the plane ‘critical radius <italic>r</italic><sub>c</sub>–density <italic>n</italic>’ (<italic>r</italic><sub>c</sub> is the critical radius where the interaction force changes its sign). We also analyzed the influence of the hardness of the ‘core’, i.e. the strength of the repulsive core part of the interaction potential, on the pattern formation.</p><p>We developed a set of criteria in order to unambiguously identify the obtained types of patterns using the following approaches: (i) the RDF <italic>g</italic>(<italic>r</italic>), (ii) the local density function <italic>I</italic> and (iii) the occupation factor <italic>A</italic>. In addition, we introduced a parameter which characterizes the degree of perfection of a lattice <italic>ϵ</italic>. Employing these approaches, we elaborated a set of criteria for the identification of the different morphologies summarized in table <xref ref-type="table" rid="nj427523t1">1</xref>.</p></sec></body><back><ack><title>Acknowledgments</title><p>We acknowledge useful discussions with Ernst Helmut Brandt, Charles Reichhardt and Cynthia Olson Reichhardt. This work was supported by the ‘Odysseus’ Program of the Flemish Government and the Flemish Science Foundation (FWO-Vl).</p></ack><ref-list content-type="numerical"><title>References</title><ref id="nj427523bib01"><label>1</label><element-citation publication-type="book"><person-group person-group-type="author"><name><surname>Vedmedenko</surname><given-names>E Y</given-names></name></person-group><year>2007</year><source>Competing Interactions and Patterns in Nanoworld</source><publisher-loc>Berlin</publisher-loc><publisher-name>Wiley-VCH</publisher-name></element-citation></ref><ref id="nj427523bib02"><label>2</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Kaganer</surname><given-names>V M</given-names></name><name><surname>Möhwald</surname><given-names>H</given-names></name><name><surname>Dutta</surname><given-names>P</given-names></name></person-group><year>1999</year><source>Rev. Mod. Phys.</source><volume>71</volume><fpage>779</fpage><lpage>819</lpage><page-range>779–819</page-range><pub-id pub-id-type="doi">10.1103/RevModPhys.71.779</pub-id></element-citation></ref><ref id="nj427523bib03"><label>3</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Suresh</surname><given-names>K A</given-names></name><name><surname>Nittmann</surname><given-names>J</given-names></name><name><surname>Rondelez</surname><given-names>F</given-names></name></person-group><year>1988</year><source>Europhys. Lett.</source><volume>6</volume><fpage>437</fpage><pub-id pub-id-type="doi">10.1209/0295-5075/6/5/011</pub-id></element-citation></ref><ref id="nj427523bib04"><label>4</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Hoffmann</surname><given-names>N</given-names></name><name><surname>Ebert</surname><given-names>F</given-names></name><name><surname>Likos</surname><given-names>C N</given-names></name><name><surname>Löwen</surname><given-names>H</given-names></name><name><surname>Maret</surname><given-names>G</given-names></name></person-group><year>2006</year><source>Phys. Rev. Lett.</source><volume>97</volume><elocation-id content-type="artnum">078301</elocation-id><pub-id pub-id-type="doi">10.1103/PhysRevLett.97.078301</pub-id></element-citation></ref><ref id="nj427523bib05"><label>5</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Lu</surname><given-names>P J</given-names></name><name><surname>Zaccarelli</surname><given-names>E</given-names></name><name><surname>Ciulla</surname><given-names>F</given-names></name><name><surname>Schofield</surname><given-names>A B</given-names></name><name><surname>Sciortino</surname><given-names>F</given-names></name><name><surname>Weitz</surname><given-names>D A</given-names></name></person-group><year>2008</year><source>Nature</source><volume>453</volume><fpage>499</fpage><pub-id pub-id-type="doi">10.1038/nature06931</pub-id></element-citation></ref><ref id="nj427523bib06"><label>6</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Campbell</surname><given-names>A I</given-names></name><name><surname>Anderson</surname><given-names>V J</given-names></name><name><surname>van Duijneveldt</surname><given-names>J S</given-names></name><name><surname>Bartlett</surname><given-names>P</given-names></name></person-group><year>2005</year><source>Phys. Rev. Lett.</source><volume>94</volume><elocation-id content-type="artnum">208301</elocation-id><pub-id pub-id-type="doi">10.1103/PhysRevLett.94.208301</pub-id></element-citation></ref><ref id="nj427523bib07"><label>7</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Klix</surname><given-names>C L</given-names></name><name><surname>Royall</surname><given-names>C P</given-names></name><name><surname>Tanaka</surname><given-names>H</given-names></name></person-group><year>2010</year><source>Phys. Rev. Lett.</source><volume>104</volume><elocation-id content-type="artnum">165702</elocation-id><pub-id pub-id-type="doi">10.1103/PhysRevLett.104.165702</pub-id></element-citation></ref><ref id="nj427523bib08"><label>8</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Munarin</surname><given-names>F F</given-names></name><name><surname>Nelissen</surname><given-names>K</given-names></name><name><surname>Ferreira</surname><given-names>W P</given-names></name><name><surname>Farias</surname><given-names>G A</given-names></name><name><surname>Peeters</surname><given-names>F M</given-names></name></person-group><year>2008</year><source>Phys. Rev. <named-content content-type="jnl-part">E</named-content></source><volume>77</volume><elocation-id content-type="artnum">031608</elocation-id><pub-id pub-id-type="other">http://pre.aps.org/abstract/PRE/v77/i3/e031608</pub-id></element-citation></ref><ref id="nj427523bib09"><label>9</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Islam</surname><given-names>M F</given-names></name><name><surname>Lin</surname><given-names>K H</given-names></name><name><surname>Lacoste</surname><given-names>D</given-names></name><name><surname>Lubensky</surname><given-names>T C</given-names></name><name><surname>Yodh</surname><given-names>A G</given-names></name></person-group><year>2003</year><source>Phys. Rev. <named-content content-type="jnl-part">E</named-content></source><volume>67</volume><elocation-id content-type="artnum">021402</elocation-id><pub-id pub-id-type="doi">10.1103/PhysRevE.67.021402</pub-id></element-citation></ref><ref id="nj427523bib10"><label>10</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Lacoste</surname><given-names>D</given-names></name><name><surname>Lubensky</surname><given-names>T C</given-names></name></person-group><year>2001</year><source>Phys. Rev. <named-content content-type="jnl-part">E</named-content></source><volume>64</volume><elocation-id content-type="artnum">041506</elocation-id><pub-id pub-id-type="doi">10.1103/PhysRevE.64.041506</pub-id></element-citation></ref><ref id="nj427523bib11"><label>11</label><element-citation publication-type="preprint"><person-group person-group-type="author"><name><surname>Rosensweig</surname><given-names>R E</given-names></name><name><surname>Zahn</surname><given-names>M</given-names></name><name><surname>Shumovich</surname><given-names>R</given-names></name></person-group><year>1983</year><source>J. Magn. Magn. Mater.</source><volume>39</volume><fpage>127</fpage><lpage>132</lpage><page-range>127–32</page-range><pub-id pub-id-type="doi">10.1016/0304-8853(83)90416-X</pub-id></element-citation></ref><ref id="nj427523bib12"><label>12</label><element-citation publication-type="book"><person-group person-group-type="author"><name><surname>Tsebers</surname><given-names>A</given-names></name><name><surname>Maiorov</surname><given-names>M</given-names></name></person-group><year>1980</year><source>Magnetohydrodynamics no 2</source><publisher-loc>New York</publisher-loc><publisher-name>Consultants Bureau</publisher-name><fpage>p 21</fpage></element-citation></ref><ref id="nj427523bib13"><label>13</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Berdiyorov</surname><given-names>G R</given-names></name><name><surname>Hernandez</surname><given-names>A D</given-names></name><name><surname>Peeters</surname><given-names>F M</given-names></name></person-group><year>2009</year><source>Phys. Rev. Lett.</source><volume>103</volume><elocation-id content-type="artnum">267002</elocation-id><pub-id pub-id-type="doi">10.1103/PhysRevLett.103.267002</pub-id></element-citation></ref><ref id="nj427523bib14"><label>14</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Horowitz</surname><given-names>C J</given-names></name><name><surname>Pérez-García</surname><given-names>M A</given-names></name><name><surname>Piekarewicz</surname><given-names>J</given-names></name></person-group><year>2004</year><source>Phys. Rev. <named-content content-type="jnl-part">C</named-content></source><volume>69</volume><elocation-id content-type="artnum">045804</elocation-id><pub-id pub-id-type="doi">10.1103/PhysRevC.69.045804</pub-id></element-citation></ref><ref id="nj427523bib15"><label>15</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Seul</surname><given-names>M</given-names></name><name><surname>Andelman</surname><given-names>D</given-names></name></person-group><year>1995</year><source>Science</source><volume>267</volume><fpage>476</fpage><lpage>483</lpage><page-range>476–83</page-range><pub-id pub-id-type="doi">10.1126/science.267.5197.476</pub-id></element-citation></ref><ref id="nj427523bib16"><label>16</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Reichhardt</surname><given-names>C</given-names></name><name><surname>Olson</surname><given-names>C J</given-names></name><name><surname>Martin</surname><given-names>I</given-names></name><name><surname>Bishop</surname><given-names>A R</given-names></name></person-group><year>2003</year><source>Europhys. Lett.</source><volume>61</volume><fpage>221</fpage><pub-id pub-id-type="doi">10.1209/epl/i2003-00222-0</pub-id></element-citation></ref><ref id="nj427523bib17"><label>17</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Reichhardt</surname><given-names>C</given-names></name><name><surname>Reichhardt</surname><given-names>C J O</given-names></name><name><surname>Martin</surname><given-names>I</given-names></name><name><surname>Bishop</surname><given-names>A R</given-names></name></person-group><year>2003</year><source>Phys. Rev. Lett.</source><volume>90</volume><elocation-id content-type="artnum">026401</elocation-id><pub-id pub-id-type="doi">10.1103/PhysRevLett.90.026401</pub-id></element-citation></ref><ref id="nj427523bib18"><label>18</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Reichhardt</surname><given-names>C</given-names></name><name><surname>Reichhardt</surname><given-names>C J O</given-names></name><name><surname>Bishop</surname><given-names>A R</given-names></name></person-group><year>2005</year><source>Europhys. Lett.</source><volume>72</volume><fpage>444</fpage><pub-id pub-id-type="doi">10.1209/epl/i2005-10250-8</pub-id></element-citation></ref><ref id="nj427523bib19"><label>19</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Sciortino</surname><given-names>F</given-names></name><name><surname>Tartaglia</surname><given-names>P</given-names></name></person-group><year>2005</year><source>Adv. Phys.</source><volume>54</volume><fpage>471</fpage><pub-id pub-id-type="doi">10.1080/00018730500414570</pub-id></element-citation></ref><ref id="nj427523bib20"><label>20</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Sciortino</surname><given-names>F</given-names></name><name><surname>Mossa</surname><given-names>S</given-names></name><name><surname>Zaccarelli</surname><given-names>E</given-names></name><name><surname>Tartaglia</surname><given-names>P</given-names></name></person-group><year>2004</year><source>Phys. Rev. Lett.</source><volume>93</volume><elocation-id content-type="artnum">055701</elocation-id><pub-id pub-id-type="doi">10.1103/PhysRevLett.93.055701</pub-id></element-citation></ref><ref id="nj427523bib21"><label>21</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Xu</surname><given-names>X B</given-names></name><name><surname>Fangohr</surname><given-names>H</given-names></name><name><surname>Ding</surname><given-names>S Y</given-names></name><name><surname>Zhou</surname><given-names>F</given-names></name><name><surname>Xu</surname><given-names>X N</given-names></name><name><surname>Wang</surname><given-names>Z H</given-names></name><name><surname>Gu</surname><given-names>M</given-names></name><name><surname>Shi</surname><given-names>D Q</given-names></name><name><surname>Dou</surname><given-names>S X</given-names></name></person-group><year>2011</year><source>Phys. Rev. <named-content content-type="jnl-part">B</named-content></source><volume>83</volume><elocation-id content-type="artnum">014501</elocation-id><pub-id pub-id-type="doi">10.1103/PhysRevB.83.014501</pub-id></element-citation></ref><ref id="nj427523bib22"><label>22</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Tarzia</surname><given-names>M</given-names></name><name><surname>Coniglio</surname><given-names>A</given-names></name></person-group><year>2006</year><source>Phys. Rev. Lett.</source><volume>96</volume><elocation-id content-type="artnum">075702</elocation-id><pub-id pub-id-type="doi">10.1103/PhysRevLett.96.075702</pub-id></element-citation></ref><ref id="nj427523bib23"><label>23</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Tarzia</surname><given-names>M</given-names></name><name><surname>Coniglio</surname><given-names>A</given-names></name></person-group><year>2007</year><source>Phys. Rev. <named-content content-type="jnl-part">E</named-content></source><volume>75</volume><elocation-id content-type="artnum">011410</elocation-id><pub-id pub-id-type="doi">10.1103/PhysRevE.75.011410</pub-id></element-citation></ref><ref id="nj427523bib24"><label>24</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Mossa</surname><given-names>S</given-names></name><name><surname>Sciortino</surname><given-names>F</given-names></name><name><surname>Tartaglia</surname><given-names>P</given-names></name><name><surname>Zaccarelli</surname><given-names>E</given-names></name></person-group><year>2004</year><source>Langmuir</source><volume>20</volume><fpage>10756</fpage><lpage>10763</lpage><page-range>10756–63</page-range><pub-id pub-id-type="doi">10.1021/la048554t</pub-id></element-citation></ref><ref id="nj427523bib25"><label>25</label><element-citation publication-type="preprint"><person-group person-group-type="author"><name><surname>de Candia</surname><given-names>A</given-names></name><name><surname>Del Gado</surname><given-names>E</given-names></name><name><surname>Fierro</surname><given-names>A</given-names></name><name><surname>Sator</surname><given-names>N</given-names></name><name><surname>Tarzia</surname><given-names>M</given-names></name><name><surname>Coniglio</surname><given-names>A</given-names></name></person-group><year>2006</year><source>Phys. Rev. <named-content content-type="jnl-part">E</named-content></source><volume>74</volume><elocation-id content-type="artnum">010403</elocation-id><pub-id pub-id-type="other">http://pre.aps.org/abstract/PRE/v74/i1/e010403</pub-id></element-citation></ref><ref id="nj427523bib26"><label>26</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Kramer</surname><given-names>L</given-names></name></person-group><year>1971</year><source>Phys. Rev. <named-content content-type="jnl-part">B</named-content></source><volume>3</volume><fpage>3821</fpage><lpage>3825</lpage><page-range>3821–5</page-range><pub-id pub-id-type="doi">10.1103/PhysRevB.3.3821</pub-id></element-citation></ref><ref id="nj427523bib27"><label>27</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Brandt</surname><given-names>E H</given-names></name><name><surname>Essmann</surname><given-names>U</given-names></name></person-group><year>1987</year><source>Phys. Status Solidi <named-content content-type="jnl-part">b</named-content></source><volume>144</volume><fpage>13</fpage><pub-id pub-id-type="doi">10.1002/(ISSN)1521-3951</pub-id></element-citation></ref><ref id="nj427523bib28"><label>28</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Brandt</surname><given-names>E</given-names></name><name><surname>Das</surname><given-names>M</given-names></name></person-group><year>2011</year><source>J. Supercond. Novel Magn.</source><volume>24</volume><fpage>57</fpage><lpage>67</lpage><page-range>57–67</page-range><pub-id pub-id-type="doi">10.1007/s10948-010-1046-8</pub-id></element-citation></ref><ref id="nj427523bib29"><label>29</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Moshchalkov</surname><given-names>V</given-names></name><name><surname>Menghini</surname><given-names>M</given-names></name><name><surname>Nishio</surname><given-names>T</given-names></name><name><surname>Chen</surname><given-names>Q H</given-names></name><name><surname>Silhanek</surname><given-names>A V</given-names></name><name><surname>Dao</surname><given-names>V H</given-names></name><name><surname>Chibotaru</surname><given-names>L F</given-names></name><name><surname>Zhigadlo</surname><given-names>N D</given-names></name><name><surname>Karpinski</surname><given-names>J</given-names></name></person-group><year>2009</year><source>Phys. Rev. Lett.</source><volume>102</volume><elocation-id content-type="artnum">117001</elocation-id><pub-id pub-id-type="doi">10.1103/PhysRevLett.102.117001</pub-id></element-citation></ref><ref id="nj427523bib30"><label>30</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Chaves</surname><given-names>A</given-names></name><name><surname>Komendová</surname><given-names>L</given-names></name><name><surname>Milošević</surname><given-names>M V</given-names></name><name><surname>Andrade</surname><given-names>J S</given-names></name><name><surname>Farias</surname><given-names>G A</given-names></name><name><surname>Peeters</surname><given-names>F M</given-names></name></person-group><year>2011</year><source>Phys. Rev. <named-content content-type="jnl-part">B</named-content></source><volume>83</volume><elocation-id content-type="artnum">214523</elocation-id><pub-id pub-id-type="doi">10.1103/PhysRevB.83.214523</pub-id></element-citation></ref><ref id="nj427523bib31"><label>31</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Komendová</surname><given-names>L</given-names></name><name><surname>Milošević</surname><given-names>M V</given-names></name><name><surname>Shanenko</surname><given-names>M V</given-names></name><name><surname>Peeters</surname><given-names>F M</given-names></name></person-group><year>2011</year><source>Phys. Rev. <named-content content-type="jnl-part">B</named-content></source><volume>84</volume><elocation-id content-type="artnum">064522</elocation-id><pub-id pub-id-type="other">http://prb.aps.org/abstract/PRB/v84/i6/e064522</pub-id></element-citation></ref><ref id="nj427523bib32"><label>32</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Babaev</surname><given-names>E</given-names></name><name><surname>Speight</surname><given-names>M</given-names></name></person-group><year>2005</year><source>Phys. Rev. <named-content content-type="jnl-part">B</named-content></source><volume>72</volume><elocation-id content-type="artnum">180502</elocation-id><pub-id pub-id-type="doi">10.1103/PhysRevB.72.180502</pub-id></element-citation></ref><ref id="nj427523bib33"><label>33</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Misko</surname><given-names>V R</given-names></name><name><surname>Bothner</surname><given-names>D</given-names></name><name><surname>Kemmler</surname><given-names>M</given-names></name><name><surname>Kleiner</surname><given-names>R</given-names></name><name><surname>Koelle</surname><given-names>D</given-names></name><name><surname>Peeters</surname><given-names>F M</given-names></name><name><surname>Nori</surname><given-names>F</given-names></name></person-group><year>2010</year><source>Phys. Rev. <named-content content-type="jnl-part">B</named-content></source><volume>82</volume><elocation-id content-type="artnum">184512</elocation-id><pub-id pub-id-type="doi">10.1103/PhysRevB.82.184512</pub-id></element-citation></ref><ref id="nj427523bib34"><label>34</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Plass</surname><given-names>R</given-names></name><name><surname>Last</surname><given-names>J A</given-names></name><name><surname>Bartelt</surname><given-names>N C</given-names></name><name><surname>Kellogg</surname><given-names>G L</given-names></name></person-group><year>2001</year><source>Nature</source><volume>412</volume><fpage>875</fpage><pub-id pub-id-type="doi">10.1038/35091143</pub-id></element-citation></ref><ref id="nj427523bib35"><label>35</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Plass</surname><given-names>R</given-names></name><name><surname>Bartelt</surname><given-names>N C</given-names></name><name><surname>Kellogg</surname><given-names>G L</given-names></name></person-group><year>2002</year><source>J. Phys.: Condens. Matter</source><volume>14</volume><fpage>4227</fpage><pub-id pub-id-type="doi">10.1088/0953-8984/14/16/313</pub-id></element-citation></ref><ref id="nj427523bib36"><label>36</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Crocker</surname><given-names>J C</given-names></name><name><surname>Grier</surname><given-names>D G</given-names></name></person-group><year>1996</year><source>J. Colloid Interface Sci.</source><volume>179</volume><fpage>298</fpage><lpage>310</lpage><page-range>298–310</page-range><pub-id pub-id-type="doi">10.1006/jcis.1996.0217</pub-id></element-citation></ref><ref id="nj427523bib37"><label>37</label><element-citation publication-type="book"><person-group person-group-type="author"><name><surname>Frenkel</surname><given-names>D</given-names></name><name><surname>Smit</surname><given-names>B</given-names></name></person-group><year>2002</year><source>Understanding Molecular Simulation: From Algorithms to Applications</source><edition>2nd edn</edition><publisher-loc>Edinburgh</publisher-loc><publisher-name>Academic</publisher-name></element-citation></ref></ref-list></back></article></istex:document></istex:metadataXml><mods version="3.6"><titleInfo><title>Analysis of pattern formation in systems with competing range interactions</title></titleInfo><titleInfo type="alternative" contentType="CDATA"><title>Analysis of pattern formation in systems with competing range interactions</title></titleInfo><name type="personal"><namePart type="given">H J</namePart><namePart type="family">Zhao</namePart><affiliation>Department of Physics, University of Antwerpen, Groenenborgerlaan 171, B-2020 Antwerpen, Belgium</affiliation><affiliation>E-mail: haijun.zhao@ua.ac.be</affiliation><role><roleTerm type="text">author</roleTerm></role></name><name type="personal"><namePart type="given">V R</namePart><namePart type="family">Misko</namePart><affiliation>Department of Physics, University of Antwerpen, Groenenborgerlaan 171, B-2020 Antwerpen, Belgium</affiliation><affiliation>E-mail: vyacheslav.misko@ua.ac.be</affiliation><role><roleTerm type="text">author</roleTerm></role></name><name type="personal"><namePart type="given">F M</namePart><namePart type="family">Peeters</namePart><affiliation>Department of Physics, University of Antwerpen, Groenenborgerlaan 171, B-2020 Antwerpen, Belgium</affiliation><affiliation>E-mail: francois.peeters@ua.ac.be</affiliation><role><roleTerm type="text">author</roleTerm></role></name><typeOfResource>text</typeOfResource><genre type="research-article" displayLabel="research-article"></genre><subject><genre>article-type</genre><topic>Paper</topic></subject><originInfo><publisher>IOP Publishing</publisher><dateIssued encoding="w3cdtf">2012-06</dateIssued><dateCreated encoding="w3cdtf">2012-06-25</dateCreated><copyrightDate encoding="w3cdtf">2012</copyrightDate></originInfo><language><languageTerm type="code" authority="iso639-2b">eng</languageTerm><languageTerm type="code" authority="rfc3066">en</languageTerm></language><physicalDescription><internetMediaType>text/html</internetMediaType></physicalDescription><abstract>We analyzed pattern formation and identified various morphologies in a system of particles interacting through a non-monotonic potential with a competing range interaction characterized by a repulsive core (rrc), using molecular-dynamics simulations. Depending on parameters, the interaction potential models the inter-particle interaction in various physical systems ranging from atoms, molecules and colloids to vortices in low type-II superconductors and in recently discovered type-1.5 superconductors. We constructed a morphology diagram in the plane critical radius rcdensity n and proposed a new approach to characterizing the different types of patterns. Namely, we elaborated a set of quantitative criteria in order to identify the different pattern types, using the radial distribution function (RDF), the local density function and the occupation factor.</abstract> <relatedItem type="host"><titleInfo><title>New Journal of Physics</title></titleInfo><genre type="journal">journal</genre><identifier type="eISSN">1367-2630</identifier><identifier type="PublisherID">nj</identifier><part><date>2012</date><detail type="volume"><caption>vol.</caption><number>14</number></detail><detail type="issue"><caption>no.</caption><number>6</number></detail><extent unit="pages"><total>21</total></extent></part></relatedItem><identifier type="istex">A36699F86F9CD28199495900E74448C936502F68</identifier><identifier type="DOI">10.1088/1367-2630/14/6/063032</identifier><identifier type="href">http://stacks.iop.org/NJP/14/063032</identifier><identifier type="ArticleID">nj427523</identifier><accessCondition type="use and reproduction" contentType="copyright">IOP Publishing and Deutsche Physikalische Gesellschaft</accessCondition><recordInfo><recordContentSource>IOP</recordContentSource></recordInfo></mods></metadata><serie></serie></istex></record>Pour manipuler ce document sous Unix (Dilib)
EXPLOR_STEP=$WICRI_ROOT/Wicri/Belgique/explor/OpenAccessBelV2/Data/Istex/Corpus
HfdSelect -h $EXPLOR_STEP/biblio.hfd -nk 001713 | SxmlIndent | more
Ou
HfdSelect -h $EXPLOR_AREA/Data/Istex/Corpus/biblio.hfd -nk 001713 | SxmlIndent | more
Pour mettre un lien sur cette page dans le réseau Wicri
{{Explor lien
|wiki= Wicri/Belgique
|area= OpenAccessBelV2
|flux= Istex
|étape= Corpus
|type= RBID
|clé= ISTEX:A36699F86F9CD28199495900E74448C936502F68
|texte= Analysis of pattern formation in systems with competing range interactions
}}
| This area was generated with Dilib version V0.6.25. |  |