We explore the chemodynamical properties of a sample of barred galaxies in the Auriga magneto-hydrodynamical cosmological zoom-in simulations, which form boxy/peanut (b/p) bulges, and compare these to the Milky Way (MW). We show that the Auriga galaxies which best reproduce the chemodynamical properties of stellar populations in the MW bulge have quiescent merger histories since redshift $z\sim3.5$: their last major merger occurs at $t_{\rm lookback}>12\,\rm Gyrs$, while subsequent mergers have a stellar mass ratio of $\leq$1:20, suggesting an upper limit of a few percent for the mass ratio of the recently proposed Gaia Sausage/Enceladus merger. These Auriga MW-analogues have a negligible fraction of ex-situ stars in the b/p region ($<1\%$), with flattened, thick disc-like metal-poor stellar populations. The average fraction of ex-situ stars in the central regions of all Auriga galaxies with b/p's is 3% -- significantly lower than in those which do not host a b/p or a bar. While the central regions of these barred galaxies contain the oldest populations, they also have stars younger than 5Gyrs (>30%) and exhibit X-shaped age and abundance distributions. Examining the discs in our sample, we find that in some cases a star-forming ring forms around the bar, which alters the metallicity of the inner regions of the galaxy. Further out in the disc, bar-induced resonances lead to metal-rich ridges in the $V_{\phi}-r$ plane -- the longest of which is due to the Outer Lindblad Resonance. Our results suggest the Milky Way has an uncommonly quiet merger history, which leads to an essentially in-situ bulge, and highlight the significant effects the bar can have on the surrounding disc.

We use magneto-hydrodynamical simulations of Milky Way-mass haloes from the Auriga project to examine the properties of surviving and destroyed dwarf galaxies that are accreted by these haloes over cosmic time. We show that the combined luminosity function of surviving and destroyed dwarfs at infall is similar in the various Auriga haloes, and is dominated by the destroyed dwarfs. There is, however, a strong dependence on infall time: destroyed dwarfs have typically early infall times, $t_{infall}<6$ Gyr, whereas the majority of dwarfs accreted at $t_{infall}>10$ Gyr have survived to the present day. Because of their late infall the surviving satellites today had higher metallicites at infall than their destroyed counterparts of similar infall mass; the difference is even more pronounced for the present-day metallicites of satellites, many of which continue to form stars after infall. In agreement with previous work, we find that a small number of relatively massive destroyed dwarf galaxies dominate the mass of the stellar haloes. However, there is a significant radial dependence: while 90 per cent of the mass in the inner regions ($<\,20\,$kpc) is contributed, on average, by only 3 massive progenitors, the outer regions ($>\,100\,$kpc) typically have $\sim8$ main progenitors of relatively lower mass. Finally, we show that a few massive progenitors dominate the metallicity distribution of accreted stars, even at the metal poor end. Contrary to common assumptions in the literature, dwarf galaxies of mass $M_{*}<10^7 \, M_{\odot}$ make up less than 10 per cent of the accreted, metal poor stars ([Fe/H] $<\,-3$) in the inner $50\,$kpc.

The circumgalactic medium (CGM) is one of the frontiers of galaxy formation and intimately connected to the galaxy via accretion of gas on to the galaxy and gaseous outflows from the galaxy. Here we analyse the magnetic field in the CGM of the Milky Way-like galaxies simulated as part of the \textsc{Auriga} project that constitutes a set of high resolution cosmological magnetohydrodynamical zoom simulations. We show that at high redshift the CGM becomes magnetised via galactic outflows that transport magnetised gas from the disk into the halo. At this time the magnetisation of the CGM closely follows its metal enrichment. We then show that at low redshift an in-situ turbulent dynamo that operates on a timescale of Gigayears further amplifies the magnetic field in the CGM and sets its strength down $z=0$. The magnetic field strength reaches a typical value of $0.1\,\mu G$ at the virial radius at $z=0$ and becomes mostly uniform within the virial radius. Its Faraday rotation signal is in excellent agreement with recent observations. For most of its evolution the magnetic field in the CGM is an unordered small scale field. Only strong coherent outflows at low redshift are able to order the magnetic field in parts of the CGM that are directly displaced by these outflows.

