Efficient numerical solution of the time fractional diffusion equation by mapping from its Brownian counterpart

The solution of a Caputo time fractional diffusion equation of order $0<\alpha<1$ is expressed in terms of the solution of a corresponding integer order diffusion equation. We demonstrate a linear time mapping between these solutions that allows for accelerated computation of the solution of the fractional order problem. In the context of an $N$-point finite difference time discretisation, the mapping allows for an improvement in time computational complexity from $O\left(N^2\right)$ to $O\left(N^\alpha\right)$, given a precomputation of $O\left(N^{1+\alpha}\ln N\right)$. The mapping is applied successfully to the least-squares fitting of a fractional advection diffusion model for the current in a time-of-flight experiment, resulting in a computational speed up in the range of one to three orders of magnitude for realistic problem sizes.Comment: 9 pages, 5 figures; added references for section

Spin-Peierls states of quantum antiferromagnets on the CaV4O9Ca V_4 O_9 lattice

We discuss the quantum paramagnetic phases of Heisenberg antiferromagnets on the 1/5-depleted square lattice found in $Ca V_4 O_9$. The possible phases of the quantum dimer model on this lattice are obtained by a mapping to a quantum-mechanical height model. In addition to the ``decoupled'' phases found earlier, we find a possible intermediate spin-Peierls phase with spontaneously-broken lattice symmetry. Experimental signatures of the different quantum paramagnetic phases are discussed.Comment: 9 pages; 2 eps figure

TEMPERATURE-DEPENDENCE OF DOMAIN-WALL COERCIVE FIELD IN MAGNETIC GARNET-FILMS

The coercive properties of magnetically uniaxial liquid-phase epitaxy garnet films were investigated between 10 K and the Neel temperature (T(N) less-than-or-equal-to 500 K). Two independent methods, the results of which are nearly identical (magnetical response of oscillating domain walls and the method of coercive loops measured in a vibrating sample magnetometer), were used. Besides the usual domain-wall coercive field, H(dw), the critical coercive pressure, p(dw), was also introduced as it describes in a direct way the interactions of the domain walls with the wall-pinning traps. Both H(dw) and p(dw) were found to increase exponentially with decreasing temperature. Three different types of wall-pinning traps were identified in the sample and their strength, their rate of change with temperature, and their temperature range of activity were determined

QUAGMIRE v1.3: a quasi-geostrophic model for investigating rotating fluids experiments

QUAGMIRE is a quasi-geostrophic numerical model for performing fast, high-resolution simulations of multi-layer rotating annulus laboratory experiments on a desktop personal computer. The model uses a hybrid finite-difference/spectral approach to numerically integrate the coupled nonlinear partial differential equations of motion in cylindrical geometry in each layer. Version 1.3 implements the special case of two fluid layers of equal resting depths. The flow is forced either by a differentially rotating lid, or by relaxation to specified streamfunction or potential vorticity fields, or both. Dissipation is achieved through Ekman layer pumping and suction at the horizontal boundaries, including the internal interface. The effects of weak interfacial tension are included, as well as the linear topographic beta-effect and the quadratic centripetal beta-effect. Stochastic forcing may optionally be activated, to represent approximately the effects of random unresolved features. A leapfrog time stepping scheme is used, with a Robert filter. Flows simulated by the model agree well with those observed in the corresponding laboratory experiments

Correlation among Cirrus Ice Content, Water Vapor and Temperature in the TTL as Observed by CALIPSO and Aura-MLS

Water vapor in the tropical tropopause layer (TTL) has a local radiative cooling effect. As a source for ice in cirrus clouds, however, it can also indirectly produce infrared heating. Using NASA A-Train satellite measurements of CALIPSO and Aura/MLS we calculated the correlation of water vapor, ice water content and temperature in the TTL. We find that temperature strongly controls water vapor (correlation r =0.94) and cirrus clouds at 100 hPa (r = 0.91). Moreover we observe that the cirrus seasonal cycle is highly (r =0.9) anticorrelated with the water vapor variation in the TTL, showing higher cloud occurrence during December-January-February. We further investigate the anticorrelation on a regional scale and find that the strong anticorrelation occurs generally in the ITCZ (Intertropical Convergence Zone). The seasonal cycle of the cirrus ice water content is also highly anticorrelated to water vapor (r = 0.91) and our results support the hypothesis that the total water at 100 hPa is roughly constant. Temperature acts as a main regulator for balancing the partition between water vapor and cirrus clouds. Thus, to a large extent, the depleting water vapor in the TTL during DJF is a manifestation of cirrus formation

