Inherent structures of phase-separating binary mixtures: Nucleation, spinodal decomposition, and pattern formation

An energy landscape view of phase separation and nonideality in binary mixtures is developed by exploring their potential energy landscape (PEL) as functions of temperature and composition. We employ molecular dynamics simulations to study a model that promotes structure breaking in the solute-solvent parent binary liquid, at low temperatures. The PEL of the system captures the potential energy distribution of the inherent structures (IS) of the system and is obtained by removing the kinetic energy (including that of intermolecular vibrations). The broader distribution of the inherent structure energy for structure breaking liquid than that of the structure making liquid demonstrates the larger role of entropy in stabilizing the parent liquid of the structure breaking type of binary mixtures. At high temperature, although the parent structure of the structure breaking binary mixture is homogenous, the corresponding inherent structure is found to be always phase separated, with a density pattern that exhibits marked correlation with the energy of its inherent structure. Over a broad range of intermediate inherent structure energy, bicontinuous phase separation prevails with interpenetrating stripes as signatures of spinodal decomposition. At low inherent structure energy, the structure is largely phase separated with one interface where as at high inherent structure energy we find nucleation type growth. Interestingly, at low temperature, the average inherent structure energy (< EIS >) exhibits a drop with temperature which signals the onset of crystallization in one of the phases while the other remains in the liquid state. The nonideal composition dependence of viscosity is anticorrelated with average inherent structure energy.

High swirl-inducing piston bowls in small diesel engines for emission reduction

Detailed three-dimensional CFD simulations involving flow and combustion chemistry are used to study the effect of swirl induced by re-entrant piston bowl geometries on pollutant emissions from a single-cylinder diesel engine. The baseline engine configuration consists of a hemispherical piston bowl and an injector with finite sac volume. The first iteration involved using a torroidal, slightly re-entrant bowl geometry, and a sac-less injector. Pollutant emission measurements indicated a reduction in emissions with this modification. Simulations on both configurations were then conducted to understand the effect of the changes. The simulation results indicate that the selected piston bowl geometry could actually be reducing the in-cylinder swirl and turbulence and the emission reduction may be entirely due to the introduction of the sac-less injector. In-cylinder air motion was then studied in a number of combustion chamber geometries, and a geometry which produced the highest in-cylinder swirl and Turbulence Kinetic Energy (TKE) around the compression top dead centre (TDC) was identified. The optimal nature of this re-entrant piston bowl geometry is confirmed by detailed combustion simulations and emission predictions. (C) 2010 Elsevier Ltd. All rights reserved.

Optical studies of SN 2009jf: a Type Ib supernova with an extremely slow decline and aspherical signature

Optical UBVRI photometry and medium-resolution spectroscopy of the Type Ib supernova SN 2009jf, during the period from similar to -15 to +250 d, with respect to the B maximum are reported. The light curves are broad, with an extremely slow decline. The early post-maximum decline rate in the V band is similar to SN 2008D; however, the late-phase decline rate is slower than other Type Ib supernovae studied. With an absolute magnitude of M-V = -17.96 +/- 0.19 at peak, SN 2009jf is a normally bright supernova. The peak bolometric luminosity and the energy deposition rate via the 56Ni -> 56Co chain indicate that similar to 0.17+0.03(-0.03) M-circle dot of 56Ni was ejected during the explosion. The He i 5876 A line is clearly identified in the first spectrum of day similar to -15, at a velocity of similar to 16 000 km s-1. The O i] 6300-6364 A line seen in the nebular spectrum has a multipeaked and asymmetric emission profile, with the blue peak being stronger. The estimated flux in this line implies that greater than or similar to 1.34 M-circle dot oxygen was ejected. The slow evolution of the light curves of SN 2009jf indicates the presence of a massive ejecta. The high expansion velocity in the early phase and broader emission lines during the nebular phase suggest it to be an explosion with a large kinetic energy. A simple qualitative estimate leads to the ejecta mass of M-ej = 4-9 M-circle dot and kinetic energy E-K = 3-8 x 1051 erg. The ejected mass estimate is indicative of an initial main-sequence mass of greater than or similar to 20-25 M-circle dot.

Optical and structural properties of reactive ion beam sputter deposited CeO2 films

