Simulation of inhomogeneous distributions of ultracold atoms in an optical lattice via a massively parallel implementation of nonequilibrium strong-coupling perturbation theory

We present a nonequilibrium strong-coupling approach to inhomogeneous systems of ultracold atoms in optical lattices. We demonstrate its application to the Mott-insulating phase of a two-dimensional Fermi-Hubbard model in the presence of a trap potential. Since the theory is formulated self-consistently, the numerical implementation relies on a massively parallel evaluation of the self-energy and the Green's function at each lattice site, employing thousands of CPUs. While the computation of the self-energy is straightforward to parallelize, the evaluation of the Green's function requires the inversion of a large sparse 10(d) x 10(d) matrix, with d > 6. As a crucial ingredient, our solution heavily relies on the smallness of the hopping as compared to the interaction strength and yields a widely scalable realization of a rapidly converging iterative algorithm which evaluates all elements of the Green's function. Results are validated by comparing with the homogeneous case via the local-density approximation. These calculations also show that the local-density approximation is valid in nonequilibrium setups without mass transport.

Free turbulent shear layer in a point vortex gas as a problem in nonequilibrium statistical mechanics

This paper attempts to unravel any relations that may exist between turbulent shear flows and statistical mechanics through a detailed numerical investigation in the simplest case where both can be well defined. The flow considered for the purpose is the two-dimensional (2D) temporal free shear layer with a velocity difference Delta U across it, statistically homogeneous in the streamwise direction (x) and evolving from a plane vortex sheet in the direction normal to it (y) in a periodic-in-x domain L x +/-infinity. Extensive computer simulations of the flow are carried out through appropriate initial-value problems for a ``vortex gas'' comprising N point vortices of the same strength (gamma = L Delta U/N) and sign. Such a vortex gas is known to provide weak solutions of the Euler equation. More than ten different initial-condition classes are investigated using simulations involving up to 32 000 vortices, with ensemble averages evaluated over up to 10(3) realizations and integration over 10(4)L/Delta U. The temporal evolution of such a system is found to exhibit three distinct regimes. In Regime I the evolution is strongly influenced by the initial condition, sometimes lasting a significant fraction of L/Delta U. Regime III is a long-time domain-dependent evolution towards a statistically stationary state, via ``violent'' and ``slow'' relaxations P.-H. Chavanis, Physica A 391, 3657 (2012)], over flow time scales of order 10(2) and 10(4)L/Delta U, respectively (for N = 400). The final state involves a single structure that stochastically samples the domain, possibly constituting a ``relative equilibrium.'' The vortex distribution within the structure follows a nonisotropic truncated form of the Lundgren-Pointin (L-P) equilibrium distribution (with negatively high temperatures; L-P parameter lambda close to -1). The central finding is that, in the intermediate Regime II, the spreading rate of the layer is universal over the wide range of cases considered here. The value (in terms of momentum thickness) is 0.0166 +/- 0.0002 times Delta U. Regime II, extensively studied in the turbulent shear flow literature as a self-similar ``equilibrium'' state, is, however, a part of the rapid nonequilibrium evolution of the vortex-gas system, which we term ``explosive'' as it lasts less than one L/Delta U. Regime II also exhibits significant values of N-independent two-vortex correlations, indicating that current kinetic theories that neglect correlations or consider them as O(1/N) cannot describe this regime. The evolution of the layer thickness in present simulations in Regimes I and II agree with the experimental observations of spatially evolving (3D Navier-Stokes) shear layers. Further, the vorticity-stream-function relations in Regime III are close to those computed in 2D Navier-Stokes temporal shear layers J. Sommeria, C. Staquet, and R. Robert, J. Fluid Mech. 233, 661 (1991)]. These findings suggest the dominance of what may be called the Kelvin-Biot-Savart mechanism in determining the growth of the free shear layer through large-scale momentum and vorticity dispersal.

Azide/thiocyanate incorporated cobalt(III)-Schiff base complexes: Characterizations and catalytic activity in aerobic epoxidation of olefins

Reaction of cobalt(II) perchlorate hexahydrate with a potentially tetradentate Schiff base ligand, HL (2-methoxy-6-(2-diethylaminoethylimino)methyl]phenol) in presence of sodium azide and sodium thiocyanate yields two complexes Co( L)( HL)(N-3)]center dot ClO4 ( 1) and Co( L)( HL)(NCS)] center dot ClO4 ( 2); both being characterized by different physicochemical methods. Crystal structure of 1 was determined by single crystal X-ray diffraction while that of 2 was reported earlier. In 1, the central cobalt(III) adopts slightly distorted octahedral geometry with same donor set to that of 2. Catalytic efficacy of the complexes towards epoxidation of different alkenes under aerobic condition were investigated in homogeneous medium which reveals that 1 is better catalyst than 2 with respect to alkene oxidation, reflected from the turn over frequencies (TOF) measured at an optimum temperature of 60 degrees C in acetonitrile. (C) 2014 Published by Elsevier B.V.