Interaction effects on 2D fermions with random hopping

We study the effects of generic short-ranged interactions on a system of 2D Dirac fermions subject to a special kind of static disorder, often referred to as ``chiral.'' The non-interacting system is a member of the disorder class BDI [M. R. Zirnbauer, J. Math. Phys. 37, 4986 (1996)]. It emerges, for example, as a low-energy description of a time-reversal invariant tight-binding model of spinless fermions on a honeycomb lattice, subject to random hopping, and possessing particle-hole symmetry. It is known that, in the absence of interactions, this disordered system is special in that it does not localize in 2D, but possesses extended states and a finite conductivity at zero energy, as well as a strongly divergent low-energy density of states. In the context of the hopping model, the short-range interactions that we consider are particle-hole symmetric density-density interactions. Using a perturbative one-loop renormalization group analysis, we show that the same mechanism responsible for the divergence of the density of states in the non-interacting system leads to an instability, in which the interactions are driven strongly relevant by the disorder. This result should be contrasted with the limit of clean Dirac fermions in 2D, which is stable against the inclusion of weak short-ranged interactions. Our work suggests a novel mechanism wherein a clean system, initially insensitive to interaction effects, can be made unstable to interactions upon the inclusion of weak static disorder.Comment: 16 pages, 10 figures; References added, figures enlarged; to be published in Phys. Rev.

Scaling and Crossover Functions for the Conductance in the Directed Network Model of Edge States

We consider the directed network (DN) of edge states on the surface of a cylinder of length L and circumference C. By mapping it to a ferromagnetic superspin chain, and using a scaling analysis, we show its equivalence to a one-dimensional supersymmetric nonlinear sigma model in the scaling limit, for any value of the ratio L/C, except for short systems where L is less than of order C^{1/2}. For the sigma model, the universal crossover functions for the conductance and its variance have been determined previously. We also show that the DN model can be mapped directly onto the random matrix (Fokker-Planck) approach to disordered quasi-one-dimensional wires, which implies that the entire distribution of the conductance is the same as in the latter system, for any value of L/C in the same scaling limit. The results of Chalker and Dohmen are explained quantitatively.Comment: 10 pages, REVTeX, 2 eps figure

Small angle neutron scattering observation of chain retraction after a large step deformation

The process of retraction in entangled linear chains after a fast nonlinear stretch was detected from time-resolved but quenched small angle neutron scattering (SANS) experiments on long, well-entangled polyisoprene chains. The statically obtained SANS data cover the relevant time regime for retraction, and they provide a direct, microscopic verification of this nonlinear process as predicted by the tube model. Clear, quantitative agreement is found with recent theories of contour length fluctuations and convective constraint release, using parameters obtained mainly from linear rheology. The theory captures the full range of scattering vectors once the crossover to fluctuations on length scales below the tube diameter is accounted for

Mesoscale theory of grains and cells: crystal plasticity and coarsening

Solids with spatial variations in the crystalline axes naturally evolve into cells or grains separated by sharp walls. Such variations are mathematically described using the Nye dislocation density tensor. At high temperatures, polycrystalline grains form from the melt and coarsen with time: the dislocations can both climb and glide. At low temperatures under shear the dislocations (which allow only glide) form into cell structures. While both the microscopic laws of dislocation motion and the macroscopic laws of coarsening and plastic deformation are well studied, we hitherto have had no simple, continuum explanation for the evolution of dislocations into sharp walls. We present here a mesoscale theory of dislocation motion. It provides a quantitative description of deformation and rotation, grounded in a microscopic order parameter field exhibiting the topologically conserved quantities. The topological current of the Nye dislocation density tensor is derived from a microscopic theory of glide driven by Peach-Koehler forces between dislocations using a simple closure approximation. The resulting theory is shown to form sharp dislocation walls in finite time, both with and without dislocation climb.Comment: 5 pages, 3 figure