Optical and structural properties of reactive ion beam sputter deposited CeO2 films as a function of oxygen partial pressures (P-O2) and substrate temperatures (T-s) have been investigated. The films deposited at ambient temperature with P-O2 of 0.01 Pa have shown a refractive index of 2.36 which increased to 2.44 at 400 degrees C. Refractive index and extinction coefficient are sensitive up to a T-s of similar to 200 degrees C. Raman spectroscopy and X-ray diffraction (XRD) have been used to characterise the structural properties. A preferential orientation of (220) was observed up to a T-s of 200 degrees C and it changed to (200) at 400 degrees C: and above. Raman line broadening, peak shift and XRD broadening indicate the formation of nanocrystalline phase for the films deposited up to a substrate temperature of 300 degrees C. However, crystallinity of the films were better for T-s values above 300 degrees C. In general both optical and structural properties were unusual compared to the films deposited by conventional electron beam evaporation, but were similar in some aspects to those deposited by ion-assisted deposition. Apart from thermal effects, this behavior is also attributed to the bombardment of backscattered ions/neutrals on the growing film as well as the higher kinetic energy of the condensing species, together resulting in increased packing density. (C) 1997 Elsevier Science S.A.

Superconductivity from repulsion: a variational view

This letter presents a new class of variational wavefunctions for Fermi systems in any dimension. These wavefunctions introduce correlations between Cooper pairs in different momentum states and the relevant correlations can be computed analytically. At half filling we have a ground state with critical superconducting correlations, that causes negligible increase of the kinetic energy. We find large enhancements in a Cooper-pair correlation function caused purely by the interplay between the uncertainty principle, repulsion and the proximity of half filling. This is surprising since there is no accompanying signature in usual charge and spin response functions, and typifies a novel kind of many-body cooperative behaviour.

Simulation and experimental study on compositional evolution of Li-Co in LiCoO(2) thin films during sputter deposition

The compositional evolution in sputter deposited LiCoO(2) thin films is influenced by process parameters involved during deposition. The electrochemical performance of these films strongly depends on their microstructure, preferential orientation and stoichiometry. The transport process of sputtered Li and Co atoms from the LiCoO(2) target to the substrate, through Ar plasma in a planar magnetron configuration, was investigated based on the Monte Carlo technique. The effect of sputtering gas pressure and the substrate-target distance (d(st)) on Li/Co ratio, as well as, energy and angular distribution of sputtered atoms on the substrate were examined. Stable Li/Co ratios have been obtained at 5 Pa pressure and d(st) in the range 5 11 cm. The kinetic energy and incident angular distribution of Li and Co atoms reaching the substrate have been found to be dependent on sputtering pressure. Simulations were extended to predict compositional variations in films prepared at various process conditions. These results were compared with the composition of films determined experimentally using x-ray photoelectron spectroscopy (XPS). Li/Co ratio calculated using XPS was in moderate agreement with that of the simulated value. The measured film thickness followed the same trend as predicted by simulation. These studies are shown to be useful in understanding the complexities in multicomponent sputtering. (C) 2011 American Institute of Physics. doi:10.1063/1.3597829]

Size and temperature dependent stability and phase transformation in single-crystal zirconium nanowire

A novel size dependent FCC (face-centered-cubic) -> HCP (hexagonally-closed-pack) phase transformation and stability of an initial FCC zirconium nanowire are studied. FCC zirconium nanowires with cross-sectional dimensions < 20 are found unstable in nature, and they undergo a FCC -> HCP phase transformation, which is driven by tensile surface stress induced high internal compressive stresses. FCC nanowire with cross-sectional dimensions > 20 , in which surface stresses are not enough to drive the phase transformation, show meta-stability. In such a case, an external kinetic energy in the form of thermal heating is required to overcome the energy barrier and achieve FCC -> HCP phase transformation. The FCC-HCP transition pathway is also studied using Nudged Elastic Band (NEB) method, to further confirm the size dependent stability/metastability of Zr nanowires. We also show size dependent critical temperature, which is required for complete phase transformation of a metastable-FCC nanowire.

Simulation and experimental study on compositional evolution of Li-Co in LiCoO2 thin films during sputter deposition

The compositional evolution in sputter deposited LiCoO2 thin films is influenced by process parameters involved during deposition. The electrochemical performance of these films strongly depends on their microstructure, preferential orientation and stoichiometry. The transport process of sputtered Li and Co atoms from the LiCoO2 target to the substrate, through Ar plasma in a planar magnetron configuration, was investigated based on the Monte Carlo technique. The effect of sputtering gas pressure and the substrate-target distance (dst) on Li/Co ratio, as well as, energy and angular distribution of sputtered atoms on the substrate were examined. Stable Li/Co ratios have been obtained at 5 Pa pressure and dst in the range 5−11 cm. The kinetic energy and incident angular distribution of Li and Co atoms reaching the substrate have been found to be dependent on sputtering pressure. Simulations were extended to predict compositional variations in films prepared at various process conditions. These results were compared with the composition of films determined experimentally using x-ray photoelectron spectroscopy (XPS). Li/Co ratio calculated using XPS was in moderate agreement with that of the simulated value. The measured film thickness followed the same trend as predicted by simulation. These studies are shown to be useful in understanding the complexities in multicomponent sputtering.

