Molecular Dynamics Simulations of Homogeneous Nucleation

Nucleation is the first step in a phase transition where small nuclei of the new phase start appearing in the metastable old phase, such as the appearance of small liquid clusters in a supersaturated vapor. Nucleation is important in various industrial and natural processes, including atmospheric new particle formation: between 20 % to 80 % of atmospheric particle concentration is due to nucleation. These atmospheric aerosol particles have a significant effect both on climate and human health. Different simulation methods are often applied when studying things that are difficult or even impossible to measure, or when trying to distinguish between the merits of various theoretical approaches. Such simulation methods include, among others, molecular dynamics and Monte Carlo simulations. In this work molecular dynamics simulations of the homogeneous nucleation of Lennard-Jones argon have been performed. Homogeneous means that the nucleation does not occur on a pre-existing surface. The simulations include runs where the starting configuration is a supersaturated vapor and the nucleation event is observed during the simulation (direct simulations), as well as simulations of a cluster in equilibrium with a surrounding vapor (indirect simulations). The latter type are a necessity when the conditions prevent the occurrence of a nucleation event in a reasonable timeframe in the direct simulations. The effect of various temperature control schemes on the nucleation rate (the rate of appearance of clusters that are equally able to grow to macroscopic sizes and to evaporate) was studied and found to be relatively small. The method to extract the nucleation rate was also found to be of minor importance. The cluster sizes from direct and indirect simulations were used in conjunction with the nucleation theorem to calculate formation free energies for the clusters in the indirect simulations. The results agreed with density functional theory, but were higher than values from Monte Carlo simulations. The formation energies were also used to calculate surface tension for the clusters. The sizes of the clusters in the direct and indirect simulations were compared, showing that the direct simulation clusters have more atoms between the liquid-like core of the cluster and the surrounding vapor. Finally, the performance of various nucleation theories in predicting simulated nucleation rates was investigated, and the results among other things highlighted once again the inadequacy of the classical nucleation theory that is commonly employed in nucleation studies.

Large-scale structures in turbulent and reverse-transitional sink flow boundary layers

Aspects of large-scale organized structures in sink flow turbulent and reverse-transitional boundary layers are studied experimentally using hot-wire anemometry. Each of the present sink flow boundary layers is in a state of 'perfect equilibrium' or 'exact self-preservation' in the sense of Townsend (The Structure of Turbulent Shear Flow, 1st and 2nd edns, 1956, 1976, Cambridge University Press) and Rotta (Progr. Aeronaut. Sci., vol. 2, 1962, pp. 1-220) and conforms to the notion of 'pure wall-flow' (Coles, J. Aerosp. Sci., vol. 24, 1957, pp. 495-506), at least for the turbulent cases. It is found that the characteristic inclination angle of the structure undergoes a systematic decrease with the increase in strength of the streamwise favourable pressure gradient. Detectable wall-normal extent of the structure is found to be typically half of the boundary layer thickness. Streamwise extent of the structure shows marked increase as the favourable pressure gradient is made progressively severe. Proposals for the typical eddy forms in sink flow turbulent and reverse-transitional flows are presented, and the possibility of structural self-organization (i.e. individual hairpin vortices forming streamwise coherent hairpin packets) in these flows is also discussed. It is further indicated that these structural ideas may be used to explain, from a structural viewpoint, the phenomenon of soft relaminarization or reverse transition of turbulent boundary layers when subjected to strong streamwise favourable pressure gradients. Taylor's 'frozen turbulence' hypothesis is experimentally shown to be valid for flows in the present study even though large streamwise accelerations are involved, the flow being even reverse transitional in some cases. Possible conditions, which are required to be satisfied for the safe use of Taylor's hypothesis in pressure-gradient-driven flows, are also outlined. Measured convection velocities are found to be fairly close to the local mean velocities (typically 90% or more) suggesting that the structure gets convected downstream almost along with the mean flow.

Optimum Design for External Seismic Stability of Geosynthetic Reinforced Soil Walls: Reliability Based Approach

In this paper, an analytical study considering the effect of uncertainties in the seismic analysis of geosynthetic-reinforced soil (GRS) walls is presented. Using limit equilibrium method and assuming sliding wedge failure mechanism, analysis is conducted to evaluate the external stability of GRS walls when subjected to earthquake loads. Target reliability based approach is used to estimate the probability of failure in three modes of failure, viz., sliding, bearing, and eccentricity failure. The properties of reinforced backfill, retained backfill, foundation soil, and geosynthetic reinforcement are treated as random variables. In addition, the uncertainties associated with horizontal seismic acceleration and surcharge load acting on the wall are considered. The optimum length of reinforcement needed to maintain the stability against three modes of failure by targeting various component and system reliability indices is obtained. Studies have also been made to study the influence of various parameters on the seismic stability in three failure modes. The results are compared with those given by first-order second moment method and Monte Carlo simulation methods. In the illustrative example, external stability of the two walls, Gould and Valencia walls, subjected to Northridge earthquake is reexamined.

