Source author record

James F. Lutsko

James F. Lutsko appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

28works
7topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

28 published item(s)

preprint2021arXiv

Classical Density Functional Theory in the Canonical Ensemble

Classical density functional theory for finite temperatures is usually formulated in the grand-canonical ensemble where arbitrary variations of the local density are possible. However, in many cases the systems of interest are closed with respect to mass, e.g. canonical systems with fixed temperature and particle number. Although the tools of standard, grand-canonical density functional theory are often used in an ad hoc manner to study closed systems, their formulation directly in the canonical ensemble has so far not been known. In this work, the fundamental theorems underlying classical DFT are revisited and carefully compared in the two ensembles showing that there are only trivial formal differences. The practicality of DFT in the canonical ensemble is then illustrated by deriving the exact Helmholtz functional for several systems: the ideal gas, certain restricted geometries in arbitrary numbers of dimensions and finally a system of two hard-spheres in one dimension (hard rods) in a small cavity. Some remarkable similarities between the ensembles are apparent even for small systems with the latter showing strong echoes of the famous exact of result of Percus in the grand-canonical ensemble.

preprint2020arXiv

Classical Density Functional Theory applied to the solid state

The standard model of classical Density Functional Theory for pair potentials consists of a hard-sphere functional plus a mean-field term accounting for long ranged attraction. However, most implementations using sophisticated Fundamental Measure hard-sphere functionals suffer from potential numerical instabilities either due to possible instabilities in the functionals themselves or due to implementations that mix real- and Fourier-space components inconsistently. Here, we present a new implementation based on a demonstrably stable hard-sphere functional that is implemented in a completely consistent manner. The present work does not depend on approximate spherical integration schemes and so is much more robust than previous algorithms. The methods are illustrated by calculating phase diagrams for the solid state using the standard Lennard-Jones potential as well as a new class of potentials recently proposed by Wang et al (Phys. Chem. Chem. Phys. 22, 10624 (2020)). The latter span the range from potentials for small molecules to those appropriate to colloidal systems simply by varying a parameter. We verify that cDFT is able to semi-quantitatively reproduce the phase diagram in all cases. We also show that for these problems computationally cheap Gaussian approximations are nearly as good as full minimization based on finite differences.

preprint2020arXiv

Explicitly stable Fundamental Measure Theory models for classical density functional theory

The derivation of the state of the art tensorial versions of Fundamental Measure Theory (a form of classical Density Functional Theory for hard spheres) are re-examined in the light of the recently introduced concept of global stability of the density functional based on its boundedness (Lutsko and Lam, Phys. Rev. E 98, 012604 (2018)). It is shown that within the present paradigm, explicitly stability of the functional can be achieved only at the cost of giving up accuracy at low densities. It is argued that this is an acceptable trade-off since the main value of DFT lies in the study of dense systems. Explicit calculations for a wide variety of systems shows that a proposed explicitly stable functional is competitive in all ways with the popular White Bear model while sharing some of its weaknesses when applied to non-close-packed solids.

preprint2020arXiv

Long-wavelength density fluctuations as nucleation precursors

Recent theories of nucleation that go beyond Classical Nucleation Theory predict that diffusion-limited nucleation of both liquid droplets and of crystals from a low-density vapor (or weak solution) begins with long-wavelength density fluctuations. This means that in the early stages of nucleation, 'clusters' can have low density but large spatial extent, which is at odds with the classical picture of arbitrarily small clusters of the condensed phase. We present the results of kinetic Monte Carlo simulations using Forward Flux Sampling to show that these predictions are confirmed: namely that on average, nucleation begins in the presence of low-amplitude, but spatially extended density fluctuations thus confirming a significant prediction of the non-classical theory.

preprint2016arXiv

Mechanism for the stabilization of protein clusters above the solubility curve

Pan, Vekilov and Lubchenko[\textit{J. Phys. Chem. B}, 2010, \textbf{114}, 7620] have proposed that dense stable protein clusters appearing in weak protein solutions above the solubility curve are composed of protein oligomers. The hypothesis is that a weak solution of oligomer species is unstable with respect to condensation causing the formation of dense, oligomer-rich droplets which are stabilized against growth by the monomer-oligomer reaction. Here, we show that such a combination of processes can be understood using a simple capillary model yielding analytic expressions for the cluster properties which can be used to interpret experimental data. We also construct a microscopic Dynamic Density Functional Theory model and show that it is consistent with the predictions of the capillary model. The viability of the mechanism is thus confirmed and it is shown how the radius of the stable clusters is related to physically interesting quantities such as the monomer-oligomer rate constants.

