Petrografia e geochimica delle vulcaniti del Monte dei Porri di Salina (Isole Eolie)

E’ stata indagata dal punto di vista stratigrafico, geochimico e petrografico la successione di vulcaniti appartenenti alle formazioni Rocce di Barcone (Tufi Grigi dei Porri o Grey Porri Tuffs – GPT), Punta Sallustro, Punta di Megna e Serro dello Sciarato (Lucchi et al., 2013a, b), riferibili all’Epoca Eruttiva 5 (~70 - 57 ka) dell’evoluzione vulcanologica dell’isola Salina (isole Eolie), durante la quale si realizza l’edificazione della porzione medio-basale dello stratocono del Monte dei Porri. Attraverso un approccio basato sull’integrazione e correlazione di metodologie stratigrafiche, petrografiche e geochimiche si è giunti ad una ricostruzione dei meccanismi eruttivi, deposizionali e petrogenetici che sottendono la successione di eventi vulcanici responsabili dell’edificazione dello stratocono suddetto. Utilizzando l’analisi di litofacies si è giunti, anche attraverso le correlazioni stratigrafiche estese sull’isola di Lipari, alla definizione di cinque unità eruttive (EU1-EU5), rappresentanti il riferimento stratigrafico per la campionatura necessaria per caratterizzarle sia petrograficamente che geochimicamente. Il significato vulcanologico delle EU1-5 ha permesso di evidenziare che il Monte dei Porri è costituito, in larga parte, da un’articolata alternanza di depositi da caduta e da correnti piroclastiche (formati sia da pomici che da scorie), alla quale si intercalano colate laviche; all’interno dell’unità i prodotti analizzati evidenziano composizione dei magmi da trachi-dacitica (EU1), a basaltica e andesitico basaltica (EU2), a francamente andesitico basaltica (EU3); quest’ultima composizione caratterizza la EU4 (che evolve anche verso termini andesitici), e la EU5. In conclusione, l’integrazione di tutti i dati (stratigrafici, petrografici e geochimici) permette di affermare che l’attività vulcanica responsabile della messa in posto delle EU1-5 cominci con una fase di apertura del condotto eruttivo (fase esplosiva pliniana) accompagnata dall’emissione dei magmi più evoluti (EU1), residenti al top di una camera magmatica zonata. Ad essa presumibilmente segue il coinvolgimento dei livelli più profondi del reservoir magmatico dove risiedono i magmi più mafici (EU2 e EU3). La ripresa dell’attività vulcanica (dopo una stasi durante il quale i magmi mafici evolvono verso composizioni mediamente evolute) vede infatti l’emissione di prodotti andesitici (EU4) seguita da magmi meno evoluti (EU5) che porta ad un progressivo svuotamento del sistema.

Rilevamento geologico-strutturale lungo la Linea di Belluno nella zona di Monte Piad (Provincia di Belluno)

Il prodotto principale del lavoro è la Carta Geologica alla scala 1:5.000 che ha portato un avanzamento delle conoscenze nella zona in esame. Dal punto di vista stratigrafico il Gruppo dei Calcari Grigi è stato suddiviso nelle formazioni già distinte nei Fogli CARG Tione, Malè, Trento e Asiago: la Formazione di Monte Zugna, Fm. Di Loppio, Fm. Di Rotzo, Fm. Di Grigno. È stato possibile anche rilevare e cartografare l’Oolite di S. Vigilio, anche detta Encrinite di Monte Agaro, con uno spessore di circa 5 m sino alla zona di Sass de Falares, cioè più ad est di quanto noto in letteratura. E’ stato osservato che la Linea di Belluno non è un piano unico, ma comprende altre due faglie: la Linea di Monte Piad che consiste in una faglia inversa che taglia in cerniera l’Anticlinale del Monte Coppolo, e poco a sud la Linea di Sasso Falares che rappresenta la faglia che ha generato la piega per propagazione di faglia del Coppolo, e proseguendo poi verso l’alto ha tagliato tutta la piega emergendo in superficie. In tal modo si forma il duplex di Sasso Falares, arrangiato in una blanda coppia anticinale-sinclinale tipica della geometria di queste strutture, che risulta delimitato a nord dalla linea omonima e a sud dalla L. di Belluno s.s. che si incontrano a quota circa 700 m. Si è potuta anche ricostruire la cinematica del piano principale del thrust di Belluno studiando l’affioramento eccezionale messo a nudo di recente. Gli assi degli sforzi, agenti in compressione, ricostruiti tramite l’analisi meso-strutturale, testimoniano una compressione orientata NNW-SSE, correlabile con l’Evento Valsuganese del Serravalliano-Tortoniano, legato ad un asse compressivo N340 che rappresenta l’evento principale nella regione sudalpina. Anche l’Evento compressivo del Cattiano-Burdigalliano con asse N30 ed una compressione orientata circa E-W correlabile con l’Evento Scledense del Messiniano-Pliocene risultano dall’analisi strutturale.

