Gravitational Lensing as a Tool on Galactic and Cosmological Scales

This work considers the reconstruction of strong gravitational lenses from their observed effects on the light distribution of background sources. After reviewing the formalism of gravitational lensing and the most common and relevant lens models, new analytical results on the elliptical power law lens are presented, including new expressions for the deflection, potential, shear and magnification, which naturally lead to a fast numerical scheme for practical calculation. The main part of the thesis investigates lens reconstruction with extended sources by means of the forward reconstruction method, in which the lenses and sources are given by parametric models. The numerical realities of the problem make it necessary to find targeted optimisations for the forward method, in order to make it feasible for general applications to modern, high resolution images. The result of these optimisations is presented in the \textsc{Lensed} algorithm. Subsequently, a number of tests for general forward reconstruction methods are created to decouple the influence of sourced from lens reconstructions, in order to objectively demonstrate the constraining power of the reconstruction. The final chapters on lens reconstruction contain two sample applications of the forward method. One is the analysis of images from a strong lensing survey. Such surveys today contain $\sim 100$ strong lenses, and much larger sample sizes are expected in the future, making it necessary to quickly and reliably analyse catalogues of lenses with a fixed model. The second application deals with the opposite situation of a single observation that is to be confronted with different lens models, where the forward method allows for natural model-building. This is demonstrated using an example reconstruction of the ``Cosmic Horseshoe''. An appendix presents an independent work on the use of weak gravitational lensing to investigate theories of modified gravity which exhibit screening in the non-linear regime of structure formation.

Properties of minimally doubled fermions

Most quark actions in lattice QCD encounter difficulties with chiral sym-rnmetry and its spontaneous breakdown. Minimally doubled fermions (MDF)rnare a category of strictly local chiral lattice fermions, whose continuum limitrnreproduces two degenerate quark flavours. The two poles of their Dirac ope-rnrator are aligned such that symmetries under charge conjugation or reflectionrnof one particular direction are explictly broken at finite lattice spacing. Pro-rnperties of MDF are scrutinised with regard to broken symmetry and mesonrnspectrum to discern their suitability for numerical studies of QCD.rnrnInteractions induce anisotropic operator mixing for MDF. Hence, resto-rnration of broken symmetries in the continuum limit requires three coun-rnterterms, one of which is power-law divergent. Counterterms and operatorrnmixing are studied perturbatively for two variants of MDF. Two indepen-rndent non-perturbative procedures for removal of the power-law divergencernare developed by means of a numerical study of hadronic observables forrnone variant of MDF in quenched approximation. Though three out of fourrnpseudoscalar mesons are affected by lattice artefacts, the spectrum’s conti-rnnuum limit is consistent with two-flavour QCD. Thus, suitability of MDF forrnnumerical studies of QCD in the quenched approximation is demonstrated.

Seasonality and diel variation in blue whale D calls recorded in the Southern California Bight

Passive acoustic data have been collected using HARPs (High-frequency Acoustic Recording Packages) and were used to assess (1) the seasonality of blue whale D calls in the Southern California Bight, (2) their interannual abundance during 2007-2012 and (3) their diel variation. This goal has been achieved running the GPL (Generalized Power-Law) automated detector. (1) Blue whale D calls were detected in the Southern California Bight from May through November with a peak in July, even though few detections were from December to April as well. A key predictor for blue whale distribution and movement in the California Current region has been identified with zooplankton aggregations, paying a particular attention to those euphausiid species, such as E. pacifica and T. spinifera, which are blue whale favorite krill. The Southern California Bight experiences seasonal upwelling, resulting in an increase of productivity and prey availability. The summer and early fall have been marked as the most favorable periods. This supports the presence of blue whales in the area at that time, supposing these marine mammals exploit the region as a feeding ground. (2) As to the interannual abundance during 2007-2012, I found a large variability. I observed a great increase of vocalizations in 2007 and 2010, whereas a decrease was shown in the other years, which is well marked in 2009. It is my belief that these fluctuations in abundance of D calls detections through the deployed period are due to the alternation of El Nino and La Nina events, which occurred in those years. (3) The assessment of the daily timing of D calls production shows that D calls are more abundant during the day than during the night with a peak at 12:00 and 13:00. Assuming that D calling is associated with feeding, the daily pattern of D calls may be linked to the prey availability. E. pacifica and T. spinifera are among those species of krill which undertake daily vertical migrations, remaining at depth during the day and slowly coming up towards the surface at night. Because of some anatomical arrangements, these euphausiids are very sensitive to the light. Given that we believe D calls have a social function, I hypothesize that blue whales may recognize the hours at the highest solar incidence as the best moment of the day in terms of prey availability, exploiting this time window to advert their conspecifics.