We use the Auriga cosmological simulations of Milky Way (MW)-mass galaxies and their surroundings to study the satellite populations of dwarf galaxies in $\Lambda$CDM. As expected from prior work, the number of satellites above a fixed stellar mass is a strong function of the mass of the primary dwarf. For galaxies as luminous as the Large Magellanic Cloud (LMC), and for halos as massive as expected for the LMC (determined by its rotation speed), the simulations predict about 3 satellites with stellar masses exceeding $M_*>10^5\, M_\odot$. If the LMC is on its first pericentric passage, then these satellites should be near the LMC and should have orbital angular momenta roughly coincident with that of the LMC. We use 3D positions and velocities from the 2nd data release of the Gaia mission to revisit which of the "classical" MW dwarf spheroidals could plausibly be LMC satellites. The new proper motions of the Fornax and Carina dwarf spheroidals place them on orbits closely aligned with the orbital plane of the Magellanic Clouds, hinting at a potential Magellanic association. Together with the Small Magellanic Cloud (SMC), this result raises to $3$ the number of LMC satellites with $M_*>10^5\, M_\odot$, as expected from simulations. This also fills the 12-mag luminosity gap between the SMC and the ultra-faints Hyi1, Car2, Hor1, and Car3, the few ultra-faint satellites confirmed to have orbits consistent with a Magellanic origin.

We investigate the formation of ultra-diffuse galaxies (UDGs) using the Auriga high-resolution cosmological magneto-hydrodynamical simulations of Milky Way-sized galaxies. We identify a sample of $92$ UDGs in the simulations that match a wide range of observables such as sizes, central surface brightness, S\'{e}rsic indices, colors, spatial distribution and abundance. Auriga UDGs have dynamical masses similar to normal dwarfs. In the field, the key to their origin is a strong correlation present in low-mass dark matter haloes between galaxy size and halo spin parameter. Field UDGs form in dark matter haloes with larger spins compared to normal dwarfs in the field, in agreement with previous semi-analytical models. Satellite UDGs, on the other hand, have two different origins: $\sim 55\%$ of them formed as field UDGs before they were accreted; the remaining $\sim 45\%$ were normal field dwarfs that subsequently turned into UDGs as a result of tidal interactions.

Galaxies are surrounded by massive gas reservoirs (i.e. the circumgalactic medium; CGM) which play a key role in their evolution. The properties of the CGM, which are dependent on a variety of internal and environmental factors, are often inferred from absorption line surveys which rely on a limited number of single lines-of-sight. In this work we present an analysis of 28 galaxy haloes selected from the Auriga project, a cosmological magneto-hydrodynamical zoom-in simulation suite of isolated Milky Way-mass galaxies, to understand the impact of CGM diversity on observational studies. Although the Auriga haloes are selected to populate a narrow range in halo mass, our work demonstrates that the CGM of L* galaxies is extremely diverse: column densities of commonly observed species span ~3-4 dex and their covering fractions range from ~5 to 90 per cent. Despite this diversity, we identify the following correlations: 1) the covering fractions (CF) of hydrogen and metals of the Auriga haloes positively correlate with stellar mass, 2) the CF of H I, C IV, and Si II anticorrelate with active galactic nucleus luminosity due to ionization effects, and 3) the CF of H I, C IV, and Si II positively correlate with galaxy disc fraction due to outflows populating the CGM with cool and dense gas. The Auriga sample demonstrates striking diversity within the CGM of L* galaxies, which poses a challenge for observations reconstructing CGM characteristics from limited samples, and also indicates that long-term merger assembly history and recent star formation are not the dominant sculptors of the CGM.

