Hydrodynamical simulations of early-type galaxies: effects of galaxy shape and stellar dynamics on hot coronae

Early-Type galaxies (ETGs) are embedded in hot (10^6-10^7 K), X-ray emitting gaseous haloes, produced mainly by stellar winds and heated by Type Ia supernovae explosions, by the thermalization of stellar motions and occasionally by the central super-massive black hole (SMBH). In particular, the thermalization of the stellar motions is due to the interaction between the stellar and the SNIa ejecta and the hot interstellar medium (ISM) already residing in the ETG. A number of different astrophysical phenomena determine the X-ray properties of the hot ISM, such as stellar population formation and evolution, galaxy structure and internal kinematics, Active Galactic Nuclei (AGN) presence, and environmental effects. With the aid of high-resolution hydrodynamical simulations performed on state-of-the-art galaxy models, in this Thesis we focus on the effects of galaxy shape, stellar kinematics and star formation on the evolution of the X-ray coronae of ETGs. Numerical simulations show that the relative importance of flattening and rotation are functions of the galaxy mass: at low galaxy masses, adding flattening and rotation induces a galactic wind, thus lowering the X-ray luminosity; at high galaxy masses the angular momentum conservation keeps the central regions of rotating galaxies at low density, whereas in non-rotating models a denser and brighter atmosphere is formed. The same dependence from the galaxy mass is present in the effects of star formation (SF): in light galaxies SF contributes to increase the spread in Lx, while at high galaxy masses the halo X-ray properties are marginally sensitive to SF effects. In every case, the star formation rate at the present epoch quite agrees with observations, and the massive, cold gaseous discs are partially or completely consumed by SF on a time-scale of few Gyr, excluding the presence of young stellar discs at the present epoch.

The electron antineutrino angular correlation coefficient a in free neutron decay: Testing the standard model with the aSPECT-spectrometer

The beta-decay of free neutrons is a strongly over-determined process in the Standard Model (SM) of Particle Physics and is described by a multitude of observables. Some of those observables are sensitive to physics beyond the SM. For example, the correlation coefficients of the involved particles belong to them. The spectrometer aSPECT was designed to measure precisely the shape of the proton energy spectrum and to extract from it the electron anti-neutrino angular correlation coefficient "a". A first test period (2005/ 2006) showed the “proof-of-principles”. The limiting influence of uncontrollable background conditions in the spectrometer made it impossible to extract a reliable value for the coefficient "a" (publication: Baessler et al., 2008, Europhys. Journ. A, 38, p.17-26). A second measurement cycle (2007/ 2008) aimed to under-run the relative accuracy of previous experiments (Stratowa et al. (1978), Byrne et al. (2002)) da/a =5%. I performed the analysis of the data taken there which is the emphasis of this doctoral thesis. A central point are background studies. The systematic impact of background on a was reduced to da/a(syst.)=0.61 %. The statistical accuracy of the analyzed measurements is da/a(stat.)=1.4 %. Besides, saturation effects of the detector electronics were investigated which were initially observed. These turned out not to be correctable on a sufficient level. An applicable idea how to avoid the saturation effects will be discussed in the last chapter.

Seismic tomographic full-waveform inversion for the Vrancea sinking lithosphere structure using the adjoint method.

The Vrancea region, at the south-eastern bend of the Carpathian Mountains in Romania, represents one of the most puzzling seismically active zones of Europe. Beside some shallow seismicity spread across the whole Romanian territory, Vrancea is the place of an intense seismicity with the presence of a cluster of intermediate-depth foci placed in a narrow nearly vertical volume. Although large-scale mantle seismic tomographic studies have revealed the presence of a narrow, almost vertical, high-velocity body in the upper mantle, the nature and the geodynamic of this deep intra-continental seismicity is still questioned. High-resolution seismic tomography could help to reveal more details in the subcrustal structure of Vrancea. Recent developments in computational seismology as well as the availability of parallel computing now allow to potentially retrieve more information out of seismic waveforms and to reach such high-resolution models. This study was aimed to evaluate the application of a full waveform inversion tomography at regional scale for the Vrancea lithosphere using data from the 1999 six months temporary local network CALIXTO. Starting from a detailed 3D Vp, Vs and density model, built on classical travel-time tomography together with gravity data, I evaluated the improvements obtained with the full waveform inversion approach. The latter proved to be highly problem dependent and highly computational expensive. The model retrieved after the first two iterations does not show large variations with respect to the initial model but remains in agreement with previous tomographic models. It presents a well-defined downgoing slab shape high velocity anomaly, composed of a N-S horizontal anomaly in the depths between 40 and 70km linked to a nearly vertical NE-SW anomaly from 70 to 180km.

