Free energy determination of phase coexistence in model C60: A comprehensive Monte Carlo study

The free energy of the solid and fluid phases of the Girifalco C60 model are determined through extensive Monte Carlo simulations. In this model the molecules interact through a spherical pair potential, characterized by a narrow and attractive well, adjacent to a harshly repulsive core. We have used the Widom test particle method and a mapping from an Einstein crystal, in order to estimate the absolute free energy in the fluid and solid phases, respectively; we have then determined the free energy along several isotherms, and the whole phase diagram, by means of standard thermodynamic integrations. We highlight how the interplay between the liquid-vapor and the liquid-solid coexistence conditions determines the existence of a narrow liquid pocket in the phase diagram, whose stability is assessed and confirmed in agreement with previous studies. In particular, the critical temperature follows closely an extended corresponding-states rule recently outlined by Noro and Frenkel [J. Chem. Phys. 113:2941 (2000)]. We discuss the emerging "energetic" properties of the system, which drive the phase behavior in systems interacting through short-range forces [A. A. Louis, Phil. Trans. R. Soc. A 359:939 (2001)], in order to explain the discrepancy between the predictions of several structural indicators and the results of full free energy calculations, to locate the fluid phase boundaries. More generally, we aim to provide extended reference data for calculations of the free energy of the C60 fullerite in the low temperature regime, as for the determination of the phase diagram of higher order fullerenes and other fullerene-related materials, whose description is based on the same model adopted in this work.Comment: RevTeX, 11 pages, 9 figure

Thermodynamic stability of fluid-fluid phase separation in binary athermal mixtures: The role of nonadditivity

We study the thermodynamic stability of fluid-fluid phase separation in binary nonadditive mixtures of hard-spheres for moderate size ratios. We are interested in elucidating the role played by small amounts of nonadditivity in determining the stability of fluid-fluid phase separation with respect to the fluid-solid phase transition. The demixing curves are built in the framework of the modified-hypernetted chain and of the Rogers-Young integral equation theories through the calculation of the Gibbs free energy. We also evaluate fluid-fluid phase equilibria within a first-order thermodynamic perturbation theory applied to an effective one-component potential obtained by integrating out the degrees of freedom of the small spheres. A qualitative agreement emerges between the two different approaches. We also address the determination of the freezing line by applying the first-order thermodynamic perturbation theory to the effective interaction between large spheres. Our results suggest that for intermediate size ratios a modest amount of nonadditivity, smaller than earlier thought, can be sufficient to drive the fluid-fluid critical point into the thermodinamically stable region of the phase diagram. These findings could be significant for rare-gas mixtures in extreme pressure and temperature conditions, where nonadditivity is expected to be rather small.Comment: 17 pages, 7 figures, to appear in J. Phys. Chem.

Structure and phase behavior of colloidal dumbbells with tunable attractive interactions

We investigate thermodynamic and structural properties of colloidal dumbbells in the framework provided by the Reference Interaction Site Model (RISM) theory of molecular fluids and Monte Carlo simulations. We consider two different models: in the first one we set identical square-well attractions on the two tangent spheres composing the molecule (SW-SW model); in the second scheme, one of square-well interactions is switched off (HS-SW model). Appreciable differences emerge between the physical properties of the two models. Specifically, the $k \to 0$ behavior of SW-SW structure factors $S(k)$ points to the presence of a gas-liquid coexistence, as confirmed by subsequent fluid phase equilibria calculations. Conversely, the HS-SW $S(k)$ develops a low-$k$ peak, signaling the presence of aggregates; such a process destabilizes the gas-liquid phase separation, promoting at low temperatures the formation of a cluster phase, whose structure depends on the system density. We further investigate such differences by studying the phase behavior of a series of intermediate models, obtained from the original SW-SW by progressively reducing the depth of one square-well interaction. RISM structural predictions positively reproduce the simulation data, including the rise of $S(k \to 0$) in the SW-SW model and the low-$k$ peak in the HS-SW structure factor. As for the phase behavior, RISM agrees with Monte Carlo simulations in predicting a gas-liquid coexistence for the SW-SW model (though the critical parameters appears overestimated by the theory) and its progressive disappearance moving toward the HS-SW model.Comment: 12 pages, 13 figures, 1 table, 78 reference

