The Two Phases of Galaxy Formation

Cosmological simulations of galaxy formation appear to show a two-phase character with a rapid early phase at z>2 during which in-situ stars are formed within the galaxy from infalling cold gas followed by an extended phase since z<3 during which ex-situ stars are primarily accreted. In the latter phase massive systems grow considerably in mass and radius by accretion of smaller satellite stellar systems formed at quite early times (z>3) outside of the virial radius of the forming central galaxy. These tentative conclusions are obtained from high resolution re-simulations of 39 individual galaxies in a full cosmological context with present-day virial halo masses ranging from 7e11 M_sun h^-1 < M_vir < 2.7e13 M_sun h^-1 and central galaxy masses between 4.5e10 M_sun h^-1 < M_* < 3.6e11 M_sun h^-1. The simulations include the effects of a uniform UV background, radiative cooling, star formation and energetic feedback from SNII. The importance of stellar accretion increases with galaxy mass and towards lower redshift. In our simulations lower mass galaxies ($M_* < 9e10 M_sun h^-1) accrete about 60 per cent of their present-day stellar mass. High mass galaxy ($M_* > 1.7e11 M_sun h^-1) assembly is dominated by accretion and merging with about 80 per cent of the stars added by the present-day. In general the simulated galaxies approximately double their mass since z=1. For massive systems this mass growth is not accompanied by significant star formation. The majority of the in-situ created stars is formed at z>2, primarily out of cold gas flows. We recover the observational result of archaeological downsizing, where the most massive galaxies harbor the oldest stars. We find that this is not in contradiction with hierarchical structure formation. Most stars in the massive galaxies are formed early on in smaller structures, the galaxies themselves are assembled late.Comment: 13 pages, 13 figures, accepted for publication in Ap

Sismicità all’Etna dal 1989 al 2010: evidenze sull’evoluzione spazio-temporale dell’attività sismica

Il Monte Etna, uno dei più attivi vulcani basaltici tra i più monitorati al mondo, è sede di una notevole attività sismica e vulcanica. Esso è ubicato in Sicilia orientale in un complesso quadro geodinamico, dove le principali strutture tettoniche regionali giocano un ruolo chiave nei processi dinamici del vulcano. La sismicità dell’Etna si manifesta con un elevato rate di terremoti di bassa e moderata energia che, a volte, a causa dell’estrema superficialità della sorgente, provocano danni ai centri abitati prossimi all’area epicentrale. Il monitoraggio sistematico dell’attività sismica etnea è effettuato sin dal 1989, mediante una rete sismica locale permanente che nel tempo è stata oggetto di importanti miglioramenti. La prima configurazione di rete era costituita da circa 10 stazioni analogiche con sensori a corto periodo gestita dall’Istituto Internazionale di Vulcanologia (IIV-CNR). Nel 1994, una rete sismica costituita da circa 40 stazioni (analogiche con sensori a corto periodo) fu installata sull’Etna nell’ambito del Progetto Poseidon. Nel 2001, le reti gestite dall’IIV-CNR e dal Progetto Poseidon confluirono nell’Istituto Nazionale di Geofisica e Vulcanologia (INGV); attualmente la rete sismica, costituita da circa 50 stazioni digitali equipaggiate con sismometri broadband a tre componenti, è gestita dalla Sezione di Catania dell’INGV. Nel periodo 1989-1999, il catalogo dei terremoti risulta costituito da circa 2000 eventi con soglia di completezza per magnitudo pari a 2.0; dal 1999 ad oggi contiene circa 6000 terremoti con soglia di completezza per magnitudo 1.5. La capacità di detezione della rete è migliorata nel tempo permettendo di registrare e localizzare anche gli eventi meno energetici (M≥1.0). In questo lavoro, vengono presentati i caratteri predominanti della sismicità etnea negli ultimi 20 anni, con un maggiore dettaglio della distribuzione spazio-temporale della sismicità verificatasi dal 1999. L’analisi della attività sismica rappresenta un utile strumento per l’interpretazione delle dinamiche che hanno contraddistinto numerose ed importanti eruzioni (2001, 2002-03, 2004, 2006, 2008-09). In particolare, la variazione del rilascio energetico della sismicità ha contribuito in maniera significativa ad identificare i probabili processi geodinamici legati alla ricarica del sistema magmatico del vulcano. La distribuzione spaziale della sismicità ha consentito di evidenziare inoltre l’esistenza di diverse aree sismogenetiche caratterizzate da un differente rate sismico, profondità focali e cinematica delle strutture associate. Infine, osservando le caratteristiche della sismicità nel lungo periodo, differenti settori del vulcano sono risultati maggiormente attivi in relazione ai più importanti recenti eventi eruttivi

Structural features of the Pernicana Fault (M. Etna, Sicily, Italy) inferred by high precise location of the microseismicity

