Search

Fundamental Measure Theory (FMT) is a successful and versatile approach for describing the properties of the hard-sphere fluid and hard-sphere mixtures within the framework of classical density functional theory (DFT). Lutsko [Phys. Rev. E 102, 062137 (2020)] introduced a version of FMT containing two free parameters, to be fixed by additional physical constraints. Whereas Lutsko focused on the stability of crystalline phases, we introduce and employ two statistical mechanical sum rules pertinent for the fluid phase, that are not automatically satisfied by FMT. By minimizing the relative deviation between different routes to calculate the excess chemical potential and the isothermal compressibility we determine the two free parameters of the theory. Our results indicate that requiring consistency with these sum rules can improve the quality of predictions of FMT for properties of the hard-sphere fluid phase. We suggest that employing these (test particle) sum rules, which apply for any interparticle pair-potential, might provide a means of testing the performance and accuracy of general DFT approximations.

We use supervised machine learning together with the concepts of classical density functional theory to investigate the effects of interparticle attraction on the pair structure, thermodynamics, bulk liquid-gas coexistence, and associated interfacial phenomena in many-body systems. Local learning of the one-body direct correlation functional is based on Monte Carlo simulations of inhomogeneous systems with randomized thermodynamic conditions, randomized planar shapes of the external potential, and randomized box sizes. Focusing on the prototypical Lennard-Jones system, we test predictions of the resulting neural attractive density functional across a broad spectrum of physical behaviour associated with liquid-gas phase coexistence in bulk and at interfaces. We analyse the bulk radial distribution function $g(r)$ obtained from automatic differentiation and the Ornstein-Zernike route and determine i) the Fisher-Widom line, i.e.\ the crossover of the asymptotic (large distance) decay of $g(r)$ from monotonic to oscillatory, ii) the (Widom) line of maximal correlation length, iii) the line of maximal isothermal compressibility and iv) the spinodal by calculating the poles of the structure factor in the complex plane. The bulk binodal and the density profile of the free liquid-gas interface are obtained from density functional minimization and the corresponding surface tension from functional line integration. We also show that the neural functional describes accurately the phenomena of drying at a hard wall and of capillary evaporation for a liquid confined in a slit pore. Our neural framework yields results that improve significantly upon standard mean-field treatments of interparticle attraction. Comparison with independent simulation results demonstrates a consistent picture of phase separation even when restricting the training to supercritical states only., Comment: 21 pages, 9 figures

Measurements of local density fluctuations are crucial to characterizing the interfacial properties of equilibrium fluids. A specific case that has been well-explored involves the heightened compressibility of water near hydrophobic entities. Commonly, a spatial profile of local fluctuation strength is constructed from measurements of the mean and variance of solvent particle number fluctuations in a set of contiguous sub-volumes of the system adjacent to the solvo/hydrophobe. An alternative measure proposed by Evans and Stewart (J. Phys.: Condens. Matter 27, 194111 (2015)) defines a local compressibility profile in terms of the chemical potential derivative of the spatial number density profile. Using Grand Canonical Monte Carlo simulation, we compare and contrast the efficacy of these two approaches for a Lennard-Jones solvent at spherical and planar solvophobic interfaces, and SPC/E water at a hydrophobic spherical solute. Our principal findings are that: (i) the local compressibility profile $\chi({\bf r})$ of Evans and Stewart is considerably more sensitive to variations in the strength of local density fluctuations than the spatial fluctuation profile $F({\bf r})$ and can resolve much more detailed structure; (ii) while the local compressibility profile is essentially independent of the choice of spatial discretization used to construct the profile, the spatial fluctuation profile exhibits strong systematic dependence on the size of the subvolumes on which the profile is defined. We clarify the origin and nature of this finite-size effect., Comment: 12 pages, 6 figures

Let $A/\mathbb{Q}$ be a Jacobian variety and let $F$ be a totally real, tamely ramified, abelian number field. Given a character $\psi$ of $F/\mathbb{Q}$, Deligne's Period Conjecture asserts the algebraicity of the suitably normalised value $\mathcal{L}(A,\psi,1)$ at $z=1$ of the Hasse-Weil-Artin $L$-function of the $\psi$-twist of $A$. We formulate a conjecture regarding the integrality properties of the family of normalised $L$-values $(\mathcal{L}(A,\psi,1))_{\psi}$, and its relation to the Tate-Shafarevich group of $A$ over $F$. We numerically investigate our conjecture through $p$-adic congruence relations between these values., Comment: 24 pages

