147 resultados para MEAN-FIELD SIMULATIONS
em Indian Institute of Science - Bangalore - Índia
Resumo:
Part I (Manjunath et al., 1994, Chem. Engng Sci. 49, 1451-1463) of this paper showed that the random particle numbers and size distributions in precipitation processes in very small drops obtained by stochastic simulation techniques deviate substantially from the predictions of conventional population balance. The foregoing problem is considered in this paper in terms of a mean field approximation obtained by applying a first-order closure to an unclosed set of mean field equations presented in Part I. The mean field approximation consists of two mutually coupled partial differential equations featuring (i) the probability distribution for residual supersaturation and (ii) the mean number density of particles for each size and supersaturation from which all average properties and fluctuations can be calculated. The mean field equations have been solved by finite difference methods for (i) crystallization and (ii) precipitation of a metal hydroxide both occurring in a single drop of specified initial supersaturation. The results for the average number of particles, average residual supersaturation, the average size distribution, and fluctuations about the average values have been compared with those obtained by stochastic simulation techniques and by population balance. This comparison shows that the mean field predictions are substantially superior to those of population balance as judged by the close proximity of results from the former to those from stochastic simulations. The agreement is excellent for broad initial supersaturations at short times but deteriorates progressively at larger times. For steep initial supersaturation distributions, predictions of the mean field theory are not satisfactory thus calling for higher-order approximations. The merit of the mean field approximation over stochastic simulation lies in its potential to reduce expensive computation times involved in simulation. More effective computational techniques could not only enhance this advantage of the mean field approximation but also make it possible to use higher-order approximations eliminating the constraints under which the stochastic dynamics of the process can be predicted accurately.
Resumo:
We present an extensive study of Mott insulator (MI) and superfluid (SF) shells in Bose-Hubbard (BH) models for bosons in optical lattices with harmonic traps. For this we apply the inhomogeneous mean-field theory developed by Sheshadri et al. Phys. Rev. Lett. 75, 4075 (1995)]. Our results for the BH model with one type of spinless bosons agree quantitatively with quantum Monte Carlo simulations. Our approach is numerically less intensive than such simulations, so we are able to perform calculations on experimentally realistic, large three-dimensional systems, explore a wide range of parameter values, and make direct contact with a variety of experimental measurements. We also extend our inhomogeneous mean-field theory to study BH models with harmonic traps and (a) two species of bosons or (b) spin-1 bosons. With two species of bosons, we obtain rich phase diagrams with a variety of SF and MI phases and associated shells when we include a quadratic confining potential. For the spin-1 BH model, we show, in a representative case, that the system can display alternating shells of polar SF and MI phases, and we make interesting predictions for experiments in such systems.
Resumo:
The near-critical behavior of the susceptibility deduced from light-scattering measurements in a ternary liquid mixture of 3-methylpyridine, water, and sodium bromide has been determined. The measurements have been performed in the one-phase region near the lower consolute points of samples with different concentrations of sodium bromide. A crossover from Ising asymptotic behavior to mean-field behavior has been observed. As the concentration of sodium bromide increases, the crossover becomes more pronounced, and the crossover temperature shifts closer to the critical temperature. The data are well described by a model that contains two independent crossover parameters. The crossover of the susceptibility critical exponent γ from its Ising value γ=1.24 to the mean-field value γ=1 is sharp and nonmonotonic. We conclude that there exists an additional length scale in the system due to the presence of the electrolyte which competes with the correlation length of the concentration fluctuations. An analogy with crossover phenomena in polymer solutions and a possible connection with multicritical phenomena is discussed.
Resumo:
A lattice-gas model of multilayer adsorption has been solved in the mean-field approximation by a different numerical method. Earlier workers obtained a single solution for all values of temperature and pressure. In the present work, multiple solutions have been obtained in certain regions of temperature and pressure which give rise to bysteresis in the adsorption isotherm. In addition, we have obtained a parameter which behaves like an order parameter for the transition. The potential-energy function shows a double minimum in the region of bysteresis and a single maximum elsewhere.
Resumo:
Several channels provided by many-body couplings — both fermion-fermion and fermion-boson — for the evolution of the chemisorption system are discussed. This provides an opportunity of a systematic study of the effects of correlations reflected through the intricate pole structure of the absorbate Green functions. The results of Newns, Anda and others in the context of chemisorption are generalized.
Resumo:
We consider the growth of an isolated precipitate when the matrix diffusivity depends on the composition. We have simulated precipitate growth using the Cahn-Hilliard model, and find good agreement between our results and those from a sharp interface theory for systems with and without a dilatational misfit. With misfit, we report (and rationalize) an interesting difference between systems with a constant diffusivity and those with a variable diffusivity in the matrix.
Resumo:
We generalize the mean-field theory for the spinless Bose-Hubbard model to account for the different types of superfluid phases that can arise in the spin-1 case. In particular, our mean-field theory can distinguish polar and ferromagnetic superfluids, Mott insulator, that arise at integer fillings at zero temperature, and normal Bose liquids into which the Mott insulators evolve at finite temperatures. We find, in contrast to the spinless case, that several of the superfluid-Mott insulator transitions are of first order at finite temperatures. Our systematic study yields rich phase diagrams that include first-order and second-order transitions and a variety of tricritical points. We discuss the possibility of realizing such phase diagrams in experimental systems.
Resumo:
The critical behavior of osmotic susceptibility in an aqueous electrolyte mixture 1-propanol (1P)+water (W)+potassium chloride is reported. This mixture exhibits re-entrant phase transitions and has a nearly parabolic critical line with its apex representing a double critical point (DCP). The behavior of the susceptibility exponent is deduced from static light-scattering measurements, on approaching the lower critical solution temperatures (TL’s) along different experimental paths (by varying t) in the one-phase region. The light-scattering data analysis substantiates the existence of a nonmonotonic crossover behavior of the susceptibility exponent in this mixture. For the TL far away from the DCP, the effective susceptibility exponent γeff as a function of t displays a nonmonotonic crossover from its single limit three-dimensional (3D)-Ising value ( ∼ 1.24) toward its mean-field value with increase in t. While for that closest to the DCP, γeff displays a sharp, nonmonotonic crossover from its nearly doubled 3D-Ising value toward its nearly doubled mean-field value with increase in t. The renormalized Ising regime extends over a relatively larger t range for the TL closest to the DCP, and a trend toward shrinkage in the renormalized Ising regime is observed as TL shifts away from the DCP. Nevertheless, the crossover to the mean-field limit extends well beyond t>10−2 for the TL’s studied. The observed crossover behavior is attributed to the presence of strong ion-induced clustering in this mixture, as revealed by various structure probing techniques. As far as the critical behavior in complex or associating mixtures with special critical points (like the DCP) is concerned, our results indicate that the influence of the DCP on the critical behavior must be taken into account not only on the renormalization of the critical exponent but also on the range of the Ising regime, which can shrink with decrease in the influence of the DCP and with the extent of structuring in the system. The utility of the field variable tUL in analyzing re-entrant phase transitions is demonstrated. The effective susceptibility exponent as a function of tUL displays a nonmonotonic crossover from its asymptotic 3D-Ising value toward a value slightly lower than its nonasymptotic mean-field value of 1. This behavior in the nonasymptotic, high tUL region is interpreted in terms of the possibility of a nonmonotonic crossover to the mean-field value from lower values, as foreseen earlier in micellar systems.
Resumo:
We report novel results obtained for the Hubbard and t-J models by various mean-field approximations.
Resumo:
We obtain metal-insulator phase diagrams at half-filling for the five-band extended Hubbard model of the square-planar CuO2 lattice treated within a Hartree-Fock mean-field approximation, allowing for spiral spin-density waves. We indicate the existence of an insulating phase (covalent insulator) characterized by strong covalency effects, not identified in the earlier Zaanen-Sawatzky-Allen phase diagram. While the insulating phase is always antiferromagnetic, we also obtain an antiferromagnetic metallic phase for a certain range of interaction parameters. Performing a nonperturbative calculation of J(eff), the in-plane antiferromagnetic interaction is presented as a function of the parameters in the model. We also calculate the band gap and magnetic moments at various sites and discuss critically the contrasting interpretation of the electronic structure of high-T(c) materials arising from photoemission and neutron-scattering experiments.
Resumo:
We present analytic results to show that the Schwinger-boson hole-fermion mean-field state exhibits non-Fermi liquid behavior due to spin-charge separation. The physical electron Green's function consists of three additive components. (a) A Fermi-liquid component associated with the bose condensate. (b) A non-Fermi liquid component which has a logarithmic peak and a long tail that gives rise to a linear density of states that is symmetric about the Fermi level and a momentum distribution function with a logarithmic discontinuity at the Fermi surface. (c) A second non-Fermi liquid component associated with the thermal bosons which leads to a constant density of states. It is shown that zero-point fluctuations associated with the spin-degrees of freedom are responsible for the logarithmic instabilities and the restoration of particle-hole symmetry close to the Fermi surface.
Resumo:
Using a lattice model for adsorption in microporous materials, pure component adsorption isotherms are obtained within a mean field approximation for methane at 300 K and xenon at 300 and 360 K in zeolite NaA. It is argued that the increased repulsive adsorbate-adsorbate interactions at high coverages must play an important role in determining the adsorption behavior. Therefore, this feature is incorporated through a "coverage-dependent interaction'' model, which introduces a free, adjustable parameter. Another important feature, the site volume reduction, has been treated in two ways: a van der Waal model and a 1D hard-rod theory [van Tassel et al., AIChE J. 40, 925 (1994)]; we have also generalized the latter to include all possible adsorbate overlap scenarios. In particular, the 1D hard-rod model, with our coverage-dependent interaction model, is shown to be in best quantitative agreement with the previous grand canonical Monte Carlo isotherms. The expressions for the isosteric heats of adsorption indicate that attractive and repulsive adsorbate-adsorbate interactions increase and decrease the heats of adsorption, respectively. It is concluded that within the mean field approximation, our simple model for repulsive interactions and the 1D hard-rod model for site volume reduction are able to capture most of the important features of adsorption in confined regions. (C) 1999 American Institute of Physics. [S0021-9606(99)70515-5].
Resumo:
We study phase transitions in the colossal-magnetoresistive manganites by using a mean-field theory both at zero and non-zero temperatures. Our Hamiltonian includes double-exchange, superexchange, and Hubbard terms with on-site and nearest-neighbour Coulomb interaction, with the parameters estimated from earlier density-functional calculations. The phase diagrams show magnetic and charge-ordered (or charge-disordered) phases as a result of the competition between the double-exchange, superexchange, and Hubbard terms, the relative effects of which are sensitively dependent on parameters such as doping, bandwidth, and temperature. In accord with the experimental observations, several important features are reproduced from our model, namely, (i) a phase transition from an insulating, charge-ordered antiferromagnetic to a metallic, charge-disordered ferromagnetic state near dopant concentration x = 1/2, (ii) the reduction of the transition temperature TAF-->F by the application of a magnetic field, (iii) melting of the charge order by a magnetic field, and (iv) phase coexistence for certain values of temperature and doping. An important feature, not reproduced in our model, is the antiferromagnetism in the electron-doped systems, e.g., La1-xCaxMnO3 over the entire range of 0.5 less than or equal to x less than or equal to 1, and we suggest that a multi-band model which includes the unoccupied t(2g) orbitals might be an important ingredient for describing this feature.
Resumo:
The statistical thermodynamics of adsorption in caged zeolites is developed by treating the zeolite as an ensemble of M identical cages or subsystems. Within each cage adsorption is assumed to occur onto a lattice of n identical sites. Expressions for the average occupancy per cage are obtained by minimizing the Helmholtz free energy in the canonical ensemble subject to the constraints of constant M and constant number of adsorbates N. Adsorbate-adsorbate interactions in the Brag-Williams or mean field approximation are treated in two ways. The local mean field approximation (LMFA) is based on the local cage occupancy and the global mean field approximation (GMFA) is based on the average coverage of the ensemble. The GMFA is shown to be equivalent in formulation to treating the zeolite as a collection of interacting single site subsystems. In contrast, the treatment in the LMFA retains the description of the zeolite as an ensemble of identical cages, whose thermodynamic properties are conveniently derived in the grand canonical ensemble. For a z coordinated lattice within the zeolite cage, with epsilon(aa) as the adsorbate-adsorbate interaction parameter, the comparisons for different values of epsilon(aa)(*)=epsilon(aa)z/2kT, and number of sites per cage, n, illustrate that for -1
Resumo:
We develop a model of the solar dynamo in which, on the one hand, we follow the Babcock-Leighton approach to include surface processes, such as the production of poloidal field from the decay of active regions, and, on the other hand, we attempt to develop a mean field theory that can be studied in quantitative detail. One of the main challenges in developing such models is to treat the buoyant rise of the toroidal field and the production of poloidal field from it near the surface. A previous paper by Choudhuri, Schüssler, & Dikpati in 1995 did not incorporate buoyancy. We extend this model by two contrasting methods. In one method, we incorporate the generation of the poloidal field near the solar surface by Durney's procedure of double-ring eruption. In the second method, the poloidal field generation is treated by a positive α-effect concentrated near the solar surface coupled with an algorithm for handling buoyancy. The two methods are found to give qualitatively similar results.