The north-eastern flank of Mt. Etna is crossed by an important and active tectonic structure, the Pernicana Fault having a mean strike WNW–ESE. It links westward to the active NE Rift and seems to have an important role in controlling instability processes affecting the eastern flank of the volcano. Recent studies suggest that Pernicana Fault is very active through sinistral, oblique-slip movements and is also characterised by frequent shallow seismicity (depth < 2 km bsl) on the uphill western segment and by remarkable creeping on the downhill eastern one. The Pernicana Fault earthquakes, which can reach magnitudes up to 4.2, sometimes with coseismic surface faulting, caused severe damages to tourist resorts and villages along or close this structure. In the last years, a strong increase of seismicity, also characterized by swarms, was recorded by INGV-CT permanent local seismic network close the Pernicana Fault. A three-step procedure was applied to calculate precise hypocentre locations. In a first step, we chose to apply cross-correlation analysis, in order to easily evaluate the similarity of waveforms useful to identify earthquakes families. In a second step, we calculate probabilistic earthquake locations using the software package NONLINLOC, which includes systematic, complete grid search and global, non-linear search methods. Subsequently, we perform relative relocation of correlated event pairs using the double-difference earthquake algorithm and the program HypoDD. The double-difference algorithm minimizes the residuals between observed and calculated travel time difference for pairs of earthquakes at common stations by iteratively adjusting the vector difference between the hypocenters. We show the recognized spatial seismic clusters identifying the most active and hazarding sectors of the structure, their geometry and depth. Finally, in order to clarify the geodynamic framework of the area, we associate these results with calculated focal mechanisms for the most energetic earthquakes

EEG alpha power is modulated by attentional changes during cognitive tasks and virtual reality immersion

Variations in alpha rhythm have a significant role in perception and attention. Recently, alpha decrease has been associated with externally directed attention, especially in the visual domain, whereas alpha increase has been related to internal processing such as mental arithmetic. However, the role of alpha oscillations and how the different components of a task (processing of external stimuli, internal manipulation/representation, and task demand) interact to affect alpha power are still unclear. Here, we investigate how alpha power is differently modulated by attentional tasks depending both on task difficulty (less/more demanding task) and direction of attention (internal/external). To this aim, we designed two experiments that differently manipulated these aspects. Experiment 1, outside Virtual Reality (VR), involved two tasks both requiring internal and external attentional components (intake of visual items for their internal manipulation) but with different internal task demands (arithmetic vs. reading). Experiment 2 took advantage of the VR (mimicking an aircraft cabin interior) to manipulate attention direction: it included a condition of VR immersion only, characterized by visual external attention, and a condition of a purely mental arithmetic task during VR immersion, requiring neglect of sensory stimuli. Results show that: (1) In line with previous studies, visual external attention caused a significant alpha decrease, especially in parieto-occipital regions; (2) Alpha decrease was significantly larger during the more demanding arithmetic task, when the task was driven by external visual stimuli; (3) Alpha dramatically increased during the purely mental task in VR immersion, whereby the external stimuli had no relation with the task. Our results suggest that alpha power is crucial to isolate a subject from the environment, and move attention from external to internal cues. Moreover, they emphasize that the emerging use of VR associated with EEG may have important implications to study brain rhythms and support the design of artificial systems

dfpk : An R-package for Bayesian dose-finding designs using Pharmacokinetics (PK) for phase I clinical trials

Background and objective Dose-finding, aiming at finding the maximum tolerated dose, and pharmacokinetics studies are the first in human studies in the development process of a new pharmacological treatment. In the literature, to date only few attempts have been made to combine pharmacokinetics and dose-finding and to our knowledge no software implementation is generally available. In previous papers, we proposed several Bayesian adaptive pharmacokinetics-based dose-finding designs in small populations. The objective of this work is to implement these dose-finding methods in an R package, called dfpk. Methods All methods were developed in a sequential Bayesian setting and Bayesian parameter estimation is carried out using the rstan package. All available pharmacokinetics and toxicity data are used to suggest the dose of the next cohort with a constraint regarding the probability of toxicity. Stopping rules are also considered for each method. The ggplot2 package is used to create summary plots of toxicities or concentration curves. Results For all implemented methods, dfpk provides a function (nextDose) to estimate the probability of efficacy and to suggest the dose to give to the next cohort, and a function to run trial simulations to design a trial (nsim). The sim.data function generates at each dose the toxicity value related to a pharmacokinetic measure of exposure, the AUC, with an underlying pharmacokinetic one compartmental model with linear absorption. It is included as an example since similar data-frames can be generated directly by the user and passed to nsim. Conclusion The developed user-friendly R package dfpk, available on the CRAN repository, supports the design of innovative dose-finding studies using PK information

