An Optimizing Code Generator for a Class of Lattice-Boltzmann Computations

The Lattice-Boltzmann method (LBM), a promising new particle-based simulation technique for complex and multiscale fluid flows, has seen tremendous adoption in recent years in computational fluid dynamics. Even with a state-of-the-art LBM solver such as Palabos, a user has to still manually write the program using library-supplied primitives. We propose an automated code generator for a class of LBM computations with the objective to achieve high performance on modern architectures. Few studies have looked at time tiling for LBM codes. We exploit a key similarity between stencils and LBM to enable polyhedral optimizations and in turn time tiling for LBM. We also characterize the performance of LBM with the Roofline performance model. Experimental results for standard LBM simulations like Lid Driven Cavity, Flow Past Cylinder, and Poiseuille Flow show that our scheme consistently outperforms Palabos-on average by up to 3x while running on 16 cores of an Intel Xeon (Sandybridge). We also obtain an improvement of 2.47x on the SPEC LBM benchmark.

Field induced domain switching as the origin of anomalous lattice strain along non-polar direction in rhombohedral BiScO3-PbTiO3 close to the morphotropic phase boundary

The lattice strain and domain switching behavior of xBiScO(3)-(1-x) PbTiO3 (x = 0.40) was investigated as a function of cyclic field and grain orientation by in situ X-ray diffraction during application of electric fields. The electric field induced 200 lattice strain was measured to be five times larger than the 111 lattice strain in pseudorhombohedral xBiScO(3)-(1-x) PbTiO3 (x = 0.40). It is shown that the anomalous 200 lattice strain is not an intrinsic phenomenon, but arises primarily due to stress associated with the reorientation of the 111 domains in dense polycrystalline ceramic. (C) 2015 AIP Publishing LLC.

Anomalous Dielectric Features and Structural Implications in Gallic Acid-lsoniazid Cocrystal Mediated by Lattice Water

An organic supramolecular ternary salt (gallic acid:isoniazid:water; GINZH) examined earlier for its proton conducting characteristics is observed to display step-like dielectric behavior across the structural phase transition mediated by loss of water of hydration at 389 K. The presence of hydration in the crystal lattice along with proton mobility between acid base pairs controls the ``ferroelectric like'' behavior until the phase transition temperature.

Renormalization group invariants and sum rules in the deflected mirage mediation supersymmetry breaking

We examine the deflected mirage mediation supersymmetry breaking (DMMSB) scenario, which combines three supersymmetry breaking scenarios, namely anomaly mediation, gravity mediation and gauge mediation using the one-loop renormalization group invariants (RGIs). We examine the effects on the RGIs at the threshold where the gauge messengers emerge, and derive the supersymmetry breaking parameters in terms of the RGIs. We further discuss whether the supersymmetry breaking mediation mechanism can be determined using a limited set of invariants, and derive sum rules valid for DMMSB below the gauge messenger scale. In addition we examine the implications of the measured Higgs mass for the DMMSB spectrum.

Effect of Super-exchange Interaction on Ground state Magnetic Properties of Spin-dependent Falicov-Kimball Model on a Triangular Lattice

Ground state magnetic properties are studied by incorporating the super-exchange interaction (J(se)) in the spin-dependent Falicov-Kimball model (FKM) between localized (f-) electrons on a triangular lattice for half filled case. Numerical diagonalization and Monte-Carlo simulation are used to study the ground state magnetic properties. We have found that the magnetic moment of (d-) and (f-) electrons strongly depend on the value of Hund's exchange (J), super-exchange interaction (J(se)) and also depends on the number of (d-) electrons (N-d). The ground state changes from antiferromagnetic (AFM) to ferromagnetic (FM) state as we decrease (N-d). Also the density of d electrons at each site depends on the value of J and J(se).

The role of intra- and inter-site exchange correlations in the extended Falicov-Kimball model on a triangular lattice

Ground state magnetic properties of the spin-dependent Falicov-Kimball model (FKM) are studied by incorporating the intrasite exchange correlation J (between itinerant d- and localized f-electrons) and intersite (superexchange) correlation J (between localized f-electrons) on a triangular lattice for two different fillings. Numerical diagonalization and Monte-Carlo techniques are used to determine the ground state magnetic properties. Transitions from antiferromagnetic to ferromagnetic and again to re-entrant antiferromagnetic phase is observed in a wide range of parameter space. The magnetic moments of d- and f-electrons are observed to depend strongly on the value off, J and also on the total number of d-electrons (N-d). (C) 2015 Elsevier Ltd. All rights reserved.

Efficient density matrix renormalization group algorithm to study Y junctions with integer and half-integer spin