Analisi e confronto di sequenze di DNA mediante modelli Markoviani

Lo scopo di questa tesi è quello di evidenziare, attraverso varie analisi statistiche ed applicazione di modelli stocastici, il comportamento strutturale e funzionale dei dinucleotidi che compongono le sequenze di DNA di diversi organismi. Gli organismi che abbiamo scelto di prendere in considerazione sono l'uomo, il topo e l'Escherichia coli. Questa scelta non è stata casuale, ma oculata, al fine di mettere in risalto alcune differenze tra organismi eucarioti, quali l'uomo e il topo, ed organismi procarioti come il batterio E.coli. Nella prima parte del nostro studio, abbiamo computato le distanze che intercorrono tra occorrenze successive dello stesso dinucleotide lungo la sequenza, usando un metodo di non sovrapposizione, ed abbiamo iterato il calcolo per tutti i 16 dinucleotidi. Dopodiché ci siamo preoccupati di graficare le distribuzioni di distanza dei 16 dinucleotidi per l'E.Coli, il topo e l'uomo; gli istogrammi evidenziano un comportamento anomalo della distribuzione di CG che accomuna gli organismi eucarioti e di cui, invece, è esente l'organismo procariote esaminato. Questo dato statistico trova una spiegazione nei processi biologici di metilazione che possono innescarsi sul dinucleotide CG nelle sequenze eucariotiche. In seguito, per determinare quanto ciascuna delle 16 distribuzioni si discosti dalle altre abbiamo usato la divergenza di Jensen-Shannon. Per quantificare le differenze sostanziali tra le distribuzioni di CG dei 3 organismi considerati abbiamo deciso di verificare quale fosse il miglior fit per tali curve tra un esponenziale ed una power-law. L'esponenziale rappresenta un buon fit per le code delle distribuzioni di CG del topo e dell'uomo; ciò rivela la presenza di una lunghezza caratteristica per entrambi gli organismi. Nella seconda parte dello studio, i risultati vengono confrontati con modelli markoviani: sequenze random generate con catene di Markov di ordine zero (basate sulle frequenze relative dei nucleotidi) e uno (basate sulle probabilità di transizione tra diversi nucleotidi). Quest'ultima riproduce abbastanza fedelmente la sequenza biologica di partenza, per cui abbiamo scelto di utilizzare la catena Markov del 1° ordine per altre analisi statistiche riguardanti le distribuzioni dei nucleotidi, dinucleotidi, ed anche dei trinucleotidi con particolare interesse per quelli in cui è contenuto CG, in modo da verificare se l'anomalia si ripercuote anche in essi. Riteniamo pertanto che metodi basati su questo approccio potrebbero essere sfruttati per confermare le peculiarità biologiche e per migliorare l'individuazione delle aree di interesse, come le isole CpG, ed eventualmente promotori e Lamina Associated Domains (LAD), nel genoma di diversi organismi.

Comportamenti di fruizione su piattaforme streaming: il caso Twitch

