Dynamical localization in a chain of hard core bosons under periodic driving

We study the dynamics of a one-dimensional lattice model of hard core bosons which is initially in a superfluid phase with a current being induced by applying a twist at the boundary. Subsequently, the twist is removed, and the system is subjected to periodic delta-function kicks in the staggered on-site potential. We present analytical expressions for the current and work done in the limit of an infinite number of kicks. Using these, we show that the current (work done) exhibits a number of dips (peaks) as a function of the driving frequency and eventually saturates to zero (a finite value) in the limit of large frequency. The vanishing of the current (and the saturation of the work done) can be attributed to a dynamic localization of the hard core bosons occurring as a consequence of the periodic driving. Remarkably, we show that for some specific values of the driving amplitude, the localization occurs for any value of the driving frequency. Moreover, starting from a half-filled lattice of hard core bosons with the particles localized in the central region, we show that the spreading of the particles occurs in a light-cone-like region with a group velocity that vanishes when the system is dynamically localized.

Majorana edge modes in the Kitaev model

We study the Majorana modes, both equilibrium and Floquet, which can appear at the edges of the Kitaev model on the honeycomb lattice. We first present the analytical solutions known for the equilibrium Majorana edge modes for both zigzag and armchair edges of a semi-infinite Kitaev model and chart the parameter regimes in which they appear. We then examine how edge modes can be generated if the Kitaev coupling on the bonds perpendicular to the edge is varied periodically in time as periodic delta-function kicks. We derive a general condition for the appearance and disappearance of the Floquet edge modes as a function of the drive frequency for a generic d-dimensional integrable system. We confirm this general condition for the Kitaev model with a finite width by mapping it to a one-dimensional model. Our numerical and analytical study of this problem shows that Floquet Majorana modes can appear on some edges in the kicked system even when the corresponding equilibrium Hamiltonian has no Majorana mode solutions on those edges. We support our analytical studies by numerics for a finite sized system which show that periodic kicks can generate modes at the edges and the corners of the lattice.

Generalized scaling of misorientation angle distributions at meso-scale in deformed materials

Scaling behaviour has been observed at mesoscopic level irrespective of crystal structure, type of boundary and operative micro-mechanisms like slip and twinning. The presence of scaling at the meso-scale accompanied with that at the nano-scale clearly demonstrates the intrinsic spanning for different deformation processes and a true universal nature of scaling. The origin of a 1/2 power law in deformation of crystalline materials in terms of misorientation proportional to square root of strain is attributed to importance of interfaces in deformation processes. It is proposed that materials existing in three dimensional Euclidean spaces accommodate plastic deformation by one dimensional dislocations and their interaction with two dimensional interfaces at different length scales. This gives rise to a 1/2 power law scaling in materials. This intrinsic relationship can be incorporated in crystal plasticity models that aim to span different length and time scales to predict the deformation response of crystalline materials accurately.

Confining Dirac electrons on a topological insulator surface using potentials and a magnetic field

We study the effects of extended and localized potentials and a magnetic field on the Dirac electrons residing at the surface of a three-dimensional topological insulator like Bi2Se3. We use a lattice model to numerically study the various states; we show how the potentials can be chosen in a way which effectively avoids the problem of fermion doubling on a lattice. We show that extended potentials of different shapes can give rise to states which propagate freely along the potential but decay exponentially away from it. For an infinitely long potential barrier, the dispersion and spin structure of these states are unusual and these can be varied continuously by changing the barrier strength. In the presence of a magnetic field applied perpendicular to the surface, these states become separated from the gapless surface states by a gap, thereby giving rise to a quasi-one-dimensional system. Similarly, a magnetic field along with a localized potential can give rise to exponentially localized states which are separated from the surface states by a gap and thereby form a zero-dimensional system. Finally, we show that a long barrier and an impurity potential can produce bound states which are localized at the impurity, and an ``L''-shaped potential can have both bound states at the corner of the L and extended states which travel along the arms of the potential. Our work opens the way to constructing wave guides for Dirac electrons.

Ising model on a random network with annealed or quenched disorder

