Mean-field Density Functional Theory Research Articles

Nonequilibrium relaxation of soft responsive colloids.

Stimuli-responsive macromolecules display large conformational changes during their dynamics, sometimes switching between states. Such a multi-stability is useful for the development of soft functional materials. Here, we introduce a mean-field dynamical density functional theory for a model of responsive colloids to study the nonequilibrium dynamics of a colloidal dispersion in time-dependent external fields, with a focus on the coupling of translational and conformational dynamics during their relaxation. Specifically, we consider soft Gaussian particles with a bimodal size distribution between two confining walls with time-dependent (switching-on and off) external gravitational and osmotic fields. We find a rich relaxation behavior of the systems in excellent agreement with particle-based Brownian dynamics computer simulations. In particular, we find time-asymmetric relaxations of integrated observables (wall pressures, mean size, and liquid center-of-mass) for activation/deactivation of external potentials, respectively, which are tunable by the ratio of translational and conformational diffusion time scales. Our work thus paves the way for studying the nonequilibrium relaxation dynamics of complex soft matter with multiple degrees of freedom and hierarchical relaxations.

External field-driven property localization in liquids of responsive macromolecules.

We explore theoretically the effects of external potentials on the spatial distribution of particle properties in a liquid of explicitly responsive macromolecules. In particular, we focus on the bistable particle size as a coarse-grained internal degree of freedom (DoF, or "property"), σ, that moves in a bimodal energy landscape, in order to model the response of a state-switching (big-to-small) macromolecular liquid to external stimuli. We employ a mean-field density functional theory (DFT) that provides the full inhomogeneous equilibrium distributions of a one-component model system of responsive colloids (RCs) interacting with a Gaussian pair potential. For systems confined between two parallel hard walls, we observe and rationalize a significant localization of the big particle state close to the walls, with pressures described by an exact RC wall theorem. Application of more complex external potentials, such as linear (gravitational), osmotic, and Hamaker potentials, promotes even stronger particle size segregation, in which macromolecules of different size are localized in different spatial regions. Importantly, we demonstrate how the degree of responsiveness of the particle size and its coupling to the external potential tune the position-dependent size distribution. The DFT predictions are corroborated by Brownian dynamics simulations. Our study highlights the fact that particle responsiveness can be used to localize liquid properties and therefore helps to control the property- and position-dependent function of macromolecules, e.g., in biomedical applications.

Collective Quantum Magnetism in Nitrogen-Doped Nanographenes.

The emergence of quantum magnetism in nanographenes provides ample opportunities to fabricate purely organic devices for spintronics and quantum information. Although heteroatom doping is a viable way to engineer the electronic properties of nanographenes, the synthesis of doped nanographenes with collective quantum magnetism remains elusive. Here, a set of nitrogen-doped nanographenes (N-NGs) with atomic precision are fabricated on Au(111) through a combination of imidazole [2+2+2]-cyclotrimerization and cyclodehydrogenation reactions. High-resolution scanning probe microscopy measurements reveal the presence of collective quantum magnetism for nanographenes with three radicals, with spectroscopic features which cannot be captured by mean-field density functional theory calculations but can be well reproduced by Heisenberg spin model calculations. In addition, the mechanism of magnetic exchange interaction of N-NGs has been revealed and compared with their counterparts with pure hydrocarbons. Our findings demonstrate the bottom-up synthesis of atomically precise N-NGs which can be utilized to fabricate low-dimensional extended graphene nanostructures for realizing ordered quantum phases.

Liquid structure of bistable responsive macromolecules using mean-field density-functional theory.

Macromolecular crowding typically applies to biomolecular and polymer-based systems in which the individual particles often feature a two-state folded/unfolded or coil-to-globule transition, such as found for proteins and peptides, DNA and RNA, or supramolecular polymers. Here, we employ a mean-field density functional theory (DFT) of a model of soft and bistable responsive colloids (RCs) in which the size of the macromolecule is explicitly resolved as a degree of freedom living in a bimodal 'Landau' energy landscape (exhibiting big and small states), thus directly responding to the crowding environment. Using this RC-DFT we study the effects of self-crowding on the liquid bulk structure and thermodynamics for different energy barriers and softnesses of the bimodal energy landscape, in conditions close to the coil-to-globule transition. We find substantial crowding effects on the internal distributions, a complex polydispersity behavior, and quasi-universal compression curves for increasing (generalized) packing fractions. Moreover, we uncover distinct signatures of bimodal versus unimodal behavior in the particle compression. Finally, the analysis of the pair structure - derived from the test particle route - reveals that the microstructure of the liquid is quite inhomogeneous due to local depletion effects, tuneable by particle softness.