From large polycyclic aromatic hydrocarbons to extended aromate-rich networks

The common ground of this study is the development of novel synthetic strategies to extended one-, two- and three-dimensional aromate-rich systems for which a number of applications are envisaged. rnThe point of departure is the synthesis and characterization of highly symmetric macrocyclic PAHs (polycyclic aromatic hydrocarbons) for which various aspects of supramolecular chemistry will be investigated. The versatility of the Yamamoto macrocyclization will be demonstrated on the basis of a set of cyclic trimers that exhibit a rich supramolecular chemistry. 1,10-phenanthroline, triphenylene and ortho-terphenyl building blocks have been successfully assembled to the corresponding macrocycles following the newly developed synthetic route. Scanning-tunneling microscopy (STM) and two-dimensional wide-angle X-ray scattering (2D-WAXS) were used to study the two- and three-dimensional self-assembly, respectively.rnSecondly, the development of chemical approaches to highly shape-anisotropic graphene nanoribbons (GNRs) and related nanographene molecules shall be discussed. Aryl-aryl coupling was used for the bottom-up fabrication of dendronized monomers, polymers and model compounds. Subsequently, these structures were converted into the final graphene material using oxidative (Scholl-type) cyclodehydrogenation. The GNRs thus obtained are characterized by an unprecedented length and lateral extension. The relevance of structural tailoring in the field of well-defined graphene materials is discussed in detail as only the chemical approach provides full geometry control. rnLastly, novel pathways towards the synthesis of extended three-dimensional networks that are dominated by nitrogen-rich motifs will be presented. If porous, these materials hold a great potential in the fields of gas and energy storage as well as for applications in catalysis. Hence, poly(aminal) networks based on melamine as crosslinking unit were synthesized and characterized with respect to the applications mentioned above. As set of conjugated poly(azomethine) networks was investigated regarding their use as a novel class of organic semiconductors for photocatalytic water splitting. The network structures described in this chapter can also be subjected to a controlled pyrolysis yielding mesoporous, nitrogen-rich carbon materials that were evaluated as active component for supercapacitors.rn

Towards the 3D attenuation imaging of active volcanoes: methods and tests on real and simulated data

The purpose of my PhD thesis has been to face the issue of retrieving a three dimensional attenuation model in volcanic areas. To this purpose, I first elaborated a robust strategy for the analysis of seismic data. This was done by performing several synthetic tests to assess the applicability of spectral ratio method to our purposes. The results of the tests allowed us to conclude that: 1) spectral ratio method gives reliable differential attenuation (dt*) measurements in smooth velocity models; 2) short signal time window has to be chosen to perform spectral analysis; 3) the frequency range over which to compute spectral ratios greatly affects dt* measurements. Furthermore, a refined approach for the application of spectral ratio method has been developed and tested. Through this procedure, the effects caused by heterogeneities of propagation medium on the seismic signals may be removed. The tested data analysis technique was applied to the real active seismic SERAPIS database. It provided a dataset of dt* measurements which was used to obtain a three dimensional attenuation model of the shallowest part of Campi Flegrei caldera. Then, a linearized, iterative, damped attenuation tomography technique has been tested and applied to the selected dataset. The tomography, with a resolution of 0.5 km in the horizontal directions and 0.25 km in the vertical direction, allowed to image important features in the off-shore part of Campi Flegrei caldera. High QP bodies are immersed in a high attenuation body (Qp=30). The latter is well correlated with low Vp and high Vp/Vs values and it is interpreted as a saturated marine and volcanic sediments layer. High Qp anomalies, instead, are interpreted as the effects either of cooled lava bodies or of a CO2 reservoir. A pseudo-circular high Qp anomaly was detected and interpreted as the buried rim of NYT caldera.

Fault detection, diagnosis and active fault tolerant control for a satellite attitude control system