An efficient density matrix renormalization group (DMRG) algorithm is presented and applied to Y junctions, systems with three arms of n sites that meet at a central site. The accuracy is comparable to DMRG of chains. As in chains, new sites are always bonded to the most recently added sites and the superblock Hamiltonian contains only new or once renormalized operators. Junctions of up to N = 3n + 1 approximate to 500 sites are studied with antiferromagnetic (AF) Heisenberg exchange J between nearest-neighbor spins S or electron transfer t between nearest neighbors in half-filled Hubbard models. Exchange or electron transfer is exclusively between sites in two sublattices with N-A not equal N-B. The ground state (GS) and spin densities rho(r) = < S-r(z)> at site r are quite different for junctions with S = 1/2, 1, 3/2, and 2. The GS has finite total spin S-G = 2S(S) for even (odd) N and for M-G = S-G in the S-G spin manifold, rho(r) > 0(< 0) at sites of the larger (smaller) sublattice. S = 1/2 junctions have delocalized states and decreasing spin densities with increasing N. S = 1 junctions have four localized S-z = 1/2 states at the end of each arm and centered on the junction, consistent with localized states in S = 1 chains with finite Haldane gap. The GS of S = 3/2 or 2 junctions of up to 500 spins is a spin density wave with increased amplitude at the ends of arms or near the junction. Quantum fluctuations completely suppress AF order in S = 1/2 or 1 junctions, as well as in half-filled Hubbard junctions, but reduce rather than suppress AF order in S = 3/2 or 2 junctions.

NH3 adsorption on PtM (Fe, Co, Ni) surfaces: Cooperating effects of charge transfer, magnetic ordering and lattice strain

Adsorption of a molecule or group with an atom which is less electronegative than oxygen (0) and directly interacting with the surface is very relevant to development of PtM (M = 3d-transition metal) catalysts with high activity. Here, we present theoretical analysis of the adsorption of NH3 molecule (N being less electronegative than 0) on (111) surfaces of PtM (Fe, Co, Ni) alloys using the first principles density functional approach. We find that, while NH3-Pt interaction is stronger than that of NH3 with the elemental M-surfaces, it is weaker than the strength of interaction of NH3 with M-site on the surface of PtM alloy. (C) 2016 Published by Elsevier B.V.

Traveling waves in trimer granular lattice I: Bifurcation structure of traveling waves in the unit-cell model

Present paper is the first one in the series devoted to the dynamics of traveling waves emerging in the uncompressed, tri-atomic granular crystals. This work is primarily concerned with the dynamics of one-dimensional periodic granular trimer (tri-atomic) chains in the state of acoustic vacuum. Each unit cell consists of three spherical particles of different masses subject to periodic boundary conditions. Hertzian interaction law governs the mutual interaction of these particles. Under the assumption of zero pre-compression, this interaction is modeled as purely nonlinear, which means the absence of linear force component. The dynamics of such chains is governed by the two system parameters that scale the mass ratios between the particles of the unit cell. Such a system supports two different classes of periodic solutions namely the traveling and standing waves. The primary objective of the present study is the numerical analysis of the bifurcation structure of these solutions with emphasis on the dynamics of traveling waves. In fact, understanding of the bifurcation structure of the traveling wave solutions emerging in the unit-cell granular trimer is rather important and can shed light on the more complex nonlinear wave phenomena emerging in semi-infinite trimer chains. (c) 2016 Elsevier B.V. All rights reserved.

Screened and unscreened phases in sedimenting suspensions

A coarse-grained stochastic hydrodynamical description of velocity and concentration fluctuations in steadily sedimenting suspensions is constructed and analyzed using self-consistent and renormalization-group methods. We find a nonequilibrium phase transition from an "unscreened" phase in which we recover the Caflisch-Luke [Phys. Fluids 28, 759 (1985)] divergence of the velocity variance to a "screened" phase where the fluctuations have a finite correlation length depending on the volume fraction phi as phi(-1/3), in agreement with Segre et al. [Phys. Rev. Lett. 79, 2574 (1997)] (if their observation of a phi-independent diffusivity is used), and the velocity variance is independent of system size.

Self-organized criticality in dynamics without branching

We demonstrate the phenomenon of self-organized criticality (SOC) in a simple random walk model described by a random walk of a myopic ant, i.e., a walker who can see only nearest neighbors. The ant acts on the underlying lattice aiming at uniform digging, i.e., reduction of the height profile of the surface but is unaffected by the underlying lattice. In one, two, and three dimensions we have explored this model and have obtained power laws in the time intervals between consecutive events of "digging." Being a simple random walk, the power laws in space translate to power laws in time. We also study the finite size scaling of asymptotic scale invariant process as well as dynamic scaling in this system. This model differs qualitatively from the cascade models of SOC.

