Tethered Monte Carlo: computing the effective potential without critical slowing down

We present Tethered Monte Carlo, a simple, general purpose method of computing the effective potential of the order parameter (Helmholtz free energy). This formalism is based on a new statistical ensemble, closely related to the micromagnetic one, but with an extended configuration space (through Creutz-like demons). Canonical averages for arbitrary values of the external magnetic field are computed without additional simulations. The method is put to work in the two-dimensional Ising model, where the existence of exact results enables us to perform high precision checks. A rather peculiar feature of our implementation, which employs a local Metropolis algorithm, is the total absence, within errors, of critical slowing down for magnetic observables. Indeed, high accuracy results are presented for lattices as large as L = 1024.

First-order transition in a three-dimensional disordered system

We present the first detailed numerical study in three dimensions of a first-order phase transition that remains first order in the presence of quenched disorder (specifically, the ferromagnetic-paramagnetic transition of the site-diluted four states Potts model). A tricritical point, which lies surprisingly near the pure-system limit and is studied by means of finite-size scaling, separates the first-order and second-order parts of the critical line. This investigation has been made possible by a new definition of the disorder average that avoids the diverging-variance probability distributions that plague the standard approach. Entropy, rather than free energy, is the basic object in this approach that exploits a recently introduced microcanonical Monte Carlo method.

First Order Phase Transition in a 3D disordered system

We present a detailed numerical study on the effects of adding quenched impurities to a three dimensional system which in the pure case undergoes a strong first order phase transition (specifically, the ferromagnetic/paramagnetic transition of the site-diluted four states Potts model). We can state that the transition remains first-order in the presence of quenched disorder (a small amount of it) but it turns out to be second order as more impurities are added. A tricritical point, which is studied by means of Finite-Size Scaling, separates the first-order and second-order parts of the critical line. The results were made possible by a new definition of the disorder average that avoids the diverging-variance probability distributions that arise using the standard methodology. We also made use of a recently proposed microcanonical Monte Carlo method in which entropy, instead of free energy, is the basic quantity.

Critical behavior of su(1|1) supersymmetric spin chains with long-range interactions

We introduce a general class of su(1|1) supersymmetric spin chains with long-range interactions which includes as particular cases the su(1|1) Inozemtsev (elliptic) and Haldane-Shastry chains, as well as the XX model. We show that this class of models can be fermionized with the help of the algebraic properties of the su(1|1) permutation operator and take advantage of this fact to analyze their quantum criticality when a chemical potential term is present in the Hamiltonian. We first study the low-energy excitations and the low-temperature behavior of the free energy, which coincides with that of a (1+1)-dimensional conformal field theory (CFT) with central charge c=1 when the chemical potential lies in the critical interval (0,E(π)), E(p) being the dispersion relation. We also analyze the von Neumann and Rényi ground state entanglement entropies, showing that they exhibit the logarithmic scaling with the size of the block of spins characteristic of a one-boson (1+1)-dimensional CFT. Our results thus show that the models under study are quantum critical when the chemical potential belongs to the critical interval, with central charge c=1. From the analysis of the fermion density at zero temperature, we also conclude that there is a quantum phase transition at both ends of the critical interval. This is further confirmed by the behavior of the fermion density at finite temperature, which is studied analytically (at low temperature), as well as numerically for the su(1|1) elliptic chain.

The equilibrium between carbon oxysulfide, carbon monoxide and sulfur ...

Thesis t.-p. pasted over the cover of the issue reprinted from the Journal of the American chemical society, vol. XXXVII, no. 9, September, 1915, which has title: The equilibrium between carbon oxysulfide, carbon monoxide and sulfur, by Gilbert N. Lewis and William N. Lacey.

Mechanical alloying of MoSi2 with ternary alloying elements. Part 2 - computer simulation

A computer model of the mechanical alloying process has been developed to simulate phase formation during the mechanical alloying of Mo and Si elemental powders with a ternary addition of Al, Mg, Ti or Zr. Using the Arhennius equation, the model balances the formation rates of the competing reactions that are observed during milling. These reactions include the formation of tetragonal C11(b) MOSi2 (t-MoSi2) by combustion, the formation of the hexagonal C40 MoSi2 polymorph (h-MoSi2), the transformation of the tetragonal to the hexagonal form, and the recovery of t-MoSi2 from h-MoSi2 and deformed t-MoSi2. The addition of the ternary additions changes the free energy of formation of the associated MoSi2 alloys, i.e. Mo(Si, Al)(2), Mo(Mg, Al)(2), (Mo, Ti)Si-2 (Mo, Zr)Si-2 and (Mo, Fe)Si-2, respectively. Variation of the energy of formation alone is sufficient for the simulation to accurately model the observed phase formation. (C) 2003 Elsevier B.V. All rights reserved.