Exact quantum Monte Carlo methods for correlated electron systems

Thema dieser Arbeit ist die Entwicklung und Kombination verschiedener numerischer Methoden, sowie deren Anwendung auf Probleme stark korrelierter Elektronensysteme. Solche Materialien zeigen viele interessante physikalische Eigenschaften, wie z.B. Supraleitung und magnetische Ordnung und spielen eine bedeutende Rolle in technischen Anwendungen. Es werden zwei verschiedene Modelle behandelt: das Hubbard-Modell und das Kondo-Gitter-Modell (KLM). In den letzten Jahrzehnten konnten bereits viele Erkenntnisse durch die numerische Lösung dieser Modelle gewonnen werden. Dennoch bleibt der physikalische Ursprung vieler Effekte verborgen. Grund dafür ist die Beschränkung aktueller Methoden auf bestimmte Parameterbereiche. Eine der stärksten Einschränkungen ist das Fehlen effizienter Algorithmen für tiefe Temperaturen.rnrnBasierend auf dem Blankenbecler-Scalapino-Sugar Quanten-Monte-Carlo (BSS-QMC) Algorithmus präsentieren wir eine numerisch exakte Methode, die das Hubbard-Modell und das KLM effizient bei sehr tiefen Temperaturen löst. Diese Methode wird auf den Mott-Übergang im zweidimensionalen Hubbard-Modell angewendet. Im Gegensatz zu früheren Studien können wir einen Mott-Übergang bei endlichen Temperaturen und endlichen Wechselwirkungen klar ausschließen.rnrnAuf der Basis dieses exakten BSS-QMC Algorithmus, haben wir einen Störstellenlöser für die dynamische Molekularfeld Theorie (DMFT) sowie ihre Cluster Erweiterungen (CDMFT) entwickelt. Die DMFT ist die vorherrschende Theorie stark korrelierter Systeme, bei denen übliche Bandstrukturrechnungen versagen. Eine Hauptlimitation ist dabei die Verfügbarkeit effizienter Störstellenlöser für das intrinsische Quantenproblem. Der in dieser Arbeit entwickelte Algorithmus hat das gleiche überlegene Skalierungsverhalten mit der inversen Temperatur wie BSS-QMC. Wir untersuchen den Mott-Übergang im Rahmen der DMFT und analysieren den Einfluss von systematischen Fehlern auf diesen Übergang.rnrnEin weiteres prominentes Thema ist die Vernachlässigung von nicht-lokalen Wechselwirkungen in der DMFT. Hierzu kombinieren wir direkte BSS-QMC Gitterrechnungen mit CDMFT für das halb gefüllte zweidimensionale anisotrope Hubbard Modell, das dotierte Hubbard Modell und das KLM. Die Ergebnisse für die verschiedenen Modelle unterscheiden sich stark: während nicht-lokale Korrelationen eine wichtige Rolle im zweidimensionalen (anisotropen) Modell spielen, ist in der paramagnetischen Phase die Impulsabhängigkeit der Selbstenergie für stark dotierte Systeme und für das KLM deutlich schwächer. Eine bemerkenswerte Erkenntnis ist, dass die Selbstenergie sich durch die nicht-wechselwirkende Dispersion parametrisieren lässt. Die spezielle Struktur der Selbstenergie im Impulsraum kann sehr nützlich für die Klassifizierung von elektronischen Korrelationseffekten sein und öffnet den Weg für die Entwicklung neuer Schemata über die Grenzen der DMFT hinaus.

Monte Carlo simulations of nucleation of colloidal crystals

