On the Luminous and Dark Matter Distribution in Early-Type Galaxies

The way mass is distributed in galaxies plays a major role in shaping their evolution across cosmic time. The galaxy's total mass is usually determined by tracing the motion of stars in its potential, which can be probed observationally by measuring stellar spectra at different distances from the galactic centre, whose kinematics is used to constrain dynamical models. A class of such models, commonly used to accurately determine the distribution of luminous and dark matter in galaxies, is that of equilibrium models. In this Thesis, a novel approach to the design of equilibrium dynamical models, in which the distribution function is an analytic function of the action integrals, is presented. Axisymmetric and rotating models are used to explain observations of a sample of nearby early-type galaxies in the Calar Alto Legacy Integral Field Area survey. Photometric and spectroscopic data for round and flattened galaxies are well fitted by the models, which are then used to get the galaxies' total mass distribution and orbital anisotropy. The time evolution of massive early-type galaxies is also investigated with numerical models. Their structural properties (mass, size, velocity dispersion) are observed to evolve, on average, with redshift. In particular, they appear to be significantly more compact at higher redshift, at fixed stellar mass, so it is interesting to investigate what drives such evolution. This Thesis focuses on the role played by dark-matter haloes: their mass-size and mass-velocity dispersion correlations evolve similarly to the analogous correlations of ellipticals; at fixed halo mass, the haloes are more compact at higher redshift, similarly to massive galaxies; a simple model, in which all the galaxy's size and velocity-dispersion evolution is due to the cosmological evolution of the underlying halo population, reproduces the observed size and velocity-dispersion of massive compact early-type galaxies up to redshift of about 2.

Charge transport in organic photovoltaic cells

Die vorliegende Dissertation dient dazu, das Verständnis des Ladungstransportes in organischen Solarzellen zu vertiefen. Mit Hilfe von Computersimulationen wird die Bewegung von Ladungsträgern in organischen Materialien rekonstruiert, und zwar ausgehend von den quantenmechanischen Prozessen auf mikroskopischer Ebene bis hin zur makroskopischen Skala, wo Ladungsträgermobilitäten quantifizierbar werden. Auf Grundlage dieses skalenübergreifenden Ansatzes werden Beziehungen zwischen der chemischen Struktur organischer Moleküle und der makroskopischen Mobilität hergestellt (Struktur-Eigenschafts-Beziehungen), die zu der Optimierung photovoltaischer Wirkungsgrade beitragen. Das Simulationsmodell beinhaltet folgende drei Schlüsselkomponenten. Erstens eine Morphologie, d. h. ein atomistisch aufgelöstes Modell der molekularen Anordnung in dem untersuchten Material. Zweitens ein Hüpfmodell des Ladungstransportes, das Ladungswanderung als eine Abfolge von Ladungstransferreaktionen zwischen einzelnen Molekülen beschreibt. Drittens ein nichtadiabatisches Modell des Ladungstransfers, das Übergangsraten durch drei Parameter ausdrückt: Reorganisationsenergien, Lageenergien und Transferintegrale. Die Ladungstransport-Simulationen richten sich auf die Materialklasse der dicyanovinyl-substituierten Oligothiophene und umfassen Morphologien von Einkristallen, Dünnschichten sowie amorphen/smektischen Mesophasen. Ein allgemeiner Befund ist, dass die molekulare Architektur, bestehend aus einer Akzeptor-Donor-Akzeptor-Sequenz und einem flexiblen Oligomergerüst, eine erhebliche Variation molekularer Dipolmomente und damit der Lageenergien bewirkt. Diese energetische Unordnung ist ungewöhnlich hoch in den Kristallen und umso höher in den Mesophasen. Für die Einkristalle wird beobachtet, dass Kristallstrukturen mit ausgeprägter π-Stapelung und entsprechend großer Transferintegrale zu verhältnismäßig niedrigen Mobilitäten führen. Dieses Verhalten wird zurückgeführt auf die Ausbildung bevorzugter Transportrichtungen, die anfällig für energetische Störungen sind. Für die Dünnschichten bestätigt sich diese Argumentation und liefert ein mikroskopisches Verständnis für experimentelle Mobilitäten. In der Tat korrelieren die Simulationsergebnisse sowohl mit gemessenen Mobilitäten als auch mit photovoltaischen Wirkungsgraden. Für die amorphen/smektischen Systeme steigt die energetische Unordnung mit der Oligomerlänge, sie führt aber auch zu einer unerwarteten Mobilitätsabnahme in dem stärker geordneten smektischen Zustand. Als Ursache dafür erweist sich, dass die smektische Schichtung der räumlichen Korrelation der energetischen Unordnung entgegensteht.

