A selfconsistent theory of current-induced switching of magnetization

A selfconsistent theory of the current-induced switching of magnetization using nonequilibrium Keldysh formalism is developed for a junction of two ferromagnets separated by a nonmagnetic spacer. It is shown that the spin-transfer torques responsible for current-induced switching of magnetization can be calculated from first principles in a steady state when the magnetization of the switching magnet is stationary. The spin-transfer torque is expressed in terms of one-electron surface Green functions for the junction cut into two independent parts by a cleavage plane immediately to the left and right of the switching magnet. The surface Green functions are calculated using a tight-binding Hamiltonian with parameters determined from a fit to an {\it ab initio} band structure.This treatment yields the spin transfer torques taking into account rigorously contributions from all the parts of the junction. To calculate the hysteresis loops of resistance versus current, and hence to determine the critical current for switching, the microscopically calculated spin-transfer torques are used as an input into the phenomenological Landau-Lifshitz equation with Gilbert damping. The present calculations for Co/Cu/Co(111) show that the critical current for switching is $\approx 10^7A/cm^2$, which is in good agreement with experiment.Comment: 23 pages, 16 figure

Statistical Mechanics of Vibration-Induced Compaction of Powders

We propose a theory which describes the density relaxation of loosely packed, cohesionless granular material under mechanical tapping. Using the compactivity concept we develope a formalism of statistical mechanics which allows us to calculate the density of a powder as a function of time and compactivity. A simple fluctuation-dissipation relation which relates compactivity to the amplitude and frequency of a tapping is proposed. Experimental data of E.R.Nowak et al. [{\it Powder Technology} 94, 79 (1997) ] show how density of initially deposited in a fluffy state powder evolves under carefully controlled tapping towards a random close packing (RCP) density. Ramping the vibration amplitude repeatedly up and back down again reveals the existence of reversible and irreversible branches in the response. In the framework of our approach the reversible branch (along which the RCP density is obtained) corresponds to the steady state solution of the Fokker-Planck equation whereas the irreversible one is represented by a superposition of "excited states" eigenfunctions. These two regimes of response are analyzed theoretically and a qualitative explanation of the hysteresis curve is offered.Comment: 11 pages, 2 figures, Latex. Revised tex

Midgap states and charge inhomogeneities in corrugated graphene

We study the changes induced by the effective gauge field due to ripples on the low energy electronic structure of graphene. We show that zero energy Landau levels will form, associated to the smooth deformation of the graphene layer, when the height corrugation, $h$, and the length of the ripple, $l$, are such that $h^2 / (l a) \gtrsim 1$, where $a$ is the lattice constant. The existence of localized levels gives rise to a large compressibility at zero energy, and to the enhancement of instabilities arising from electron-electron interactions including electronic phase separation. The combined effect of the ripples and an external magnetic field breaks the valley symmetry of graphene leading to the possibility of valley selection

The Data Processing Pipeline for the Herschel-HIFI Instrument

The HIFI data processing pipeline was developed to systematically process diagnostic, calibration and astronomical observations taken with the HIFI science instrumentas part of the Herschel mission. The HIFI pipeline processed data from all HIFI observing modes within the Herschel automated processing environment, as well as, within an interactive environment. A common software framework was developed to best support the use cases required by the instrument teams and by the general astronomers. The HIFI pipeline was built on top of that and was designed with a high degree of modularity. This modular design provided the necessary flexibility and extensibility to deal with the complexity of batch-processing eighteen different observing modes, to support the astronomers in the interactive analysis and to cope with adjustments necessary to improve the pipeline and the quality of the end-products. This approach to the software development and data processing effort was arrived at by coalescing the lessons learned from similar research based projects with the understanding that a degree of foresight was required given the overall length of the project. In this article, both the successes and challenges of the HIFI software development process are presented. To support future similar projects and retain experience gained lessons learned are extracted.Comment: 18 pages, 5 figure

Mathematics of random growing interfaces

We establish a thermodynamic limit and Gaussian fluctuations for the height and surface width of the random interface formed by the deposition of particles on surfaces. The results hold for the standard ballistic deposition model as well as the surface relaxation model in the off-lattice setting. The results are proved with the aid of general limit theorems for stabilizing functionals of marked Poisson point processes.Comment: 12 page

Statistical comparison of pooled nitrogen washout data of various altitude decompression response groups

This analysis was done to determine whether various decompression response groups could be characterized by the pooled nitrogen (N2) washout profiles of the group members, pooling individual washout profiles provided a smooth time dependent function of means representative of the decompression response group. No statistically significant differences were detected. The statistical comparisons of the profiles were performed by means of univariate weighted t-test at each 5 minute profile point, and with levels of significance of 5 and 10 percent. The estimated powers of the tests (i.e., probabilities) to detect the observed differences in the pooled profiles were of the order of 8 to 30 percent

A mean field description of jamming in non-cohesive frictionless particulate systems

A theory for kinetic arrest in isotropic systems of repulsive, radially-interacting particles is presented that predicts exponents for the scaling of various macroscopic quantities near the rigidity transition that are in agreement with simulations, including the non-trivial shear exponent. Both statics and dynamics are treated in a simplified, one-particle level description, and coupled via the assumption that kinetic arrest occurs on the boundary between mechanically stable and unstable regions of the static parameter diagram. This suggests the arrested states observed in simulations are at (or near) an elastic buckling transition. Some additional numerical evidence to confirm the scaling of microscopic quantities is also provided.Comment: 9 pages, 3 figs; additional clarification of different elastic moduli exponents, plus typo fix. To appear in PR

Examples of mathematical modeling tales from the crypt

Mathematical modeling is being increasingly recognized within the biomedical sciences as an important tool that can aid the understanding of biological systems. The heavily regulated cell renewal cycle in the colonic crypt provides a good example of how modeling can be used to find out key features of the system kinetics, and help to explain both the breakdown of homeostasis and the initiation of tumorigenesis. We use the cell population model by Johnston et al. (2007) Proc. Natl. Acad. Sci. USA 104, 4008-4013, to illustrate the power of mathematical modeling by considering two key questions about the cell population dynamics in the colonic crypt. We ask: how can a model describe both homeostasis and unregulated growth in tumorigenesis; and to which parameters in the system is the model most sensitive? In order to address these questions, we discuss what type of modeling approach is most appropriate in the crypt. We use the model to argue why tumorigenesis is observed to occur in stages with long lag phases between periods of rapid growth, and we identify the key parameters

On the proportion of cancer stem cells in a tumour

It is now generally accepted that cancers contain a sub-population, the cancer stem cells (CSCs), which initiate and drive a tumour’s growth. At least until recently it has been widely assumed that only a small proportion of the cells in a tumour are CSCs. Here we use a mathematical model, supported by experimental evidence, to show that such an assumption is unwarranted. We show that CSCs may comprise any possible proportion of the tumour, and that the higher the proportion the more aggressive the tumour is likely to be