Modern control systems are becoming more and more complex and control algorithms more and more sophisticated. Consequently, Fault Detection and Diagnosis (FDD) and Fault Tolerant Control (FTC) have gained central importance over the past decades, due to the increasing requirements of availability, cost efficiency, reliability and operating safety. This thesis deals with the FDD and FTC problems in a spacecraft Attitude Determination and Control System (ADCS). Firstly, the detailed nonlinear models of the spacecraft attitude dynamics and kinematics are described, along with the dynamic models of the actuators and main external disturbance sources. The considered ADCS is composed of an array of four redundant reaction wheels. A set of sensors provides satellite angular velocity, attitude and flywheel spin rate information. Then, general overviews of the Fault Detection and Isolation (FDI), Fault Estimation (FE) and Fault Tolerant Control (FTC) problems are presented, and the design and implementation of a novel diagnosis system is described. The system consists of a FDI module composed of properly organized model-based residual filters, exploiting the available input and output information for the detection and localization of an occurred fault. A proper fault mapping procedure and the nonlinear geometric approach are exploited to design residual filters explicitly decoupled from the external aerodynamic disturbance and sensitive to specific sets of faults. The subsequent use of suitable adaptive FE algorithms, based on the exploitation of radial basis function neural networks, allows to obtain accurate fault estimations. Finally, this estimation is actively exploited in a FTC scheme to achieve a suitable fault accommodation and guarantee the desired control performances. A standard sliding mode controller is implemented for attitude stabilization and control. Several simulation results are given to highlight the performances of the overall designed system in case of different types of faults affecting the ADCS actuators and sensors.

Constraining mass and shape of galaxy clusters through large scale structures

The mass estimation of galaxy clusters is a crucial point for modern cosmology, and can be obtained by several different techniques. In this work we discuss a new method to measure the mass of galaxy clusters connecting the gravitational potential of the cluster with the kinematical properties of its surroundings. We explore the dynamics of the structures located in the region outside virialized cluster, We identify groups of galaxies, as sheets or filaments, in the cluster outer region, and model how the cluster gravitational potential perturbs the motion of these structures from the Hubble fow. This identification is done in the redshift space where we look for overdensities with a filamentary shape. Then we use a radial mean velocity profile that has been found as a quite universal trend in simulations, and we fit the radial infall velocity profile of the overdensities found. The method has been tested on several cluster-size haloes from cosmological N-body simulations giving results in very good agreement with the true values of virial masses of the haloes and orientation of the sheets. We then applied the method to the Coma cluster and even in this case we found a good correspondence with previous. It is possible to notice a mass discrepancy between sheets with different alignments respect to the center of the cluster. This difference can be used to reproduce the shape of the cluster, and to demonstrate that the spherical symmetry is not always a valid assumption. In fact, if the cluster is not spherical, sheets oriented along different axes should feel a slightly different gravitational potential, and so give different masses as result of the analysis described before. Even this estimation has been tested on cosmological simulations and then applied to Coma, showing the actual non-sphericity of this cluster.

Towards antibody-mediated targeting of dendritic cells for cancer immunotherapy with multivalent polymer-antigen conjugates

