Two-dimensional flow past circular cylinders using finite volume lattice Boltzmann formulation

In this article, an extension to the total variation diminishing finite volume formulation of the lattice Boltzmann equation method on unstructured meshes was presented. The quadratic least squares procedure is used for the estimation of first-order and second-order spatial gradients of the particle distribution functions. The distribution functions were extrapolated quadratically to the virtual upwind node. The time integration was performed using the fourth-order RungeKutta procedure. A grid convergence study was performed in order to demonstrate the order of accuracy of the present scheme. The formulation was validated for the benchmark two-dimensional, laminar, and unsteady flow past a single circular cylinder. These computations were then investigated for the low Mach number simulations. Further validation was performed for flow past two circular cylinders arranged in tandem and side-by-side. Results of these simulations were extensively compared with the previous numerical data. Copyright (C) 2011 John Wiley & Sons, Ltd.

Search on a fractal lattice using a quantum random walk

The spatial search problem on regular lattice structures in integer number of dimensions d >= 2 has been studied extensively, using both coined and coinless quantum walks. The relativistic Dirac operator has been a crucial ingredient in these studies. Here, we investigate the spatial search problem on fractals of noninteger dimensions. Although the Dirac operator cannot be defined on a fractal, we construct the quantum walk on a fractal using the flip-flop operator that incorporates a Klein-Gordon mode. We find that the scaling behavior of the spatial search is determined by the spectral (and not the fractal) dimension. Our numerical results have been obtained on the well-known Sierpinski gaskets in two and three dimensions.

Activation of Lattice Oxygen of TiO2 by Pd2+ Ion: Correlation of Low-Temperature CO and Hydrocarbon Oxidation with Structure of Ti1-xPdxO2-x (x=0.01-0.03)

Lattice oxygen of TiO2 is activated by the substitution of Pd ion in its lattice. Ti1-xPdxO2-x (x = 0.01-0.03) have been synthesized by solution combustion method crystallizing in anatase TiO2 structure. Pd is in +2 oxidation state and Ti is in +4 oxidation state in the catalyst. Pd is more ionic in TiO2 lattice compared to Pd in PdO. Oxygen storage capacity defined by ``amount of oxygen that is used reversibly to oxidize CO'' is as high as 5100 mu mol/g of Ti0.97Pd0.03O1.97. Oxygen is extracted by CO to CO2 in absence of feed oxygen even at room temperature which is more than 20 times compared to pure TiO2. Rate of CO oxidation is 2.75 mu mol g(-1) s(-1) at 60 degrees C over Ti0.97Pd0.03O1.97 and C2H2 gets oxidized to CO2 and H2O at room temperature. Catalyst is not poisoned on long time operation of the reactor. Such high catalytic activity is due to activated lattice oxygen created by the substitution of Pd ion as seen from first-principles density functional theory (DFT) calculations with 96 atom supercells of Ti32O64, Ti31Pd1O63, Ti30Pd2O62, and Ti29Pd3O61. The compounds crystallize in anatase TiO2 structure with Pd2+ ion in nearly square planar geometry and TiO6 octahedra are distorted by the creation of weakly bound oxygens. Structural analysis of Ti31Pd1O63 which is close to 3% Pd ion substituted TiO2 shows that oxygens associated with both Ti and Pd ions in the lattice show bond valence sum of 1.87, a low value characteristic of weak oxygen in the lattice compared to oxygens with valence 2 and above in the same lattice. Exact positions of activated oxygens have been identified in the lattice from DFT calculations.

Correlation between enhanced lattice polarizability and high piezoelectric response in BiScO3-PbTiO3

Piezoelectric and ex situ electric-field induced structural studies were carried out on closely spaced compositions in the morphotropic phase boundary region of (1 - x) PbTiO3-(x)BiScO3. While the common approach of zero field structural analysis failed to provide a unique relationship between the anomalous piezoresponse of x = 0.3725 and structural factor(s), ex situ study of electric-field induced structural changes revealed that the composition exhibiting the highest piezoelectric response is the one which also exhibits significantly enhanced polarizability of the lattices of both coexisting (monoclinic and tetragonal) phases. The enhanced lattice polarizability manifests as a significant fraction of the monoclinic phase transforming irreversibly to the tetragonal phase after electric poling. DOI: 10.1103/PhysRevB.87.064106

Chapman-Enskog analysis for finite-volume formulation of lattice Boltzmann equation