Regional flood frequency analysis using kernel- based fuzzy clustering approach

Regionalization approaches are widely used in water resources engineering to identify hydrologically homogeneous groups of watersheds that are referred to as regions. Pooled information from sites (depicting watersheds) in a region forms the basis to estimate quantiles associated with hydrological extreme events at ungauged/sparsely gauged sites in the region. Conventional regionalization approaches can be effective when watersheds (data points) corresponding to different regions can be separated using straight lines or linear planes in the space of watershed related attributes. In this paper, a kernel-based Fuzzy c-means (KFCM) clustering approach is presented for use in situations where such linear separation of regions cannot be accomplished. The approach uses kernel-based functions to map the data points from the attribute space to a higher-dimensional space where they can be separated into regions by linear planes. A procedure to determine optimal number of regions with the KFCM approach is suggested. Further, formulations to estimate flood quantiles at ungauged sites with the approach are developed. Effectiveness of the approach is demonstrated through Monte-Carlo simulation experiments and a case study on watersheds in United States. Comparison of results with those based on conventional Fuzzy c-means clustering, Region-of-influence approach and a prior study indicate that KFCM approach outperforms the other approaches in forming regions that are closer to being statistically homogeneous and in estimating flood quantiles at ungauged sites. Key Points Kernel-based regionalization approach is presented for flood frequency analysis Kernel procedure to estimate flood quantiles at ungauged sites is developed A set of fuzzy regions is delineated in Ohio, USA

From (Au5Sn + AuSn) physical mixture to phase pure AuSn and Au5Sn intermetallic nanocrystals with tailored morphology: digestive ripening assisted approach

Here we present digestive ripening facilitated interatomic diffusion for the phase controlled synthesis of homogeneous intermetallic nanocrystals of Au-Sn system. Au and Sn metal nanoparticles synthesized by a solvated metal atom dispersion (SMAD) method are employed as precursors for the fabrication of AuSn and Au5Sn which are Au-rich Au-Sn intermetallic nanocrystals. By optimizing the stoichiometry of Au and Sn in the reaction mixture, and by employing growth directing agents, the formation of phase pure intermetallic AuSn and Au5Sn nanocrystals could be realized. The as-prepared Au and Sn colloidal nanoparticles and the resulting intermetallic nanocrystals are thoroughly characterized by powder X-ray diffraction, transmission electron microscopy (TEM and STEM-EDS), and optical spectroscopy. The results obtained here demonstrate the potential of solution chemistry which allows synthesizing phase pure Au-Sn intermetallics with tailored morphology.

Modeling of Bubble-Size Distribution in Water and Freon Co-Blown Free Rise Polyurethane Foams

A model has been developed to simulate the foam characteristics obtained, when chemical (water) and physical (Freon) blowing agents are used together for the formation of polyurethane foams. The model considers the rate of reaction, the consequent rise in temperature of the reaction mixture, nucleation of bubbles, and mass transfer of CO2 and Freon to them till the time of gelation. The model is able to explain the experimental results available in literature. It further predicts that the nucleation period gets reduced with increase in water (at constant Freon content), whereas with increase in Freon (at constant water) concentration nucleation period decreases marginally leading to narrower bubble-size distribution. By the use of uniform sized nuclei added initially, the model predicts that the bubble-size distribution can be made independent of the rate of homogeneous nucleation and can, thus, offer an extra parameter for its control. (C) 2014 Wiley Periodicals, Inc.

Flag structure for operators in the Cowen-Douglas class

The explicit description of homogeneous operators and localization of a Hilbert module naturally leads to the definition of a class of Cowen-Douglas operators possessing a flag structure. These operators are irreducible. We show that the flag structure is rigid in the sense that the unitary equivalence class of the operator and the flag structure determine each other. We obtain a complete set of unitary invariants which are somewhat more tractable than those of an arbitrary operator in the Cowen-Douglas class. (C) 2014 Academie des sciences. Published by Elsevier Masson SAS. All rights reserved.

Statistics of the inverse-cascade regime in two-dimensional magnetohydrodynamic turbulence