A crystal nucleus in a finite volume may exhibit phase coexistence with a surrounding fluid. The thermodynamic properties of the coexisting fluid (pressure and chemical potential) are enhanced relative to their coexistence values. This enhancement is uniquely related to the surface excess free energy. rnA model for weakly attractive soft colloidal particles is investigated, the so called Asakura-Oosawa model. In simulations, this model allows for the calculation of the pressure in the liquid using the virial formula directly. The phase coexistence pressure in the thermodynamic limit is obtained from the interface velocity method. We introduce a method by which the chemical potential in dense liquids can be measured. There is neither a need to locate the interface nor to compute the anisotropic interfacial tension to obtain nucleation barriers. Therefore, our analysis is appropriate for nuclei of arbitrary shape. Monte Carlo simulations over a wide range of nucleus volumes yield to nucleation barriers independent from the total system volume. The interfacial tension is determined via the ensemble-switch method, hence a detailed test of classical nucleation theory is possible. The anisotropy of the interfacial tension and the resulting non-spherical shape has only a minor effect on the barrier for the Asakura-Oosawa model.

Verifica radioprotezionistica con tecniche Monte Carlo del bunker per radioterapia dell' ASMN-IRCCS di Reggio Emilia

Conoscere con dettaglio il campo di radiazione che si genera nell'utilizzo di un acceleratore lineare di elettroni durante una seduta di radioterapia è essenziale sia per i pazienti sia per gli operatori. L'utilizzo del codice Monte Carlo MCNPX 2.7.0 permette di stimare dati dosimetrici dettagliati in zone dove può essere complicato effettuare misurazioni.Lo scopo di questo lavoro è indagare il comportamento del fascio fotonico prodotto nel bunker di radioterapia dell'ASMN-IRCCS di Reggio Emilia, valutando con precisione in particolare la produzione di fotoneutroni secondari. L'obiettivo è la verifica dell'efficacia delle barriere offerte dalla struttura tenendo in considerazione anche il canale di penetrazione degli impianti di servizio che costituisce un punto di fuga per le radiazioni.

"Descrizione della circolazione dei fluidi e della sequenza di deformazione nella Formazione di Monte Venere"

Il presente lavoro ha come obiettivo quello di descrivere, alla scala dell'affioramento, 1) la sequenza di deformazione fragile: fratture di modo I (joints e vene) e di modo II-III (faglie) e 2) le modalità di circolazione dei fluidi in relazione alle strutture stesse. Da una parte i joint hanno agevolato la circolazione generando dei particolari fronti reattivi (diffusione e ossidazione). Dall'altra vene e faglie sembrano aver agito da barriera, disturbando la circolazione dei fluidi.

Analisi petrografica di inclusi in depositi piroclastici del Monte Vesuvio

Il lavoro svolto mira alla classificazione dei litici presi in carica dal magma durante le eruzioni del Monte Somma-Vesuvio

Realizzazione e convalida del modello Monte Carlo di un calibratore di attività per radiofarmaci

I calibratori di attività sono strumenti molto importanti per la pratica, diagnostica e terapeutica, in medicina nucleare, perché permettono di associare ad un radiofarmaco una misura accurata dell’attività dell’isotopo in esso contenuto; questo è fondamentale in quanto l’attività della sorgente esprime la quantità di farmaco somministrata al paziente. In questo lavoro è stato sviluppato il modello Monte Carlo di un calibratore di attività ampiamente diffuso nei laboratori di radiofarmacia (Capintec CRC-15), utilizzando il codice Monte Carlo FLUKA. Per realizzare il modello si è posta estrema attenzione nel riprodurre al meglio tutti i dettagli delle componenti geometriche della camera e dei campioni delle sorgenti radioattive utilizzati. A tale scopo, la camera di ionizzazione di un calibratore è stata studiata mediante imaging TAC. Un’analisi preliminare è stata eseguita valutando il confronto tra l’andamento sperimentale dell’efficienza della camera in funzione dell’energia dei fotoni incidenti e quello ottenuto in simulazione. In seguito si è proceduto con la validazione: si sono studiati a questo proposito la risposta del calibratore in funzione dell’altezza della sorgente e i confronti tra i fattori relativi (rispetto ad una sorgente certificata di 137Cs) e le misure di confronto sono state eseguite con diverse sorgenti certificate di 133Ba, 68Ge-68Ga, 177Lu ed uno standard tarato internamente di 99mTc. In tale modo, si è ricoperto l'intero campo di interesse dei principali radionuclidi impiegati nelle applicazioni diagnostiche e terapeutiche di Medicina Nucleare. Il modello sviluppato rappresenta un importante risultato per l’eventuale determinazione di nuovi fattori di calibrazione o per un futuro studio relativo all’ottimizzazione della risposta del calibratore.

Ricostruzione di segnali multipli del rivelatore a tempo di volo (ToF) di ALICE: confronto fra dati e Monte Carlo

