Nonmonotonic Magnetic Susceptibility of Dipolar Hard-Spheres at Low Temperature and Density

We investigate, via numerical simulations, mean field, and density functional theories, the magnetic response of a dipolar hard sphere fluid at low temperatures and densities, in the region of strong association. The proposed parameter-free theory is able to capture both the density and temperature dependence of the ring-chain equilibrium and the contribution to the susceptibility of a chain of generic length. The theory predicts a nonmonotonic temperature dependence of the initial (zero field) magnetic susceptibility, arising from the competition between magnetically inert particle rings and magnetically active chains. Monte Carlo simulation results closely agree with the theoretical findings. DOI: 10.1103/PhysRevLett.110.148306

Effective isotropic potential for dipolar hard spheres

A new effective isotropic potential is proposed for the dipolar hard-sphere fluid, on the basis of recent results by others for its angle-averaged radial distribution function. The new effective potential is shown to exhibit oscillations even for moderately high densities and moderately strong dipole moments, which are absent from earlier effective isotropic potentials. The validity and significance of this result are briefly discussed.

Low temperature structural transitions in dipolar hard spheres: the influence on magnetic properties

We investigate the structural chain-to-ring transition at low temperature in a gas of dipolar hard spheres (DRS). Due to the weakening of entropic contribution, ring formation becomes noticeable when the effective dipole-dipole magnetic interaction increases, It results in the redistribution of particles from usually observed flexible chains into flexible rings. The concentration (rho) of DI-IS plays a crucial part in this transition: at a very low rho only chains and rings are observed, whereas even a slight increase of the volume fraction leads to the formation of branched or defect structures. As a result, the fraction of DHS aggregated in defect-free rings turns out to be a non-monotonic function of rho. The average ring size is found to be a slower increasing function of rho when compared Lo that of chains. Both theory and computer simulations confirm the dramatic influence of the ring formation on the rho-dependence of the initial magnetic susceptibility (chi) when the temperature decreases. The rings clue to their zero total dipole moment are irresponsive to a weak magnetic field and drive to the strong decrease of the initial magnetic susceptibility. (C) 2014 Elsevier B.V. All rights reserved.

Branching points in the low-temperature dipolar hard sphere fluid

In this contribution, we investigate the low-temperature, low-density behaviour of dipolar hard-sphere (DHS) particles, i.e., hard spheres with dipoles embedded in their centre. We aim at describing the DHS fluid in terms of a network of chains and rings (the fundamental clusters) held together by branching points (defects) of different nature. We first introduce a systematic way of classifying inter-cluster connections according to their topology, and then employ this classification to analyse the geometric and thermodynamic properties of each class of defects, as extracted from state-of-the-art equilibrium Monte Carlo simulations. By computing the average density and energetic cost of each defect class, we find that the relevant contribution to inter-cluster interactions is indeed provided by (rare) three-way junctions and by four-way junctions arising from parallel or anti-parallel locally linear aggregates. All other (numerous) defects are either intra-cluster or associated to low cluster-cluster interaction energies, suggesting that these defects do not play a significant part in the thermodynamic description of the self-assembly processes of dipolar hard spheres. (C) 2013 AIP Publishing LLC.

Structure and anomalous solubility for hard spheres in an associating lattice gas model