We present a detailed direct numerical simulation of statistically steady, homogeneous, isotropic, two-dimensional magnetohydrodynamic turbulence. Our study concentrates on the inverse cascade of the magnetic vector potential. We examine the dependence of the statistical properties of such turbulence on dissipation and friction coefficients. We extend earlier work significantly by calculating fluid and magnetic spectra, probability distribution functions (PDFs) of the velocity, magnetic, vorticity, current, stream-function, and magnetic-vector-potential fields, and their increments. We quantify the deviations of these PDFs from Gaussian ones by computing their flatnesses and hyperflatnesses. We also present PDFs of the Okubo-Weiss parameter, which distinguishes between vortical and extensional flow regions, and its magnetic analog. We show that the hyperflatnesses of PDFs of the increments of the stream function and the magnetic vector potential exhibit significant scale dependence and we examine the implication of this for the multiscaling of structure functions. We compare our results with those of earlier studies.

Maximal Contractive Tuples

Maximality of a contractive tuple of operators is considered. A characterization for a contractive tuple to be maximal is obtained. The notion of maximality for a submodule of the Drury-Arveson module on the -dimensional unit ball is defined. For , it is shown that every submodule of the Hardy module over the unit disc is maximal. But for we prove that any homogeneous submodule or submodule generated by polynomials is not maximal. A characterization of maximal submodules is obtained.

Octahedral distortion induced magnetic anomalies in LaMn0.5Co0.5O3 single crystals

Single crystals of LaMn0.5Co0.5O3 belonging to the ferromagnetic-insulator and distorted perovskite class were grown using a four-mirror optical float zone furnace. The as-grown crystal crystallizes into an orthorhombic Pbnm structure. The spatially resolved 2D Raman scan reveals a strain-induced distribution of transition metal (TM)-oxygen (O) octahedral deformation in the as-grown crystal. A rigorous annealing process releases the strain, thereby generating homogeneous octahedral distortion. The octahedra tilt by reducing the bond angle TM-O-TM, resulting in a decline of the exchange energy in the annealed crystal. The critical behavior is investigated from the bulk magnetization. It is found that the ground state magnetic behavior assigned to the strain-free LaMn0.5Co0.5O3 crystal is of the 3D Heisenberg kind. Strain induces mean field-like interaction in some sites, and consequently, the critical exponents deviate from the 3D Heisenberg class in the as-grown crystal. The temperature-dependent Raman scattering study reveals strong spin-phonon coupling and the existence of two magnetic ground states in the same crystal. (C) 2014 AIP Publishing LLC.

A framework for post-silicon realization of arbitrary instruction extensions on reconfigurable data-paths

In this paper we present a framework for realizing arbitrary instruction set extensions (IE) that are identified post-silicon. The proposed framework has two components viz., an IE synthesis methodology and the architecture of a reconfigurable data-path for realization of the such IEs. The IE synthesis methodology ensures maximal utilization of resources on the reconfigurable data-path. In this context we present the techniques used to realize IEs for applications that demand high throughput or those that must process data streams. The reconfigurable hardware called HyperCell comprises a reconfigurable execution fabric. The fabric is a collection of interconnected compute units. A typical use case of HyperCell is where it acts as a co-processor with a host and accelerates execution of IEs that are defined post-silicon. We demonstrate the effectiveness of our approach by evaluating the performance of some well-known integer kernels that are realized as IEs on HyperCell. Our methodology for realizing IEs through HyperCells permits overlapping of potentially all memory transactions with computations. We show significant improvement in performance for streaming applications over general purpose processor based solutions, by fully pipelining the data-path. (C) 2014 Elsevier B.V. All rights reserved.

Thermoelectric properties of Zn doped Cu2SnSe3

Zn doped ternary compounds Cu2ZnxSn1-xSe3 (x = 0, 0.025, 0.05, 0.075) were prepared by solid state synthesis. The undoped compound showed a monoclinic crystal structure as a major phase, while the doped compounds showed a cubic crystal structure confirmed by powder XRD (X-Ray Diffraction). The surface morphology and elemental composition analysis for all the samples were studied by SEM (Scanning Electron Microscopy) and EPMA (Electron Probe Micro Analyzer), respectively. SEM micrographs of the hot pressed samples showed the presence of continuous and homogeneous grains confirming sufficient densification. Elemental composition of all the samples revealed an off-stoichiometry, which was determined by EPMA. Transport properties were measured between 324 K and 773 K. The electrical resistivity decreased up to the samples with Zn content x = 0.05 in Cu2ZnxSn1-xSe3, and slightly increased in the sample Cu2Zn0.075Sn0.925Se3. This behavior is consistent with the changes in the carrier concentration confirmed by room temperature Hall coefficient data. Temperature dependent electrical resistivity of all samples showed heavily doped semiconductor behavior. All the samples exhibit positive Seebeck coefficient (S) and Hall coefficient indicating that the majority of the carriers are holes. A linear increase in Seebeck coefficient with increase in temperature indicates the degenerate semiconductor behavior. The total thermal conductivity of the doped samples increased with a higher amount of doping, due to the increase in the carrier contribution. The total and lattice thermal conductivity of all samples showed 1/1 dependence, which points toward the dominance of phonon scattering at high temperatures. The maximum 1/TZF = 0.48 at 773 K was obtained for the sample Cu2SnSe3 due to a low thermal conductivity compared to the doped samples. (C) 2014 Elsevier B.V. All rights reserved.