Simulations of water near extended hydrophobic spherical solutes have revealed the presence of a region of depleted density and accompanying enhanced density fluctuations.The physical origin of both phenomena has remained somewhat obscure. We investigate these effects employing a mesoscopic binding potential analysis, classical density functional theory (DFT) calculations for a simple Lennard-Jones (LJ) solvent and Grand Canonical Monte Carlo (GCMC) simulations of a monatomic water (mw) model. We argue that the density depletion and enhanced fluctuations are near-critical phenomena. Specifically, we show that they can be viewed as remnants of the critical drying surface phase transition that occurs at bulk liquid-vapor coexistence in the macroscopic planar limit, i.e.~as the solute radius $R_s\to\infty$. Focusing on the radial density profile $\rho(r)$ and a sensitive spatial measure of fluctuations, the local compressibility profile $\chi(r)$, our binding potential analysis provides explicit predictions for the manner in which the key features of $\rho(r)$ and $\chi(r)$ scale with $R_s$, the strength of solute-water attraction $\varepsilon_{sf}$, and the deviation from liquid-vapor coexistence of the chemical potential, $\delta\mu$. These scaling predictions are confirmed by our DFT calculations and GCMC simulations. As such our theory provides a firm basis for understanding the physics of hydrophobic solvation., Comment: 18 pages

We investigate density fluctuations in a liquid close to a solvophobic substrate at which a surface critical drying transition occurs. Using classical density functional theory, we determine three measures of the spatial extent and strength of the fluctuations, i.e., the local compressibility $\chi_{\mu}(z)$, the local thermal susceptibility $\chi_T(z)$ and the reduced density $\chi_*(z)$; $z$ is the distance from the substrate. Whilst the first measure is frequently used, the second and third were introduced very recently by Eckert et.al., Phys. Rev. Lett. 125, 268004 (2020). For state points in the critical drying regime, all three measures, each scaled by its bulk value, exhibit very similar forms and the ratio of $\chi_T(z)$ to $\chi_{\mu}(z)$, for fixed $z$ in the vapour-liquid transition region, is constant. Using a scaling treatment of surface thermodynamics we show that such behaviour is to be expected on general grounds., Comment: 7 pages, 3 figures. Accepted by Physical Review E

We investigate the origin of the density depletion and enhanced density fluctuations that occur in water in the vicinity of an extended hydrophobic solute. We argue that both phenomena are remnants of the critical drying surface phase transition that occurs at liquid-vapor coexistence in the macroscopic planar limit, ie. as the solute radius $R_s\to\infty$. Focusing on the density profile $\rho(r)$ and a sensitive spatial measure of fluctuations, the local compressibility profile $\chi(r)$, we develop a scaling theory which expresses the extent of the density depletion and enhancement in compressibility in terms of $R_s$, the strength of solute-water attraction $\varepsilon_s$, and the deviation from liquid-vapor coexistence $\delta\mu$. Testing the predictions against results of classical density functional theory for a simple solvent and Grand Canonical Monte Carlo simulations of a popular water model, we find that the theory provides a firm physical basis for understanding how water behaves at a hydrophobe., Comment: 5 pages, 4 figures. To appear in Phys. Rev. Lett

Density functional calculations can fail for want of an accurate exchange-correlation approximation. The energy can instead be extracted from a sequence of density functional calculations of conditional probabilities (CP-DFT). Simple CP approximations yield usefully accurate results for two-electron ions, the hydrogen dimer, and the uniform gas at all temperatures. CP-DFT has no self-interaction error for one electron, and correctly dissociates H2, both major challenges. For warm dense matter, classical CP-DFT calculations can overcome the convergence problems of Kohn-Sham DFT., Comment: 6 pages, 3 figures, 2 tables

Clarifying the factors that control the contact angle of a liquid on a solid substrate is a long-standing scientific problem pertinent across physics, chemistry and materials science. Progress has been hampered by the lack of a comprehensive and unified understanding of the physics of wetting and drying phase transitions. Using various theoretical and simulational techniques applied to realistic fluid models, we elucidate how the character of these transitions depends sensitively on both the range of fluid-fluid and substrate-fluid interactions and the temperature. Our calculations uncover previously unrecognised classes of surface phase diagram which differ from that established for simple lattice models and often assumed to be universal. The differences relate both to the topology of the phase diagram and to the nature of the transitions, with a remarkable feature being a difference between drying and wetting transitions which persists even in the approach to the bulk critical point. Most experimental and simulational studies of liquids at a substrate belong to one of these previously unrecognised classes. We predict that while there appears to be nothing particularly special about water with regard to its wetting and drying behavior, superhydrophobic behavior should be more readily observable in experiments conducted at high temperatures than at room temperature., Comment: 10 pages plus SI