Internet e il Web rappresentano una risorsa informativa immensa, utilizzabile anche per compiere indagini di mercato, soprattutto a fronte di nuovi fenomeni di consumo. Un caso emergente è rappresentato dai servizi di video streaming. La piattaforma presa in considerazione dal mio studio di analisi è Twitch, leader nel segmento di mercato dei video-giochi. Attraverso di essa milioni di follower seguono in tempo reali i comportamenti di gioco di migliaia di player, fra i quali i principali campioni mondiali di e-sports. La tesi ha esaminato, attraverso lo studio di fonti di letteratura e di mercato, le principali caratteristiche di relazione degli utenti. Su tale base è stato definito un campione attraverso uso delle API del sito, oggetto di studio finalizzato alla comprensione dei comportamenti di scelta da parte dei follower. Fra i vari esiti, emerge confermata la presenza e tutti i livelli di analisi della distribuzione power-law, tipica anche di Internet nel suoi complesso.

Strong Dynamical Heterogeneity and Universal Scaling in Driven Granular Fluids

Large-scale simulations of two-dimensional bidisperse granular fluids allow us to determine spatial correlations of slow particles via the four-point structure factor S-4 (q, t). Both cases, elastic (epsilon = 1) and inelastic (epsilon < 1) collisions, are studied. As the fluid approaches structural arrest, i.e., for packing fractions in the range 0.6 <= phi <= 0.805, scaling is shown to hold: S-4 (q, t)/chi(4)(t) = s(q xi(t)). Both the dynamic susceptibility chi(4)(tau(alpha)) and the dynamic correlation length xi(tau(alpha)) evaluated at the alpha relaxation time tau(alpha) can be fitted to a power law divergence at a critical packing fraction. The measured xi(tau(alpha)) widely exceeds the largest one previously observed for three-dimensional (3d) hard sphere fluids. The number of particles in a slow cluster and the correlation length are related by a robust power law, chi(4)(tau(alpha)) approximate to xi(d-p) (tau(alpha)), with an exponent d - p approximate to 1.6. This scaling is remarkably independent of epsilon, even though the strength of the dynamical heterogeneity at constant volume fraction depends strongly on epsilon.

Multiple time scales of temporal response in pyramidal and fast spiking cortical neurons

Neural dynamic processes correlated over several time scales are found in vivo, in stimulus-evoked as well as spontaneous activity, and are thought to affect the way sensory stimulation is processed. Despite their potential computational consequences, a systematic description of the presence of multiple time scales in single cortical neurons is lacking. In this study, we injected fast spiking and pyramidal (PYR) neurons in vitro with long-lasting episodes of step-like and noisy, in-vivo-like current. Several processes shaped the time course of the instantaneous spike frequency, which could be reduced to a small number (1-4) of phenomenological mechanisms, either reducing (adapting) or increasing (facilitating) the neuron's firing rate over time. The different adaptation/facilitation processes cover a wide range of time scales, ranging from initial adaptation (<10 ms, PYR neurons only), to fast adaptation (<300 ms), early facilitation (0.5-1 s, PYR only), and slow (or late) adaptation (order of seconds). These processes are characterized by broad distributions of their magnitudes and time constants across cells, showing that multiple time scales are at play in cortical neurons, even in response to stationary stimuli and in the presence of input fluctuations. These processes might be part of a cascade of processes responsible for the power-law behavior of adaptation observed in several preparations, and may have far-reaching computational consequences that have been recently described.

Labeling the human skeleton with 41Ca to assess changes in bone calcium metabolism