N-body simulations make unambiguous predictions for the abundance of substructures within dark matter halos. However, the inclusion of baryons in the simulations changes the picture because processes associated with the presence of a large galaxy in the halo can destroy subhalos and substantially alter the mass function and velocity distribution of subhalos. We compare the effect of galaxy formation on subhalo populations in two state-of-the-art sets of hydrodynamical CDM simulations of Milky Way mass halos, APOSTLE and AURIGA. We introduce a new method for tracking the orbits of subhalos between simulation snapshots that gives accurate results down to a few kiloparsecs from the centre of the halo. Relative to a dark matter-only simulation, the abundance of subhalos in APOSTLE is reduced by 50% near the centre and by 10% within r200. In AURIGA the corresponding numbers are 80% and 40%. The velocity distributions of subhalos are also affected by the presence of the galaxy, much more so in AURIGA than in APOSTLE . The differences on subhalo properties in the two simulations can be traced back to the mass of the central galaxies, which in AURIGA are typically twice as massive as those in APOSTLE . We show that some of the results from previous studies are inaccurate due to systematic errors in the modelling of subhalo orbits near the centre of halos.

Measurements of the rotation curves of dwarf galaxies are often interpreted as requiring a constant density core at the centre, at odds with the "cuspy" inner profiles predicted by $N$-body simulations of cold dark matter (CDM) haloes. It has been suggested that this conflict could be resolved by fluctuations in the inner gravitational potential caused by the periodic removal of gas following bursts of star formation. Earlier work has suggested that core formation requires a bursty and extended star formation history (SFH). Here we investigate the structure of CDM haloes of dwarf galaxies ($M_{{\rm DM}} \sim 10^9-5\times10^{10}\,{\rm M}_\odot$) formed in the APOSTLE ('A Project of Simulating the Local Environment') and AURIGA cosmological hydrodynamic simulations. Our simulations have comparable or better resolution than others that make cores ($M_{{\rm gas}} \sim 10^4\,{\rm M}_\odot$, gravitational softening $\sim 150$ pc). Yet, we do not find evidence of core formation at {\it any} mass or any correlation between the inner slope of the DM density profile and temporal variations in the SFH. APOSTLE and AURIGA dwarfs display a similar diversity in their cumulative SFHs to available data for Local Group dwarfs. Dwarfs in both simulations are DM-dominated on all resolved scales at all times, likely limiting the ability of gas outflows to alter significantly the central density profiles of their haloes. We conclude that recurrent bursts of star formation are not sufficient to cause the formation of cores, and that other conditions must also be met for baryons to be able to modify the central DM cusp.

We use the APOSTLE and Auriga cosmological simulations to study the star formation histories (SFHs) of field and satellite dwarf galaxies. Despite sizeable galaxy-to-galaxy scatter, the SFHs of APOSTLE and Auriga dwarfs exhibit robust average trends with galaxy stellar mass: faint field dwarfs ($10^5<M_{\rm star}/M_\odot<10^{6.5}$) have, on average, steadily declining SFHs, whereas brighter dwarfs ($10^{7.5}<M_{\rm star}/M_\odot<10^{9}$) show the opposite trend. Intermediate-mass dwarfs have roughly constant SFHs. Satellites exhibit similar average trends, but with substantially suppressed star formation in the most recent $\sim 5$ Gyr, likely as a result of gas loss due to tidal and ram-pressure stripping after entering the haloes of their primaries. These simple mass and environmental trends are in good agreement with the derived SFHs of Local Group (LG) dwarfs whose photometry reaches the oldest main sequence turnoff. SFHs of galaxies with less deep data show deviations from these trends, but this may be explained, at least in part, by the large galaxy-to-galaxy scatter, the limited sample size, and the large uncertainties of the inferred SFHs. Confirming the predicted mass and environmental trends will require deeper photometric data than currently available, especially for isolated dwarfs.