Der Fokus dieser Arbeit lag auf der definierten Synthese multifunktioneller Polymer-Konjugate zur Anwendung in der Krebs-Immunotherapie. Durch gezielte Variation der Kon-jugationsbedingungen wurde Zusammensetzung, Größe und Aggregationsverhalten in Zell-medium sowie in humanem Serum untersucht. Nach definierter physikalisch-chemischer Charakterisierung wurde dann die induzierte Antigen-Präsentation zur Aktivierung der T-Zellproliferation analysiert.rnDafür wurden zwei verschiedene polymere Carrier-Systeme gewählt, lineares Poly-L-lysin und eine Polylysinbürste (PLL-Bürste). Es wird vermutet, dass die PLL-Bürste aufgrund der anisotropen Form eine bessere Verteilung im Körper und eine verlängerte Zirkulationsdauer zeigen wird. Die zu konjugierenden biologisch aktiven Komponenten waren der antiDEC205-Antikörper (aDEC205) für die gezielte Adressierung CD8-positiver dendritischer Zellen (DC), ein Ovalbumin (OVA)-spezifisches Antigen mit der Kernsequenz SIINFEKL für die Spezifität der Immunantwort gegen Krebszellen, die dieses Antigen tragen, und ein immunaktivieren-der TLR9-Ligand, CpG1826. Die Effizienz dieses Konjugates dendritische Zellen zu aktivieren, welche wiederum eine Immunantwort gegen OVA-exprimierende Krebszellen induzieren, wurde durch die Konjugation aller Komponenten am identischen Trägermolekül deutlich höher erwartet.rnLineares Poly-L-lysin diente als Modellsystem um die Konjugationschemie zu etablieren und dann auf die zylindrische Polylysinbürste zu übertragen. Anhand dieser polymeren Träger wurde das Verhalten der verschiedenen Topologien des Knäuels und der Bürste im Hinblick auf den Einfluss struktureller Unterschiede sowohl auf Konjugationsreaktionen als auch auf das in situ und in vitro Verhalten untersucht.rnFluoreszenzmarkiertes Antigen und der CpG Aktivator konnten jeweils aufgrund einer Thiol-Modifizierung an die Thiol-reaktive Maleimidgruppe des heterobifunktionellen Linkers Sulfo-SMCC an PLL-AlexaFluor48 konjugiert werden. Anschließend wurde aDEC205-AlexaFluor647 an PLL gekoppelt, entweder durch Schiff Base-Reaktion des oxidierten Antikörpers mit PLL und anschließender Reduzierung oder durch Click-Reaktion des PEG-Azids modifizierten An-tikörpers mit Dicyclobenzylcyclooctin (DIBO)-funktionalisiertem PLL. Die Konjugation der biologisch aktiven Komponenten wurde mit Durchflusszytometrie (FACS) und konfokaler Laser Scanning Mikroskopie (CLSM) untersucht und die Zusammensetzung des Konjugatesrnmittels UV/Vis-Spektroskopie bestimmt. Die PLL-Bürste alleine zeigte eine hohe Zytotoxizität bei HeLa und JAWS II Zelllinien, wohingegen lineares PLL und PLL-Konjugate sowie die PLL Bürsten-Konjugate keine ausgeprägte Zytotoxizität aufwiesen. Die Polymer-Konjugate wie-sen keine Aggregation in Zellmedium oder humanem Serum auf, was mittels winkelabhängi-ger dynamischer Lichtstreuung bestimmt wurde. CLSM Aufnahmen zeigten Kolokalisation der an die einzelnen Komponenten gebundenen Fluoreszenzfarbstoffe in dendritischen Zel-len, was die erfolgreiche Konjugation und Internalisierung der Konjugate in die Zellen bele-gen konnte. FACS Messungen ergaben eine geringfügig erhöhte Aufnahme des adressierten PLL-Antigen-Antikörper-Konjugates verglichen mit dem PLL-Antigen-Konjugat. Experimente mit dem „Specific Hybridization Internalization Sensor“ (SHIP) zeigten jedoch nur Aufnahme der PLL-Konjugate in CD8+ unreife DC, nicht in reife DC, die nicht mehr unspezifisch, sondern nur noch über Rezeptoren internalisieren. Dies bewies die unspezifische Aufnahme des Kon-jugates, da Antikörper-Konjugation keine Rezeptor-vermittelte Endozytose in reife DC indu-zieren konnte. T-Zell-Proliferationsassays ergaben eine Aktivierung von CD8+ T-Zellen indu-ziert durch Antigen-tragende Konjugate, wohingegen Konjugate ohne Antigen als Negativ-kontrollen dienten und keine T-Zell-Proliferation erzielten. Es konnte jedoch kein Unter-schied zwischen adressierten und nicht adressierten Konjugaten aufgrund der unspezifischen Aufnahme durch das Polymer beobachtet werden. Lösliches SIINFEKL alleine bewirkte schon bei geringeren Konzentrationen eine T-Zell-Proliferation.rnEs war somit möglich, drei biologischen Komponenten an einen polymeren Träger zu konju-gieren und diese Konjugate im Hinblick auf Zusammensetzung, Größe, Internalisierung in dendritische Zellen und Aktivierung der T-Zell-Proliferation zu untersuchen. Außerdem wur-de die Konjugationschemie erfolgreich von dem Modellsystem des linearen PLL auf die PLL-Bürste übertragen. Die Polymer-Konjugate werde unspezifisch in DC aufgenommen und in-duzieren T-Zellproliferation, die mit Antigen-Präsentationsassays nachgewiesen wird. Es konnte jedoch durch Konjugation des Antikörpers keine Rezeptor-vermittelte Aufnahme in CD8+ DC erzielt werden.rnDiese Studien stellen einen erfolgsversprechenden ersten Schritt zur Entwicklung neuer Na-nomaterialien für die Anwendung in Krebs-Immuntherapie dar.

Computer simulation methods to study interfacial tensions : from the Ising model to colloidal crystals