The classical Chapman-Enskog expansion is performed for the recently proposed finite-volume formulation of lattice Boltzmann equation (LBE) method D.V. Patil, K.N. Lakshmisha, Finite volume TVD formulation of lattice Boltzmann simulation on unstructured mesh, J. Comput. Phys. 228 (2009) 5262-5279]. First, a modified partial differential equation is derived from a numerical approximation of the discrete Boltzmann equation. Then, the multi-scale, small parameter expansion is followed to recover the continuity and the Navier-Stokes (NS) equations with additional error terms. The expression for apparent value of the kinematic viscosity is derived for finite-volume formulation under certain assumptions. The attenuation of a shear wave, Taylor-Green vortex flow and driven channel flow are studied to analyze the apparent viscosity relation.

Cl-35 NQR frequency and spin lattice relaxation time in 3,4-dichlorophenol as a function of pressure and temperature

The pressure dependences of Cl-35 nuclear quadrupole resonance (NQR) frequency, temperature and pressure variation of spin lattice relaxation time (T-1) were investigated in 3,4-dichlorophenol. T-1 was measured in the temperature range 77-300 K. Furthermore, the NQR frequency and T-1 for these compounds were measured as a function of pressure up to 5 kbar at 300 K. The temperature dependence of the average torsional lifetimes of the molecules and the transition probabilities W-1 and W-2 for the Delta m = +/- 1 and Delta m = +/- 2 transitions were also obtained. A nonlinear variation of NQR frequency with pressure has been observed and the pressure coefficients were observed to be positive. A thermodynamic analysis of the data was carried out to determine the constant volume temperature coefficients of the NQR frequency. An attempt is made to compare the torsional frequencies evaluated from NQR data with those obtained by IR spectra. On selecting the appropriate mode from IR spectra, a good agreement with torsional frequency obtained from NQR data is observed. The previously mentioned approach is a good illustration of the supplementary nature of the data from IR studies, in relation to NQR studies of compounds in solid state.

Strong-coupling expansion for ultracold bosons in an optical lattice at finite temperatures in the presence of superfluidity

We develop a strong-coupling (t << U) expansion technique for calculating the density profile for bosonic atoms trapped in an optical lattice with an overall harmonic trap at finite temperature and finite on-site interaction in the presence of superfluid regions. Our results match well with quantum Monte Carlo simulations at finite temperature. We also show that the superfluid order parameter never vanishes in the trap due to the proximity effect. Our calculations for the scaled density in the vacuum-to-superfluid transition agree well with the experimental data for appropriate temperatures. We present calculations for the entropy per particle as a function of temperature which can be used to calibrate the temperature in experiments. We also discuss issues connected with the demonstration of universal quantum critical scaling in the experiments.

Atomic Structure of Quantum Gold Nanowires: Quantification of the Lattice Strain

Theoretical studies exist to compute the atomic arrangement in gold nanowires and the influence on their electronic behavior with decreasing diameter. Experimental studies, e.g., by transmission electron microscopy, on chemically synthesized ultrafine wires are however lacking owing to the unavailability of suitable protocols for sample preparation and the stability of the wires under electron beam irradiation. In this work, we present an atomic scale structural investigation on quantum single crystalline gold nanowires of 2 nm diameter, chemically prepared on a carbon film grid. Using low dose aberration-corrected high resolution (S)TEM, we observe an inhomogeneous strain distribution in the crystal, largely concentrated at the twin boundaries and the surface along with the presence of facets and surface steps leading to a noncircular cross section of the wires. These structural aspects are critical inputs needed to determine their unique electronic character and their potential as a suitable catalyst material. Furthermore, electron-beam-induced structural changes at the atomic scale, having implications on their mechanical behavior and their suitability as interconnects, are discussed.

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.

Electrostatic modulation of periodic potentials in a two-dimensional electron gas: from antidot lattice to quantum dot lattice

We use a dual gated device structure to introduce a gate-tuneable periodic potential in a GaAs/AlGaAs two dimensional electron gas (2DEG). Using only a suitable choice of gate voltages we can controllably alter the potential landscape of the bare 2DEG, inducing either a periodic array of antidots or quantum dots. Antidots are artificial scattering centers, and therefore allow for a study of electron dynamics. In particular, we show that the thermovoltage of an antidot lattice is particularly sensitive to the relative positions of the Fermi level and the antidot potential. A quantum dot lattice, on the other hand, provides the opportunity to study correlated electron physics. We find that its current-voltage characteristics display a voltage threshold, as well as a power law scaling, indicative of collective Coulomb blockade in a disordered background.