The highly r-process enhanced (r-II) metal-poor halo stars we observe today could play a key role in understanding early ultra-faint dwarf galaxies, the smallest building blocks of the Milky Way. If a significant fraction of metal-poor r-II halo stars originated in the ultra-faint dwarf galaxies that merged to help form the Milky Way, observations of r-II stars could help us study these now-destroyed systems and probe the formation history of our Galaxy. To conduct our initial investigation into this possible connection, we use high-resolution cosmological simulations of Milky-Way-mass galaxies from the Caterpillar suite in combination with a simple, empirically motivated treatment of r-process enrichment. We determine the fraction of metal-poor halo stars that could have formed from highly r-process enhanced gas in now-destroyed low-mass ultra-faint dwarf galaxies, the simulated r-II fraction, and compare it to the "as observed" r-II fraction. We find that the simulated fraction, f_{r-II,sim} ~ 1-2%, can account for around half of the "as observed" fraction, f_{r-II,obs} ~ 2-4%. The "as observed" fraction likely overrepresents the fraction of r-II stars due to incomplete sampling, though, meaning f_{r-II,sim} likely accounts for more than half of the true f_{r-II,obs}. Further considering some parameter variations and scatter between individual simulations, the simulated fraction can account for around 20-80% of the "as observed" fraction.

We present the physical characterization of FSR 1758, a new large, massive object very recently discovered in the Galactic Bulge. The combination of optical data from the 2nd Gaia Data Release (GDR2) and the DECam Plane Survey (DECaPS), and near-IR data from the VISTA Variables in the V\'{\i}a L\'actea Extended Survey (VVVX) led to a clean sample of likely members. Based on this integrated dataset, position, distance, reddening, size, metallicity, absolute magnitude, and proper motion of this object are measured. We estimate the following parameters: $\alpha=17:31:12$, $\delta=-39:48:30$ (J2000), $D=11.5 \pm 1.0$ kpc, $E(J-Ks)=0.20 \pm 0.03$ mag, $R_c=10$ pc, $R_t=150$ pc, $[Fe/H]=-1.5 \pm 0.3$ dex, $M_i < -8.6 \pm 1.0$, $\mu_{\alpha} = -2.85$ mas yr$^{-1}$, and $\mu_{\delta} = 2.55$ mas yr$^{-1}$. The nature of this object is discussed. If FRS 1758 is a genuine globular cluster, it is one of the largest in the Milky Way, with a size comparable or even larger than that of $\omega$ Cen, being also an extreme outlier in the size vs. Galactocentric distance diagram. The presence of a concentration of long-period RR Lyrae variable stars and blue horizontal branch stars suggests that it is a typical metal-poor globular cluster of Oosterhoff type II. Further exploration of a larger surrounding field reveals common proper motion stars, suggesting either tidal debris or that FRS\,1758 is actually the central part of a larger extended structure such as a new dwarf galaxy, tentatively named as Scorpius. In either case, this object is remarkable, and its discovery graphically illustrates the possibility to find other large objects hidden in the Galactic Bulge using future surveys.

Using the astrometry from the ESA's Gaia mission, previous works have shown that the Milky Way stellar halo is dominated by metal-rich stars on highly eccentric orbits. To shed light on the nature of this prominent halo component, we have analysed 28 Galaxy analogues in the Auriga suite of cosmological hydrodynamics zoom-in simulations. Some three quarters of the Auriga galaxies contain significant components with high radial velocity anisotropy, beta > 0.6. However, only in one third of the hosts do the high-beta stars contribute significantly to the accreted stellar halo overall, similar to what is observed in the Milky Way. For this particular subset we reveal the origin of the dominant stellar halo component with high metallicity, [Fe/H]~-1, and high orbital anisotropy, beta>0.8, by tracing their stars back to the epoch of accretion. It appears that, typically, these stars come from a single dwarf galaxy with a stellar mass of order of 10^9-10^10 Msol that merged around 6-10 Gyr ago, causing a sharp increase in the halo mass. Our study therefore establishes a firm link between the excess of radially anisotropic stellar debris in the Milky Way halo and an ancient head-on collision between the young Milky Way and a massive dwarf galaxy