Spin populations in one-dimensional alternant spin systems: A theoretical approach

Spin-density maps, deduced from polarized neutron diffraction experiments, for both the pair and chain compounds of the system Mn2+Cu2+ have been reported recently. These results have motivated us to investigate theoretically the spin populations in such alternant mixed-spin systems. In this paper, we report our studies on the one-dimensional ferrimagnetic systems (S-A,S-B)(N) where hi is the number of AB pairs. We have considered all cases in which the spin Sri takes on allowed values in the range I to 7/2 while the spin S-B is held fixed at 1/2. The theoretical studies have been carried out on the isotropic Heisenberg model, using the density matrix renormalization group method. The effect of the magnitude of the larger spin SA On the quantum fluctuations in both A and B sublattices has been studied as a function of the system size N. We have investigated systems with both periodic and open boundary conditions, the latter with a view to understanding end-of-chain effects. The spin populations have been followed as a function of temperature as well as an applied magnetic field. High-magnetic fields are found to lead to interesting re-entrant behavior. The ratio of spin populations P-A-P-B is not sensitive to temperature at low temperatures.

Exact ground state and kink-like excitations of a two-dimensional Heisenberg antiferromagnet

A rare example of a two-dimensional Heisenberg model with an exact dimerized ground state is presented. This model, which can be regarded as a variation on the kagome' lattice, has several features of interest: it has a highly (but not macroscopically) degenerate ground state; it is closely related to spin chains studied by earlier authors; in particular, it exhibits domain-wall-like "kink" excitations normally associated only with one-dimensional systems. In some limits it decouples into noninteracting chains; unusually, this happens in the limit of strong, rather than weak, interchain coupling. [S0163-1829(99)50338-X].

Surface chemical studies on mica and galena using dextrin

The interactions of dextrin with biotite mica and galena have been investigated through adsorption, flotation, and electrokinetic measurements. The adsorption densities of dextrin onto mica continuously increase with increase of pH, while those onto galena show a maximum at pH 11.5. It is observed that the adsorption density of dextrin onto galena is quite high compared to that on mica. Both the adsorption isotherms exhibit Langmuirian behavior. Electrokinetic measurements portray conformational rearrangements of macromolecules with the loading, resulting in a shift of the shear plane, further away from the interface. Dissolution experiments indicate release of the lattice metal ions from mica and galena. Coprecipitation tests confirm polymer-metal ion interaction in the bulk solution. Dextrin does not exhibit any depressant action toward mica, whereas, with galena, the flotation recovery is decreased with an increase in pH beyond 9, in the presence of dextrin, complementing the adsorption results. Differential flotation results on a synthetic mixture of mica and galena show that mica can be selectively separated from galena using dextrin as a depressant for galena above pH 10. Possible mechanisms of interaction between dextrin and mica/galena are discussed.

Synthesis and properties of main chain discotic liquid crystalline polyethers based on rufigallol

Two series of thermotropic main chain discotic liquid crystalline polyethers, PR4m-n, based on rufigallol were prepared starting from the symmetric tetraethers of rufigallol, R4m; m and n represent the number of carbon atoms in the side chain and spacer segment, respectively. The symmetric tetraethers were in turn readily prepared by selective alkylation of rufigallol under controlled phase-transfer conditions. GPC analysis of the polymers suggested that they were all of moderate molecular weights, with M-n varying between 5400 and 17 000. The length of the spacer segment n in these polyethers was systematically varied, and its effect on the phase transition temperatures and the mesophase structure was examined using DSC, polarized light microscopy, and X-ray diffraction. It is noticed that when the spacer lengths are relatively long(n greater than or equal to 2m), the isotropization temperature (TD-i) decreases as the spacer length n increases, an observation that is in accordance with those previously made. However, when the spacer lengths are relatively small (n < 2m), the dependence of TD-i is quite the opposite; TD-i actually increases with an increase in spacer length. Furthermore, X-ray diffraction studies indicate that, in the discotic columnar mesophases that are formed, the columns pack in a hexagonal manner when n greater than or equal to 2m, while they do so in a rectangular lattice when n < 2m, leading to the formation of Dh and Dr mesophases, respectively. Finally, comparison of the discotic polyethers with their low molar mass analogues confirms the role of polymerization in stabilizing the mesophase; while all the polymers exhibit columnar mesophases, some of their low molar mass analogues are not liquid crystalline.