In condensed matter systems, the interfacial tension plays a central role for a multitude of phenomena. It is the driving force for nucleation processes, determines the shape and structure of crystalline structures and is important for industrial applications. Despite its importance, the interfacial tension is hard to determine in experiments and also in computer simulations. While for liquid-vapor interfacial tensions there exist sophisticated simulation methods to compute the interfacial tension, current methods for solid-liquid interfaces produce unsatisfactory results.rnrnAs a first approach to this topic, the influence of the interfacial tension on nuclei is studied within the three-dimensional Ising model. This model is well suited because despite its simplicity, one can learn much about nucleation of crystalline nuclei. Below the so-called roughening temperature, nuclei in the Ising model are not spherical anymore but become cubic because of the anisotropy of the interfacial tension. This is similar to crystalline nuclei, which are in general not spherical but more like a convex polyhedron with flat facets on the surface. In this context, the problem of distinguishing between the two bulk phases in the vicinity of the diffuse droplet surface is addressed. A new definition is found which correctly determines the volume of a droplet in a given configuration if compared to the volume predicted by simple macroscopic assumptions.rnrnTo compute the interfacial tension of solid-liquid interfaces, a new Monte Carlo method called ensemble switch method'' is presented which allows to compute the interfacial tension of liquid-vapor interfaces as well as solid-liquid interfaces with great accuracy. In the past, the dependence of the interfacial tension on the finite size and shape of the simulation box has often been neglected although there is a nontrivial dependence on the box dimensions. As a consequence, one needs to systematically increase the box size and extrapolate to infinite volume in order to accurately predict the interfacial tension. Therefore, a thorough finite-size scaling analysis is established in this thesis. Logarithmic corrections to the finite-size scaling are motivated and identified, which are of leading order and therefore must not be neglected. The astounding feature of these logarithmic corrections is that they do not depend at all on the model under consideration. Using the ensemble switch method, the validity of a finite-size scaling ansatz containing the aforementioned logarithmic corrections is carefully tested and confirmed. Combining the finite-size scaling theory with the ensemble switch method, the interfacial tension of several model systems, ranging from the Ising model to colloidal systems, is computed with great accuracy.

Efficiency limits for silicon/perovskite tandem solar cells: a theoretical model

The primary goal of this work is related to the extension of an analytic electro-optical model. It will be used to describe single-junction crystalline silicon solar cells and a silicon/perovskite tandem solar cell in the presence of light-trapping in order to calculate efficiency limits for such a device. In particular, our tandem system is composed by crystalline silicon and a perovskite structure material: metilammoniumleadtriiodide (MALI). Perovskite are among the most convenient materials for photovoltaics thanks to their reduced cost and increasing efficiencies. Solar cell efficiencies of devices using these materials increased from 3.8% in 2009 to a certified 20.1% in 2014 making this the fastest-advancing solar technology to date. Moreover, texturization increases the amount of light which can be absorbed through an active layer. Using Green’s formalism it is possible to calculate the photogeneration rate of a single-layer structure with Lambertian light trapping analytically. In this work we go further: we study the optical coupling between the two cells in our tandem system in order to calculate the photogeneration rate of the whole structure. We also model the electronic part of such a device by considering the perovskite top cell as an ideal diode and solving the drift-diffusion equation with appropriate boundary conditions for the silicon bottom cell. We have a four terminal structure, so our tandem system is totally unconstrained. Then we calculate the efficiency limits of our tandem including several recombination mechanisms such as Auger, SRH and surface recombination. We focus also on the dependence of the results on the band gap of the perovskite and we calculare an optimal band gap to optimize the tandem efficiency. The whole work has been continuously supported by a numerical validation of out analytic model against Silvaco ATLAS which solves drift-diffusion equations using a finite elements method. Our goal is to develop a simpler and cheaper, but accurate model to study such devices.

Statistical shape modeling of pathological scoliotic vertebrae: A comparative analysis

Statistical shape models (SSMs) have been used widely as a basis for segmenting and interpreting complex anatomical structures. The robustness of these models are sensitive to the registration procedures, i.e., establishment of a dense correspondence across a training data set. In this work, two SSMs based on the same training data set of scoliotic vertebrae, and registration procedures were compared. The first model was constructed based on the original binary masks without applying any image pre- and post-processing, and the second was obtained by means of a feature preserving smoothing method applied to the original training data set, followed by a standard rasterization algorithm. The accuracies of the correspondences were assessed quantitatively by means of the maximum of the mean minimum distance (MMMD) and Hausdorf distance (H(D)). Anatomical validity of the models were quantified by means of three different criteria, i.e., compactness, specificity, and model generalization ability. The objective of this study was to compare quasi-identical models based on standard metrics. Preliminary results suggest that the MMMD distance and eigenvalues are not sensitive metrics for evaluating the performance and robustness of SSMs.

Numerical model of the human cornea based on stromal microstructure: identification of mechanical properties for two age groups