Bone research is limited by the methods available for detecting changes in bone metabolism. While dual X-ray absorptiometry is rather insensitive, biochemical markers are subject to significant intra-individual variation. In the study presented here, we evaluated the isotopic labeling of bone using 41Ca, a long-lived radiotracer, as an alternative approach. After successful labeling of the skeleton, changes in the systematics of urinary 41Ca excretion are expected to directly reflect changes in bone Ca metabolism. A minute amount of 41Ca (100 nCi) was administered orally to 22 postmenopausal women. Kinetics of tracer excretion were assessed by monitoring changes in urinary 41Ca/40Ca isotope ratios up to 700 days post-dosing using accelerator mass spectrometry and resonance ionization mass spectrometry. Isotopic labeling of the skeleton was evaluated by two different approaches: (i) urinary 41Ca data were fitted to an established function consisting of an exponential term and a power law term for each individual; (ii) 41Ca data were analyzed by population pharmacokinetic (NONMEM) analysis to identify a compartmental model that describes urinary 41Ca tracer kinetics. A linear three-compartment model with a central compartment and two sequential peripheral compartments was found to best fit the 41Ca data. Fits based on the use of the combined exponential/power law function describing urinary tracer excretion showed substantially higher deviations between predicted and measured values than fits based on the compartmental modeling approach. By establishing the urinary 41Ca excretion pattern using data points up to day 500 and extrapolating these curves up to day 700, it was found that the calculated 41Ca/40Ca isotope ratios in urine were significantly lower than the observed 41Ca/40Ca isotope ratios for both techniques. Compartmental analysis can overcome this limitation. By identifying relative changes in transfer rates between compartments in response to an intervention, inaccuracies in the underlying model cancel out. Changes in tracer distribution between compartments were modeled based on identified kinetic parameters. While changes in bone formation and resorption can, in principle, be assessed by monitoring urinary 41Ca excretion over the first few weeks post-dosing, assessment of an intervention effect is more reliable approximately 150 days post-dosing when excreted tracer originates mainly from bone.

Mechanisms of intrinsic beating variability in cardiac cell cultures and model pacemaker networks

Heart rate variability (HRV) exhibits fluctuations characterized by a power law behavior of its power spectrum. The interpretation of this nonlinear HRV behavior, resulting from interactions between extracardiac regulatory mechanisms, could be clinically useful. However, the involvement of intrinsic variations of pacemaker rate in HRV has scarcely been investigated. We examined beating variability in spontaneously active incubating cultures of neonatal rat ventricular myocytes using microelectrode arrays. In networks of mathematical model pacemaker cells, we evaluated the variability induced by the stochastic gating of transmembrane currents and of calcium release channels and by the dynamic turnover of ion channels. In the cultures, spontaneous activity originated from a mobile focus. Both the beat-to-beat movement of the focus and beat rate variability exhibited a power law behavior. In the model networks, stochastic fluctuations in transmembrane currents and stochastic gating of calcium release channels did not reproduce the spatiotemporal patterns observed in vitro. In contrast, long-term correlations produced by the turnover of ion channels induced variability patterns with a power law behavior similar to those observed experimentally. Therefore, phenomena leading to long-term correlated variations in pacemaker cellular function may, in conjunction with extracardiac regulatory mechanisms, contribute to the nonlinear characteristics of HRV.

Measuring energy spectra of TeV gamma-ray emission from the Cygnus region of our galaxy with Milagro

High energy gamma rays can provide fundamental clues to the origins of cosmic rays. In this thesis, TeV gamma-ray emission from the Cygnus region is studied. Previously the Milagro experiment detected five TeV gamma-ray sources in this region and a significant excess of TeV gamma rays whose origin is still unclear. To better understand the diffuse excess the separation of sources and diffuse emission is studied using the latest and most sensitive data set of the Milagro experiment. In addition, a newly developed technique is applied that allows the energy spectrum of the TeV gamma rays to be reconstructed using Milagro data. No conclusive statement can be made about the spectrum of the diffuse emission from the Cygnus region because of its low significance of 2.2 σ above the background in the studied data sample. The entire Cygnus region emission is best fit with a power law with a spectral index of α=2.40 (68% confidence interval: 1.35-2.92) and a exponential cutoff energy of 31.6 TeV (10.0-251.2 TeV). In the case of a simple power law assumption without a cutoff energy the best fit yields a spectral index of α=2.97 (68% confidence interval: 2.83-3.10). Neither of these best fits are in good agreement with the data. The best spectral fit to the TeV emission from MGRO J2019+37, the brightest source in the Cygnus region, yields a spectral index of α=2.30 (68% confidence interval: 1.40-2.70) with a cutoff energy of 50.1 TeV (68% confidence interval: 17.8-251.2 TeV) and a spectral index of α=2.75 (68% confidence interval: 2.65-2.85) when no exponential cutoff energy is assumed. According to the present analysis, MGRO J2019+37 contributes 25% to the differential flux from the entire Cygnus at 15 TeV.