We study the equilibrium properties of an Ising model on a disordered random network where the disorder can be quenched or annealed. The network consists of fourfold coordinated sites connected via variable length one-dimensional chains. Our emphasis is on nonuniversal properties and we consider the transition temperature and other equilibrium thermodynamic properties, including those associated with one-dimensional fluctuations arising from the chains. We use analytic methods in the annealed case, and a Monte Carlo simulation for the quenched disorder. Our objective is to study the difference between quenched and annealed results with a broad random distribution of interaction parameters. The former represents a situation where the time scale associated with the randomness is very long and the corresponding degrees of freedom can be viewed as frozen, while the annealed case models the situation where this is not so. We find that the transition temperature and the entropy associated with one-dimensional fluctuations are always higher for quenched disorder than in the annealed case. These differences increase with the strength of the disorder up to a saturating value. We discuss our results in connection to physical systems where a broad distribution of interaction strengths is present.

Atom Lithography

We review the two kinds of forces that near-resonant light exerts on atoms the spontaneous force that is used for laser cooling, and the stimulated force that is used for coherent manipulation of atoms. We will discuss an experiment where laser cooling is used to collimate an atomic beam of sodium atoms, and the stimulated force within one period of a one-dimensional standing wave is used as a lens to focus the atoms to a narrow line about 20 nm wide. This kind of atom lithography is an example of the general field of atom optics in which light is used to manipulate atoms.

Organic solar cell by using vertically aligned nanostructured ZnO nanorods

A low temperature solution approach was employed to grow zinc oxide (ZnO) nanorods with various aspect ratios. Various sizes (diameter-10-25nm) of the nanorods were grown by changing the concentrations of the growth solution. The length (50nm-500nm) of nanorods was controlled using growth times. These one-dimensional (1D) nanostructures with direct paths for a charge transport with high surface area for light harvesting, are promising candidates for organic photovoltaics (OPV). The structural and optical properties of the prepared ZnO nanorods have been studied using SEM, XRD and UV-Vis absorption spectroscopy. Using as-grown ZnO inverted OPV was fabricated. ZnO nanorods were subjected to various doses of UV-ozone irradiation which led to improvement in transmission and hence enhanced device performance.

Isotropic finite volume discretization

Finite volume methods traditionally employ dimension by dimension extension of the one-dimensional reconstruction and averaging procedures to achieve spatial discretization of the governing partial differential equations on a structured Cartesian mesh in multiple dimensions. This simple approach based on tensor product stencils introduces an undesirable grid orientation dependence in the computed solution. The resulting anisotropic errors lead to a disparity in the calculations that is most prominent between directions parallel and diagonal to the grid lines. In this work we develop isotropic finite volume discretization schemes which minimize such grid orientation effects in multidimensional calculations by eliminating the directional bias in the lowest order term in the truncation error. Explicit isotropic expressions that relate the cell face averaged line and surface integrals of a function and its derivatives to the given cell area and volume averages are derived in two and three dimensions, respectively. It is found that a family of isotropic approximations with a free parameter can be derived by combining isotropic schemes based on next-nearest and next-next-nearest neighbors in three dimensions. Use of these isotropic expressions alone in a standard finite volume framework, however, is found to be insufficient in enforcing rotational invariance when the flux vector is nonlinear and/or spatially non-uniform. The rotationally invariant terms which lead to a loss of isotropy in such cases are explicitly identified and recast in a differential form. Various forms of flux correction terms which allow for a full recovery of rotational invariance in the lowest order truncation error terms, while preserving the formal order of accuracy and discrete conservation of the original finite volume method, are developed. Numerical tests in two and three dimensions attest the superior directional attributes of the proposed isotropic finite volume method. Prominent anisotropic errors, such as spurious asymmetric distortions on a circular reaction-diffusion wave that feature in the conventional finite volume implementation are effectively suppressed through isotropic finite volume discretization. Furthermore, for a given spatial resolution, a striking improvement in the prediction of kinetic energy decay rate corresponding to a general two-dimensional incompressible flow field is observed with the use of an isotropic finite volume method instead of the conventional discretization. (C) 2014 Elsevier Inc. All rights reserved.

Finite size scaling in crossover among different random matrix ensembles in microscopic lattice models

