The hard-particle model for a dense granular flow down an inclined plane

Perfectly hard particles are those which experience an infinite repulsive force when they overlap, and no force when they do not overlap. In the hard-particle model, the only static state is the isostatic state where the forces between particles are statically determinate. In the flowing state, the interactions between particles are instantaneous because the time of contact approaches zero in the limit of infinite particle stiffness. Here, we discuss the development of a hard particle model for a realistic granular flow down an inclined plane, and examine its utility for predicting the salient features both qualitatively and quantitatively. We first discuss Discrete Element simulations, that even very dense flows of sand or glass beads with volume fraction between 0.5 and 0.58 are in the rapid flow regime, due to the very high particle stiffness. An important length scale in the shear flow of inelastic particles is the `conduction length' delta = (d/(1 - e(2))(1/2)), where d is the particle diameter and e is the coefficient of restitution. When the macroscopic scale h (height of the flowing layer) is larger than the conduction length, the rates of shear production and inelastic dissipation are nearly equal in the bulk of the flow, while the rate of conduction of energy is O((delta/h)(2)) smaller than the rate of dissipation of energy. Energy conduction is important in boundary layers of thickness delta at the top and bottom. The flow in the boundary layer at the top and bottom is examined using asymptotic analysis. We derive an exact relationship showing that the a boundary layer solution exists only if the volume fraction in the bulk decreases as the angle of inclination is increased. In the opposite case, where the volume fraction increases as the angle of inclination is increased, there is no boundary layer solution. The boundary layer theory also provides us with a way of understanding the cessation of flow when at a given angle of inclination when the height of the layer is decreased below a value h(stop), which is a function of the angle of inclination. There is dissipation of energy due to particle collisions in the flow as well as due to particle collisions with the base, and the fraction of energy dissipation in the base increases as the thickness decreases. When the shear production in the flow cannot compensate for the additional energy drawn out of the flow due to the wall collisions, the temperature decreases to zero and the flow stops. Scaling relations can be derived for h(stop) as a function of angle of inclination.

Heisenberg-like critical properties in ferromagnetic Nd1-xPbxMnO3 single crystals

Static magnetization for single crystals of insulating Nd0.85Pb0.15MnO3 and marginally conducting Nd0.70Pb0.30MnO3 has been studied around the ferromagnetic to paramagnetic transition temperature T-C. Results of measurements carried out in the critical range vertical bar(T - T-C)/T-C vertical bar <= 0.1 are reported. Critical exponents beta and gamma for the thermal behaviour of magnetization and susceptibility have been obtained both by modified Arrott plots and the Kouvel-Fisher method. The exponent delta independently obtained from the critical isotherm was found to satisfy the Widom scaling relation delta = gamma/beta + 1. For both compositions the values of exponents are consistent with those expected for isotropic magnets belonging to the Heisenberg universality class with short-range exchange in three dimensions. Correspondingly, the specific heat displays only a cusp-like anomaly at the critical temperature of these crystals which is consistent with an exponent alpha < 0. The results show that the ferromagnetic ordering transition in Nd1-xPbxMnO3 in the composition range 0.15 <= x <= 0.40 is continuous. This mixed-valent manganite displays the conventional properties of a Heisenberg-like ferromagnet, irrespective of the differing transport properties and in spite of low ordering temperatures T-C = 109 and 147.2 K for x = 0.15 and 0.30, respectively.

Spatial survival probability for one-dimensional fluctuating interfaces in the steady state

We report numerical and analytic results for the spatial survival probability for fluctuating one-dimensional interfaces with Edwards-Wilkinson or Kardar-Parisi-Zhang dynamics in the steady state. Our numerical results are obtained from analysis of steady-state profiles generated by integrating a spatially discretized form of the Edwards-Wilkinson equation to long times. We show that the survival probability exhibits scaling behavior in its dependence on the system size and the "sampling interval" used in the measurement for both "steady-state" and "finite" initial conditions. Analytic results for the scaling functions are obtained from a path-integral treatment of a formulation of the problem in terms of one-dimensional Brownian motion. A "deterministic approximation" is used to obtain closed-form expressions for survival probabilities from the formally exact analytic treatment. The resulting approximate analytic results provide a fairly good description of the numerical data.

Memory effect in Dy0.5Sr0.5MnO3 single crystals

We have performed a series of magnetic aging experiments on single crystals of Dy0.5Sr0.5MnO3. The results demonstrate striking memory and chaos-like effects in this insulating half-doped perovskite manganite and suggest the existence of strong magnetic relaxation mechanisms of a clustered magnetic state. The spin-glass-like state established below a temperature T-sg approximate to 34 K originates from quenched disorder arising due to the ionic-radii mismatch at the rare earth site. However, deviations from the typical behavior seen in canonical spin glass materials are observed which indicate that the glassy magnetic properties are due to cooperative and frustrated dynamics in a heterogeneous or clustered magnetic state. In particular, the microscopic spin flip time obtained from dynamical scaling near the spin glass freezing temperature is four orders of magnitude larger than microscopic times found in atomic spin glasses. The magnetic viscosity deduced from the time dependence of the zero-field-cooled magnetization exhibits a peak at a temperature T < T-sg and displays a marked dependence on waiting time in zero field.

Infrared Spectra of Dimethylquinolines in the Gas Phase: Experiment and Theory

Infrared spectra of atmospherically important dimethylquinolines (DMQs), namely 2,4-DMQ, 2,6-DMQ, 2,7-DMQ, and 2,8-DMQ in the gas phase at 80 degrees C were recorded using a long variable path-length cell. DFT calculations were carried out to assign the bands in the experimentally observed spectra at the B3LYP/6-31G* level of theory. The spectral assignments particularly for the C-H stretching modes could not be made unambiguously using calculated anharmonic or scaled harmonic frequencies. To resolve this problem, a scaled force field method of assignment was used. Assignment of fundamental modes was confirmed by potential energy distributions (PEDs) of the normal modes derived by the scaled force fields using a modified version of the UMAT program in the QCPE package. We demonstrate that for large molecules such as the DMQs, the scaling of the force field is more effective in arriving at the correct assignment of the fundamentals for a quantitative vibrational analysis. An error analysis of the mean deviation of the calculated harmonic, anharmonic, and force field fitted frequencies from the observed frequency provides strong evidence for the correctness of the assignment.

Generalized asymptotic expansions for coupled wavenumbers in fluid-filled cylindrical shells

Analytical expressions are found for the coupled wavenumbers in an infinite fluid-filled cylindrical shell using the asymptotic methods. These expressions are valid for any general circumferential order (n).The shallow shell theory (which is more accurate at higher frequencies)is used to model the cylinder. Initially, the in vacua shell is dealt with and asymptotic expressions are derived for the shell wavenumbers in the high-and the low-frequency regimes. Next, the fluid-filled shell is considered. Defining a relevant fluid-loading parameter p, we find solutions for the limiting cases of small and large p. Wherever relevant, a frequency scaling parameter along with some ingenuity is used to arrive at an elegant asymptotic expression. In all cases.Poisson's ratio v is used as an expansion variable. The asymptotic results are compared with numerical solutions of the dispersion equation and the dispersion relation obtained by using the more general Donnell-Mushtari shell theory (in vacuo and fluid-filled). A good match is obtained. Hence, the contribution of this work lies in the extension of the existing literature to include arbitrary circumferential orders(n). (C) 2010 Elsevier Ltd. All rights reserved.

Theoretical Estimation of Length Dependent In-Plane Stiffness of Single Walled Carbon Nanotubes Using the Nonlocal Elasticity Theory

In the present paper, Eringen's nonlocal elasticity theory is employed to evaluate the length dependent in-plane stiffness of single-walled carbon nanotubes (SWCNTs). The SWCNT is modeled as an Euler-Bernoulli beam and is analyzed for various boundary conditions to evaluate the length dependent in-plane stiffness. It has been found that the nonlocal scaling parameter has a significant effect on the length dependent in-plane stiffness of SWCNTs. It has been observed that as the nonlocal scale parameter increases the stiffness ratio of SWCNT decreases. In nonlocality, the cantilever SWCNT has high in-plane stiffness as compared to the simply-supported and the clamped cases.

Functional Surfaces for Microfluidics in Proteomic Analysis

Microchips for use in biomolecular analysis show a lot of promise for medical diagnostics and biomedical basic research. Among the potential advantages are more sensitive and faster analyses as well as reduced cost and sample consumption. Due to scaling laws, the surface are to volume ratios of microfluidic chips is very high. Because of this, tailoring the surface properties and surface functionalization are very important technical issues for microchip development. This thesis studies two different types of functional surfaces, surfaces for open surface capillary microfluidics and surfaces for surface assisted laser desorption ionization mass spectrometry, and combinations thereof. Open surface capillary microfluidics can be used to transport and control liquid samples on easily accessible open surfaces simply based on surface forces, without any connections to pumps or electrical power sources. Capillary filling of open partially wetting grooves is shown to be possible with certain geometries, aspect ratios and contact angles, and a theoretical model is developed to identify complete channel filling domains, as well as partial filling domains. On the other hand, partially wetting surfaces with triangular microstructures can be used for achieving directional wetting, where the water droplets do not spread isotropically, but instead only spread to a predetermined sector. Furthermore, by patterning completely wetting and superhydrophobic areas on the same surface, complex droplet shapes are achieved, as the water stretches to make contact with the wetting surface, but does not enter into the superhydrophobic domains. Surfaces for surface assisted laser desorption ionization mass spectrometry are developed by applying various active thin film coatings on multiple substrates, in order to separate surface and bulk effects. Clear differences are observed between both surface and substrate layers. The best performance surfaces consisted of amorphous silicon coating and an inorganic-organic hybrid substrate, with nanopillars and nanopores. These surfaces are used for matrix-free ionization of drugs, peptides and proteins, and for some analytes, the detection limits were in the high attomoles. Microfluidics and laser desorption ionization surfaces are combined on a functionalized drying platforms, where the surface is used to control the shape of the deposited analyte droplet, and the shape of the initial analyte droplet affects the dried droplet solute deposition pattern. The deposited droplets can then directly detected by mass spectrometry. Utilizing this approach, results of analyte concentration, splitting and separation are demonstrated.

Phase-shift randomness in one-dimensional disordered conductors in the quasi-metallic domain

An invariant imbedding method yields exact analytical results for the distribution of the phase theta (L) of the reflection amplitude and for low-order resistance moments (pn) for a disordered conductor of length L in the quasi-metallic regime L<scaling of resistance in one dimension is thus obtained.

Investigation of the effect of nonlocal scale on ultrasonic wave dispersion characteristics of a monolayer graphene

In this article, an ultrasonic wave propagation in graphene sheet is studied using nonlocal elasticity theory incorporating small scale effects. The graphene sheet is modeled as an isotropic plate of one-atom thick. For this model, the nonlocal governing differential equations of motion are derived from the minimization of the total potential energy of the entire system. An ultrasonic type of wave propagation model is also derived for the graphene sheet. The nonlocal scale parameter introduces certain band gap region in in-plane and flexural wave modes where no wave propagation occurs. This is manifested in the wavenumber plots as the region where the wavenumber tends to infinite or wave speed tends to zero. The frequency at which this phenomenon occurs is called the escape frequency. The explicit expressions for cutoff frequencies and escape frequencies are derived. The escape frequencies are mainly introduced because of the nonlocal elasticity. Obviously these frequencies are function of nonlocal scaling parameter. It has also been obtained that these frequencies are independent of y-directional wavenumber. It means that for any type of nanostructure, the escape frequencies are purely a function of nonlocal scaling parameter only. It is also independent of the geometry of the structure. It has been found that the cutoff frequencies are function of nonlocal scaling parameter (e(0)a) and the y-directional wavenumber (k(y)). For a given nanostructure, nonlocal small scale coefficient can be obtained by matching the results from molecular dynamics (MD) simulations and the nonlocal elasticity calculations. At that value of the nonlocal scale coefficient, the waves will propagate in the nanostructure at that cut-off frequency. In the present paper, different values of e(o)a are used. One can get the exact e(0)a for a given graphene sheet by matching the MD simulation results of graphene with the results presented in this paper. (C) 2010 Elsevier B.V. All rights reserved.

Spatially dependent correlation functions in the Anderson model

We report results from a first principles calculation of spatially dependent correlation functions around a magnetic impurity in metals described by the nondegenerate Anderson model. Our computations are based on a combination of perturbative scaling theory and numerical renormalization group methods. Results for the conduction election charge density around the impurity and correlation functions involving the conduction electron and impurity charge and spin densities will be presented. The behavior in various regimes including the mixed valent regime will be explored.

Root system traits of Norway spruce, Scots pine, and silver birch in mixed boreal forests: an analysis of root architecture, morphology, and anatomy

The aim of this thesis was to unravel the functional-structural characteristics of root systems of Betula pendula Roth., Picea abies (L.) Karst., and Pinus sylvestris L. in mixed boreal forest stands differing in their developmental stage and site fertility. The root systems of these species had similar structural regularities: horizontally-oriented shallow roots defined the horizontal area of influence, and within this area, each species placed fine roots in the uppermost soil layers, while sinker roots defined the maximum rooting depth. Large radial spread and high ramification of coarse roots, and the high specific root length (SRL) and root length density (RLD) of fine roots indicated the high belowground competitiveness and root plasticity of B. pendula. Smaller radial root spread and sparser branching of coarse roots, and low SRL and RLD of fine roots of the conifers could indicate their more conservative resource use and high association with and dependence on ectomycorrhiza-forming fungi. The vertical fine root distributions of the species were mostly overlapping, implying the possibility for intense belowground competition for nutrients. In each species, conduits tapered and their frequency increased from distal roots to the stem, from the stem to the branches, and to leaf petioles in B. pendula. Conduit tapering was organ-specific in each species violating the assumptions of the general vascular scaling model (WBE). This reflects the hierarchical organization of a tree and differences between organs in the relative importance of transport, safety, and mechanical demands. The applied root model was capable of depicting the mass, length and spread of coarse roots of B. pendula and P. abies, and to the lesser extent in P. sylvestris. The roots did not follow self-similar fractal branching, because the parameter values varied within the root systems. Model parameters indicate differences in rooting behavior, and therefore different ecophysiological adaptations between species.

Search on a hypercubic lattice using a quantum random walk. I. d > 2

Random walks describe diffusion processes, where movement at every time step is restricted to only the neighboring locations. We construct a quantum random walk algorithm, based on discretization of the Dirac evolution operator inspired by staggered lattice fermions. We use it to investigate the spatial search problem, that is, to find a marked vertex on a d-dimensional hypercubic lattice. The restriction on movement hardly matters for d > 2, and scaling behavior close to Grover's optimal algorithm (which has no restriction on movement) can be achieved. Using numerical simulations, we optimize the proportionality constants of the scaling behavior, and demonstrate the approach to that for Grover's algorithm (equivalent to the mean-field theory or the d -> infinity limit). In particular, the scaling behavior for d = 3 is only about 25% higher than the optimal d -> infinity value.

Search on a hypercubic lattice using a quantum random walk. II. d = 2

We investigate the spatial search problem on the two-dimensional square lattice, using the Dirac evolution operator discretized according to the staggered lattice fermion formalism. d = 2 is the critical dimension for the spatial search problem, where infrared divergence of the evolution operator leads to logarithmic factors in the scaling behavior. As a result, the construction used in our accompanying article A. Patel and M. A. Rahaman, Phys. Rev. A 82, 032330 (2010)] provides an O(root N ln N) algorithm, which is not optimal. The scaling behavior can be improved to O(root N ln N) by cleverly controlling the massless Dirac evolution operator by an ancilla qubit, as proposed by Tulsi Phys. Rev. A 78, 012310 (2008)]. We reinterpret the ancilla control as introduction of an effective mass at the marked vertex, and optimize the proportionality constants of the scaling behavior of the algorithm by numerically tuning the parameters.

Experimental studies and phase field modeling of microstructure evolution during solidification with electromagnetic stirring

Thixocasting requires manufacturing of billets with non-dendritic microstructure. Aluminum alloy A356 billets were produced by rheocasting in a mould placed inside a linear electromagnetic stirrer. Subsequent heat treatment was used to produce a transition from rosette to globular microstructure. The current and the duration of stirring were explored as control parameters. Simultaneous induction heating of the billet during stirring was quantified using experimentally determined thermal profiles. The effect of processing parameters on the dendrite fragmentation was discussed. Corresponding computational modeling of the process was performed using phase-field modeling of alloy solidification in order to gain insight into the process of morphological changes of a solid during this process. A non-isothermal alloy solidification model was used for simulations. The morphological evolution under such imposed thermal cycles was simulated and compared with experimentally determined one. Suitable scaling using the thermosolutal diffusion distances was used to overcome computational difficulties in quantitative comparison at system scale. The results were interpreted in the light of existing theories of microstructure refinement and globularisation.