Dynamics of settling charged particles in turbulence : theory and experiments

It has been proposed that inertial clustering may lead to an increased collision rate of water droplets in clouds. Atmospheric clouds and electrosprays contain electrically charged particles embedded in turbulent flows, often under the influence of an externally imposed, approximately uniform gravitational or electric force. In this thesis, we present the investigation of charged inertial particles embedded in turbulence. We have developed a theoretical description for the dynamics of such systems of charged, sedimenting particles in turbulence, allowing radial distribution functions to be predicted for both monodisperse and bidisperse particle size distributions. The governing parameters are the particle Stokes number (particle inertial time scale relative to turbulence dissipation time scale), the Coulomb-turbulence parameter (ratio of Coulomb ’terminalar speed to turbulence dissipation velocity scale), and the settling parameter (the ratio of the gravitational terminal speed to turbulence dissipation velocity scale). For the monodispersion particles, The peak in the radial distribution function is well predicted by the balance between the particle terminal velocity under Coulomb repulsion and a time-averaged ’drift’ velocity obtained from the nonuniform sampling of fluid strain and rotation due to finite particle inertia. The theory is compared to measured radial distribution functions for water particles in homogeneous, isotropic air turbulence. The radial distribution functions are obtained from particle positions measured in three dimensions using digital holography. The measurements support the general theoretical expression, consisting of a power law increase in particle clustering due to particle response to dissipative turbulent eddies, modulated by an exponential electrostatic interaction term. Both terms are modified as a result of the gravitational diffusion-like term, and the role of ’gravity’ is explored by imposing a macroscopic uniform electric field to create an enhanced, effective gravity. The relation between the radial distribution functions and inward mean radial relative velocity is established for charged particles.

Lunar energetic neutral atom (ENA) spectra measured by the interstellar boundary explorer (IBEX)

The solar wind continuously flows out from the Sun, filling interplanetary space and directly interacting with the surfaces of small planetary bodies and other objects throughout the solar system. A significant fraction of these ions backscatter from the surface as energetic neutral atoms (ENAs). The first observations of these ENA emissions from the Moon were recently reported from the Interstellar Boundary Explorer (IBEX). These observations yielded a lunar ENA albedo of ˜10% and showed that the Moon reflects ˜150 metric tons of neutral hydrogen per year. More recently, a survey of the first 2.5 years of IBEX observations of lunar ENAs was conducted for times when the Moon was in the solar wind. Here, we present the first IBEX ENA observations when the Moon is inside the terrestrial magnetosheath and compare them with observations when the Moon is in the solar wind. Our analysis shows that: (1) the ENA intensities are on average higher when the Moon is in the magnetosheath, (2) the energy spectra are similar above ~0.6* solar wind energy but below there are large differences of the order of a factor of 10, (3) the energy spectra resemble a power law with a "hump" at ˜0.6 * solar wind energy, and (4) this "hump" is broader when the Moon is in the magnetosheath. We explore potential scenarios to explain the differences, namely the effects of the topography of the lunar surface and the consequences of a very different Mach number in the solar wind versus in the magnetosheath.

Multi-scale micro-structure generation strategy for up-scaling transport in clays