Application of density functional theory to equilibrium adsorption of argon and nitrogen on amorphous silica surface

We present a new version of non-local density functional theory (NL-DFT) adapted to description of vapor adsorption isotherms on amorphous materials like non-porous silica. The novel feature of this approach is that it accounts for the roughness of adsorbent surface. The solid–fluid interaction is described in the same framework as in the case of fluid–fluid interactions, using the Weeks–Chandler–Andersen (WCA) scheme and the Carnahan–Starling (CS) equation for attractive and repulsive parts of the Helmholtz free energy, respectively. Application to nitrogen and argon adsorption isotherms on non-porous silica LiChrospher Si-1000 at their boiling points, recently published by Jaroniec and co-workers, has shown an excellent correlative ability of our approach over the complete range of pressures, which suggests that the surface roughness is mostly the reason for the observed behavior of adsorption isotherms. From the analysis of these data, we found that in the case of nitrogen adsorption short-range interactions between oxygen atoms on the silica surface and quadrupole of nitrogen molecules play an important role. The approach presented in this paper may be further used in quantitative analysis of adsorption and desorption isotherms in cylindrical pores such as MCM-41 and carbon nanotubes.

Adsorption of N-2, Ar, acetone, chloroform and acetone-chloroform mixture on carbon black in the framework of molecular layer structure theory (MLST)

Adsorption of pure nitrogen, argon, acetone, chloroform and acetone-chloroform mixture on graphitized thermal carbon black is considered at sub-critical conditions by means of molecular layer structure theory (MLST). In the present version of the MLST an adsorbed fluid is considered as a sequence of 2D molecular layers, whose Helmholtz free energies are obtained directly from the analysis of experimental adsorption isotherm of pure components. The interaction of the nearest layers is accounted for in the framework of mean field approximation. This approach allows quantitative correlating of experimental nitrogen and argon adsorption isotherm both in the monolayer region and in the range of multi-layer coverage up to 10 molecular layers. In the case of acetone and chloroform the approach also leads to excellent quantitative correlation of adsorption isotherms, while molecular approaches such as the non-local density functional theory (NLDFT) fail to describe those isotherms. We extend our new method to calculate the Helmholtz free energy of an adsorbed mixture using a simple mixing rule, and this allows us to predict mixture adsorption isotherms from pure component adsorption isotherms. The approach, which accounts for the difference in composition in different molecular layers, is tested against the experimental data of acetone-chloroform mixture (non-ideal mixture) adsorption on graphitized thermal carbon black at 50 degrees C. (C) 2005 Elsevier Ltd. All rights reserved.

The effects of parabens on the mechanosensitive channels of E. coli

Parabens are alkyl esters of p-hydroxybenzoic acid used as preservatives in a wide range of food, pharmaceutical, and cosmetic products (Soni et al. Food Chem. Toxicol. 39:513-532, 2001). Despite their common use for over 50 years, their mechanism of action is still unclear. In this study we examined the effects of ethyl and propyl paraben, on gating of the E. coli mechanosensitive channel of large conductance (MscL) reconstituted into azolectin liposomes. We found that propyl and ethyl paraben spontaneously activate MscL. Moreover, the addition of propyl paraben caused an increase in MscL activity and the lowering Of P-1/2, the pressure at which the MscL was opened 50% of the time, the AGO, the free energy required to open the MscL, and the parameter a, which describes the channel sensitivity to pressure. In addition, in silico studies showed that propyl paraben binds to the channel gate of the MscL. The mechanosensitive channel of small conductance was also found to be spontaneously activated by parabens. In summary, our study indicates that one of the previously unidentified mechanisms of action of parabens as antimicrobial agents is via an interaction with the mechanosensitive channels to upset the osmotic gradients in bacteria.

Modeling of adsorption and nucleation in infinite cylindrical pores by two-dimensional density functional theory

In this paper, we present an analysis of argon adsorption in cylindrical pores having amorphous silica structure by means of a nonlocal density functional theory (NLDFT). In the modeling, we account for the radial and longitudinal density distributions, which allow us to consider the interface between the liquidlike and vaporlike fluids separated by a hemispherical meniscus in the canonical ensemble. The Helmholtz free energy of the meniscus was determined as a function of pore diameter. The canonical NLDFT simulations show the details of density rearrangement at the vaporlike and liquidlike spinodal points. The limits of stability of the smallest bridge and the smallest bubble were also determined with the canonical NLDFT. The energy of nucleation as a function of the bulk pressure and the pore diameter was determined with the grand canonical NLDFT using an additional external potential field. It was shown that the experimentally observed reversibility of argon adsorption isotherms at its boiling point up to the pore diameter of 4 nm is possible if the potential barrier of 22kT is overcome due to density fluctuations.

Adsorption in slit pores and pore-size distribution: A molecular layer structure theory