Theory and simulation of short-range models of globular protein solutions

We report theoretical and simulation studies of phase coexistence in model globular protein solutions, based on short-range, central, pair potential representations of the interaction among macro-particles. After reviewing our previous investigations of hard-core Yukawa and generalised Lennard-Jones potentials, we report more recent results obtained within a DLVO-like description of lysozyme solutions in water and added salt. We show that a one-parameter fit of this model based on Static Light Scattering and Self-Interaction Chromatography data in the dilute protein regime, yields demixing and crystallization curves in good agreement with experimental protein-rich/protein-poor and solubility envelopes. The dependence of cloud and solubility points temperature of the model on the ionic strength is also investigated. Our findings highlight the minimal assumptions on the properties of the microscopic interaction sufficient for a satisfactory reproduction of the phase diagram topology of globular protein solutions.Comment: 17 pages, 8 figures, Proc. of Conference "Structural Arrest Transitions in Colloidal Systems with Short-Range Attractions", Messina (ITALY) 17-20 December 200

Theoretical description of phase coexistence in model C60

We have investigated the phase diagram of the Girifalco model of C60 fullerene in the framework provided by the MHNC and the SCOZA liquid state theories, and by a Perturbation Theory (PT), for the free energy of the solid phase. We present an extended assessment of such theories as set against a recent Monte Carlo study of the same model [D. Costa et al, J. Chem. Phys. 118:304 (2003)]. We have compared the theoretical predictions with the corresponding simulation results for several thermodynamic properties. Then we have determined the phase diagram of the model, by using either the SCOZA, or the MHNC, or the PT predictions for one of the coexisting phases, and the simulation data for the other phase, in order to separately ascertain the accuracy of each theory. It turns out that the overall appearance of the phase portrait is reproduced fairly well by all theories, with remarkable accuracy as for the melting line and the solid-vapor equilibrium. The MHNC and SCOZA results for the liquid-vapor coexistence, as well as for the corresponding critical points, are quite accurate. All results are discussed in terms of the basic assumptions underlying each theory. We have selected the MHNC for the fluid and the first-order PT for the solid phase, as the most accurate tools to investigate the phase behavior of the model in terms of purely theoretical approaches. The overall results appear as a robust benchmark for further theoretical investigations on higher order C(n>60) fullerenes, as well as on other fullerene-related materials, whose description can be based on a modelization similar to that adopted in this work.Comment: RevTeX4, 15 pages, 7 figures; submitted to Phys. Rev.

Using the Process Digital Twin as a tool for companies to evaluate the Return on Investment of manufacturing automation

The fourth industrial revolution is gaining momentum, but still lacks full realization. Several studies suggest that many companies around the world have begun the digital transformation undertaking, but most are still far from full adoption and yet fail to see the full economic potential, being stuck in what has been called "pilot purgatory”. Digitalization is largely recognized as an accelerator and enabler for full automation in manufacturing, but companies are still struggling to assess the return on investment and the impact on operational performance indicators. Therefore, companies, especially SMEs characterized by dynamic, high-value, high-mix, and low-volume contexts, are reluctant to invest further. By incorporating simulation, data analytics and behavioral models, digital twins may also be used to support automation solutions ramp-up, demonstrate their impact evaluation, usage scenarios, eliminating the need for physical prototypes, reducing development time, and improving quality. Few forward-thinking companies are pursuing the digital transformation path, while the majority are clipping the wings of a transformation that is essential for a sustainable manufacturing. This paper describes a theoretical approach to exploit the digital twin technology to gather insights towards a realistic economical assessment of full automation solutions, to back and encourage investments to realize the potential of the digital manufacturing transformation. The approach is being tested under the European Union’s Horizon 2020 research and innovation program under grant agreement No. 958363, which provides an opportunity to assess how the various components of the method are constructed, how complex they are, and what level of effort is required, using a practical example