Jump at the onset of saltation

We reveal a discontinuous transition in the saturated flux for aeolian saltation by simulating explicitly particle motion in turbulent flow. The discontinuity is followed by a coexistence interval with two metastable solutions. The modification of the wind profile due to momentum exchange exhibits a second maximum at high shear strength. The saturated flux depends on the strength of the wind as $q_s=q_0+A(u_*-u_t)(u_*^2+u_t^2)$

On the analogy between streamlined magnetic and solid obstacles

Analogies are elaborated in the qualitative description of two systems: the magnetohydrodynamic (MHD) flow moving through a region where an external local magnetic field (magnetic obstacle) is applied, and the ordinary hydrodynamic flow around a solid obstacle. The former problem is of interest both practically and theoretically, and the latter one is a classical problem being well understood in ordinary hydrodynamics. The first analogy is the formation in the MHD flow of an impenetrable region -- core of the magnetic obstacle -- as the interaction parameter $N$, i.e. strength of the applied magnetic field, increases significantly. The core of the magnetic obstacle is streamlined both by the upstream flow and by the induced cross stream electric currents, like a foreign insulated insertion placed inside the ordinary hydrodynamic flow. In the core, closed streamlines of the mass flow resemble contour lines of electric potential, while closed streamlines of the electric current resemble contour lines of pressure. The second analogy is the breaking away of attached vortices from the recirculation pattern produced by the magnetic obstacle when the Reynolds number $Re$, i.e. velocity of the upstream flow, is larger than a critical value. This breaking away of vortices from the magnetic obstacle is similar to that occurring past a real solid obstacle. Depending on the inlet and/or initial conditions, the observed vortex shedding can be either symmetric or asymmetric.Comment: minor changes, accepted for PoF, 26 pages, 7 figure

Regularity for eigenfunctions of Schr\"odinger operators

We prove a regularity result in weighted Sobolev spaces (or Babuska--Kondratiev spaces) for the eigenfunctions of a Schr\"odinger operator. More precisely, let K_{a}^{m}(\mathbb{R}^{3N}) be the weighted Sobolev space obtained by blowing up the set of singular points of the Coulomb type potential V(x) = \sum_{1 \le j \le N} \frac{b_j}{|x_j|} + \sum_{1 \le i < j \le N} \frac{c_{ij}}{|x_i-x_j|}, x in \mathbb{R}^{3N}, b_j, c_{ij} in \mathbb{R}. If u in L^2(\mathbb{R}^{3N}) satisfies (-\Delta + V) u = \lambda u in distribution sense, then u belongs to K_{a}^{m} for all m \in \mathbb{Z}_+ and all a \le 0. Our result extends to the case when b_j and c_{ij} are suitable bounded functions on the blown-up space. In the single-electron, multi-nuclei case, we obtain the same result for all a<3/2.Comment: to appear in Lett. Math. Phy

The brain decade in debate: II. Panic or anxiety? From animal models to a neurobiological basis

This article is a transcription of an electronic symposium sponsored by the Brazilian Society of Neuroscience and Behavior (SBNeC). Invited researchers from the European Union, North America and Brazil discussed two issues on anxiety, namely whether panic is a very intense anxiety or something else, and what aspects of clinical anxiety are reproduced by animal models. Concerning the first issue, most participants agreed that generalized anxiety and panic disorder are different on the basis of clinical manifestations, drug response and animal models. Also, underlying brain structures, neurotransmitter modulation and hormonal changes seem to involve important differences. It is also common knowledge that existing animal models generate different types of fear/anxiety. A challenge for future research is to establish a good correlation between animal models and nosological classification.Universidade Federal do Paraná Departamento de Farmacologia Laboratório de Fisiologia e Farmacologia do Sistema Nervoso CentralUniversity of Hawaii Department of NeurobiologyUniversity of Hawaii Department of PsychologyUniversidade de São Paulo Faculdade de Filosofia Ciências e Letras de Ribeirão Preto Departamento de PsicobiologiaUniversidade de São Paulo Faculdade de Medicina de Ribeirão Preto Departamento de FisiologiaUniversidade de São Paulo Faculdade de Medicina de Ribeirão Preto Departamento de NeuropsiquiatriaUniversidade Federal de Santa Catarina Departamento de FarmacologiaCentral Nervous System Research Department Sanofi SynthelaboAston University Institute of Pharmaceutical SciencesHoffmann-La Roche Ltd.Universidade Federal de São Paulo (UNIFESP) Escola Paulista de Medicina Departamento de PsicologiaUniversity of Leeds Department of Psychology Ethopharmacology LaboratoryUniversidade Federal do Espírito Santo Centro de Biomedicina Departamento de Ciências FisiológicasUNIFESP, EPM, Depto. de PsicologiaSciEL