Faraday rotation is one of the most widely used observables to infer the strength and configuration of the magnetic field in the ionised gas of the Milky Way and nearby spiral galaxies. Here we compute synthetic Faraday rotation maps at $z=0$ for a set of disk galaxies from the Auriga high-resolution cosmological simulations, for different observer positions within and outside the galaxy. We find that the strength of the Faraday rotation of our simulated galaxies for a hypothetic observer at the solar circle is broadly consistent with the Faraday rotation seen for the Milky Way. The same holds for an observer outside the galaxy and the observed signal of the nearby spiral galaxy M51. However, we also find that the structure and angular power spectra of the synthetic all-sky Faraday rotation maps vary strongly with azimuthal position along the solar circle. We argue that this variation is a result of the structure of the magnetic field of the galaxy that is dominated by an azimuthal magnetic field ordered scales of several kpc, but has radial and vertical magnetic field components that are only ordered on scales of 1-2 kpc. Because the magnetic field strength decreases exponentially with height above the disk, the Faraday rotation for an observer at the solar circle is dominated by the local environment. This represents a severe obstacle for attempts to reconstruct the global magnetic field of the Milky Way from Faraday rotation maps alone without including additional observables.

We combine Gaia data release 1 astrometry with Sloan Digital Sky Survey (SDSS) images taken some ~10-15 years earlier, to measure proper motions of stars in the halo of our Galaxy. The SDSS-Gaia proper motions have typical statistical errors of 2 mas/yr down to r ~ 20 mag, and are robust to variations with magnitude and colour. Armed with this exquisite set of halo proper motions, we identify RR Lyrae, blue horizontal branch (BHB), and K giant stars in the halo, and measure their net rotation with respect to the Galactic disc. We find evidence for a gently rotating prograde signal (< $V_\phi$ > ~ 5-25 km/s) in the halo stars, which shows little variation with Galactocentric radius out to 50 kpc. The average rotation signal for the three populations is < $V_\phi$ > = 14 +/- 2 +/- 10 (syst.) km/s. There is also tentative evidence for a kinematic correlation with metallicity, whereby the metal richer BHB and K giant stars have slightly stronger prograde rotation than the metal poorer stars. Using the Auriga simulation suite we find that the old (T >10 Gyr) stars in the simulated halos exhibit mild prograde rotation, with little dependence on radius or metallicity, in general agreement with the observations. The weak halo rotation suggests that the Milky Way has a minor in situ halo component, and has undergone a relatively quiet accretion history.

The magnetic fields observed in the Milky~Way and nearby galaxies appear to be in equipartition with the turbulent, thermal, and cosmic ray energy densities, and hence are expected to be dynamically important. However, the origin of these strong magnetic fields is still unclear, and most previous attempts to simulate galaxy formation from cosmological initial conditions have ignored them altogether. Here, we analyse the magnetic fields predicted by the simulations of the Auriga Project, a set of 30 high-resolution cosmological zoom simulations of Milky~Way-like galaxies, carried out with a moving-mesh magneto-hydrodynamics code and a detailed galaxy formation physics model. We find that the magnetic fields grow exponentially at early times owing to a small-scale dynamo with an e-folding time of roughly $100\,\rm{Myr}$ in the center of halos until saturation occurs around $z=2-3$, when the magnetic energy density reaches about $10\%$ of the turbulent energy density with a typical strength of $10-50\,\rm{\mu G}$. In the galactic centers the ratio between magnetic and turbulent energy remains nearly constant until $z=0$. At larger radii, differential rotation in the disks leads to linear amplification that typically saturates around $z=0.5$ to $z=0$. The final radial and vertical variations of the magnetic field strength can be well described by two joint exponential profiles, and are in good agreement with observational constraints. Overall, the magnetic fields have only little effect on the global evolution of the galaxies as it takes too long to reach equipartition. We also demonstrate that our results are well converged with numerical resolution.

A recent observational study of haloes of nearby Milky Way-like galaxies shows that only half (four out of eight) of the current sample exhibits strong negative metallicity ([Fe/H]) gradients. This is at odds with predictions from hydrodynamical simulations where such gradients are ubiquitous. In this Letter, we use high resolution cosmological hydrodynamical simulations to study the [Fe/H] distribution of galactic haloes. We find that kinematically selected stellar haloes, including both in-situ and accreted particles, have an oblate [Fe/H] distribution. Spherical [Fe/H] radial profiles show strong negative gradients within 100 kpc, in agreement with previous numerical results. However, the projected median [Fe/H] profiles along the galactic disc minor axis, typically obtained in observations, are significantly flatter. The median [Fe/H] values at a given radius are larger for the spherical profiles than for the minor axis profiles by as much as 0.4 dex within the inner 50 kpc. Similar results are obtained if only the accreted stellar component is considered indicating that the differences between spherical and minor axis profiles are not purely driven by heated disc star particles formed in situ. Our study highlights the importance of performing careful comparisons between models and observations of halo [Fe/H] distributions.