Smoothed Functional Algorithms for Stochastic Optimization Using q-Gaussian Distributions

Smoothed functional (SF) schemes for gradient estimation are known to be efficient in stochastic optimization algorithms, especially when the objective is to improve the performance of a stochastic system However, the performance of these methods depends on several parameters, such as the choice of a suitable smoothing kernel. Different kernels have been studied in the literature, which include Gaussian, Cauchy, and uniform distributions, among others. This article studies a new class of kernels based on the q-Gaussian distribution, which has gained popularity in statistical physics over the last decade. Though the importance of this family of distributions is attributed to its ability to generalize the Gaussian distribution, we observe that this class encompasses almost all existing smoothing kernels. This motivates us to study SF schemes for gradient estimation using the q-Gaussian distribution. Using the derived gradient estimates, we propose two-timescale algorithms for optimization of a stochastic objective function in a constrained setting with a projected gradient search approach. We prove the convergence of our algorithms to the set of stationary points of an associated ODE. We also demonstrate their performance numerically through simulations on a queuing model.

Mixed finite elements for electromagnetic analysis

The occurrence of spurious solutions is a well-known limitation of the standard nodal finite element method when applied to electromagnetic problems. The two commonly used remedies that are used to address this problem are (i) The addition of a penalty term with the penalty factor based on the local dielectric constant, and which reduces to a Helmholtz form on homogeneous domains (regularized formulation); (ii) A formulation based on a vector and a scalar potential. Both these strategies have some shortcomings. The penalty method does not completely get rid of the spurious modes, and both methods are incapable of predicting singular eigenvalues in non-convex domains. Some non-zero spurious eigenvalues are also predicted by these methods on non-convex domains. In this work, we develop mixed finite element formulations which predict the eigenfrequencies (including their multiplicities) accurately, even for nonconvex domains. The main feature of the proposed mixed finite element formulation is that no ad-hoc terms are added to the formulation as in the penalty formulation, and the improvement is achieved purely by an appropriate choice of finite element spaces for the different variables. We show that the formulation works even for inhomogeneous domains where `double noding' is used to enforce the appropriate continuity requirements at an interface. For two-dimensional problems, the shape of the domain can be arbitrary, while for the three-dimensional ones, with our current formulation, only regular domains (which can be nonconvex) can be modeled. Since eigenfrequencies are modeled accurately, these elements also yield accurate results for driven problems. (C) 2014 Elsevier Ltd. All rights reserved.

An unusual demixing behavior in PS-PVME blends in the presence of nanoparticles

The effect of silver nanoparticles (sNP) on the demixing and the evolution of morphology in off-critical blends of 90/10 and 10/90 (wt/wt) PS/PVME polystyrene/poly(vinyl methyl ether)] was probed here using shear rheology and optical microscopy. The faster component (PVME) has a higher molecular weight (80 kDa) than the slower component (PS, 35 kDa), which makes this system quite interesting to study with respect to the evolving morphology, as the blends transit through the binodal and the spinodal envelopes. An unusual demixing behavior was observed in both PVME rich and PS rich blends. Temperature modulated differential scanning calorimetry measurements showed that the T-g value for the blends with sNP was slightly lower than that of the neat blends. A decreased volume of cooperativity at T-g suggests confined segmental dynamics in the presence of sNP. Although, the addition of sNP had no influence on the thermodynamic demixing temperature, it significantly altered the elasticity of the minor component during the transition of the blend from the homogeneous to the heterogeneous state. This is manifested from energetically driven localization of the sNP in the PVME phase during demixing. As a direct consequence of this, the formation of the microstructures upon demixing was observed to be delayed in the presence of sNP. Interestingly, in the intermediate quench depth, the higher viscoelastic phase evolved as an interconnected network, which subsequently coarsened into discrete droplets in the late stages for the 90/10 PS/PVME blends. Similar observations were made for 10/90 PS/PVME blends where threads of PVME appeared at deeper quench depths in the presence of sNP. The interconnected network formation of the minor phase (here PVME), which is also the faster component in the blend, was different from the usual demixing behavior.