A new approach is developed to analyze the thermodynamic properties of a sub-critical fluid adsorbed in a slit pore of activated carbon. The approach is based on a representation that an adsorbed fluid forms an ordered structure close to a smoothed solid surface. This ordered structure is modelled as a collection of parallel molecular layers. Such a structure allows us to express the Helmholtz free energy of a molecular layer as the sum of the intrinsic Helmholtz free energy specific to that layer and the potential energy of interaction of that layer with all other layers and the solid surface. The intrinsic Helmholtz free energy of a molecular layer is a function (at given temperature) of its two-dimensional density and it can be readily obtained from bulk-phase properties, while the interlayer potential energy interaction is determined by using the 10-4 Lennard-Jones potential. The positions of all layers close to the graphite surface or in a slit pore are considered to correspond to the minimum of the potential energy of the system. This model has led to accurate predictions of nitrogen and argon adsorption on carbon black at their normal boiling points. In the case of adsorption in slit pores, local isotherms are determined from the minimization of the grand potential. The model provides a reasonable description of the 0-1 monolayer transition, phase transition and packing effect. The adsorption of nitrogen at 77.35 K and argon at 87.29 K on activated carbons is analyzed to illustrate the potential of this theory, and the derived pore-size distribution is compared favourably with that obtained by the Density Functional Theory (DFT). The model is less time-consuming than methods such as the DFT and Monte-Carlo simulation, and most importantly it can be readily extended to the adsorption of mixtures and capillary condensation phenomena.

C-terminal charged cluster of MscL, RKKEE, functions as a pH sensor

The RKKEE cluster of charged residues located within the cytoplasmic helix of the bacterial mechanosensitive channel, MscL, is essential for the channel function. The structure of MscL determined by x-ray crystallography and electron paramagnetic resonance spectroscopy has revealed discrepancies toward the C-terminus suggesting that the structure of the C-terminal helical bundle differs depending on the pH of the cytoplasm. In this study we examined the effect of pH as well as charge reversal and residue substitution within the RKKEE cluster on the mechanosensitivity of Escherichia coli MscL reconstituted into liposomes using the patch-clamp technique. Protonation of either positively or negatively charged residues within the cluster, achieved by changing the experimental pH or residue substitution within the RKKEE cluster, significantly increased the free energy of activation for the MscL channel due to an increase in activation pressure. Our data suggest that the orientation of the C-terminal helices relative to the aqueous medium is pH dependent, indicating that the RKKEE cluster functions as a proton sensor by adjusting the channel sensitivity to membrane tension in a pH-dependent fashion. A possible implication of our results for the physiology of bacterial cells is briefly discussed.

Effect of adsorption deformation on thermodynamic characteristics of a fluid in slit pores at sub-critical conditions

Solvation. pressure due to adsorption of fluids in porous materials is the cause of elastic deformation of an adsorbent, which is accessible to direct experimental measurements. Such a deformation contributes to the Helmholtz free energy of the whole adsorbent-adsorbate system due to accumulation of compression or tension energy by the solid. It means that in the general case the solid has to be considered as not solely a source of the external potential field for the fluid confined in the pore volume, but also as thermodynamically nonmert component of the solid-fluid system. We present analysis of nitrogen adsorption isotherms and heat of adsorption in slit graphitic pores accounting for the adsorption deformation by means of nonlocal density functional theory. (c) 2006 Elsevier Ltd. All rights reserved.

Equation of state of a superfluid Fermi gas in the BCS-BEC crossover

We present a theory for a superfluid Fermi gas near the BCS-BEC crossover, including pairing fluctuation contributions to the free energy similar to that considered by Nozieres and Schmitt-Rink for the normal phase. In the strong coupling limit, our theory is able to recover the Bogoliubov theory of a weakly interacting Bose gas with a molecular scattering length very close to the known exact result. We compare our results with recent Quantum Monte Carlo simulations both for the ground state and at finite temperature. Excellent agreement is found for all interaction strengths where simulation results are available.

Sviluppo e Applicazione di Metodologie per l'Ottimizzazione di Sistemi Energetici