Two-Phase Thermodynamic Model for Efficient and Accurate Absolute Entropy of Water from Molecular Dynamics Simulations

Presented here is the two-phase thermodynamic (2PT) model for the calculation of energy and entropy of molecular fluids from the trajectory of molecular dynamics (MD) simulations. In this method, the density of state (DoS) functions (including the normal modes of translation, rotation, and intramolecular vibration motions) are determined from the Fourier transform of the corresponding velocity autocorrelation functions. A fluidicity parameter (f), extracted from the thermodynamic state of the system derived from the same MD, is used to partition the translation and rotation modes into a diffusive, gas-like component (with 3Nf degrees of freedom) and a nondiffusive, solid-like component. The thermodynamic properties, including the absolute value of entropy, are then obtained by applying quantum statistics to the solid component and applying hard sphere/rigid rotor thermodynamics to the gas component. The 2PT method produces exact thermodynamic properties of the system in two limiting states: the nondiffusive solid state (where the fluidicity is zero) and the ideal gas state (where the fluidicity becomes unity). We examine the 2PT entropy for various water models (F3C, SPC, SPC/E, TIP3P, and TIP4P-Ew) at ambient conditions and find good agreement with literature results obtained based on other simulation techniques. We also validate the entropy of water in the liquid and vapor phases along the vapor-liquid equilibrium curve from the triple point to the critical point. We show that this method produces converged liquid phase entropy in tens of picoseconds, making it an efficient means for extracting thermodynamic properties from MD simulations.

Reliability Analysis of Earth Dams

In the present study, results of reliability analyses of four selected rehabilitated earth dam sections, i.e., Chang, Tapar, Rudramata, and Kaswati, under pseudostatic loading conditions, are presented. Using the response surface methodology, in combination with first order reliability method and numerical analysis, the reliability index (beta) values are obtained and results are interpreted in conjunction with conventional factor of safety values. The influence of considering variability in the input soil shear strength parameters, horizontal seismic coefficient (alpha(h)), and location of reservoir full level on the stability assessment of the earth dam sections is discussed in the probabilistic framework. A comparison of results with those obtained from other method of reliability analysis, viz., Monte Carlo simulations combined with limit equilibrium approach, provided a basis for discussing the stability of earth dams in probabilistic terms, and the results of the analysis suggest that the considered earth dam sections are reliable and are expected to perform satisfactorily.

Optimum Design for External Seismic Stability of Geosynthetic Reinforced Soil Walls: Reliability Based Approach

In this paper, an analytical study considering the effect of uncertainties in the seismic analysis of geosynthetic-reinforced soil (GRS) walls is presented. Using limit equilibrium method and assuming sliding wedge failure mechanism, analysis is conducted to evaluate the external stability of GRS walls when subjected to earthquake loads. Target reliability based approach is used to estimate the probability of failure in three modes of failure, viz., sliding, bearing, and eccentricity failure. The properties of reinforced backfill, retained backfill, foundation soil, and geosynthetic reinforcement are treated as random variables. In addition, the uncertainties associated with horizontal seismic acceleration and surcharge load acting on the wall are considered. The optimum length of reinforcement needed to maintain the stability against three modes of failure by targeting various component and system reliability indices is obtained. Studies have also been made to study the influence of various parameters on the seismic stability in three failure modes. The results are compared with those given by first-order second moment method and Monte Carlo simulation methods. In the illustrative example, external stability of the two walls, Gould and Valencia walls, subjected to Northridge earthquake is reexamined.

Studies on interaction of colloidal Ag nanoparticles with Bovine Serum Albumin (BSA)

