Dynamical cluster approximation: Nonlocal dynamics of correlated electron systems

We recently introduced the dynamical cluster approximation (DCA), a technique that includes short-ranged dynamical correlations in addition to the local dynamics of the dynamical mean-field approximation while preserving causality. The technique is based on an iterative self-consistency scheme on a finite-size periodic cluster. The dynamical mean-field approximation (exact result) is obtained by taking the cluster to a single site (the thermodynamic limit). Here, we provide details of our method, explicitly show that it is causal, systematic, Phi derivable, and that it becomes conserving as the cluster size increases. We demonstrate the DCA by applying it to a quantum Monte Carlo and exact enumeration study of the two-dimensional Falicov-Kimball model. The resulting spectral functions preserve causality, and the spectra and the charge-density-wave transition temperature converge quickly and systematically to the thermodynamic limit as the cluster size increases.

Polarization Characteristics of Porous Electrode Systems with Adsorbed Intermediates Participating in the Electrode Reaction

Current-potential relationships are derived for porous electrode systems following a homogeneous model and whenadsorbed intermediates participate in the electrode reaction. Limiting Tafel slopes were deduced and compared with thecorresponding behavior on planar electrode systems. The theoretical results showed doubling of Tafel slopes when theslow-step is a charge-transfer reaction and a nonlogarithmic current-voltage behavior when the slow-step is a chemical reaction.Comparison of the experimental results with theory for the case of oxygen reduction on carbon surfaces in alkalinemedia indicates that a slow chemical reaction following the initial charge-transfer reaction to be the likely rate-controllingstep. Theoretical relationships are utilized to determine the exchange current density and the surface coverage by the adsorbedintermediates during the course of oxygen reduction from alkaline solutions on "carbon." Tafel slope measurementson planar and porous electrodes for the same reaction are suggested as one of the diagnostic criteria for elucidatingthe mechanistic pathways of electrochemical reactions.

Simulation of magnetovolume effects in Fe65Ni35 nanoparticles

We compare magnetovolume effects in bulk and nanoparticles by performing Monte Carlo simulations of a spin-analogous model with coupled spatial and magnetic degrees of freedom and chemical disorder. We find that correlations between surface and bulk atoms lead with decreasing particle size to a substantial modification of the magnetic and elastic behavior at low temperatures.

Specification-driven approach for protocol design of distributed computing systems

An important issue in the design of a distributed computing system (DCS) is the development of a suitable protocol. This paper presents an effort to systematize the protocol design procedure for a DCS. Protocol design and development can be divided into six phases: specification of the DCS, specification of protocol requirements, protocol design, specification and validation of the designed protocol, performance evaluation, and hardware/software implementation. This paper describes techniques for the second and third phases, while the first phase has been considered by the authors in their earlier work. Matrix and set theoretic based approaches are used for specification of a DCS and for specification of the protocol requirements. These two formal specification techniques form the basis of the development of a simple and straightforward procedure for the design of the protocol. The applicability of the above design procedure has been illustrated by considering an example of a computing system encountered on board a spacecraft. A Petri-net based approach has been adopted to model the protocol. The methodology developed in this paper can be used in other DCS applications.

Direct numerical simulation of autoignition in a non-premixed, turbulent medium

In this paper, direct numerical simulation of autoignition in an initially non-premixed medium under isotropic, homogeneous, and decaying turbulence is presented. The pressure-based method developed herein is a spectral implementation of the sequential steps followed in the predictor-corrector type of algorithms; it includes the effects of density fluctuations caused by spatial inhomogeneities ill temperature and species. The velocity and pressure field are solved in the spectral space while the scalars and density field are solved in the physical space. The presented results reveal that the autoignition spots originate and evolve at locations where (1) the composition corresponds to a small range around a specific mixture fraction, and (2) the conditional scaler dissipation rate is low. A careful examination of the data obtained indicates that the autoignition spots originate in the vortex cores, and the hot gases travel outward as combustion progresses. Hence, the applicability of the transient laminar flamelet model for this problem is questioned. The dependence of autoignition characteristics on parameters such as (1) die initial eddy-turnover time and (2) the initial ratio of length scale of scalars to that of velocities are investigated. Certain implications of new results on the conditional moment closure modeling are discussed.

Computer simulation of protein—carbohydrate complexes: application to arabinose-binding protein and pea lectin