preprint2016arXiv

Mechanism for the stabilization of protein clusters above the solubility curve: the role of non-ideal chemical reactions

Dense protein clusters are known to play an important role in nucleation of protein crystals from dilute solutions. While these have generally been thought to be formed from a metastable phase, the observation of similar, if not identical, clusters above the critical point for the dilute-solution/strong-solution phase transition has thrown this into doubt. Furthermore, the observed clusters are stable for relatively long times. Because protein aggregation plays an important role in some pathologies, understanding the nature of such clusters is an important problem. One mechanism for the stabilization of such structures was proposed by Pan, Vekilov and Lubchenko and was investigated using a DDFT model which confirmed the viability of the model. Here, we revisit that model and incorporate additional physics in the form of state-dependent reaction rates. We show by a combination of numerical results and general arguments that the state-dependent rates disrupt the stability mechanism. Finally, we argue that the state-depedent reactions correct unphysical aspects of the model with ideal (state-independent) reactions and that this necessarily leads to the failure of the proposed mechanism.

preprint2016arXiv

Temporal Fokker-Planck Equations

The temporal Fokker-Plank equation [{\it J. Stat. Phys.}, {\bf 3/4}, 527 (2003)] or propagation-dispersion equation was derived to describe diffusive processes with temporal dispersion rather than spatial dispersion as in classical diffusion. %\cite{boon-grosfils-lutsko}. We present two generalizations of the temporal Fokker-Plank equation for the first passage distribution function $f_j(r,t)$ of a particle moving on a substrate with time delays $τ_j$. Both generalizations follow from the first visit master equation. In the first case, the time delays depend on the local concentration, that is the time delay probability $P_j$ is a functional of the particle distribution function and we show that when the functional dependence is of the power law type, $P_j \propto f_j^{ν- 1}$, the generalized Fokker-Plank equation exhibits a structure similar to that of the nonlinear spatial diffusion equation where the roles of space and time are reversed. In the second case, we consider the situation where the time delays are distributed according to a power law, $P_j \propto τ_j^{-1-α}$ (with $0 < α< 2$), in which case we obtain a fractional propagation-dispersion equation which is the temporal analog of the fractional spatial diffusion equation (with space and time interchanged).

preprint2015arXiv

Unification of classical nucleation theories via unified Itô-Stratonovich stochastic equation

Classical nucleation theory (CNT) is the most widely used framework to describe the early stage of first-order phase transitions. Unfortunately the different points of view adopted to derive it yield different kinetic equations for the probability density function, e.g. Zeldovich-Frenkel or Becker-Döring-Tunitskii equations. Starting from a phenomenological stochastic differential equation a unified equation is obtained in this work. In other words, CNT expressions are recovered by selecting one or another stochastic calculus. Moreover, it is shown that the unified CNT thus obtained produces the same Fokker-Planck equation as that from a recent update of CNT [J.F. Lutsko and M.A. Durán-Olivencia, J. Chem. Phys., 2013, 138, 244908] when mass transport is governed by diffusion. Finally, we derive a general induction-time expression along with specific approximations of it to be used under different scenarios. In particular, when the mass-transport mechanism is governed by direct impingement, volume diffusion, surface diffusion or interface transfer.

preprint2014arXiv

A Two-parameter Extension of Classical Nucleation Theory

A two-variable stochastic model for diffusion-limited nucleation is developed using a formalism derived from fluctuating hydrodynamics. The model is a direct generalization of the standard Classical Nucleation Theory. The nucleation rate and pathway are calculated in the weak-noise approximation and are shown to be in good agreement with direct numerical simulations for the weak-solution/strong-solution transition in globular proteins. We find that Classical Nucleation Theory underestimates the time needed for the formation of a critical cluster by two orders of magnitude and that this discrepancy is due to the more complex dynamics of the two variable model and not, as often is assumed, a result of errors in the estimation of the free energy barrier.

preprint2014arXiv

Molecular theory of anomalous diffusion - Application to Fluorescence Correlation Spectroscopy