La fisica delle collisioni ad alte energie è, ad oggi, uno dei campi di ricerca più interessante per la verifica di modelli teorici che spieghino la nascita e la formazione dell'universo in cui viviamo. In quest'ottica lavorano esperimenti presso i più importanti acceleratori di particelle: tra questi anche l'esperimento ALICE, presso il Large Hadron Collider LHC del CERN di Ginevra. Il suo scopo principale è quello di verificare e ampliare l'insieme delle prove sperimentali alla base sull'esistenza di un nuovo stato della materia: il Quark Gluon Plasma. La presenza della transizione di fase della materia adronica ordinaria a QGP era stata teorizzata da diversi studi termodinamici e da calcoli di QCD su reticolo: in particolare si prevedeva l'esistenza di uno stato della materia in cui i quark sono deconfinati. Il QGP è dunque un plasma colorato e densissimo di quark e gluoni, liberi di interagire tra loro. Queste condizioni sarebbero state quelle dell'universo primordiale, circa 1µs dopo il Big Bang: a seguito di una transizione di fase, si sarebbe poi formata tutta la materia adronica ordinaria. Per riprodurre le condizioni necessarie alla formazione del QGP occorrono collisioni ad energie ultrarelativistiche come quelle prodotte, negli ultimi anni, dall'acceleratore LHC. Uno dei principali rivelatori dedicati all'identificazione di particelle in ALICE è il sistema a tempo di volo TOF. Nonostante un attento processo di ottimizzazione della risoluzione delle sue componenti, persistono residui errori strumentali che limitano la qualità (già ottima) del segnale, tutt'oggi oggetto di studio. L'elaborato presentato in questa tesi è suddiviso in tre capitoli: nel primo ripercorriamo gli aspetti teorici del Modello Standard e del Quark Gluon Plasma. Nel secondo descriviamo la struttura di rivelazione di ALICE, analizzando il funzionamento delle singole componenti. Nel terzo, infine, verifichiamo le principali correzioni al TOF ad oggi note, confrontando i dati a nostra disposizione con delle simulazioni Monte Carlo: questo ci permette da un lato di approfondirne la conoscenza, dall'altro di cercarne di migliorare la descrizione del comportamento del rivelatore.

Studio idrogeologico per la progettazione e l'esecuzione di un tracciamento multiplo dell'acquifero carsico del Massiccio del Cansiglio-Monte Cavallo (alta valle Livenza, Pordenone)

Il Massiccio del Cansiglio-Monte Cavallo costituisce uno dei maggiori complessi carsici italiani; con questo lavoro si cerca di fornirne una caratterizzazione idrogeologica sulla quale poter basare una futura prova di tracciamento multiplo, al fine di definirne con precisione la configurazione interna. Lo studio è incentrato sul monitoraggio in continuo di portate, temperature e conducibilità elettrica compensata a 25 °C delle tre sorgenti maggiori del fiume Livenza, che sgorgano alla base del Massiccio, e su un monitoraggio idrochimico discontinuo di nove sorgenti (tra cui figurano anche le tre principali) e di quattro grotte rappresentanti le aree di infiltrazione concentrata. L’enorme volume d’acqua contenuto in questo grande acquifero carsico sembra riversarsi quasi completamente nelle aree sorgive situate al margine occidentale della pianura friulana, che cinge il lato orientale del Massiccio. In linea generale, queste acque, tipicamente carbonatiche, non hanno tempi di residenza lunghi all’interno del bacino di provenienza e fluiscono con portate elevate (> 1 m3/s) da ognuna delle tre sorgenti maggiori; nel periodo estivo (luglio-agosto) si registra una fase di svuotamento, assai accentuata al Gorgazzo, sorgente, tra le tre, posta più a nord. Molinetto, situata più a sud, pare differenziarsi maggiormente dalle altre due grandi sorgenti, Santissima e Gorgazzo, mentre le sorgenti minori sembrano, al di là di alcune caratteristiche peculiari, abbastanza simili tra loro. In particolare, il Molinetto mostra, in risposta agli eventi meteorici intensi, un incremento di portata meno netto, accompagnato da un evidente fenomeno di pistonaggio, almeno durante il periodo di morbida. Al contrario, il Gorgazzo risponde immediatamente alle forti precipitazioni piovose e le acque di pioggia viaggiano rapidamente dalla zona di infiltrazione alla sorgente stessa. La Santissima, sorgente con valori medi di portata maggiori e valori di conducibilità elettrica e temperatura mediamente inferiori alle altre, sembra dipendere da un sistema con caratteristiche intermedie tra quelle delle due sorgenti adiacenti. Emerge la complessità di questo vasto sistema carsico che sembra mostrare un grado di carsificazione in aumento da sud-ovest a nord-est. Sulla base dei dati raccolti si suppone che al suo interno possano svilupparsi due sottosistemi: il primo sembra che riversi la sua riserva idrica nelle sorgenti Molinetto e Santissima, il secondo accomuna invece la Santissima con il Gorgazzo, il quale ne costituisce un “troppo pieno”.