Theoretical description of coherent doublon creation via lattice modulation spectroscopy

Using a recently developed strong-coupling method, we present a comprehensive theory for doublon production processes in modulation spectroscopy of a three-dimensional system of ultracold fermionic atoms in an optical lattice with a trap. The theoretical predictions compare well to the experimental time traces of doublon production. For experimentally feasible conditions, we provide a quantitative prediction for the presence of a nonlinear ``two-photon'' excitation at strong modulation amplitudes.

Lattice Reduction Aided Detection in Large-MIMO Systems

Lattice reduction (LR) aided detection algorithms are known to achieve the same diversity order as that of maximum-likelihood (ML) detection at low complexity. However, they suffer SNR loss compared to ML performance. The SNR loss is mainly due to imperfect orthogonalization and imperfect nearest neighbor quantization. In this paper, we propose an improved LR-aided (ILR) detection algorithm, where we specifically target to reduce the effects of both imperfect orthogonalization and imperfect nearest neighbor quantization. The proposed ILR detection algorithm is shown to achieve near-ML performance in large-MIMO systems and outperform other LR-aided detection algorithms in the literature. Specifically, the SNR loss incurred by the proposed ILR algorithm compared to ML performance is just 0.1 dB for 4-QAM and < 0.5 dB for 16-QAM in 16 x 16 V-BLAST MIMO system. This performance is superior compared to those of other LR-aided detection algorithms, whose SNR losses are in the 2 dB to 9 dB range.

Fermionic superfluid from a bilayer band insulator in an optical lattice

We propose a model to realize a fermionic superfluid state in an optical lattice circumventing the cooling problem. Our proposal exploits the idea of tuning the interaction in a characteristically low-entropy state, a band insulator in an optical bilayer system, to obtain a superfluid. By performing a detailed analysis of the model including fluctuations and augmented by a variational quantum Monte Carlo calculation of the ground state, we show that the superfluid state obtained has a high transition temperature of the order of the hopping energy. Our system is designed to suppress other competing orders such as a charge density wave. We suggest a laboratory realization of this model via an orthogonally shaken optical lattice bilayer.

Study of ground state phases for spin-1/2 Falicov-Kimball model on a triangular lattice

The spin dependent Falicov-Kimball model (FKM) is studied on a triangular lattice using numerical diagonalization technique and Monte-Carlo simulation algorithm. Magnetic properties have been explored for different values of parameters: on-site Coulomb correlation U, exchange interaction J and filling of electrons. We have found that the ground state configurations exhibit long range Neel order, ferromagnetism or a mixture of both as J is varied. The magnetic moments of itinerant (d) and localized U) electrons are also studied. For the one-fourth filling case we found no magnetic moment from d- and f-electrons for U less than a critical value. `.2014 Elsevier Ltd. All rights reserved.

Optimal sequential wireless relay placement on a random lattice path

Our work is motivated by impromptu (or ``as-you-go'') deployment of wireless relay nodes along a path, a need that arises in many situations. In this paper, the path is modeled as starting at the origin (where there is the data sink, e.g., the control center), and evolving randomly over a lattice in the positive quadrant. A person walks along the path deploying relay nodes as he goes. At each step, the path can, randomly, either continue in the same direction or take a turn, or come to an end, at which point a data source (e.g., a sensor) has to be placed, that will send packets to the data sink. A decision has to be made at each step whether or not to place a wireless relay node. Assuming that the packet generation rate by the source is very low, and simple link-by-link scheduling, we consider the problem of sequential relay placement so as to minimize the expectation of an end-to-end cost metric (a linear combination of the sum of convex hop costs and the number of relays placed). This impromptu relay placement problem is formulated as a total cost Markov decision process. First, we derive the optimal policy in terms of an optimal placement set and show that this set is characterized by a boundary (with respect to the position of the last placed relay) beyond which it is optimal to place the next relay. Next, based on a simpler one-step-look-ahead characterization of the optimal policy, we propose an algorithm which is proved to converge to the optimal placement set in a finite number of steps and which is faster than value iteration. We show by simulations that the distance threshold based heuristic, usually assumed in the literature, is close to the optimal, provided that the threshold distance is carefully chosen. (C) 2014 Elsevier B.V. All rights reserved.