Hunting for warped extra dimensions via loop-induced processes

This thesis is on loop-induced processes in theories with warped extra dimensions where the fermions and gauge bosons are allowed to propagate in the bulk, while the Higgs sector is localized on or near the infra-red brane. These so-called Randall-Sundrum (RS) models have the potential to simultaneously explain the hierarchy problem and address the question of what causes the large hierarchies in the fermion sector of the Standard Model (SM). The Kaluza-Klein (KK) excitations of the bulk fields can significantly affect the loop-level processes considered in this thesis and, hence, could indirectly indicate the existence of warped extra dimensions. The analytical part of this thesis deals with the detailed calculation of three loop-induced processes in the RS models in question: the Higgs production process via gluon fusion, the Higgs decay into two photons, and the flavor-changing neutral current b → sγ. A comprehensive, five-dimensional (5D) analysis will show that the amplitudes of the Higgs processes can be expressed in terms of integrals over 5D propagators with the Higgs-boson profile along the extra dimension, which can be used for arbitrary models with a compact extra dimension. To this end, both the boson and fermion propagators in a warped 5D background are derived. It will be shown that the seemingly contradictory results for the gluon fusion amplitude in the literature can be traced back to two distinguishable, not smoothly-connected incarnations of the RS model. The investigation of the b → sγ transition is performed in the KK decomposed theory. It will be argued that summing up the entire KK tower leads to a finite result, which can be well approximated by a closed, analytical expression.rnIn the phenomenological part of this thesis, the analytic results of all relevant Higgs couplings in the RS models in question are compared with current and in particular future sensitivities of the Large Hadron Collider (LHC) and the planned International Linear Collider. The latest LHC Higgs data is then used to exclude significant portions of the parameter space of each RS scenario. The analysis will demonstrate that especially the loop-induced Higgs couplings are sensitive to KK particles of the custodial RS model with masses in the multi tera-electronvolt range. Finally, the effect of the RS model on three flavor observables associated with the b → sγ transition are examined. In particular, we study the branching ratio of the inclusive decay B → X_s γ

Applications of elliptic functions to solve differential equations

In questa tesi si mostrano alcune applicazioni degli integrali ellittici nella meccanica Hamiltoniana, allo scopo di risolvere i sistemi integrabili. Vengono descritte le funzioni ellittiche, in particolare la funzione ellittica di Weierstrass, ed elenchiamo i tipi di integrali ellittici costruendoli dalle funzioni di Weierstrass. Dopo aver considerato le basi della meccanica Hamiltoniana ed il teorema di Arnold Liouville, studiamo un esempio preso dal libro di Moser-Integrable Hamiltonian Systems and Spectral Theory, dove si prendono in considerazione i sistemi integrabili lungo la geodetica di un'ellissoide, e il sistema di Von Neumann. In particolare vediamo che nel caso n=2 abbiamo un integrale ellittico.

Whither PQL?

Generalized linear mixed models (GLMM) are generalized linear models with normally distributed random effects in the linear predictor. Penalized quasi-likelihood (PQL), an approximate method of inference in GLMMs, involves repeated fitting of linear mixed models with “working” dependent variables and iterative weights that depend on parameter estimates from the previous cycle of iteration. The generality of PQL, and its implementation in commercially available software, has encouraged the application of GLMMs in many scientific fields. Caution is needed, however, since PQL may sometimes yield badly biased estimates of variance components, especially with binary outcomes. Recent developments in numerical integration, including adaptive Gaussian quadrature, higher order Laplace expansions, stochastic integration and Markov chain Monte Carlo (MCMC) algorithms, provide attractive alternatives to PQL for approximate likelihood inference in GLMMs. Analyses of some well known datasets, and simulations based on these analyses, suggest that PQL still performs remarkably well in comparison with more elaborate procedures in many practical situations. Adaptive Gaussian quadrature is a viable alternative for nested designs where the numerical integration is limited to a small number of dimensions. Higher order Laplace approximations hold the promise of accurate inference more generally. MCMC is likely the method of choice for the most complex problems that involve high dimensional integrals.