Using numerical diagonalization we study the crossover among different random matrix ensembles (Poissonian, Gaussian orthogonal ensemble (GOE), Gaussian unitary ensemble (GUE) and Gaussian symplectic ensemble (GSE)) realized in two different microscopic models. The specific diagnostic tool used to study the crossovers is the level spacing distribution. The first model is a one-dimensional lattice model of interacting hard-core bosons (or equivalently spin 1/2 objects) and the other a higher dimensional model of non-interacting particles with disorder and spin-orbit coupling. We find that the perturbation causing the crossover among the different ensembles scales to zero with system size as a power law with an exponent that depends on the ensembles between which the crossover takes place. This exponent is independent of microscopic details of the perturbation. We also find that the crossover from the Poissonian ensemble to the other three is dominated by the Poissonian to GOE crossover which introduces level repulsion while the crossover from GOE to GUE or GOE to GSE associated with symmetry breaking introduces a subdominant contribution. We also conjecture that the exponent is dependent on whether the system contains interactions among the elementary degrees of freedom or not and is independent of the dimensionality of the system.

Hetero-tri-spin 2p-3d-4f] Chain Compounds Based on Nitronyl Nitroxide Lanthanide Metallo-Ligands: Synthesis, Structure, and Magnetic Properties

Employing nitronyl nitroxide lanthanide(III) complexes as metallo-ligands allowed the efficient and highly selective preparation of three series of unprecedented heterotri-spin (Cu Ln-radical) one-dimensional compounds. These 2p-3d-4f spin systems, namely Ln(3)Cu(hfac)II(NitPhOAII)41 (Ln(III)=Gd 1(Gd), Tb 1(Tb), Dy 1(Dy); NitPhOAII=2-(4'-allyloxyphenyl)-4,4,5,5-tetramethylimidazoline-1-oxyl-3- oxide), Ln(3)Cu(hfac)II(NitPhOPO4] (1-nrn=Gd 2Gd, Tb 2Tb, Dy 2(Dy), Ho 2HOf Yb 2yb; NitPhOPr= 2-(4'-propoxyphenyI)-4,4,5,5-tetramethyl-imidazoline-1-oxyl-3-oxide) and Ln3Cu(hfac)II(NitPhOB441 (LnIm=Gd 3Gd, Tb 3Tb, Dy 3(Dy); NitPhOBz=2-(4'-benzyloxy- phenyl)-4,4,5,5-tetramethyl-imidazoline-1-oxyl-3-oxide) involve O-bound nitronyl nitroxide radicals as bridging ligands in chain structures with a Cu-Nit-Ln-Nit-Ln-Nit-Ln-Nit] repeating unit. The dc magnetic studies show that ferromagnetic metal radical interactions take place in these heterotri-spin chain complexes, these and the next-neighbor interactions have been quantified for the Gd derivatives. Complexes 1Tb and 2Tb exhibit frequency dependence of ac magnetic susceptibilities, indicating single-chain magnet behavior.

In a hot bubble: why does superbubble feedback work, but isolated supernovae do not?

Using idealized one-dimensional Eulerian hydrodynamic simulations, we contrast the behaviour of isolated supernovae with the superbubbles driven by multiple, collocated supernovae. Continuous energy injection via successive supernovae exploding within the hot/dilute bubble maintains a strong termination shock. This strong shock keeps the superbubble over-pressured and drives the outer shock well after it becomes radiative. Isolated supernovae, in contrast, with no further energy injection, become radiative quite early (less than or similar to 0.1Myr, tens of pc), and stall at scales less than or similar to 100 pc. We show that isolated supernovae lose almost all of their mechanical energy by 1 Myr, but superbubbles can retain up to similar to 40 per cent of the input energy in the form of mechanical energy over the lifetime of the star cluster (a few tens of Myr). These conclusions hold even in the presence of realistic magnetic fields and thermal conduction. We also compare various methods for implementing supernova feedback in numerical simulations. For various feedback prescriptions, we derive the spatial scale below which the energy needs to be deposited in order for it to couple to the interstellar medium. We show that a steady thermal wind within the superbubble appears only for a large number (greater than or similar to 10(4)) of supernovae. For smaller clusters, we expect multiple internal shocks instead of a smooth, dense thermalized wind.

Self-Assembled, Aligned ZnO Nanorod Buffer Layers for High-Current-Density, Inverted Organic Photovoltaics

Two different soft-chemical, self-assembly-based solution approaches are employed to grow zinc oxide (ZnO) nanorods with controlled texture. The methods used involve seeding and growth on a substrate. Nanorods with various aspect ratios (1-5) and diameters (15-65 nm) are grown. Obtaining highly oriented rods is determined by the way the substrate is mounted within the chemical bath. Furthermore, a preheat and centrifugation step is essential for the optimization of the growth solution. In the best samples, we obtain ZnO nanorods that are almost entirely oriented in the (002) direction; this is desirable since electron mobility of ZnO is highest along this crystallographic axis. When used as the buffer layer of inverted organic photovoltaics (I-OPVs), these one-dimensional (1D) nanostructures offer: (a) direct paths for charge transport and (b) high interfacial area for electron collection. The morphological, structural, and optical properties of ZnO nanorods are studied using scanning electron microscopy, X-ray diffraction, and ultraviolet-visible light (UV-vis) absorption spectroscopy. Furthermore, the surface chemical features of ZnO films are studied using X-ray photoelectron spectroscopy and contact angle measurements. Using as-grown ZnO, inverted OPVs are fabricated and characterized. For improving device performance, the ZnO nanorods are subjected to UV-ozone irradiation. UV-ozone treated ZnO nanorods show: (i) improvement in optical transmission, (ii) increased wetting of active organic components, and (iii) increased concentration of Zn-O surface bonds. These observations correlate well with improved device performance. The devices fabricated using these optimized buffer layers have an efficiency of similar to 3.2% and a fill factor of 0.50; this is comparable to the best I-OPVs reported that use a P3HT-PCBM active layer.

Modelling of transient heat conduction with diffuse interface methods

We present a survey on different numerical interpolation schemes used for two-phase transient heat conduction problems in the context of interface capturing phase-field methods. Examples are general transport problems in the context of diffuse interface methods with a non-equal heat conductivity in normal and tangential directions to the interface. We extend the tonsorial approach recently published by Nicoli M et al (2011 Phys. Rev. E 84 1-6) to the general three-dimensional (3D) transient evolution equations. Validations for one-dimensional, two-dimensional and 3D transient test cases are provided, and the results are in good agreement with analytical and numerical reference solutions.

Riemann shock tube: 1D normal shocks in air, simulations and experiments

This paper presents numerical simulation of the evolution of one-dimensional normal shocks, their propagation, reflection and interaction in air using a single diaphragm Riemann shock tube and validate them using experimental results. Mathematical model is derived for one-dimensional compressible flow of viscous and conducting medium. Dimensionless form of the mathematical model is used to construct space-time finite element processes based on minimization of the space-time residual functional. The space-time local approximation functions for space-time p-version hierarchical finite elements are considered in higher order GRAPHICS] spaces that permit desired order of global differentiability of local approximations in space and time. The resulting algebraic systems from this approach yield unconditionally positive-definite coefficient matrices, hence ensure unique numerical solution. The evolution is computed for a space-time strip corresponding to a time increment Delta t and then time march to obtain the evolution up to any desired value of time. Numerical studies are designed using recently invented hand-driven shock tube (Reddy tube) parameters, high/low side density and pressure values, high- and low-pressure side shock tube lengths, so that numerically computed results can be compared with actual experimental measurements.

Appraisal of observational method for consolidation analysis

Two of the aims of laboratory one-dimensional consolidation tests are prediction of the end of primary settlement, and determination of the coefficient of consolidation of soils required for the time rate of consolidation analysis from time-compression data. Of the many methods documented in the literature to achieve these aims, Asaoka's method is a simple and useful tool, and yet the most neglected one since its inception in the geotechnical engineering literature more than three decades ago. This paper appraises Asaoka's method, originally proposed for the field prediction of ultimate settlement, from the perspective of laboratory consolidation analysis along with recent developments. It is shown through experimental illustrations that Asaoka's method is simpler than the conventional and popular methods, and makes a satisfactory prediction of both the end of primary compression and the coefficient of consolidation from laboratory one-dimensional consolidation test data.