Density functional thermodynamic description of spin, phonon and displacement degrees of freedom in antiferromagnetic-to-paramagnetic phase transition in YNiO3

This work demonstrates a direct density functional description of the finite-temperature thermodynamic properties of solids exhibiting phase transitions through positional and spin symmetry breaking degrees of freedom. A classic example addressed here is the rare-earth (R) nickelates RNiO3 where the ground state is characterized by crystallographic and magnetic (e.g., antiferromagnetic) long-range order (LRO), whereas the higher temperature paramagnetic phase manifests a range of local spin and positional symmetry breaking motifs with short-range order (SRO). Unlike time-dependent simulations of spin and positional degrees of freedom, in the present work, phases are described via a superposition of static configurations constructed by populating a periodic base lattice supercell allowing for the formation of energy lowing distribution of positional and spin local motifs. The thermal populations of the configurations in such a superposition phase are obtained from the energy-minimized Density Functional Theory (DFT)-calculated partition functions at different temperatures. This approach offers flexible inclusion of different physical contributions to the free energy, such as elastic, electronic and phonon free energies, all obtained from the same underlying DFT total energy calculations of periodic structures. The thermodynamic and magnetic properties of both LRO and SRO crystallographic and spin phases, including antiferromagnetic (AFM) to paramagnetic (PM) Néel phase transition in YNiO3 are studied. Including spin and phonon contributions, we find a DFT-calculated Néel temperature to be 144 K in satisfactory agreement with the experimental value of 145 K; whereas omitting the phonon contribution, one obtains a Néel temperature of 81 K. We present phonon contributions to the DFT-calculated temperature-dependent SRO, heat capacities, and the polymorphous distribution of nonzero local magnetic moments in the PM phase. This approach thus extends to finite temperatures the symmetry-broken DFT description of both the AFM and PM phases, demonstrating that a thermodynamic superposition approach based on symmetry broken configurations evaluated by a mean-field like DFT is sufficient to obtain a consistent description of the thermal physics of the AFM, PM phases and their interconversion in 3d oxides illustrated by YNiO3.

Analytical treatment of the structure for systems interacting via core-softened potentials

We consider a double ramp (DR) model of the family of core-softened potentials which possess a liquid–liquid critical point when the energy parameter h is negative, as in the Jagla potential [E. A. Jagla, Phys. Rev. E 63, 061501–5 (2001)]. Our model offers the possibility of gradually passing from an attractive potential to a purely repulsive potential by varying the energy parameter h. Using a mean field density functional theory and the MC simulation, we determine the structure factor of the DR model for different values of the parameters. We show that the DR model belongs to the class of adaptable potentials continuously evolving between waterlike and hard sphere (HS) systems. The HS structure factor necessary for the analytical treatment is defined in terms of the effective HS diameter σe which is a function of temperature and density, and which makes it possible to predict the structure as a function of temperature with fairly good precision.

Nanostructures of Ionic Liquids Confined in Pores of SBA-15: Insights from Experimental Studies and Mean-Field Density Functional Theory

Nanostructures of Ionic Liquids Confined in Pores of SBA-15: Insights from Experimental Studies and Mean-Field Density Functional Theory

Prediction of chemical ordering in refractory high-entropy superalloys

Refractory high-entropy superalloys (RHESs) exhibit impressive nanostructure-property relationships and have promise in next-generation high-temperature structural applications, which has motivated extensive research into these materials. The design space, however, is compositionally vast and complex due to the presence of multiple phases that differ in the composition and chemical order. To address these obstacles in a computationally efficient manner, an advanced approach combining the mean-field density functional theory with parameters determined using machine learning tools has been developed. This approach was implemented here to investigate AlMo0.5NbTa0.5TiZr, which exhibits a nanostructure consisting of cuboidal BCC precipitates coherently embedded within the B2 matrix. It was found that Al and Zr were responsible for the formation of the B2 matrix. In addition, the matrix and the precipitate were found to have very different elastic characteristics. The matrix has a small elastic moduli and large anisotropy while the precipitate is elastically stiff and nearly isotropic. Beyond the current findings, the parameters for the mean field approach are given in the supplementary material and these can be used in future efforts to predict chemical orders, phase partitioning, and elastic properties of RHESs as a function of chemical composition.

