The Excursion Set Theory of Halo Mass Functions, Halo Clustering, and Halo Growth

I review the excursion set theory (EST) of dark matter halo formation and clustering. I recount the Press-Schechter argument for the mass function of bound objects and review the derivation of the Press-Schechter mass function in EST. The EST formalism is powerful and can be applied to numerous problems. I review the EST of halo bias and the properties of void regions. I spend considerable time reviewing halo growth in the EST. This section culminates with descriptions of two Monte Carlo methods for generating halo mass accretion histories. In the final section, I emphasize that the standard EST approach is the result of several simplifying assumptions. Dropping these assumptions can lead to more faithful predictions and a more versatile formalism. One such assumption is the constant height of the barrier for nonlinear collapse. I review implementations of the excursion set approach with arbitrary barrier shapes. An application of this is the now well-known improvement to standard EST that follows from the ellipsoidal-collapse barrier. Additionally, I emphasize that the statement that halo accretion histories are independent of halo environments is a simplifying assumption, rather than a prediction of the theory. I review the method for constructing correlated random walks of the density field in more general cases. I construct a simple toy model with correlated walks and I show that excursion set theory makes a qualitatively simple and general prediction for the relation between halo accretion histories and halo environments: regions of high density preferentially contain late-forming halos and conversely for regions of low density. I conclude with a brief discussion of this prediction in the context of recent numerical studies of the environmental dependence of halo properties. (Abridged)Comment: 62 pages, 19 figures. Review article based on lectures given at the Sixth Summer School of the Helmholtz Institute for Supercomputational Physics. Accepted for Publication in IJMPD. Comments Welcom

Obliquity histories of Earth and Mars: Influence of inertial and dissipative core-mantle coupling

For both the Earth and Mars, secular variations in the angular separation of the spin axis from the orbit normal are suspected of driving major climatic changes. There is considerable interest in determining the amplitude and timing of these obliquity variations. If the orientation of the orbital plane were inertially fixed, and the planet were to act as a rigid body in it response to precessional torques, the spin axis would simply precess around the orbit at a fixed obliquity and at a uniform angular rate. The precession rate parameter depends on the principal moments of inertia and rotation rate of the perturbed body, and on the gravitational masses and semiminor axes of the perturbing bodies. For Mars, the precession rate is not well known, but probably lies in the interval 8 to 10 arcsec/year. Gravitational interactions between the planets lead to secular motions of the orbit planes. In the rigid body case, the spin axis still attempts to precess about the instantaneous orbit normal, but now the obliquity varies. The hydrostatic figure of a planet represents a compromise between gravitation, which attempts to attain spherical symmetry, and rotation, which prefers cylindrical symmetry. Due to their higher mean densities the cores of the Earth and Mars will be more nearly spherical than the outer layers of these planets. On short time scales it is appropriate to consider the core to be an inviscid fluid constrained to move with the ellipsoidal region bounded by the rigid mantle. The inertial coupling provided by this mechanism is effective whenever the ellipticicy of the container exceeds the ratio of precessional to rotational rates. If the mantle were actually rigid, this would be an extremely effective type of coupling. However, on sufficiently long time scales, the mantle will deform viscously and can accommodate the motions of the core fluid. A fundamentally different type of coupling is provided by electromagnetic or viscous torques. This type of coupling is likely to be most important on longer time scales. In each case, the mantle exerts an equal and opposite torque on the core

An efficient Jeans modelling of axisymmetric galaxies with multiple stellar components

Dynamical models of stellar systems represent a powerful tool to study their internal structure and dynamics, to interpret the observed morphological and kinematical fields, and also to support numerical simulations of their evolution. We present a method especially designed to build axisymmetric Jeans models of galaxies, assumed as stationary and collisionless stellar systems. The aim is the development of a rigorous and flexible modelling procedure of multicomponent galaxies, composed of different stellar and dark matter distributions, and a central supermassive black hole. The stellar components, in particular, are intended to represent different galaxy structures, such as discs, bulges, halos, and can then have different structural (density profile, flattening, mass, scale-length), dynamical (rotation, velocity dispersion anisotropy), and population (age, metallicity, initial mass function, mass-to-light ratio) properties. The theoretical framework supporting the modelling procedure is presented, with the introduction of a suitable nomenclature, and its numerical implementation is discussed, with particular reference to the numerical code JASMINE2, developed for this purpose. We propose an approach for efficiently scaling the contributions in mass, luminosity, and rotational support, of the different matter components, allowing for fast and flexible explorations of the model parameter space. We also offer different methods of the computation of the gravitational potentials associated of the density components, especially convenient for their easier numerical tractability. A few galaxy models are studied, showing internal, and projected, structural and dynamical properties of multicomponent galaxies, with a focus on axisymmetric early-type galaxies with complex kinematical morphologies. The application of galaxy models to the study of initial conditions for hydro-dynamical and $N$-body simulations of galaxy evolution is also addressed, allowing in particular to investigate the large number of interesting combinations of the parameters which determine the structure and dynamics of complex multicomponent stellar systems

Visualization and analysis of diffusion tensor fields