Numerical Computations with H(div)-Finite Elements for the Brinkman Problem

The H(div)-conforming approach for the Brinkman equation is studied numerically, verifying the theoretical a priori and a posteriori analysis in previous work of the authors. Furthermore, the results are extended to cover a non-constant permeability. A hybridization technique for the problem is presented, complete with a convergence analysis and numerical verification. Finally, the numerical convergence studies are complemented with numerical examples of applications to domain decomposition and adaptive mesh refinement.Comment: Minor clarifications, added references. Reordering of some figures. To appear in Computational Geosciences, final article available at http://www.springerlink.co

Regularized Linear Inversion with Randomized Singular Value Decomposition

In this work, we develop efficient solvers for linear inverse problems based on randomized singular value decomposition (RSVD). This is achieved by combining RSVD with classical regularization methods, e.g., truncated singular value decomposition, Tikhonov regularization, and general Tikhonov regularization with a smoothness penalty. One distinct feature of the proposed approach is that it explicitly preserves the structure of the regularized solution in the sense that it always lies in the range of a certain adjoint operator. We provide error estimates between the approximation and the exact solution under canonical source condition, and interpret the approach in the lens of convex duality. Extensive numerical experiments are provided to illustrate the efficiency and accuracy of the approach.Comment: 20 pages, 4 figure

Tips for implementing multigrid methods on domains containing holes

As part of our development of a computer code to perform 3D `constrained evolution' of Einstein's equations in 3+1 form, we discuss issues regarding the efficient solution of elliptic equations on domains containing holes (i.e., excised regions), via the multigrid method. We consider as a test case the Poisson equation with a nonlinear term added, as a means of illustrating the principles involved, and move to a "real world" 3-dimensional problem which is the solution of the conformally flat Hamiltonian constraint with Dirichlet and Robin boundary conditions. Using our vertex-centered multigrid code, we demonstrate globally second-order-accurate solutions of elliptic equations over domains containing holes, in two and three spatial dimensions. Keys to the success of this method are the choice of the restriction operator near the holes and definition of the location of the inner boundary. In some cases (e.g. two holes in two dimensions), more and more smoothing may be required as the mesh spacing decreases to zero; however for the resolutions currently of interest to many numerical relativists, it is feasible to maintain second order convergence by concentrating smoothing (spatially) where it is needed most. This paper, and our publicly available source code, are intended to serve as semi-pedagogical guides for those who may wish to implement similar schemes.Comment: 18 pages, 11 figures, LaTeX. Added clarifications and references re. scope of paper, mathematical foundations, relevance of work. Accepted for publication in Classical & Quantum Gravit

Pharmacological Alterations of Anxious Behaviour in Mice Depending on Both Strain and the Behavioural Situation

A previous study comparing non-emotive mice from the strain C57BL/6/ByJ with ABP/Le mice showed ABP/Le to be more anxious in an open-field situation. In the present study, several compounds affecting anxiety were assayed on ABP/Le and C57BL/6/ByJ mice using three behavioural models of anxiety: the elevated plus-maze, the light-dark discrimination test and the free exploratory paradigm. The compounds used were the full benzodiazepine receptor agonist, chlordiazepoxide, and the antagonist, flumazenil, the GABAA antagonist, bicuculline, the full 5-HT1A agonist 8-OH-DPAT, and the mixed 5-HT1A/5-HT1B agonist, RU 24969. Results showed the effect of the compounds to be dependent on both the strain and the behavioural task. Several compounds found to be anxiolytic in ABP/Le mice had an anxiogenic effect on C57BL/6/ByJ mice. More behavioural changes were observed for ABP/Le in the elevated plus-maze, but the clearest findings for C57BL/6/ByJ mice were observed in the light-dark discrimination apparatus. These data demonstrate that anxious behaviour is a complex phenomenon which cannot be described by a single behavioural task nor by the action of a single compound

Towards a unified theory of Sobolev inequalities

We discuss our work on pointwise inequalities for the gradient which are connected with the isoperimetric profile associated to a given geometry. We show how they can be used to unify certain aspects of the theory of Sobolev inequalities. In particular, we discuss our recent papers on fractional order inequalities, Coulhon type inequalities, transference and dimensionless inequalities and our forthcoming work on sharp higher order Sobolev inequalities that can be obtained by iteration.Comment: 39 pages, made some changes to section 1