The nonlinear theory of anomalous diffusion is based on particle interactions giving an explicit microscopic description of diffusive processes leading to sub-, normal, or super-diffusion as a result competitive effects between attractive and repulsive interactions. We present the explicit analytical solution to the nonlinear diffusion equation which we then use to compute the correlation function which is experimentally measured by correlation spectroscopy. The theoretical results are applicable in particular to the analysis of fluorescence correlation spectroscopy of marked molecules in biological systems. More specifically we consider the case of fluorescently labeled lipids and we find that the nonlinear correlation spectrum reproduces very well the experimental data indicating sub-diffusive molecular motion of lipid molecules in the cell membrane.

preprint2014arXiv

New comments on "Possible divergences in Tsallis' thermostatistics''

In a recent letter ({\it{EPL}}, {\bf{104}} (2013) 60003; see also {\it {arXiv:1309.5645}}), Plastino and Rocca suggest that the divergences inherent to the formulation of nonextensive statistical mechanics can be eliminated {\it {via}} the use of $q$-Laplace transformation which is illustrated for the case of a kinetic Hamiltonian system, the harmonic oscillator. The suggested new formulation raises questions which are discussed in the present comment.

preprint2013arXiv

Molecular theory of anomalous diffusion

We present a Master Equation formulation based on a Markovian random walk model that exhibits sub-diffusion, classical diffusion and super-diffusion as a function of a single parameter. The non-classical diffusive behavior is generated by allowing for interactions between a population of walkers. At the macroscopic level, this gives rise to a nonlinear Fokker-Planck equation. The diffusive behavior is reflected not only in the mean-squared displacement ($<r^2(t)>\sim t^γ$ with $0 <γ\leq 1.5$) but also in the existence of self-similar scaling solutions of the Fokker-Planck equation. We give a physical interpretation of sub- and super-diffusion in terms of the attractive and repulsive interactions between the diffusing particles and we discuss analytically the limiting values of the exponent $γ$. Simulations based on the Master Equation are shown to be in agreement with the analytical solutions of the nonlinear Fokker-Planck equation in all three diffusion regimes.

preprint2012arXiv

Classical nucleation theory from a dynamical approach to nucleation

It is shown that diffusion-limited classical nucleation theory (CNT) can be recovered as a simple limit of the recently proposed dynamical theory of nucleation based on fluctuating hydrodynamics (Lutsko, JCP 136, 034509 (2012)). The same framework is also used to construct a more realistic theory in which clusters have finite interfacial width. When applied to the dilute solution/dense solution transition in globular proteins, it is found that the extension gives corrections to the the nucleation rate even for the case of small supersaturations due to changes in the monomer distribution function and to the excess free energy. It is also found that the monomer attachement/detachment picture breaks down at high supersaturations corresponding to clusters smaller than about 100 molecules. The results also confirm the usual assumption that most important corrections to CNT can be acheived by means of improved estimates of the free energy barrier. The theory also illustrates two topics that have received considerable attention in the recent literature on nucleation: the importance sub-dominant corrections to the capillary model for the free energy and of the correct choice of the reaction coordinate.

preprint2012arXiv

Nucleation of colloids and macromolecules in a finite volume

A recently formulated description of homogeneous nucleation for Brownian particles in the over-damped limit based on fluctuating hydrodynamics is used to determine the nucleation pathway, characterized as the most likely path (MLP), for the nucleation of a dense-concentration droplet of globular protein from a dilute solution in a small, finite container. The calculations are performed by directly discretizing the equations for the MLP and it is found that they confirm previous results obtained for infinite systems: the process of homogeneous nucleation begins with a long-wavelength, spatially-extended concentration fluctuation that it condenses to form the pre-critical cluster. This is followed by a classical growth processes. The calculations show that the post-critical growth involves the formation of a depletion zone around the cluster whereas no such depletion is observed in the pre-critical cluster. The approach therefore captures dynamical effects not found in classical Density Functional Theory studies while consistently describing the formation of the pre-critical cluster.

preprint2012arXiv

Nucleation of colloids and macromolecules: does the nucleation pathway matter?

A recent description of diffusion-limited nucleation based on fluctuating hydrodynamics that extends classical nucleation theory predicts a very non-classical two-step scenario whereby nucleation is most likely to occur in spatially-extended, low-amplitude density fluctuations. In this paper, it is shown how the formalism can be used to determine the maximum probability of observing \emph{any} proposed nucleation pathway, thus allowing one to address the question as to their relative likelihood, including of the newly proposed pathway compared to classical scenarios. Calculations are presented for the nucleation of high-concentration bubbles in a low-concentration solution of globular proteins and it is found that the relative probabilities (new theory compared to classical result) for reaching a critical nucleus containing $N_c$ molecules scales as $e^{-N_c/3}$ thus indicating that for all but the smallest nuclei, the classical scenario is extremely unlikely.