Clays and claystones are used as backfill and barrier materials in the design of waste repositories, because they act as hydraulic barriers and retain contaminants. Transport through such barriers occurs mainly by molecular diffusion. There is thus an interest to relate the diffusion properties of clays to their structural properties. In previous work, we have developed a concept for up-scaling pore-scale molecular diffusion coefficients using a grid-based model for the sample pore structure. Here we present an operational algorithm which can generate such model pore structures of polymineral materials. The obtained pore maps match the rock’s mineralogical components and its macroscopic properties such as porosity, grain and pore size distributions. Representative ensembles of grains in 2D or 3D are created by a lattice Monte Carlo (MC) method, which minimizes the interfacial energy of grains starting from an initial grain distribution. Pores are generated at grain boundaries and/or within grains. The method is general and allows to generate anisotropic structures with grains of approximately predetermined shapes, or with mixtures of different grain types. A specific focus of this study was on the simulation of clay-like materials. The generated clay pore maps were then used to derive upscaled effective diffusion coefficients for non-sorbing tracers using a homogenization technique. The large number of generated maps allowed to check the relations between micro-structural features of clays and their effective transport parameters, as is required to explain and extrapolate experimental diffusion results. As examples, we present a set of 2D and 3D simulations and investigated the effects of nanopores within particles (interlayer pores) and micropores between particles. Archie’s simple power law is followed in systems with only micropores. When nanopores are present, additional parameters are required; the data reveal that effective diffusion coefficients could be described by a sum of two power functions, related to the micro- and nanoporosity. We further used the model to investigate the relationships between particle orientation and effective transport properties of the sample.

How common is the lipid body-containing interstitial cell in the mammalian lung?

Pulmonary lipofibroblasts are thought to be involved in lung development, regeneration, vitamin A storage, and surfactant synthesis. Most of the evidence for these important functions relies on mouse or rat studies. Therefore, the present study was designed to investigate the presence of lipofibroblasts in a variety of early postnatal and adult mammalian species (including humans) to evaluate the ability to generalize functions of this cell type for other species. For this purpose, lung samples from 14 adult mammalian species as well as from postnatal mice, rats, and humans were investigated using light and electron microscopic stereology to obtain the volume fraction and the total volume of lipid bodies. In adult animals, lipid bodies were observed only, but not in all rodents. In all other species, no lipofibroblasts were observed. In rodents, lipid body volume scaled with body mass with an exponent b = 0.73 in the power law equation. Lipid bodies were not observed in postnatal human lungs but showed a characteristic postnatal increase in mice and rats and persisted at a lower level in the adult animals. Among 14 mammalian species, lipofibroblasts were only observed in rodents. The great increase in lipid body volume during early postnatal development of the mouse lung confirms the special role of lipofibroblasts during rodent lung development. It is evident that the cellular functions of pulmonary lipofibroblasts cannot be transferred easily from rodents to other species, in particular humans.

Imaging the Heliosphere using neutral Atoms from Solar Wind Energy down to 15 eV

We study the spatial and temporal distribution of hydrogen energetic neutral atoms (ENAs) from the heliosheath observed with the IBEX-Lo sensor of the Interstellar Boundary EXplorer (IBEX) from solar wind energies down to the lowest available energy (15 eV). All available IBEX-Lo data from 2009 January until 2013 June were included. The sky regions imaged when the spacecraft was outside of Earth's magnetosphere and when the Earth was moving toward the direction of observation offer a sufficient signal-to-noise ratio even at very low energies. We find that the ENA ribbon—a 20° wide region of high ENA intensities—is most prominent at solar wind energies whereas it fades at lower energies. The maximum emission in the ribbon is located near the poles for 2 keV and closer to the ecliptic plane for energies below 1 keV. This shift is an evidence that the ENA ribbon originates from the solar wind. Below 0.1 keV, the ribbon can no longer be identified against the globally distributed ENA signal. The ENA measurements in the downwind direction are affected by magnetospheric contamination below 0.5 keV, but a region of very low ENA intensities can be identified from 0.1 keV to 2 keV. The energy spectra of heliospheric ENAs follow a uniform power law down to 0.1 keV. Below this energy, they seem to become flatter, which is consistent with predictions. Due to the subtraction of local background, the ENA intensities measured with IBEX agree with the upper limit derived from Lyα observations.