La riduzione dei consumi di combustibili fossili e lo sviluppo di tecnologie per il risparmio energetico sono una questione di centrale importanza sia per l’industria che per la ricerca, a causa dei drastici effetti che le emissioni di inquinanti antropogenici stanno avendo sull’ambiente. Mentre un crescente numero di normative e regolamenti vengono emessi per far fronte a questi problemi, la necessità di sviluppare tecnologie a basse emissioni sta guidando la ricerca in numerosi settori industriali. Nonostante la realizzazione di fonti energetiche rinnovabili sia vista come la soluzione più promettente nel lungo periodo, un’efficace e completa integrazione di tali tecnologie risulta ad oggi impraticabile, a causa sia di vincoli tecnici che della vastità della quota di energia prodotta, attualmente soddisfatta da fonti fossili, che le tecnologie alternative dovrebbero andare a coprire. L’ottimizzazione della produzione e della gestione energetica d’altra parte, associata allo sviluppo di tecnologie per la riduzione dei consumi energetici, rappresenta una soluzione adeguata al problema, che può al contempo essere integrata all’interno di orizzonti temporali più brevi. L’obiettivo della presente tesi è quello di investigare, sviluppare ed applicare un insieme di strumenti numerici per ottimizzare la progettazione e la gestione di processi energetici che possa essere usato per ottenere una riduzione dei consumi di combustibile ed un’ottimizzazione dell’efficienza energetica. La metodologia sviluppata si appoggia su un approccio basato sulla modellazione numerica dei sistemi, che sfrutta le capacità predittive, derivanti da una rappresentazione matematica dei processi, per sviluppare delle strategie di ottimizzazione degli stessi, a fronte di condizioni di impiego realistiche. Nello sviluppo di queste procedure, particolare enfasi viene data alla necessità di derivare delle corrette strategie di gestione, che tengano conto delle dinamiche degli impianti analizzati, per poter ottenere le migliori prestazioni durante l’effettiva fase operativa. Durante lo sviluppo della tesi il problema dell’ottimizzazione energetica è stato affrontato in riferimento a tre diverse applicazioni tecnologiche. Nella prima di queste è stato considerato un impianto multi-fonte per la soddisfazione della domanda energetica di un edificio ad uso commerciale. Poiché tale sistema utilizza una serie di molteplici tecnologie per la produzione dell’energia termica ed elettrica richiesta dalle utenze, è necessario identificare la corretta strategia di ripartizione dei carichi, in grado di garantire la massima efficienza energetica dell’impianto. Basandosi su un modello semplificato dell’impianto, il problema è stato risolto applicando un algoritmo di Programmazione Dinamica deterministico, e i risultati ottenuti sono stati comparati con quelli derivanti dall’adozione di una più semplice strategia a regole, provando in tal modo i vantaggi connessi all’adozione di una strategia di controllo ottimale. Nella seconda applicazione è stata investigata la progettazione di una soluzione ibrida per il recupero energetico da uno scavatore idraulico. Poiché diversi layout tecnologici per implementare questa soluzione possono essere concepiti e l’introduzione di componenti aggiuntivi necessita di un corretto dimensionamento, è necessario lo sviluppo di una metodologia che permetta di valutare le massime prestazioni ottenibili da ognuna di tali soluzioni alternative. Il confronto fra i diversi layout è stato perciò condotto sulla base delle prestazioni energetiche del macchinario durante un ciclo di scavo standardizzato, stimate grazie all’ausilio di un dettagliato modello dell’impianto. Poiché l’aggiunta di dispositivi per il recupero energetico introduce gradi di libertà addizionali nel sistema, è stato inoltre necessario determinare la strategia di controllo ottimale dei medesimi, al fine di poter valutare le massime prestazioni ottenibili da ciascun layout. Tale problema è stato di nuovo risolto grazie all’ausilio di un algoritmo di Programmazione Dinamica, che sfrutta un modello semplificato del sistema, ideato per lo scopo. Una volta che le prestazioni ottimali per ogni soluzione progettuale sono state determinate, è stato possibile effettuare un equo confronto fra le diverse alternative. Nella terza ed ultima applicazione è stato analizzato un impianto a ciclo Rankine organico (ORC) per il recupero di cascami termici dai gas di scarico di autovetture. Nonostante gli impianti ORC siano potenzialmente in grado di produrre rilevanti incrementi nel risparmio di combustibile di un veicolo, è necessario per il loro corretto funzionamento lo sviluppo di complesse strategie di controllo, che siano in grado di far fronte alla variabilità della fonte di calore per il processo; inoltre, contemporaneamente alla massimizzazione dei risparmi di combustibile, il sistema deve essere mantenuto in condizioni di funzionamento sicure. Per far fronte al problema, un robusto ed efficace modello dell’impianto è stato realizzato, basandosi sulla Moving Boundary Methodology, per la simulazione delle dinamiche di cambio di fase del fluido organico e la stima delle prestazioni dell’impianto. Tale modello è stato in seguito utilizzato per progettare un controllore predittivo (MPC) in grado di stimare i parametri di controllo ottimali per la gestione del sistema durante il funzionamento transitorio. Per la soluzione del corrispondente problema di ottimizzazione dinamica non lineare, un algoritmo basato sulla Particle Swarm Optimization è stato sviluppato. I risultati ottenuti con l’adozione di tale controllore sono stati confrontati con quelli ottenibili da un classico controllore proporzionale integrale (PI), mostrando nuovamente i vantaggi, da un punto di vista energetico, derivanti dall’adozione di una strategia di controllo ottima.