We re-visit the competition between attractive and repulsive interparticle forces in simple fluids and how this governs and connects the macroscopic phase behavior and structural properties as manifest in pair correlation functions. We focus on the asymptotic decay of the total correlation function $h(r)$ which is, in turn, controlled by the form of the pair direct correlation function $c(r)$. The decay of $r h(r)$ to zero can be either exponential (monotonic) if attraction dominates repulsion and exponentially damped oscillatory otherwise. The Fisher-Widom (FW) line separates the phase diagram into two regions characterized by the two different types of asymptotic decay. We show that there is a new and physically intuitive thermodynamic criterion which approximates well the actual FW line. This new criterion defines a line where the isothermal compressibility takes its ideal gas value $\chi_T=\chi_T^\text{id}$. We test our hypothesis by considering four commonly used models for simple fluids. In all cases the new criterion yields a line in the phase diagram that is close to the actual FW line for the thermodynamic state points that are most relevant. We also investigate (Widom) lines of maximal correlation length, emphasizing the importance of distinguishing between the true and Ornstein-Zernike correlation lengths

This is an investigation into the possible existence and consequences of a Birch-Swinnerton-Dyer-type formula for L-functions of elliptic curves twisted by Artin representations. We translate expected properties of L-functions into purely arithmetic predictions for elliptic curves, and show that these force some peculiar properties of the Tate-Shafarevich group, which do not appear to be tractable by traditional Selmer group techniques. In particular we exhibit settings where the different p-primary components of the Tate-Shafarevich group do not behave independently of one another. We also give examples of "arithmetically identical" settings for elliptic curves twisted by Artin representations, where the associated L-values can nonetheless differ, in contrast to the classical Birch-Swinnerton-Dyer conjecture., Comment: 27 pages, new version

Understanding the phase behaviour of mixtures continues to pose challenges, even for systems that might be considered "simple". Here we consider a very simple mixture of two colloidal and one non-adsorbing polymer species which can be simplified even further to a size-asymmetrical binary mixture, in which the effective colloid-colloid interactions depend on the polymer concentration. We show that this basic system exhibits surprisingly rich phase behaviour. In particular, we enquire whether such a system features only a liquid-vapor phase separation (as in one-component colloid-polymer mixtures) or whether, additionally, liquid-liquid demixing of two colloidal phases can occur. Particle-resolved experiments show demixing-like behaviour, but when combined with bespoke Monte Carlo simulations, this proves illusory, and we reveal that only a single liquid-vapor transition occurs. Progressive migration of the small particles to the liquid phase as the polymer concentration increases gives rise to composition inversion - a maximum in the large particle concentration in the liquid phase. Near criticality the density fluctuations are found to be dominated by the larger colloids., Comment: submitted to J. Chem. Phys

We describe and then elaborate the model of trading zones first presented in Collins et al 2007, Trading Zones and Interactional Expertise. We believe this expanded version of the model includes some very important but previously overlooked ways for separate language communities to communicate. We illustrate the argument with examples., Comment: 9 pages, 1 Figure

In classical density functional theory (DFT) the part of the Helmholtz free energy functional arising from attractive inter-particle interactions is often treated in a mean-field or van der Waals approximation. On the face of it, this is a somewhat crude treatment as the resulting functional generates the simple random phase approximation (RPA) for the bulk fluid pair direct correlation function. We explain why using standard mean-field DFT to describe inhomogeneous fluid structure and thermodynamics is more accurate than one might expect based on this observation. By considering the pair correlation function $g(x)$ and structure factor $S(k)$ of a one-dimensional model fluid, for which exact results are available, we show that the mean-field DFT, employed within the test-particle procedure, yields results much superior to those from the RPA closure of the bulk Ornstein-Zernike equation. We argue that one should not judge the quality of a DFT based solely on the approximation it generates for the bulk pair direct correlation function., Comment: 9 pages, 3 figures