Biofunctionalization of noble metal nanoparticles like Ag, Au is essential to obtain biocompatibility for specific biomedical applications. Silver nanciparticles are being increasingly used in bio-sensing applications owing to excellent optoelectronic properties. Among the serum albumins, the most abundant proteins in plasma, a wide range of physiological functions of Bovine Serum Albumin (BSA) has made it a model system for biofunctionalization. In absence of adequate prior reports, this study aims to investigate the interaction between silver nanoparticles and BSA. The interaction of BSA [0.05-0.85% concentrations] with Ag nanoparticles [50 ppm concentration] in aqueous dispersion was Studied through UV-vis spectral changes, morphological and surface structural changes. At pH 7, which is More than the isoelectric point of BSA, a decrease in absorbance at plasmon peak of uninteracted nanciparticles (425 mn) was noted till 0.45% BSA, beyond that a blue shift towards 410 urn was observed. The blue shift may be attributed to enhanced electron density on the particle surfaces. Increasing pH to 12 enhanced the blue shift further to 400 rim. The conformational changes in BSA at alkaline pH ranges and consequent hydrophobic interactions also played an important role. The equilibrium adsorption data fitted better to Freundlich isotherm compared to Langmuir Curve. The X-ray diffraction study revealed complete coverage of Ag nanoparticles by BSA. The scanning electron microscopic study of the interacted nanoparticles was also carried Out to decipher morphological changes. This study established that tailoring the concentration of BSA and pH of the interaction it was possible to reduce aggregation of nanoparticles. Biofunctionalized Ag nanoparticles with reduced aggregation will be more amenable towards bio-sensing applications. (C) 2009 Elsevier B.V. All rights reserved.

Structure and electron transport anomalies in InSe-In6Se 7 composite

Rapid solidification of an equiatomic In-Se alloy resulted in the formation of an equilibrium InSe-In6Se7 phase mixture. The InSe phase was found to be polytypic and exhibited the structural variants 2H, 3H, and 4H. The 4H polytype was found to be in considerably higher proportion compared to 2H and 3H types. The In6Se7 phase was found to be hexagonal with a=0.8919 nm and c=1.4273 nm. Both In6Se 7 and the polytypes of InSe could be identified with the space group P61. The conductivity σ variation with temperature was found to be similar to that observed in disordered semiconducting materials. For temperatures >200 K, ln σ decreased linearly with T-1, phonon-assisted carrier excitation. For temperatures <200 K, ln σ decrease followed T-1/3 behavior, representative of variable-range hopping conduction of electrons.

Solution concepts in two-person multicriteria games

In this paper, we propose new solution concepts for multicriteria games and compare them with existing ones. The general setting is that of two-person finite games in normal form (matrix games) with pure and mixed strategy sets for the players. The notions of efficiency (Pareto optimality), security levels, and response strategies have all been used in defining solutions ranging from equilibrium points to Pareto saddle points. Methods for obtaining strategies that yield Pareto security levels to the players or Pareto saddle points to the game, when they exist, are presented. Finally, we study games with more than two qualitative outcomes such as combat games. Using the notion of guaranteed outcomes, we obtain saddle-point solutions in mixed strategies for a number of cases. Examples illustrating the concepts, methods, and solutions are included.

A chemical approach to an understanding of the fast ion conduction in silver iodide-silver oxysalt glasses

A number of AgI based fast ion conducting glasses, with a general formula AgI---Ag2O---MxOy (MxOy=MoO3, SeO3, WO3, V2O5, P2O5, GeO2, B2O3, As2O3, CrO3) have been studied. A chemical approach is made to investigate the origin of fast ion conduction in these glasses. An index known as Image tructural Image npinning Image umber, SUN (S), has been defined for the purpose, based on the unscreened nuclear charge of silver ions and the equilibrium lectronegativities of the halide-oxyanion matrix in these glasses. The variation of the glass transition temperature, Tg, conductivity, σ, and the energy of activation, Ea, with the concentration of AgI are discussed in the light of the structural unpinning number. Conductivities increase uniformly in any given glass series as a smooth function of S and level off at very high values. The entire range of conductivity appears to vary as ln Image , where ln σ0 corresponds roughly to the conductivity of the hypothetical AgI glass and “a” is a constant which could be obtained as the slope in the graph of ln Ea versus S. It is suggested that the increase in the concentration of AgI beyond 75–80 mole% in the glass is not advantageous from the conductivity point of view.

Thermal behaviours of sulphur-selenium mixtures

The thermal behaviours of sulphur, selenium and their mixtures have been studied over the range 40–450Dagger. It has been shown that the polymerization threshold temperature of sulphur,T Ø, decreases with increasing selenium content and follows the equilibrium copolymerization model proposed by Tobolsky and Owen. The formation of octa-atomic species Se8–xSx, where 8 >x > 4, takes place only after sulphur is in the liquid state. The rate of polymerization is enhanced by the addition of increasing amounts of selenium and this is reflected in the higher polymerization peak temperatures. The X-ray powder diffractograms show that all the sulphur-selenium melts belong to the same phase as that of SeS, though the constituent atoms are randomly distributed.

Chemical potentials of oxygen for fayalite-quartz-lron and fayalite-quartz-magnetite equilibria