Studio mineralogico e chimico della serie della "Vena del Gesso" della cava di Monte Tondo (Borgo Rivola, Riolo Terme - RA)

Questo elaborato consiste nello studio mineralogico e chimico di una parte della successione della Vena del Gesso, coincidente con il sito estrattivo di Monte Tondo. Oltre all'inquadramento geologico della zona oggetto di studio e alcuni accenni sulle attività della cava, il presente elaborato contiene tabelle con le analisi di fluorescenza XRF per tutti i campioni prelevati, grafici dello spettro XRD per due campioni, completi di tabelle con i valori in allegato, diagrammi di correlazione per alcuni elementi, una analisi termogravimetrica ed una termodifferenziale, e i grafici con le rette di regressione per gli elementi maggiori e in tracce, con l'esclusione di certi elementi che hanno dato concentrazioni sotto il limite di rilevabilità.

Monte Carlo dose calculation improvements for low energy electron beams using eMC

The electron Monte Carlo (eMC) dose calculation algorithm in Eclipse (Varian Medical Systems) is based on the macro MC method and is able to predict dose distributions for high energy electron beams with high accuracy. However, there are limitations for low energy electron beams. This work aims to improve the accuracy of the dose calculation using eMC for 4 and 6 MeV electron beams of Varian linear accelerators. Improvements implemented into the eMC include (1) improved determination of the initial electron energy spectrum by increased resolution of mono-energetic depth dose curves used during beam configuration; (2) inclusion of all the scrapers of the applicator in the beam model; (3) reduction of the maximum size of the sphere to be selected within the macro MC transport when the energy of the incident electron is below certain thresholds. The impact of these changes in eMC is investigated by comparing calculated dose distributions for 4 and 6 MeV electron beams at source to surface distance (SSD) of 100 and 110 cm with applicators ranging from 6 x 6 to 25 x 25 cm(2) of a Varian Clinac 2300C/D with the corresponding measurements. Dose differences between calculated and measured absolute depth dose curves are reduced from 6% to less than 1.5% for both energies and all applicators considered at SSD of 100 cm. Using the original eMC implementation, absolute dose profiles at depths of 1 cm, d(max) and R50 in water lead to dose differences of up to 8% for applicators larger than 15 x 15 cm(2) at SSD 100 cm. Those differences are now reduced to less than 2% for all dose profiles investigated when the improved version of eMC is used. At SSD of 110 cm the dose difference for the original eMC version is even more pronounced and can be larger than 10%. Those differences are reduced to within 2% or 2 mm with the improved version of eMC. In this work several enhancements were made in the eMC algorithm leading to significant improvements in the accuracy of the dose calculation for 4 and 6 MeV electron beams of Varian linear accelerators.

Comparison of monte carlo collimator transport methods for photon treatment planning in radiotherapy

The aim of this work was a Monte Carlo (MC) based investigation of the impact of different radiation transport methods in collimators of a linear accelerator on photon beam characteristics, dose distributions, and efficiency. Thereby it is investigated if it is possible to use different simplifications in the radiation transport for some clinical situations in order to save calculation time.

Monte Carlo Dose Calculation on Deforming Anatomy

This article presents the implementation and validation of a dose calculation approach for deforming anatomical objects. Deformation is represented by deformation vector fields leading to deformed voxel grids representing the different deformation scenarios. Particle transport in the resulting deformed voxels is handled through the approximation of voxel surfaces by triangles in the geometry implementation of the Swiss Monte Carlo Plan framework. The focus lies on the validation methodology which uses computational phantoms representing the same physical object through regular and irregular voxel grids. These phantoms are chosen such that the new implementation for a deformed voxel grid can be compared directly with an established dose calculation algorithm for regular grids. Furthermore, separate validation of the aspects voxel geometry and the density changes resulting from deformation is achieved through suitable design of the validation phantom. We show that equivalent results are obtained with the proposed method and that no statistically significant errors are introduced through the implementation for irregular voxel geometries. This enables the use of the presented and validated implementation for further investigations of dose calculation on deforming anatomy.