The optical quality of the human eye mainly depends on the refractive performance of the cornea. The shape of the cornea is a mechanical balance between intraocular pressure and tissue intrinsic stiffness. Several surgical procedures in ophthalmology alter the biomechanics of the cornea to provoke local or global curvature changes for vision correction. Legitimated by the large number of surgical interventions performed every day, the demand for a deeper understanding of corneal biomechanics is rising to improve the safety of procedures and medical devices. The aim of our work is to propose a numerical model of corneal biomechanics, based on the stromal microstructure. Our novel anisotropic constitutive material law features a probabilistic weighting approach to model collagen fiber distribution as observed on human cornea by Xray scattering analysis (Aghamohammadzadeh et. al., Structure, February 2004). Furthermore, collagen cross-linking was explicitly included in the strain energy function. Results showed that the proposed model is able to successfully reproduce both inflation and extensiometry experimental data (Elsheikh et. al., Curr Eye Res, 2007; Elsheikh et. al., Exp Eye Res, May 2008). In addition, the mechanical properties calculated for patients of different age groups (Group A: 65-79 years; Group B: 80-95 years) demonstrate an increased collagen cross-linking, and a decrease in collagen fiber elasticity from younger to older specimen. These findings correspond to what is known about maturing fibrous biological tissue. Since the presented model can handle different loading situations and includes the anisotropic distribution of collagen fibers, it has the potential to simulate clinical procedures involving nonsymmetrical tissue interventions. In the future, such mechanical model can be used to improve surgical planning and the design of next generation ophthalmic devices.

PSGL-1 is dispensible for the development of active experimental autoimmune encephalomyelitis in SJL/J mice

The adhesion molecule P-selectin glycoprotein ligand (PSGL)-1 has been suggested to be involved in the immunopathogenesis of multiple sclerosis (MS). However, in C57BL/6 mice PSGL-1 was found to be dispensible for the development of MOG(aa35-55)-induced experimental autoimmune encephalomyelitis (EAE), an animal model for MS. To study, if involvement of PSGL-1 to EAE pathogenesis can be observed in another common mouse model, we backcrossed PSGL-1(-/-) mice for at least 12 generations into the SJL/J background and compared PLP(aa139-151) induced EAE in PSGL-1(-/-) SJL/J mice versus wild-type SJL/J mice. Here, we demonstrate that PSGL-1(-/-) SJL/J mice exhibited EAE pathogenesis indistinguishable from wild-type SJL/J mice. Our present study underscores and emphasizes previous observations that PSGL-1 is dispensible for EAE pathogenesis.

Automated cement segmentation in vertebroplasty

Vertebroplasty is a minimally invasive procedure with many benefits; however, the procedure is not without risks and potential complications, of which leakage of the cement out of the vertebral body and into the surrounding tissues is one of the most serious. Cement can leak into the spinal canal, venous system, soft tissues, lungs and intradiscal space, causing serious neurological complications, tissue necrosis or pulmonary embolism. We present a method for automatic segmentation and tracking of bone cement during vertebroplasty procedures, as a first step towards developing a warning system to avoid cement leakage outside the vertebral body. We show that by using active contours based on level sets the shape of the injected cement can be accurately detected. The model has been improved for segmentation as proposed in our previous work by including a term that restricts the level set function to the vertebral body. The method has been applied to a set of real intra-operative X-ray images and the results show that the algorithm can successfully detect different shapes with blurred and not well-defined boundaries, where the classical active contours segmentation is not applicable. The method has been positively evaluated by physicians.

Mesh-based vs. Image-based Statistical Appearance Model of the Human Femur: a Preliminary Comparison Study for the Creation of Finite Element Meshes

Statistical models have been recently introduced in computational orthopaedics to investigate the bone mechanical properties across several populations. A fundamental aspect for the construction of statistical models concerns the establishment of accurate anatomical correspondences among the objects of the training dataset. Various methods have been proposed to solve this problem such as mesh morphing or image registration algorithms. The objective of this study is to compare a mesh-based and an image-based statistical appearance model approaches for the creation of nite element(FE) meshes. A computer tomography (CT) dataset of 157 human left femurs was used for the comparison. For each approach, 30 finite element meshes were generated with the models. The quality of the obtained FE meshes was evaluated in terms of volume, size and shape of the elements. Results showed that the quality of the meshes obtained with the image-based approach was higher than the quality of the mesh-based approach. Future studies are required to evaluate the impact of this finding on the final mechanical simulations.