Mass enhancement in 3d and s−p perovskites from symmetry breaking

In some d-electron oxides the measured effective mass m(exptl) has long been known to be significantly larger than the model effective mass m(model) deduced from mean-field band theory, i.e., m(exptl) = beta * m(model), where beta > 1 is the "mass enhancement", or "mass renormalization" factor. Previous applications of density functional theory (DFT), based on the smallest number of possible magnetic, orbital, and structural degrees of freedom, missed such mass enhancement, a fact that was taken as evidence of strong electronic correlation being the exclusive enabling physics. The current paper reports that known modalities of energy-lowering symmetry-broken spin and structural effects included in mean-field DFT show mass enhancement for both electrons and holes in a range of d-electron perovskites SrVO3, SrTiO3, BaTiO3, and LaMnO3 as well as p-electron perovskites CsPbI3 and SrBiO3, including metals (SrVO3) and insulators (the rest). This is revealed only when enlarged unit cells of the same parent global symmetry, which are large enough to allow for symmetry breaking distortions and concomitant variations in spin order, are explored for their ability to lower the total energy. The paper analyzes the contributions of different symmetry-broken modalities to mass enhancement, finds common effects in the range of d- as well as p-electron perovskites. Thus, symmetry breaking in mean-field theory could describe effects that were previously attributed exclusively to complex correlated treatments. Indeed, broken symmetries are often experimentally observed, e.g., octahedral tilting, Jahn-Teller distortions, or bond disproportionation, and often provide intuitive explanations in terms of degeneracy removal due to lowering of structural symmetry, or as shifting of band energies in explainable directions.

Freezing of a soft-core fluid in a one-dimensional potential: appearance of a locked smectic phase

ABSTRACT We investigate the phase behaviour of a two-dimensional colloidal model system of ultra-soft particles on a substrate which varies periodically along one spatial direction. Our calculations are based on mean-field density functional theory for a system of particles interacting via an ultra-soft potential, that is, the generalised exponential model with index four (GEM-4). For suitable substrate periodicities (with commensurability parameter p = 2), we find a succession of phase transitions from a modulated liquid to a locked smectic and then to a locked floating solid phase. The appearance of a locked smectic phase is consistent with earlier theoretical predictions and experiments for freezing of more repulsive systems on structured surfaces (with p = 2). However, the present ultra-soft system does not display re-entrant melting. We here investigate the details of the density distributions of the different phases, thereby supplementing earlier work on GEM-4 systems with p = 1 [Phys. Rev. E 101, 012609 (2020)]. Interestingly, the observed succession of phase transitions can be triggered through different paths along which physical control parameters are changed.

Structural transitions at electrodes, immersed in simple ionic liquid models.

We used a recently developed classical Density Functional Theory (DFT) method to study the structures, phase transitions, and electrochemical behaviours of two coarse-grained ionic fluid models, in the presence of a perfectly conducting model electrode. Common to both is that the charge of the cationic component is able to approach the electrode interface more closely than the anion charge. This means that the cations are specifically attracted to the electrode, due to surface polarization effects. Hence, for a positively charged electrode, there is competition at the surface between cations and anions, where the latter are attracted by the positive electrode charge. This generates demixing, for a range of positive voltages, where the two phases are structurally quite different. The surface charge density is also different between the two phases, even at the same potential. The DFT formulation contains an approximate treatment of ion correlations, and surface polarization, where the latter is modelled via screened image interactions. Using a mean-field DFT, where ion correlations are neglected, causes the phase transition to vanish for both models, but there is still a dramatic drop in the differential capacitance as proximal cations are replaced by anions, for increasing surface potentials. While these findings were obtained for relatively crude coarse-grained models, we argue that the findings can also be relevant in "real" systems, where we note that many ionic liquids are composed of a spherically symmetric anion, and a cation that is asymmetric both from a steric and a charge distribution point of view.

A theory of ordering of elongated and curved proteins on membranes driven by density and curvature.

Cell membranes interact with a myriad of curvature-active proteins that control membrane morphology and are responsible for mechanosensation and mechanotransduction. Some of these proteins, such as those containing BAR domains, are curved and elongated, and hence may adopt different states of orientational order, from isotropic to maximize entropy to nematic as a result of crowding or to adapt to the curvature of the underlying membrane. Here, extending the classical work of Onsager for ordering in hard particle systems and that of [E. S. Nascimento et al., Phys. Rev. E, 2017, 96, 022704], we develop a mean-field density functional theory to predict the orientational order and evaluate the free energy of ensembles of elongated and curved objects on curved membranes. This theory depends on the microscopic properties of the particles and explains how a density-dependent isotropic-to-nematic transition is modified by anisotropic curvature. We also examine the coexistence of isotropic and nematic phases. This theory predicts how ordering depends on geometry but we assume here that the geometry is fixed. It also lays the ground to understand the interplay between membrane reshaping by BAR proteins and molecular order, examined by [Le Roux et al., submitted, 2020].