Phase-plane techniques for the solution of transient-stability problems

The paper presents a graphical-numerical method for determining the transient stability limits of a two-machine system under the usual assumptions of constant input, no damping and constant voltage behind transient reactance. The method presented is based on the phase-plane criterion,1, 2 in contrast to the usual step-by-step and equal-area methods. For the transient stability limit of a two-machine system, under the assumptions stated, the sum of the kinetic energy and the potential energy, at the instant of fault clearing, should just be equal to the maximum value of the potential energy which the machines can accommodate with the fault cleared. The assumption of constant voltage behind transient reactance is then discarded in favour of the more accurate assumption of constant field flux linkages. Finally, the method is extended to include the effect of field decrement and damping. A number of examples corresponding to each case are worked out, and the results obtained by the proposed method are compared with those obtained by the usual methods.

On the athermal nature of the beta to omega transformation

One characteristic feature of the athermal beta -> omega transformation is the short time scale of the transformation. So far, no clear understanding of this issue exists. Here we construct a model that includes contributions from a Landau sixth-order free energy density, kinetic energy due to displacement, and the Rayleigh dissipation function to account for the dissipation arising from the rapid movement of the parent product interface during rapid nucleation. We also include the contribution from omega-like fluctuations to local stress. The model shows that the transformation is complete on a time scale comparable to the velocity of sound. The estimated nucleation rate is several orders higher than that for diffusion-controlled transformations. The model predicts that the athermal omega phase is limited to a certain range of alloying composition. The estimated nucleation rate and the size of ``isothermal'' particles beyond 17% Nb are also consistent with experimental results. The model provides an explanation for the reprecipitation process of the omega particles in the ``cleared'' channels formed during deformation of omega-forming alloys. The model also predicts that acoustic emission should be detectable during the formation of the athermal phase. (C) 2011 Acta Materialia Inc. Published by Elsevier Ltd. All rights reserved.

The Nature of Bond Critical Points in Dinuclear Copper(I) Complexes

Closed-shell contacts between two copper(I) ions are expected to be repulsive. However, such contacts are quite frequent and are well documented. Crystallographic characterization of such contacts in unsupported and bridged multinuclear copper(I) complexes has repeatedly invited debates on the existence of cuprophilicity. Recent developments in the application of Baders theory of atoms-in-molecules (AIM) to systems in which weak hydrogen bonds are involved suggests that the copper(I)copper(I) contacts would benefit from a similar analysis. Thus the nature of electron-density distributions in copper(I) dimers that are unsupported, and those that are bridged, have been examined. A comparison of complexes that are dimers of symmetrical monomers and those that are dimers of two copper(I) monomers with different coordination spheres has also been made. AIM analysis shows that a bond critical point (BCP) between two Cu atoms is present in most cases. The nature of the BCP in terms of the electron density, ?, and its Laplacian is quite similar to the nature of critical points observed in hydrogen bonds in the same systems. The ? is inversely correlated to Cu?Cu distance. It is higher in asymmetrical systems than what is observed in corresponding symmetrical systems. By examining the ratio of the local electron potential-energy density (Vc) to the kinetic energy density (Gc), |Vc|/Gc at the critical point suggests that these interactions are not perfectly ionic but have some shared nature. Thus an analysis of critical points by using AIM theory points to the presence of an attractive metallophilic interaction similar to other well-documented weak interactions like hydrogen bonding.

Segregation into Chiral Enantiomeric Conformations of an Achiral Molecule by Concomitant Polymorphism

An analogue of the green fluorescent protein (GFP) luminophore crystallizes from a methanol solution impregnated with dichloromethane, into a pair of chiral crystals. Thermal analysis, fluorescence emission studies, and crystal packing analysis show that the two crystals are different materials. The two polymorphs arise from the rotation of a monosubstituted benzene ring about a C-N bond which results in the formation of two strong bifurcated C-H center dot center dot center dot O intermolecular bonds to oxygen O(6). The color difference has been ascribed to a difference in the packing of the two crystal forms. Theoretical studies supported by low temperature NMR show low kinetic energy barriers (similar to 10 kJ mol(-1)) separating the asymmetric units of the two crystal structures, suggesting that the driving force for the polymorphism could be the result of packing of two different asymmetric units.