preprint2011arXiv

A dynamical theory of homogeneous nucleation for colloids and macromolecules

Homogeneous nucleation is formulated within the context of fluctuating hydrodynamics. It is shown that for a colloidal or macromolecular system in the strong damping limit the most likely path for nucleation can be determined by gradient descent in density space governed by a nontrivial metric fixed by the dynamics. The theory provides a justification and extension of more heuristic equilibrium approaches based solely on the free energy. It is illustrated by application to liquid-vapor nucleation where it is shown that, in contrast to most free energy-based studies, the smallest clusters correspond to long wavelength, small amplitude perturbations.

preprint2011arXiv

A dynamical theory of nucleation for colloids and macromolecules

A general theory of nucleation for colloids and macromolecules in solution is formulated within the context of fluctuating hydrodynamics. A formalism for the determination of nucleation pathways is developed and stochastic differential equations for the evolution of order parameters are given. The conditions under which the elements of Classical Nucleation Theory are recovered are determined. The theory provides a justification and extension of more heuristic equilibrium approaches based solely on the free energy. It is illustrated by application to the low-concentration/high-concentration transition in globular proteins where a novel two-step mechanism is identified where the first step involves the formation of long-wavelength density fluctuations and the second step is the actual nucleation event occuring within the fluctuation.

preprint2011arXiv

A microscopic approach to nonlinear Reaction-Diffusion: the case of morphogen gradient formation

We develop a microscopic theory for reaction-difusion (R-D) processes based on a generalization of Einstein's master equation with a reactive term and we show how the mean field formulation leads to a generalized R-D equation with non-classical solutions. For the $n$-th order annihilation reaction $A+A+A+...+A\rightarrow 0$, we obtain a nonlinear reaction-diffusion equation for which we discuss scaling and non-scaling formulations. We find steady states with either solutions exhibiting long range power law behavior (for $n>α$) showing the relative dominance of sub-diffusion over reaction effects in constrained systems, or conversely solutions (for $n<α<n+1$) with finite support of the concentration distribution describing situations where diffusion is slow and extinction is fast. Theoretical results are compared with experimental data for morphogen gradient formation.

preprint2011arXiv

Density Functional Theory of Inhomogeneous Liquids IV. Squared-gradient approximation and Classical Nucleation Theory

The Squared-Gradient approximation to the Modified-Core Van der Waals density functional theory model is developed. A simple, explicit expression for the SGA coefficient involving only the bulk equation of state and the interaction potential is given. The model is solved for planar interfaces and spherical clusters and is shown to be quantitatively accurate in comparisons to computer simulations. An approximate technique for solving the SGA based on piecewise-linear density profiles is introduced and is shown to give reasonable zeroth-order approximations to the numerical solution of the model. The piecewise-linear models of spherical clusters are shown to be a natural extension of Classical Nucleation Theory and serve to clarify some of the non-classical effects previously observed in liquid-vapor nucleation. Nucleation pathways are investigated using both constrained energy-minimization and steepest-descent techniques.

preprint2011arXiv

On the role of metastable intermediate states in the homogeneous nucleation of solids from solution

The role of metastable liquid phases in vapor-crystal nucleation is studied using Density Functional Theory(DFT). The model gives a semi-quantitatively accurate description of both the vapor-liquid-solid phase diagram for both simple fluids (Lennard-Jones interactions) and of the low-density/high-density/crystal phase diagram for model globular proteins (ten Wolde-Frenkel interaction). The density profile is characterized by two local order parameters, the average density and the crystallinity. The bulk free energy model is supplemented by squared-derivative terms in these order parameters to account for inhomogeneities thus producing a model similar in spirit to phase-field theory. It is shown that for both interaction models, the vapor-crystal part of the phase-diagram can be separated into regions for which metastable liquid phases are more or less stable than the vapor, but always less stable than the solid. The former case allows for the possibility of \emph{double nucleation} whereby liquid droplets nucleate from the vapor followed by a separate nucleation of the solid phase within the liquid droplets. Whether or not this actually occurs depends on the relative free energy barriers for vapor-solid and vapor-liquid nucleation and it is shown that for simple fluids, double nucleation is indeed favored at sufficiently large supersaturation. Finally, by studying the minimum free energy path from the vapor to the solid, the separate possibility of \emph{transient nucleation}, where the vapor-solid transition involves a single nucleation event but the sub-critical clusters tend to be liquid-like while small, is shown to also be possible when the metastable liquid exists, but the supersaturation is too small for double nucleation.