Toward a density-functional theory for the Jagla fluid.

The so-called Jagla fluid is well known to exhibit, in addition to the usual gas-liquid critical point, also a liquid-liquid critical point, as well as a density anomaly. This makes it an interesting toy model for water, for which a liquid-liquid critical point is considered to exist but so far eludes experimental verification due to crystallization occurring in the corresponding metastable, deeply supercooled state. With the Jagla fluid being understood quite well in bulk-mostly via simulation studies-the focus of the present study is to describe the spatially inhomogeneous fluid in terms of classical density-functional theory (DFT) with the aim to be able to control its phase behavior on changing the shape or the nature of the confinement of the fluid. This information might contribute to guide potential experimental tests of the liquid-liquid critical point of actual water. We first determine the bulk phase diagram for the Jagla fluid by using thermodynamical perturbation theory. In doing so we explain why the perturbation theories of Barker and Henderson as well as of Weeks, Chandler, and Anderson fail to describe the Jagla fluid. We then continue to construct a perturbative DFT based on our bulk model, which shows significant improvement over the standard mean-field DFT valid at high temperatures. But ultimately the perturbative DFT breaks down at state points close to the binodal line and at low temperatures. This prevents us from achieving the original aim to study a highly confined, inhomogeneous Jagla fluid close to its liquid-liquid binodal.

Modeling the Impact of Mesoporous Silica Microstructures on the Adsorption Hysteresis Loop

We present a mean-field density functional theory (MFDFT) study of adsorption and desorption for nitrogen at 77 K in three-dimensional (3D) geometrical models of ordered and random mesoporous silic...

Mott gapping in 3dABO3 perovskites without Mott-Hubbard interelectronic repulsion energy U

The existence of band gaps in Mott insulators such as perovskite oxides with partially filled 3d shells has been traditionally explained in terms of strong, dynamic inter-electronic repulsion codified by the on-site repulsion energy U in the Hubbard Hamiltonian. The success of the "DFT+U approach" where an empirical on-site potential term U is added to the exchange-and correlation Density Functional Theory (DFT) raised questions on whether U in DFT+U represents interelectronic correlation in the same way as it does in the Hubbard Hamiltonian, and if empiricism in selecting U can be avoided. Here we illustrate that ab-initio DFT without any U is able to predict gapping trends and structural symmetry breaking (octahedra rotations, Jahn-Teller modes, bond disproportionation) for all ABO3 3d perovskites from titanates to nickelates in both spin-ordered and spin disordered paramagnetic phases. We describe the paramagnetic phases as a supercell where individual sites can have different local environments thereby allowing DFT to develop finite moments on different sites as long as the total cell has zero moment. We use a recently developed exchange and correlation functional ("SCAN") that is sanctioned by the usual single-determinant, mean-field DFT paradigm with static correlations, but has a more precise rendering of self-interaction cancelation. Our results suggest that strong dynamic electronic correlations are not playing a universal role in gapping of 3d ABO3 Mott insulators, and opens the way for future applications of DFT for studying a plethora of complexity effects that depend on the existence of gaps, such as doping, defects, and band alignment in ABO3 oxides.

Coarse-Grained Modelling and Temperature Effect on the Morphology of PS-b-PI Copolymer.

Spontaneous spatial organization behavior and the aggregate morphology of polystyrene-block-polyisoprene (PS-b-PI) copolymer were investigated. Molecular dynamic (MD) and mesoscopic simulations using the dynamic of mean field density functional theory (DDF) were adopted to investigate the morphology changes exhibited by this block copolymer (BCP). In the mesoscopic simulations, several atoms in repeating units were grouped together into a bead representing styrene or isoprene segments as a coarse-grained model. Inter-bead interactions and essential parameters for mesoscopic models were optimized from MD simulations. Study indicated that morphology alternations can be induced in this system at annealing temperature of 393, 493, and 533 K. From our simulations, lamellar, bicontinuous, and hexagonally packed cylindrical equilibrium morphologies were achieved. Our simulated morphologies agree well with the reported experimental evidence at the selected temperature. The process of aggregate formation and morphology evolution were concretely clarified.