The oxygen potentials corresponding to fayalite-quartz-iron (FQI) and fayalite-quartz-magnetite (FQM) equilibria have been determined using solid-state galvanic cells: Pt,Fe + Fe2SiO4 + SiO2/(Y2O3)ZrO2/Fe + \r"FeO,\l"Pt and Pt, Fe3O4 + Fe2SiO4 + SiO2/(Y2O3)ZrO2/Ni + NiO, Pt in the temperature ranges 900 to 1400 K and 1080 to 1340 K, respectively. The cells are written such that the right-hand electrodes are positive. Silica used in this study had the quartz structure. The emf of both cells was found to be reversible and to vary linearly with temperature. From the emf, Gibbs energy changes were deduced for the reactions: 0.106Fe (s) + 2Fe0.947O (r.s.) + SiO2 (qz) → Fe2SiO4 (ol) δG‡= -39,140+ 15.59T(± 150) J mol-1 and 3Fe2SiO4 (ol) + O2 (g) → 2Fe3O4 (sp) + 3SiO2 (qz) δG‡ = -471,750 + 160.06 T±} 1100) J mol-1 The “third-law≓ analysis of fayalite-quartz-wustite and fayalite-quartz-magnetite equilibria gives value for δH‡298 as -35.22 (±0.1) and -528.10 (±0.1) kJ mol-1, respectively, independent of temperature. The Gibbs energy of formation of the spinel form of Fe2SiO4 is derived by com-bining the present results on FQI equilibrium with the high-pressure data on olivine to spinel transformation of Fe2SiO4.

Protein sequence design on the basis of topology optimization techniques -Using continuous modeling of discrete amino acid types

The notion of optimization is inherent in protein design. A long linear chain of twenty types of amino acid residues are known to fold to a 3-D conformation that minimizes the combined inter-residue energy interactions. There are two distinct protein design problems, viz. predicting the folded structure from a given sequence of amino acid monomers (folding problem) and determining a sequence for a given folded structure (inverse folding problem). These two problems have much similarity to engineering structural analysis and structural optimization problems respectively. In the folding problem, a protein chain with a given sequence folds to a conformation, called a native state, which has a unique global minimum energy value when compared to all other unfolded conformations. This involves a search in the conformation space. This is somewhat akin to the principle of minimum potential energy that determines the deformed static equilibrium configuration of an elastic structure of given topology, shape, and size that is subjected to certain boundary conditions. In the inverse-folding problem, one has to design a sequence with some objectives (having a specific feature of the folded structure, docking with another protein, etc.) and constraints (sequence being fixed in some portion, a particular composition of amino acid types, etc.) while obtaining a sequence that would fold to the desired conformation satisfying the criteria of folding. This requires a search in the sequence space. This is similar to structural optimization in the design-variable space wherein a certain feature of structural response is optimized subject to some constraints while satisfying the governing static or dynamic equilibrium equations. Based on this similarity, in this work we apply the topology optimization methods to protein design, discuss modeling issues and present some initial results.

A one point shock capturing kinetic scheme for hyperbolic conservation laws

We have presented a new low dissipative kinetic scheme based on a modified Courant Splitting of the molecular velocity through a parameter φ. Conditions for the split fluxes derived based on equilibrium determine φ for a one point shock. It turns out that φ is a function of the Left and Right states to the shock and that these states should satisfy the Rankine-Hugoniot Jump condition. Hence φ is utilized in regions where the gradients are sufficiently high, and is switched to unity in smooth regions. Numerical results confirm a discrete shock structure with a single interior point when the shock is aligned with the grid.

Thermodynamics of GdMnO3 and GdMn2O5 phases determined by the EMF method

Using solid oxide galvanic cells of the type: MnO + Gd2O3 + GdMnO3/O-2/Ni + NiO and Mn3O4 + GdMnO3 + GdMn2O5/O-2/air the equilibrium oxygen pressure for the following reactions :MnO + 1/2Gd(2)O(3) + 1/4O(2) = GdMnO3 1/3Mn(3)O(4) + GdMnO3 + 1/3O(2) = GdMn2O5 was determined in the temperature range from 1073 to 1450 K. From the determined equilibrium oxygen partial pressure the corresponding G i b b s free energy change for these reactions was derived: Delta G(f,GdMnO3)(0) (+/- 425J) 132721(+/ - 2240) +51.91(+/ - 0.81)T Delta G(f,GdMn2O5)(0)(+/- 670J) 121858(+/ - 6176) + 79.52(+/ - 4.83)T From these data, standard G i b b s energies, enthalpies and entropies of formation of GdMnO3 and GdMn2O5 from component oxides and from the elements are derived. Thermodynamic data tables for the two ternary phases are compiled from 298.15 to 1400 K.