technical reportThe power of medical imaging modalities to measure and characterize biological tissue is amplified by visualization and analysis methods that help researchers to see and understand the structures within their data. Diffusion tensor magnetic resonance imaging can measure microstructural properties of biological tissue, such as the coherent linear organization of white matter of the central nervous system, or the fibrous texture of muscle tissue. This dissertation describes new methods for visualizing and analyzing the salient structure of diffusion tensor datasets. Glyphs from superquadric surfaces and textures from reactiondiffusion systems facilitate inspection of data properties and trends. Fiber tractography based on vector-tensor multiplication allows major white matter pathways to be visualized. The generalization of direct volume rendering to tensor data allows large-scale structures to be shaded and rendered. Finally, a mathematical framework for analyzing the derivatives of tensor values, in terms of shape and orientation change, enables analytical shading in volume renderings, and a method of feature detection important for feature-preserving filtering of tensor fields. Together, the combination of methods enhances the ability of diffusion tensor imaging to provide insight into the local and global structure of biological tissue

Self-organization of polymers in bulk and at interfaces

Fully atomistic analysis of polymeric systems is computationally very demanding because the time and length scales involved span over several orders of magnitude. At the same time many properties of polymers are universal in the sense that they do not depend on the chemical nature of the comprising monomers. This makes coarse-grained methods, such as self-consistent field (SCF) modeling, an ideal tool for studying them. In this thesis we employ SCF modeling to study intra- and intermolecular self-organization organization of polymers and ordering of polymers near interfaces. Where possible, the results are compared to experiments and predictions of analytical theories

Spectral Analysis and Parameter Estimation in Fibre Levitated Optomechanics

In levitated optomechanics, nano-scale objects are optically trapped so that their motion can be studied. These trapped nanoparticles are held in a 3D quadratic potential and act as damped harmonic oscillators; they are thermally and mechanically decoupled from the apparatus and their position is measured interfer-ometrically to picometre accuracy. These systems are well suited to sensing and metrology applications, as any external disturbance of the particle can be observed using the scattered trapping light.When examining the motion of a levitated nanoparticle, it’s position is recorded and used to estimate a power spectral density (PSD), from which state parameters can be estimated. In this thesis an experi-mental setup is presented, optimised for maximum collection of particle position information in 1D, using a fibre-based parabolic mirror trap and heterodyne measurement system in order to produce spectra with minimal noise and unwanted artefacts.A novel application of the Middleton expansion from RF engineering is used to generate a complete power spectrum that depends on the physical parameters of the system. This method treats the particle as a stochastic harmonic oscillator, phase modulated by a Gaussian random process with known PSD. We reproduce the PSD of intensity at a detector, a quantity that is sinusoidally dependent on particle posi-tion. This technique generates a single, full PSD using modified Bessel functions, and does not depend on assumptions about the relative phases of the interfered fields, highlighting the non-linear dependence of measured signal on position. Theoretical spectra are fitted to a measured PSD and the phase modulation depth is extracted; this is used to calculate the particle oscillation amplitude and, by an equipartition ar-gument, the centre of mass temperature to mass ratio. State parameters are tracked as environmental conditions change and an increase in centre of mass temperature as a function of decreasing background gas pressure is observed

Whisker Movements Reveal Spatial Attention: A Unified Computational Model of Active Sensing Control in the Rat

Spatial attention is most often investigated in the visual modality through measurement of eye movements, with primates, including humans, a widely-studied model. Its study in laboratory rodents, such as mice and rats, requires different techniques, owing to the lack of a visual fovea and the particular ethological relevance of orienting movements of the snout and the whiskers in these animals. In recent years, several reliable relationships have been observed between environmental and behavioural variables and movements of the whiskers, but the function of these responses, as well as how they integrate, remains unclear. Here, we propose a unifying abstract model of whisker movement control that has as its key variable the region of space that is the animal's current focus of attention, and demonstrate, using computer-simulated behavioral experiments, that the model is consistent with a broad range of experimental observations. A core hypothesis is that the rat explicitly decodes the location in space of whisker contacts and that this representation is used to regulate whisker drive signals. This proposition stands in contrast to earlier proposals that the modulation of whisker movement during exploration is mediated primarily by reflex loops. We go on to argue that the superior colliculus is a candidate neural substrate for the siting of a head-centred map guiding whisker movement, in analogy to current models of visual attention. The proposed model has the potential to offer a more complete understanding of whisker control as well as to highlight the potential of the rodent and its whiskers as a tool for the study of mammalian attention

Scientific visualization of stress tensor information with applications to stress determination by X-ray and neutron diffraction

Includes bibliographical references (leaves 232-249).The visual analysis of mechanical stress facilitates physical understanding of the tensor quantity which is concealed in scalar and vector methods. In this study, the principles and techniques of scientific visualization are used to develop a visual analysis of mechanical stresses. Scientific visualization is not only applied to the final tensorial quantity obtained from the diffraction measurements, but the visual methods are developed from, and integrated into current residual stress analysis practices by relating the newly developed visual techniques to the conventional techniques, highlighting its advantages. This study consists of the mathematical analysis of the tensor character of mechanical stresses, discussion of the principles and techniques of scientific visualization (visual data analysis) in physical research, and tensor determination, visual analysis and presentation of residual stresses obtained from diffraction measurements