The CCEM method (Contact Criteria and Energy Minimisation) has been developed and applied to study protein-carbohydrate interactions. The method uses available X-ray data even on the native protein at low resolution (above 2.4 Å) to generate realistic models of a variety of proteins with various ligands.The two examples discussed in this paper are arabinose-binding protein (ABP) and pea lectin. The X-ray crystal structure data reported on ABP-β-l-arabinose complex at 2.8, 2.4 and 1.7 Å resolution differ drastically in predicting the nature of the interactions between the protein and ligand. It is shown that, using the data at 2.4 Å resolution, the CCEM method generates complexes which are as good as the higher (1.7 Å) resolution data. The CCEM method predicts some of the important hydrogen bonds between the ligand and the protein which are missing in the interpretation of the X-ray data at 2.4 Å resolution. The theoretically predicted hydrogen bonds are in good agreement with those reported at 1.7 Å resolution. Pea lectin has been solved only in the native form at 3 Å resolution. Application of the CCEM method also enables us to generate complexes of pea lectin with methyl-α-d-glucopyranoside and methyl-2,3-dimethyl-α-d-glucopyranoside which explain well the available experimental data in solution.

Trends in the stability of ternary oxides: Systems M-Pb-O (M = Ca, Sr, Ba)

Recent experimental investigations of phase equilibria and thermodynamic properties of the systems M-Pb-O, where M = Ca, Sr or Ba, indicate a regular increase in thermodynamic stability of ternary oxides, MPbO3 and M2PbO4, with increasing basicity of the oxide of the alkaline-earth metal. Number of stable interoxide compounds at 1100 K in the systems M-Pb-O (M = Mg, Ca, Sr, Ba) increases in unit increments from Mg to Ba. In this paper, experimentally determined standard Gibbs energies of formation of M2PbO4 (M = Ca, Sr, Ba) and MPbO3 (M = Sr, Ba) from their component binary monoxides and oxygen gas are combined with an estimated value for CaPbO3 to delineate systematic trends in thermodynamic stability of the ternary oxides. The trends are interpreted using concepts of tolerance factor and acid-base interactions. All the ternary oxides in these systems contain lead in the tetravalent state. The small Pb4+ ions polarize the surrounding oxygen ions and cause the formation of oxyanions which are acidic in character. Hence, the higher oxidation state of lead is stabilized in the presence of basic oxides of alkaline-earth group. A schematic subsolidus temperature-composition phase diagram is presented for the system BaO-PbO-O-2 to illustrate the change in oxidation states in binary and ternary oxides with temperature.

Sequential design of decentralized control systems

Our main result is a new sequential method for the design of decentralized control systems. Controller synthesis is conducted on a loop-by-loop basis, and at each step the designer obtains an explicit characterization of the class C of all compensators for the loop being closed that results in closed-loop system poles being in a specified closed region D of the s-plane, instead of merely stabilizing the closed-loop system. Since one of the primary goals of control system design is to satisfy basic performance requirements that are often directly related to closed-loop pole location (bandwidth, percentage overshoot, rise time, settling time), this approach immediately allows the designer to focus on other concerns such as robustness and sensitivity. By considering only compensators from class C and seeking the optimum member of that set with respect to sensitivity or robustness, the designer has a clearly-defined limited optimization problem to solve without concern for loss of performance. A solution to the decentralized tracking problem is also provided. This design approach has the attractive features of expandability, the use of only 'local models' for controller synthesis, and fault tolerance with respect to certain types of failure.

Numerical simulation of solidification of liquid aluminum alloy flowing on cooling slope

Preparation of semisolid slurry using a cooling slope is increasingly becoming popular, primarily because of the simplicity in design and ease control of the process. In this process, liquid alloy is poured down an inclined surface which is cooled from underneath. The cooling enables partial solidification and the incline provides the necessary shear for producing semisolid slurry. However, the final microstructure of the ingot depends on several process parameters such as cooling rate, incline angle of the cooling slope, length of the slope and initial melt superheat. In this work, a CFD model using volume of fluid (VOF) method for simulating flow along the cooling slope was presented. Equations for conservation of mass, momentum, energy and species were solved to predict hydrodynamic and thermal behavior, in addition to predicting solid fraction distribution and macrosegregation. Solidification was modeled using an enthalpy approach and a volume averaged technique for the different phases. The mushy region was modeled as a multi-layered porous medium consisting of fixed columnar dendrites and mobile equiaxed/fragmented grains. The alloy chosen for the study was aluminum alloy A356, for which adequate experimental data were available in the literature. The effects of two key process parameters, namely the slope angle and the pouring temperature, on temperature distribution, velocity distribution and macrosegregation were also studied.