We report a theoretical and simulation study of the drying and wetting phase transitions of a truncated Lennard-Jones fluid at a flat structureless wall. Binding potential calculations predict that the nature of these transitions depends on whether the wall-fluid attraction has a long ranged (LR) power law decay, or is instead truncated, rendering it short ranged (SR). Using grand canonical Monte Carlo simulation and classical density functional theory we examine both cases in detail. We find that for the LR case wetting is first order, while drying is continuous (critical) and occurs exactly at zero attractive wall strength, ie. in the limit of a hard wall. In the SR case, drying is also critical but the order of the wetting transition depends on the truncation range of the wall-fluid potential. We characterize the approach to critical drying and wetting in terms of the density and local compressibility profiles and via the finite-size scaling properties of the probability distribution of the overall density. For the LR case, where the drying point is known exactly, this analysis allows us to estimate the exponent $\nu_\parallel$ which controls the parallel correlation length, i.e. the extent of vapor bubbles at the wall. Surprisingly, the value we obtain is over twice that predicted by mean field and renormalization group calculations, despite the fact that our three dimensional system is at the upper critical dimension where mean field theory for critical exponents is expected to hold. We suggest reasons for this discrepancy., Comment: 22 pages, 19 figures

Using classical density functional theory (DFT) we calculate the density profile $\rho({\mathbf r})$ and local compressibility $\chi({\mathbf r})$ of a simple liquid solvent in which a pair of blocks with (microscopic) rectangular cross-section are immersed. We consider blocks that are solvophobic, solvophilic and also ones that have both solvophobic and solvophilic patches. Large values of $\chi({\mathbf r})$ correspond to regions in space where the liquid density is fluctuating most strongly. We seek to elucidate how enhanced density fluctuations correlate with the solvent mediated force between the blocks, as the distance between the blocks and the chemical potential of the liquid reservoir vary. For sufficiently solvophobic blocks, at small block separations and small deviations from bulk gas-liquid coexistence, we observe a strongly attractive (near constant) force, stemming from capillary evaporation to form a low density gas-like intrusion between the blocks. The accompanying $\chi({\mathbf r})$ exhibits structure which reflects the incipient gas-liquid interfaces that develop. We argue that our model system provides a means to understanding the basic physics of solvent mediated interactions between nanostructures, and between objects such as proteins in water, that possess hydrophobic and hydrophilic patches., Comment: 19 pages, 21 figures

The notion of interactional expertise is explained starting with its origins and discussing its many applications. Interactional expertise is the ability to understand a technical area purely be deeply immersed in its 'practice-language' without actually practising. One of its many applications is to explain how large sciences are managed., Comment: 39 pages, one table

We report a detailed simulation and classical density functional theory study of the drying transition in a realistic model fluid at a smooth substrate. This transition (in which the contact angle $\theta\to 180^\circ$) is shown to be critical for both short ranged and long-ranged substrate-fluid interaction potentials. In the latter case critical drying occurs at exactly zero attractive substrate strength. This observation permits the accurate elucidation of the character of the transition via a finite-size scaling analysis of the density probability function. We find that the critical exponent $\nu_\parallel$ that controls the parallel correlation length, i.e. the extent of vapor bubbles at the wall, is over twice as large as predicted by mean field and renormalization group calculations. We suggest a reason for the discrepancy. Our findings shed new light on fluctuation phenomena in fluids near hydrophobic and solvophobic interfaces., Comment: 5 pages plus supplementary information. To appear in Physical Review Letters

For binary fluid mixtures of spherical particles in which the two species are sufficiently different in size, the dominant wavelength of oscillations of the pair correlation functions is predicted to change from roughly the diameter of the large species to that of the small species along a sharp crossover line in the phase diagram [C. Grodon, M. Dijkstra, R. Evans & R. Roth, J.Chem.Phys. 121, 7869 (2004)]. Using particle-resolved colloid experiments in 3d we demonstrate that crossover exists and that its location in the phase diagram is in quantitative agreement with the results of both theory and our Monte-Carlo simulations. In contrast with previous work [J. Baumgartl, R. Dullens, M. Dijkstra, R. Roth & C. Bechinger, Phys.Rev.Lett. 98, 198303 (2007)], where a correspondence was drawn between crossover and percolation of both species, in our 3d study we find that structural crossover is unrelated to percolation.

Employing smart Monte Carlo sampling techniques within the grand canonical ensemble, we investigate the properties of water at a model hydrophobic substrate. By reducing the strength of substrate-water attraction we find that fluctuations in the local number density, quantified by a rigorous definition of the local compressibility $\chi(z)$, increase rapidly for distances $z$ within $1$ or $2$ molecular diameters from the substrate as the degree of hydrophobicity, measured by the macroscopic contact angle $\theta$, increases. Our simulations provide evidence for a continuous (critical) drying transition as the substrate-water interaction becomes very weak: $\cos(\theta)\to -1$. We speculate that the existence of such a transition might account for earlier simulation observations of strongly enhanced density fluctuations.