Arterial blood pressure during early sepsis and outcome

OBJECTIVE: To evaluate the association between arterial blood pressure (ABP) during the first 24 h and mortality in sepsis. DESIGN: Retrospective cohort study. SETTING: Multidisciplinary intensive care unit (ICU). PATIENTS AND PARTICIPANTS: A total of 274 septic patients. INTERVENTIONS: None. MEASUREMENTS AND RESULTS: Hemodynamic, and laboratory parameters were extracted from a PDMS database. The hourly time integral of ABP drops below clinically relevant systolic arterial pressure (SAP), mean arterial pressure (MAP), and mean perfusion pressure (MPP = MAP - central venous pressure) levels was calculated for the first 24 h after ICU admission and compared with 28-day-mortality. Binary and linear regression models (adjusted for SAPS II as a measure of disease severity), and a receiver operating characteristic (ROC) analysis were applied. The areas under the ROC curve were largest for the hourly time integrals of ABP drops below MAP 60 mmHg (0.779 vs. 0.764 for ABP drops below MAP 55 mmHg; P < or = 0.01) and MPP 45 mmHg. No association between the hourly time integrals of ABP drops below certain SAP levels and mortality was detected. One or more episodes of MAP < 60 mmHg increased the risk of death by 2.96 (CI 95%, 1.06-10.36, P = 0.04). The area under the ROC curve to predict the need for renal replacement therapy was highest for the hourly time integral of ABP drops below MAP 75 mmHg. CONCLUSIONS: A MAP level > or = 60 mmHg may be as safe as higher MAP levels during the first 24 h of ICU therapy in septic patients. A higher MAP may be required to maintain kidney function.

Pulmonary venous flow velocity patterns in 404 individuals without cardiovascular disease

OBJECTIVE To determine the pulmonary venous flow velocity (PVFV) values in a large normal population. DESIGN Prospective study in consecutive individuals. SETTING University hospital. METHODS Among 404 normal individuals, the flow velocity pattern in the right upper pulmonary vein was recorded in 315 subjects using transthoracic echocardiography, and in both upper pulmonary veins in 100 subjects using transoesophageal echocardiography. Subjects were divided into five age groups. The PVFV values were compared between transthoracic and transoesophageal echocardiography within the age groups, and intraindividually between the right and left upper pulmonary veins in transoesophageal echocardiography. RESULTS Normal PVFV values for the right upper pulmonary vein in transthoracic and transoesophageal echocardiography are presented. The duration of flow reversal at atrial contraction was overestimated using transthoracic echocardiography (mean (SD): 96 (21) ms in transoesophageal echocardiography, 120 (28) ms in transthoracic echocardiography, p < 0.0001). Systolic to diastolic peak flow velocity ratio (S:D) increased earlier with advancing age with transoesophageal echocardiography than with transthoracic echocardiography. Similar results were found for the corresponding time-velocity integrals. Data from the left and right upper pulmonary veins differed with respect to onset and deceleration of flow velocities, but not for flow durations or peak velocities. CONCLUSIONS Normal PVFV values generally show a wide range. The data presented will be of value in assessing left ventricular diastolic function and mitral regurgitation using the PVFV pattern.

A Stieltjes approach to static hedges

Static hedging of complicated payoff structures by standard instruments becomes increasingly popular in finance. The classical approach is developed for quite regular functions, while for less regular cases, generalized functions and approximation arguments are used. In this note, we discuss the regularity conditions in the classical decomposition formula due to P. Carr and D. Madan (in Jarrow ed, Volatility, pp. 417–427, Risk Publ., London, 1998) if the integrals in this formula are interpreted as Lebesgue integrals with respect to the Lebesgue measure. Furthermore, we show that if we replace these integrals by Lebesgue–Stieltjes integrals, the family of representable functions can be extended considerably with a direct approach.

A scrutiny of hard pion chiral perturbation theory

In the recently proposed framework of hard pion chiral perturbation theory, the leading chiral logarithms are predicted to factorize with respect to the energy dependence in the chiral limit. We have scrutinized this assumption in the case of vector and scalar pion form factors FV;S(s) by means of standard chiral perturbation theory and dispersion relations. We show that this factorization property is valid for the elastic contribution to the dispersion integrals for FV;S(s) but it is violated starting at three loops when the inelastic four-pion contributions arise.

Factorization of leading chiral logarithms in the pion form factors