The Large Volume Detector (LVD) has been continuously taking data since 1992 at the INFN Gran Sasso National Laboratory. LVD is sensitive to neutrino bursts from gravitational stellar collapses with full detection probability over the Galaxy. We have searched for neutrino bursts in LVD data taken in 7335 days of operation. No evidence of neutrino signals has been found between June 1992 and December 2013. The 90% C.L. upper limit on the rate of core-collapse and failed supernova explosions out to distances of 25 kpc is found to be 0.114/y.

To asses stability against 1/f noise, the Low Frequency Instrument (LFI) onboard the Planck mission will acquire data at a rate much higher than the data rate allowed by its telemetry bandwith of 35.5 kbps. The data are processed by an onboard pipeline, followed onground by a reversing step. This paper illustrates the LFI scientific onboard processing to fit the allowed datarate. This is a lossy process tuned by using a set of 5 parameters Naver, r1, r2, q, O for each of the 44 LFI detectors. The paper quantifies the level of distortion introduced by the onboard processing, EpsilonQ, as a function of these parameters. It describes the method of optimizing the onboard processing chain. The tuning procedure is based on a optimization algorithm applied to unprocessed and uncompressed raw data provided either by simulations, prelaunch tests or data taken from LFI operating in diagnostic mode. All the needed optimization steps are performed by an automated tool, OCA2, which ends with optimized parameters and produces a set of statistical indicators, among them the compression rate Cr and EpsilonQ. For Planck/LFI the requirements are Cr = 2.4 and EpsilonQ <= 10% of the rms of the instrumental white noise. To speedup the process an analytical model is developed that is able to extract most of the relevant information on EpsilonQ and Cr as a function of the signal statistics and the processing parameters. This model will be of interest for the instrument data analysis. The method was applied during ground tests when the instrument was operating in conditions representative of flight. Optimized parameters were obtained and the performance has been verified, the required data rate of 35.5 Kbps has been achieved while keeping EpsilonQ at a level of 3.8% of white noise rms well within the requirements.

We give the calibration and scientific performance parameters of the Planck Low Frequency Instrument (LFI) measured during the ground cryogenic test campaign. These parameters characterise the instrument response and constitute our best pre-launch knowledge of the LFI scientific performance. The LFI shows excellent $1/f$ stability and rejection of instrumental systematic effects; measured noise performance shows that LFI is the most sensitive instrument of its kind. The set of measured calibration parameters will be updated during flight operations through the end of the mission.

In this paper we present the Low Frequency Instrument (LFI), designed and developed as part of the Planck space mission, the ESA program dedicated to precision imaging of the cosmic microwave background (CMB). Planck-LFI will observe the full sky in intensity and polarisation in three frequency bands centred at 30, 44 and 70 GHz, while higher frequencies (100-850 GHz) will be covered by the HFI instrument. The LFI is an array of microwave radiometers based on state-of-the-art Indium Phosphide cryogenic HEMT amplifiers implemented in a differential system using blackbody loads as reference signals. The front-end is cooled to 20K for optimal sensitivity and the reference loads are cooled to 4K to minimise low frequency noise. We provide an overview of the LFI, discuss the leading scientific requirements and describe the design solutions adopted for the various hardware subsystems. The main drivers of the radiometric, optical and thermal design are discussed, including the stringent requirements on sensitivity, stability, and rejection of systematic effects. Further details on the key instrument units and the results of ground calibration are provided in a set of companion papers.