Stabilized SPH-based simulations of impact dynamics using acceleration-corrected artificial viscosity

Artificial viscosity in SPH-based computations of impact dynamics is a numerical artifice that helps stabilize spurious oscillations near the shock fronts and requires certain user-defined parameters. Improper choice of these parameters may lead to spurious entropy generation within the discretized system and make it over-dissipative. This is of particular concern in impact mechanics problems wherein the transient structural response may depend sensitively on the transfer of momentum and kinetic energy due to impact. In order to address this difficulty, an acceleration correction algorithm was proposed in Shaw and Reid (''Heuristic acceleration correction algorithm for use in SPH computations in impact mechanics'', Comput. Methods Appl. Mech. Engrg., 198, 3962-3974) and further rationalized in Shaw et al. (An Optimally Corrected Form of Acceleration Correction Algorithm within SPH-based Simulations of Solid Mechanics, submitted to Comput. Methods Appl. Mech. Engrg). It was shown that the acceleration correction algorithm removes spurious high frequency oscillations in the computed response whilst retaining the stabilizing characteristics of the artificial viscosity in the presence of shocks and layers with sharp gradients. In this paper, we aim at gathering further insights into the acceleration correction algorithm by further exploring its application to problems related to impact dynamics. The numerical evidence in this work thus establishes that, together with the acceleration correction algorithm, SPH can be used as an accurate and efficient tool in dynamic, inelastic structural mechanics. (C) 2011 Elsevier Ltd. All rights reserved.

Quasiuniversal Transient Behavior of a Nonequilibrium Mott Insulator Driven by an Electric Field

We use a self-consistent strong-coupling expansion for the self-energy (perturbation theory in the hopping) to describe the nonequilibrium dynamics of strongly correlated lattice fermions. We study the three-dimensional homogeneous Fermi-Hubbard model driven by an external electric field showing that the damping of the ensuing Bloch oscillations depends on the direction of the field and that for a broad range of field strengths a long-lived transient prethermalized state emerges. This long-lived transient regime implies that thermal equilibrium may be out of reach of the time scales accessible in present cold atom experiments but shows that an interesting new quasiuniversal transient state exists in nonequilibrium governed by a thermalized kinetic energy but not a thermalized potential energy. In addition, when the field strength is equal in magnitude to the interaction between atoms, the system undergoes a rapid thermalization, characterized by a different quasiuniversal behavior of the current and spectral function for different values of the hopping. DOI: 10.1103/PhysRevLett.109.260402

Dynamics of the Indian-Ocean oxygen minimum zones

In the Indian Ocean, mid-depth oxygen minimum zones (OMZs) occur in the Arabian Sea and the Bay of Bengal. The lower part of the Arabian-Sea OMZ (ASOMZ; below 400 m) intensifies northward across the basin; in contrast, its upper part (above 400 m) is located in the central/eastern basin, well east of the most productive regions along the western boundary. The Bay-of-Bengal OMZ (BBOMZ), although strong, is weaker than the ASOMZ. To investigate the processes that maintain the Indian-Ocean OMZs, we obtain a suite of solutions to a coupled biological/physical model. Its physical component is a variable-density, 6 1/2-layer model, in which each layer corresponds to a distinct dynamical regime or water-mass type. Its biological component has six compartments: nutrients, phytoplankton, zooplankton, two size classes of detritus, and oxygen. Because the model grid is non-eddy resolving (0.5 degrees), the biological model also includes a parameterization of enhanced mixing based on the eddy kinetic energy derived from satellite observations. To explore further the impact of local processes on OMZs, we also obtain analytic solutions to a one-dimensional, simplified version of the biological model. Our control run is able to simulate basic features of the oxygen, nutrient, and phytoplankton fields throughout the Indian Ocean. The model OMZs result from a balance, or lack thereof, between a sink of oxygen by remineralization and subsurface oxygen sources due primarily to northward spreading of oxygenated water from the Southern Hemisphere, with a contribution from Persian-Gulf water in the northern Arabian Sea. The northward intensification of the lower ASOMZ results mostly from horizontal mixing since advection is weak in its depth range. The eastward shift of the upper ASOMZ is due primarily to enhanced advection and vertical eddy mixing in the western Arabian Sea, which spread oxygenated waters both horizontally and vertically. Advection carries small detritus from the western boundary into the central/eastern Arabian Sea, where it provides an additional source of remineralization that drives the ASOMZ to suboxic levels. The model BBOMZ is weaker than the ASOMZ because the Bay lacks a remote source of detritus from the western boundary. Although detritus has a prominent annual cycle, the model OMZs do not because there is not enough time for significant remineralization to occur.