In the recently proposed framework of hard pion chiral perturbation theory, the leading chiral logarithms are predicted to factorize with respect to the energy dependence in the chiral limit. We have scrutinized this assumption in the case of vector and scalar pion form factors FV;S(s) by means of standard chiral perturbation theory and dispersion relations. We show that this factorization property is valid for the elastic contribution to the dispersion integrals for FV;S(s) but it is violated starting at three loops when the inelastic four-pion contributions arise.

Effect of surface restoring on subsurface variability in a climate model during 1949–2005

Initializing the ocean for decadal predictability studies is a challenge, as it requires reconstructing the little observed subsurface trajectory of ocean variability. In this study we explore to what extent surface nudging using well-observed sea surface temperature (SST) can reconstruct the deeper ocean variations for the 1949–2005 period. An ensemble made with a nudged version of the IPSLCM5A model and compared to ocean reanalyses and reconstructed datasets. The SST is restored to observations using a physically-based relaxation coefficient, in contrast to earlier studies, which use a much larger value. The assessment is restricted to the regions where the ocean reanalyses agree, i.e. in the upper 500 m of the ocean, although this can be latitude and basin dependent. Significant reconstruction of the subsurface is achieved in specific regions, namely region of subduction in the subtropical Atlantic, below the thermocline in the equatorial Pacific and, in some cases, in the North Atlantic deep convection regions. Beyond the mean correlations, ocean integrals are used to explore the time evolution of the correlation over 20-year windows. Classical fixed depth heat content diagnostics do not exhibit any significant reconstruction between the different existing bservation-based references and can therefore not be used to assess global average time-varying correlations in the nudged simulations. Using the physically based average temperature above an isotherm (14°C) alleviates this issue in the tropics and subtropics and shows significant reconstruction of these quantities in the nudged simulations for several decades. This skill is attributed to the wind stress reconstruction in the tropics, as already demonstrated in a perfect model study using the same model. Thus, we also show here the robustness of this result in an historical and observational context.

Steiner’s formula in the Heisenberg group

Steiner’s tube formula states that the volume of an ϵ-neighborhood of a smooth regular domain in Rn is a polynomial of degree n in the variable ϵ whose coefficients are curvature integrals (also called quermassintegrals). We prove a similar result in the sub-Riemannian setting of the first Heisenberg group. In contrast to the Euclidean setting, we find that the volume of an ϵ-neighborhood with respect to the Heisenberg metric is an analytic function of ϵ that is generally not a polynomial. The coefficients of the series expansion can be explicitly written in terms of integrals of iteratively defined canonical polynomials of just five curvature terms.

Phytoplankton daily production in the Rhode River, Maryland (USA), 1990-2009

Parameters in the photosynthesis-irradiance (P-E) relationship of phytoplankton were used to calculate daily production at weekly to bi-weekly intervals for 20 years at 6 stations on the Rhode River, Maryland (USA). The objectives of this work were to determine the patterns and controls on the P-E parameters and primary production of phytoplankton in a shallow eutrophic estuary. Additional measurements that are components of calculated daily rates of primary productivity are given: the light-saturation irradiance, photoperiod, maximal noon incident irradiance, optical depth, dimensionless depth integrals, and a correction for spectral selectivity of light absorption. P-E parameters and chlorophyll a concentrations were given in a related dataset, Gallegos (2012, doi:10.1594/PANGAEA.816494).

Theoretical Analysis of the No-Slip Boundary Condition Enforcement in SPH Methods

The aim of the present work is to provide an in-depth analysis of the most representative mirroring techniques used in SPH to enforce boundary conditions (BC) along solid profiles. We specifically refer to dummy particles, ghost particles, and Takeda et al. [Prog. Theor. Phys. 92 (1994), 939] boundary integrals. The analysis has been carried out by studying the convergence of the first- and second-order differential operators as the smoothing length (that is, the characteristic length on which relies the SPH interpolation) decreases. These differential operators are of fundamental importance for the computation of the viscous drag and the viscous/diffusive terms in the momentum and energy equations. It has been proved that close to the boundaries some of the mirroring techniques leads to intrinsic inaccuracies in the convergence of the differential operators. A consistent formulation has been derived starting from Takeda et al. boundary integrals (see the above reference). This original formulation allows implementing no-slip boundary conditions consistently in many practical applications as viscous flows and diffusion problems.