In this paper we investigate the solubility of a hard-sphere gas in a solvent modeled as an associating lattice gas. The solution phase diagram for solute at 5% is compared with the phase diagram of the original solute free model. Model properties are investigated both through Monte Carlo simulations and a cluster approximation. The model solubility is computed via simulations and is shown to exhibit a minimum as a function of temperature. The line of minimum solubility (TmS) coincides with the line of maximum density (TMD) for different solvent chemical potentials, in accordance with the literature on continuous realistic models and on the "cavity" picture. (C) 2012 American Institute of Physics. [http://dx.doi.org/10.1063/1.4743635]

Equilibrium fluid-solid coexistence of hard spheres

We present a tethered Monte Carlo simulation of the crystallization of hard spheres. Our method boosts the traditional umbrella sampling to the point of making practical the study of constrained Gibbs’ free energies depending on several crystalline order parameters. We obtain high-accuracy estimates of the fluid-crystal coexistence pressure for up to 2916 particles (enough to accommodate fluid-solid interfaces). We are able to extrapolate to infinite volume the coexistence pressure [p_(co) = 11.5727(10)k_(B)T/σ^(3)] and the interfacial free energy [γ_({100}) = 0.636(11)k_(B)T/σ^(2)].

Re-entrant phase behaviour of network fluids: A patchy particle model with temperature-dependent valence

We study a model consisting of particles with dissimilar bonding sites ("patches"), which exhibits self-assembly into chains connected by Y-junctions, and investigate its phase behaviour by both simulations and theory. We show that, as the energy cost epsilon(j) of forming Y-junctions increases, the extent of the liquid-vapour coexistence region at lower temperatures and densities is reduced. The phase diagram thus acquires a characteristic "pinched" shape in which the liquid branch density decreases as the temperature is lowered. To our knowledge, this is the first model in which the predicted topological phase transition between a fluid composed of short chains and a fluid rich in Y-junctions is actually observed. Above a certain threshold for epsilon(j), condensation ceases to exist because the entropy gain of forming Y-junctions can no longer offset their energy cost. We also show that the properties of these phase diagrams can be understood in terms of a temperature-dependent effective valence of the patchy particles. (C) 2011 American Institute of Physics. [doi: 10.1063/1.3605703]

Reentrant Phase Diagram of Network Fluids

We introduce a microscopic model for particles with dissimilar patches which displays an unconventional "pinched'' phase diagram, similar to the one predicted by Tlusty and Safran in the context of dipolar fluids [Science 290, 1328 (2000)]. The model-based on two types of patch interactions, which account, respectively, for chaining and branching of the self-assembled networks-is studied both numerically via Monte Carlo simulations and theoretically via first-order perturbation theory. The dense phase is rich in junctions, while the less-dense phase is rich in chain ends. The model provides a reference system for a deep understanding of the competition between condensation and self-assembly into equilibrium-polymer chains.

Three-dimensional patchy lattice model for empty fluids

The phase diagram of a simple model with two patches of type A and ten patches of type B (2A10B) on the face centred cubic lattice has been calculated by simulations and theory. Assuming that there is no interaction between the B patches the behavior of the system can be described in terms of the ratio of the AB and AA interactions, r. Our results show that, similarly to what happens for related off-lattice and two-dimensional lattice models, the liquid-vapor phase equilibria exhibit reentrant behavior for some values of the interaction parameters. However, for the model studied here the liquid-vapor phase equilibria occur for values of r lower than 1/3, a threshold value which was previously thought to be universal for 2AnB models. In addition, the theory predicts that below r = 1/3 (and above a new condensation threshold which is < 1/3) the reentrant liquid-vapor equilibria are so extreme that it exhibits a closed loop with a lower critical point, a very unusual behavior in single-component systems. An order-disorder transition is also observed at higher densities than the liquid-vapor equilibria, which shows that the liquid-vapor reentrancy occurs in an equilibrium region of the phase diagram. These findings may have implications in the understanding of the condensation of dipolar hard spheres given the analogy between that system and the 2AnB models considered here. (C) 2012 American Institute of Physics. [http://dx.doi.org/10.1063/1.4771591]

Self-assembly in chains, rings, and branches: a single component system with two critical points

Agências financiadoras: FCT - PEstOE/FIS/UI0618/2011; PTDC/FIS/098254/2008 ERC-PATCHYCOLLOIDS e MIUR-PRIN

Three-Dimensional patchy lattice model for empty fluids

The phase diagram of a simple model with two patches of type A and ten patches of type B (2A10B) on the face centred cubic lattice has been calculated by simulations and theory. Assuming that there is no interaction between the B patches the behavior of the system can be described in terms of the ratio of the AB and AA interactions, r. Our results show that, similarly to what happens for related off-lattice and two-dimensional lattice models, the liquid-vapor phase equilibria exhibit reentrant behavior for some values of the interaction parameters. However, for the model studied here the liquid-vapor phase equilibria occur for values of r lower than 1/3, a threshold value which was previously thought to be universal for 2AnB models. In addition, the theory predicts that below r = 1/3 (and above a new condensation threshold which is < 1/3) the reentrant liquid-vapor equilibria are so extreme that it exhibits a closed loop with a lower critical point, a very unusual behavior in single-component systems. An order-disorder transition is also observed at higher densities than the liquid-vapor equilibria, which shows that the liquid-vapor reentrancy occurs in an equilibrium region of the phase diagram. These findings may have implications in the understanding of the condensation of dipolar hard spheres given the analogy between that system and the 2AnB models considered here.

Temperature-induced structural transitions in self-assembling magnetic nanocolloids

With the help of a unique combination of density functional theory and computer simulations, we discover two possible scenarios, depending on concentration, for the hierarchical self-assembly of magnetic nanoparticles on cooling. We show that typically considered low temperature clusters, i.e. defect-free chains and rings, merge into more complex branched structures through only three types of defects: four-way X junctions, three-way Y junctions and two-way Z junctions. Our accurate calculations reveal the predominance of weakly magnetically responsive rings cross-linked by X defects at the lowest temperatures. We thus provide a strategy to fine-tune magnetic and thermodynamic responses of magnetic nanocolloids to be used in medical and microfluidics applications.

Dynamik unterkühlter molekularer Flüssigkeiten in Modenkopplungsapproximation - die Rolle der Orientierungsfreiheitsgrade für den Bosonenpeak

Die von Schilling et al. 1997 aufgestellten molekularenModenkopplungsgleichungen sind für zwei Systeme, dipolareharte Kugeln undharte Rotationsellipsoide, numerisch gelöst worden.In den Kompressibilitätsspektrenund den entsprechenden Korrelationsspektren der beidenSysteme fiel im Vergleich mit den Spektren einfacherFlüssigkeiten einzusätzlicher Peak (bp) auf.Die Position des zusätzlichen Peaks ist etwa eine Dekade unterhalb der Position derHochfrequenzpeaks, also im mesoskopischen Frequenzbereich.In molekularen Flüssigkeiten sind Phasenübergänge derOrientierung derglqq Molekülegrqq möglich.Bei dipolaren harten Kugeln kann ferroelektrische Ordnung,bei hartenEllipsoiden nematische Ordnung auftreten. Als Vorläuferphänomen der Orientierungsordnung bilden sich in den molekularenFlüssigkeiten orientierungsgeordnete Domänen. Die Kurzzeitdynamik dieser Domänen findetauf mesoskopischer Zeitskala statt. In der unterkühlten Flüssigkeit ist die Dynamik aller Modenan die Dynamikdieser Orientierungsmoden gekoppelt. Unsere Untersuchungen zeigen starke Hinweise darauf, dassdiese Kopplung denbp erzeugt.Demnach ist die Streuung nichtlokalisierter longitudinalerakustischer Modenan lokalisierten Orientierungsmoden die Ursache für den bp.hat nahezu keine q-Abhängigkeit und zeigt einenIsotopeneffekt. Weiterhin verglichen wir die Eigenschaften des bp mit denEigenschaften des Bosonenpeaks (BP).Der BP befindet sich in Spektren realer Systeme ebenfallseine Dekadeunterhalb der Hochfrequenzpeaks.Dieser Vergleich und der Vergleich mit den Ergebnissen einerMD-Simulation desHE-Systems ließ uns vermuten, dass es sich bei dem bp und bei dem BP in Spektrenmolekularer Systeme um dasselbe Phänomen handelt.In einem schematischen Modenkopplungsmodell konnte dasVerhalten des bpzusätzlich untersucht werden.

Patching up dipoles: can dipolar particles be viewed as patchy colloids?

We investigate whether the liquid-vapour phase transition of strongly dipolar fluids can be understood using a model of patchy colloids. These consist of hard spherical particles with three short-ranged attractive sites (patches) on their surfaces. Two of the patches are of type A and one is of type B. Patches A on a particle may bond either to a patch A or to a patch B on another particle. Formation of an AA (AB) bond lowers the energy by epsilon AA (epsilon AB). In the limit [image omitted], this patchy model exhibits condensation driven by AB-bonds (Y-junctions). Y-junctions are also present in low-density, strongly dipolar fluids, and have been conjectured to play a key role in determining their critical behaviour. We map the dipolar Yukawa hard-sphere (DYHS) fluid onto this 2A + 1B patchy model by requiring that the latter reproduce the correct DYHS critical point as a function of the isotropic interaction strength epsilon Y. This is achieved for sensible values of epsilon AB and the bond volumes. Results for the internal energy and the particle coordination number are in qualitative agreement with simulations of DYHSs. Finally, by taking the limit [image omitted], we arrive at a new estimate for the critical point of the dipolar hard-sphere fluid, which agrees with extrapolations from simulation.