This paper provides an overview of the Low Frequency Instrument (LFI) programme within the ESA Planck mission. The LFI instrument has been developed to produce high precision maps of the microwave sky at frequencies in the range 27-77 GHz, below the peak of the cosmic microwave background (CMB) radiation spectrum. The scientific goals are described, ranging from fundamental cosmology to Galactic and extragalactic astrophysics. The instrument design and development are outlined, together with the model philosophy and testing strategy. The instrument is presented in the context of the Planck mission. The LFI approach to ground and inflight calibration is described. We also describe the LFI ground segment. We present the results of a number of tests demonstrating the capability of the LFI data processing centre (DPC) to properly reduce and analyse LFI flight data, from telemetry information to calibrated and cleaned time ordered data, sky maps at each frequency (in temperature and polarization), component emission maps (CMB and diffuse foregrounds), catalogs for various classes of sources (the Early Release Compact Source Catalogue and the Final Compact Source Catalogue). The organization of the LFI consortium is briefly presented as well as the role of the core team in data analysis and scientific exploitation. All tests carried out on the LFI flight model demonstrate the excellent performance of the instrument and its various subunits. The data analysis pipeline has been tested and its main steps verified. In the first three months after launch, the commissioning, calibration, performance, and verification phases will be completed, after which Planck will begin its operational life, in which LFI will have an integral part.

We study the evolution of satellite debris to establish the most suitable space to identify past merger events. We confirm that the space of orbital frequencies is very promising in this respect. In frequency space individual streams can be easily identified, and their separation provides a direct measurement of the time of accretion. We are able to show for a few idealised gravitational potentials that these features are preserved also in systems that have evolved strongly in time. Furthermore, this time evolution is imprinted in the distribution of streams in frequency space. We have also tested the power of the orbital frequencies in a fully self-consistent (live) N-body simulation of the merger between a disk galaxy and a massive satellite. Even in this case streams can be easily identified and the time of accretion of the satellite can be accurately estimated.

We study the orbital properties of stars in four (published) simulations of thick disks formed by: i) accretion from disrupted satellites, ii) heating of a pre-existing thin disk by a minor merger, iii) radial migration and iv) gas rich mergers. We find that the distribution of orbital eccentricities are predicted to be different for each model: a prominent peak at low eccentricity is expected for the heating, migration and gas-rich merging scenarios, while the eccentricity distribution is broader and shifted towards higher values for the accretion model. These differences can be traced back to whether the bulk of the stars in each case is formed 'in-situ' or is 'accreted', and are robust to the peculiarities of each model. A simple test based on the eccentricity distribution of nearby thick disk stars may thus help elucidate the dominant formation mechanism of the Galactic thick disk.

We study the phase-space behaviour of nearby trajectories in integrable potentials. We show that the separation of nearby orbits initially diverges very fast, mimicking a nearly exponential behaviour, while at late times it grows linearly. This initial exponential phase, known as Miller's instability, is commonly found in N-body simulations, and has been attributed to short-term (microscopic) N-body chaos. However we show here analytically that the initial divergence is simply due to the shape of an orbit in phase-space. This result confirms previous suspicions that this transient phenomenon is not related to an instability in the sense of non-integrable behaviour in the dynamics of N-body systems.

In this paper we show the capabilities of the Large Volume Detector (INFN Gran Sasso National Laboratory) to identify a neutrino burst associated to a supernova explosion, in the absence of an "external trigger", e.g., an optical observation. We describe how the detector trigger and event selection have been optimized for this purpose, and we detail the algorithm used for the on-line burst recognition. The on-line sensitivity of the detector is defined and discussed in terms of supernova distance and electron anti-neutrino intensity at the source.

We report on a search for low-energy neutrino (antineutrino) bursts in correlation with the 8 time coincident events observed by the gravitational waves detectors EXPLORER and NAUTILUS (GWD) during the year 2001. The search, conducted with the LVD detector (INFN Gran Sasso National Laboratory, Italy), has considered several neutrino reactions, corresponding to different neutrino species, and a wide range of time intervals around the (GWD) observed events. No evidence for statistically significant correlated signals in LVD has been found. Assuming two different origins for neutrino emission, the cooling of a neutron star from a core-collapse supernova or from coalescing neutron stars and the accretion of shocked matter, and taking into account neutrino oscillations, we derive limits to the total energy emitted in neutrinos and to the amount of accreting mass, respectively.