Caratterizzazione sismica del sistema strutturale Pernicana - Provenzana (settore NE dell'Etna) attraverso l'utilizzo di differenti tecniche di rilocalizzazione.

Il fianco nord-orientale dell’Etna è interessato da un noto sistema strutturale denominato Pernicana-Provenzana, che ha un andamento WNW–ESE. Esso è collegato ad ovest ad un altro importante elemento strutturale, il Rift di Nord-Est, che mostra avere un ruolo importante nel controllo dei fenomeni di instabilità del fianco orientale del vulcano. La sismicità associata a questo sistema strutturale è di tipo superficiale (max 2-3 km b.s.l.) e rilevanti fenomeni di creeping sono rilevabili sul suo segmento orientale. I terremoti associati a questo sistema di faglie, che possono raggiungere magnitudo sino a 4.3, qualche volta con fenomeni di fagliazione superficiale, hanno provocato danni importanti alle principali strutture alberghiere ed ai paesi ubicati in prossimità della struttura. Nel presente lavoro, sono riportati i risultati di uno studio di dettaglio della sismicità localizzata lungo tale sistema strutturale, nel periodo 1999-2009. I terremoti registrati dalla rete sismica permanente dell’Istituto Nazionale di Geofisica e Vulcanologia – Sezione di Catania, localizzati con un modello 1D utilizzando l’algoritmo Hypoellipse (Gruppo Analisi Dati Sismici, 2010), sono stati rilocalizzati applicando due differenti tecniche di localizzazione: NonLinLoc sviluppato da Lomax et al. (2000) e HypoDD proposto da Waldhauser & Ellsworth (2000). La prima metodologia è basata su un processo di ricerca globale, nello spazio 3D, dei parametri di localizzazione che possono essere ottenuti utilizzando diversi algoritmi. Il metodo HypoDD, che non prevede l’utilizzo di un modello 3D, è invece basato sull’algoritmo della doppia differenza che minimizza i residui tra le differenze dei traveltime osservati e calcolati per coppie di terremoti a stazioni comuni. L’applicazione di tali tecniche ha permesso di ottenere localizzazioni ipocentrali di migliore qualità, fondamentali per la caratterizzazione sismica della struttura. L’applicazione di queste differenti metodologie ha permesso di evidenziare che il sistema strutturale Pernicana- Provenzana risulta composto da segmenti caratterizzati da differenti rilasci di energia sismica. Sono stati individuati due cluster principali di terremoti, la cui distribuzione spaziale ha evidenziato un differente verso nell’immersione dei piani di faglia collegabili a questa sismicità. Infine, l’applicazione di tecniche di cross-correlazione delle forme d’onda registrate nel periodo indagato ha consentito di individuare “famiglie” di terremoti. L’analisi spazio – temporale delle famiglie individuate ha evidenziato per alcune di esse, una ricorrenza temporale ed ha permesso di ipotizzare che l’applicazione di un campo di stress sul sistema Pernicana-Provenzana potrebbe essere capace di attivare le stesse sorgenti sismiche in differenti periodi

Changes in brain rhythms and connectivity tracking fear acquisition and reversal

Fear conditioning is used to investigate the neural bases of threat and anxiety, and to understand their flexible modifications when the environment changes. This study aims to examine the temporal evolution of brain rhythms using electroencephalographic signals recorded in healthy volunteers during a protocol of Pavlovian fear conditioning and reversal. Power changes and Granger connectivity in theta, alpha, and gamma bands are investigated from neuroelectrical activity reconstructed on the cortex. Results show a significant increase in theta power in the left (contralateral to electrical shock) portion of the midcingulate cortex during fear acquisition, and a significant decrease in alpha power in a broad network over the left posterior-frontal and parietal cortex. These changes occur since the initial trials for theta power, but require more trials (3/4) to develop for alpha, and are also present during reversal, despite being less pronounced. In both bands, relevant changes in connectivity are mainly evident in the last block of reversal, just when power differences attenuate. No significant changes in the gamma band were detected. We conclude that the increased theta rhythm in the cingulate cortex subserves fear acquisition and is transmitted to other cortical regions via increased functional connectivity allowing a fast theta synchronization, whereas the decrease in alpha power can represent a partial activation of motor and somatosensory areas contralateral to the shock side in the presence of a dangerous stimulus. In addition, connectivity changes at the end of reversal may reflect long-term alterations in synapses necessary to reverse the previously acquired contingencies

Intracranial pressure dynamics in patients with acute brain damage

The time pattern of intracranial pressure (ICP) during pressure-volume index (PVI) tests was analyzed in 20 patients with severe acute brain damage by means of a simple mathematical model. In most cases, a satisfactory fitting between model response and patient data was achieved by adjusting only four parameters: the cerebrospinal fluid (CSF) outflow resistance, the intracranial elastance coefficient, and the gain and time constant of cerebral autoregulation. The correlation between the parameter estimates was also analyzed to elucidate the main mechanisms responsible for ICP changes in each patient. Starting from information on the estimated parameter values and their correlation, the patients were classified into two main classes: those with weak autoregulation (8 of 20 patients) and those with strong autoregulation (12 of 20 patients). In the first group of patients, ICP mainly reflects CSF circulation and passive cerebral blood volume changes. In the second group, ICP exhibits paradoxical responses attributable to active changes in cerebral blood volume. Moreover, in two patients of the second group, the time constant of autoregulation is significantly increased (>40 s). The correlation between the parameter estimates was significantly different in the two groups of patients, suggesting the existence of different mechanisms responsible for ICP changes. Moreover, analysis of the correlation between the parameter estimates might give information on the directions of parameter changes that have a greater impact on ICP