preprint2010arXiv

Kinetic roughening-like transition with finite nucleation barrier

Recent observations of the growth of protein crystals have identified two different growth regimes. At low supersaturation, the surface of the crystal is smooth and increasing in size due to the nucleation of steps at defects and the subsequent growth of the steps. At high supersaturation, nucleation occurs at many places simultaneously, the crystal surface becomes rough and the growth velocity increases more rapidly with increasing supersaturation than in the smooth regime. Kinetic roughening transitions are typically assumed to be due to the vanishing of the barrier for two-dimension nucleation on the surface of the crystal. We show here, by means of both analytic mean field models and kinetic Monte Carlo simulations that a transition between different growth modes reminiscent of kinetic roughening can also arise as a kinetic effect occurring at finite nucleation barriers.

preprint2010arXiv

Kinetics of intermediate-mediated self-assembly in nano-sized materials: a generic model

We propose in this paper a generic model of a non-standard aggregation mechanism for self-assembly processes of a class of materials involving the mediation of intermediates consisting of a polydisperse population of nano-sized particles. The model accounts for a long induction period in the process. The proposed mechanism also gives insight on future experiments aiming at a more comprehensive picture of the role of self-organization in self-assembly processes.

preprint2010arXiv

The low-density/high-density liquid phase transition for model globular proteins

The effect of molecule size (excluded volume) and the range of interaction on the surface tension, phase diagram and nucleation properties of a model globular protein is investigated using a combinations of Monte Carlo simulations and finite temperature classical Density Functional Theory calculations. We use a parametrized potential that can vary smoothly from the standard Lennard-Jones interaction characteristic of simple fluids, to the ten Wolde-Frenkel model for the effective interaction of globular proteins in solution. We find that the large excluded volume characteristic of large macromolecules such as proteins is the dominant effect in determining the liquid-vapor surface tension and nucleation properties. The variation of the range of the potential only appears important in the case of small excluded volumes such as for simple fluids. The DFT calculations are then used to study homogeneous nucleation of the high-density phase from the low-density phase including the nucleation barriers, nucleation pathways and the rate. It is found that the nucleation barriers are typically only a few $k_{B}T$ and that the nucleation rates substantially higher than would be predicted by Classical Nucleation Theory.

preprint2009arXiv

Is the Tsallis entropy stable?

The question of whether the Tsallis entropy is Lesche-stable is revisited. It is argued that when physical averages are computed with the escort probabilities, the correct application of the concept of Lesche-stability requires use of the escort probabilities. As a consequence, as shown here, the Tsallis entropy is unstable but the thermodynamic averages are stable. We further show that Lesche stability as well as thermodynamic stability can be obtained if the homogeneous entropy is used as the basis of the formulation of non-extensive thermodynamics. In this approach, the escort distribution arises naturally as a secondary structure.

preprint2009arXiv

Phase behavior of a confined nano-droplet in the grand-canonical ensemble: the reverse liquid-vapor transition

The equilibrium density distribution and thermodynamic properties of a Lennard-Jones fluid confined to nano-sized spherical cavities at constant chemical potential was determined using Monte Carlo simulations. The results describe both a single cavity with semipermeable walls as well as a collection of closed cavities formed at constant chemical potential. The results are compared to calculations using classical Density Functional Theory (DFT). It is found that the DFT calculations give a quantitatively accurate description of the pressure and structure of the fluid. Both theory and simulation show the presence of a ``reverse'' liquid-vapor transition whereby the equilibrium state is a liquid at large volumes but becomes a vapor at small volumes.

preprint2005arXiv

Theoretical evidence for a dense fluid precursor to crystallization

We present classical density functional theory calculations of the free energy landscape for fluids below their triple point as a function of density and crystallinity. We find that for both a model globular protein and for a simple atomic fluid modeled with a Lennard-Jones interaction, it is free-energetically easier to crystallize by passing through a metastable dense fluid in accord with the Ostwald rule of stages but in contrast to the alternative of ordering and densifying at once as assumed in the classical picture of crystallization.