Random phase approximation for the non-uniform Yukawa fluid

Mean-field density functional theory can be used to estimate the free energy of non-uniform fluids. The second functional derivative with respect to density of the free energy is related to the direct correlation function of the fluid and, in principle, this can be inverted to find an improved approximation for the pair correlation function and hence the free energy, the so-called ‘random phase approximation’. If the repulsive molecular interaction is approximated by the local density approximation and the attractive interaction is assumed to be of the Yukawa form, the problem reduces to that of finding the eigenvalues of Schrödinger-like equations, which, for certain models (such as the ‘Φ4 model’), can be done analytically in the planar case. The relationship between this approach and field theoretical treatment of the vapour-liquid interface is discussed. The ultraviolet divergence of the expression can be eliminated by separating the first term in the expansion, although quantitative results still depend on the behaviour of the attractive potential in the repulsive core. In the case of a spherical droplet of radius R, correction terms to the free energy involving lnR appear due to (i) cluster translational invariance, (ii) the unstable mode corresponding to droplet growth, and (iii) capillary waves. The net effect of these terms is to modify the classical expression for the nucleation rate by a factor proportional to R4/3.

Adsorption and interfacial phenomena of a Lennard-Jones fluid adsorbed in slit pores: DFT and GCMC simulations

ABSTRACTConfinement of fluids in porous media leads to the presence of solid–fluid (SF) interfaces that play a key role in many different fields. The experimental characterisation of SF interfacial properties, in particular the surface tension, is challenging or not accessible. In this work, we apply mean-field density functional theory (DFT) to determine the surface tension and also density profile of a Lennard-Jones fluid in slit-shaped pores for realistic amounts of adsorbed molecules. We consider the pore walls to interact with fluid molecules through the well-known 10-4-3 Steele potential. The results are compared with those obtained from Monte Carlo simulations in the Grand Canonical Ensemble (GCMC) using the test-area method. We analyse the effect on the adsorption and interfacial phenomena of volume and energy factors, in particular, the pore diameter and the ratio between SF and fluid–fluid dispersive energy parameters, respectively. Results from DFT and GCMC simulations were found to be comparable, which points to their reliability.

Mean-field density functional theory of a nanoconfined classical, three-dimensional Heisenberg fluid. II. The interplay between molecular packing and orientational order.

As in Paper I of this series of papers [S. M. Cattes et al., J. Chem. Phys. 144, 194704 (2016)], we study a Heisenberg fluid confined to a nanoscopic slit pore with smooth walls. The pore walls can either energetically discriminate specific orientations of the molecules next to them or are indifferent to molecular orientations. Unlike in Paper I, we employ a version of classical density functional theory that allows us to explicitly account for the stratification of the fluid (i.e., the formation of molecular layers) as a consequence of the symmetry-breaking presence of the pore walls. We treat this stratification within the White Bear version (Mark I) of fundamental measure theory. Thus, in this work, we focus on the interplay between local packing of the molecules and orientational features. In particular, we demonstrate why a critical end point can only exist if the pore walls are not energetically discriminating specific molecular orientations. We analyze in detail the positional and orientational order of the confined fluid and show that reorienting molecules across the pore space can be a two-dimensional process. Last but not least, we propose an algorithm based upon a series expansion of Bessel functions of the first kind with which we can solve certain types of integrals in a very efficient manner.

Morphology of Fluids Confined in Physically Reconstructed Mesoporous Silica: Experiment and Mean Field Density Functional Theory.

Three-dimensional physical reconstruction of the random mesopore network in a hierarchically structured, macroporous-mesoporous silica monolith via electron tomography has been used to generate a lattice model of amorphous, mesoporous silica. This geometrical model has subsequently been employed in mean field density functional theory (MFDFT) calculations of adsorption and desorption. Comparison of the results with experimental sorption isotherms for nitrogen at 77 K shows a good qualitative agreement, with both experiment and theory producing isotherms characterized by type H2 hysteresis. In addition to the isotherms, MFDFT provides the three-dimensional density distribution for the fluid in the porous material for each state studied. We use this information to map the phase distribution in the mesopore network in the hysteresis region. Phase distributions on the desorption boundary curve are compared to those on the adsorption boundary curve for both fixed pressure and fixed density, revealing insights into the relationship between phase distribution and hysteresis.