Using the NOvA neutrino detectors, a broad search has been performed for any signal coincident with 28 gravitational wave events detected by the LIGO/Virgo Collaboration between September 2015 and July 2019. For all of these events, NOvA is sensitive to possible arrival of neutrinos and cosmic rays of GeV and higher energies. For five (seven) events in the NOvA Far (Near) Detector, timely public alerts from the LIGO/Virgo Collaboration allowed recording of MeV-scale events. No signal candidates were found.

We present a grid-based non-parametric approach to obtain a triaxial three-dimensional luminosity density from its surface brightness distribution. Triaxial deprojection is highly degenerate and our approach illustrates the resulting difficulties. Fortunately, for massive elliptical galaxies, many deprojections for a particular line of sight can be discarded, because their projection along other lines of sight does not resemble elliptical galaxies. The near-elliptical isophotes of these objects imply near ellipsoidal intrinsic shapes. In fact, deprojection is unique for densities distributed on ellipsoidal shells. The constrained non-parametric deprojection method we present here relaxes this constraint and assumes that the contours of the luminosity density are boxy/discy ellipsoids with radially varying axis ratios. With this approach we are able to reconstruct the intrinsic triaxial densities of our test models, including one drawn from an $N$-body simulation. The method also allows to compare the relative likelihood of deprojections at different viewing angles. We show that the viewing orientations of individual galaxies with nearly ellipsoidal isophotal shapes can be constrained from photometric data alone.

We study the association between active galactic nuclei (AGN) and environment at scales of $0.01-1\ h^{-1}$Mpc in the IllustrisTNG (TNG100) simulated universe. We identify supermassive black hole (BH) pairs and multiples within scales of 0.01, 0.1, & 1 $h^{-1}$Mpc and examine their AGN activity in relation to randomly-selected pairs and multiples. The number density of BHs in TNG100 is $n=0.06\,h^3$Mpc$^{-3}$ at $z\lesssim1.5$ ($n=0.02\,h^3$ Mpc$^{-3}$ at $z=3$). About $\sim10$% and $\sim1$% of them live in pairs and multiples, respectively, within 0.1 $h^{-1}$Mpc scales. We find that BH systems have enhanced likelihood (up to factors of 3-6) of containing high Eddington ratio ($\eta\gtrsim0.7$) AGN compared to random pairs and multiples. Conversely, the likelihood of an AGN to live in 0.1$h^{-1}$Mpc scale systems is also higher (by factors $\sim4$ for $\eta\gtrsim0.7$) compared to random pairs and multiples. We also estimate that $\sim10$% of ultra-hard X-ray selected AGN in TNG100 have detectable 2-10 keV AGN companions on $0.1\ h^{-1}$Mpc scales, in agreement with observations. On larger spatial scales ($\sim 1$ $h^{-1}$Mpc), however, no significant enhancement is associated with BH pairs and multiples, even at high Eddington ratios. The enhancement of AGN activity in rich, small-scale ($\lesssim0.1$ $h^{-1}$Mpc) environments is therefore likely to be driven by galaxy interactions and mergers. Nonetheless, the overall percentage of AGN that live in $\lesssim0.1$ $h^{-1}$Mpc scale multiples is still subdominant (at most $\sim40$% for the highest Eddington ratio AGN). Furthermore, the enhancement in Eddington ratios of BH systems(as well as merging BHs) is only up to factors of $\sim2-3$. Thus, our results support the existence of a merger-AGN connection, but they also suggest that mergers and interactions play a relatively minor role in fueling the AGN population as a whole.

Massive black hole (MBH) binary inspiral time scales are uncertain, and their spins are even more poorly constrained. Spin misalignment, along with unequal mass ratios and spin magnitudes, introduces asymmetry in the gravitational radiation, which imparts a recoil kick to the merged MBH. Understanding how MBH binary spins evolve is crucial for determining their recoil velocities, their gravitational wave (GW) waveforms detectable with LISA, as well as their post-merger retention rate in galaxies and thus their subsequent merger rate. Here we present a novel study that introduces a sub-resolution model for gas- and GW-driven MBH binary spin evolution using a population of accreting MBHs from the Illustris cosmological hydrodynamics simulations. We also model sub-resolution binary inspiral via dynamical friction, stellar scattering, viscous gas drag, and GW emission. Our model assumes differential accretion, which causes greater alignment of the secondary MBH spin in unequal-mass mergers. We find that 47% of the MBHs in our population merge by $z=0$. Of these, 19% have misaligned primaries and 10% have misaligned secondaries at the time of merger in our (conservative) fiducial model. The MBH misalignment fraction depends strongly on the accretion disc parameters, however. Reducing accretion rates by a factor of 100, in a thicker disc, yields 79% and 42% misalignment for primaries and secondaries, respectively. Even in the fiducial model, more than 12% of binaries experience recoils of $>500$ km s$^{-1}$, which could displace them at least temporarily from galactic nuclei. We additionally find that a significant number of systems experience strong precession.

We use high spatial resolution stellar velocity maps from the Gemini GMOS integral-field spectrograph (IFS) and wide-field velocity maps from the McDonald Mitchell IFS to study the stellar velocity profiles and kinematic misalignments from $\sim 200$ pc to $\sim 20$ kpc in 20 early-type galaxies with stellar mass $M_* > 10^{11.7} M_\odot$ in the MASSIVE survey. While 80% of the galaxies have low spins ($\lambda < 0.1$) and low rotational velocities ($< 50$ km/s) in both the central region and the main body, we find a diverse range of velocity features and misalignment angles. For the 18 galaxies with measurable central kinematic axes, 10 have well aligned kinematic axis and photometric major axis, and the other 8 galaxies have misalignment angles that are distributed quite evenly from $15^\circ$ to the maximal value of $90^\circ$. There is a strong correlation between central kinematic misalignment and galaxy spin, where all 4 galaxies with significant spins have well aligned kinematic and photometric axes, but only 43% of the low-spin galaxies are well aligned. The central and main-body kinematic axes within a galaxy are not always aligned. When the two kinematic axes are aligned ($\sim 60$% of the cases), they are either also aligned with the photometric major axis or orthogonal to it. We find 13 galaxies to also exhibit noticeable local kinematic twists, and one galaxy to have a counter-rotating core. A diverse assembly history consisting of multiple gas-poor mergers of a variety of progenitor mass ratios and orbits is likely to be needed to account for the predominance of low spins and the wide range of central and main-body velocity features reported here for local massive ETGs.

Holm 15A, the brightest cluster galaxy (BCG) of the galaxy cluster Abell 85, has an ultra-diffuse central region, 2 mag fainter than the faintest depleted core of any early-type galaxy (ETG) that has been dynamically modelled in detail. We use orbit-based, axisymmetric Schwarzschild models to analyse the stellar kinematics of Holm 15A from new high-resolution, wide-field spectral observations obtained with MUSE at the VLT. We find a supermassive black hole (SMBH) with a mass of (4.0 +- 0.80) x 10^10 solar masses at the center of Holm 15A. This is the most massive black hole with a direct dynamical detection in the local universe. We find that the distribution of stellar orbits is increasingly biased towards tangential motions inside the core. However, the tangential bias is less than in other cored elliptical galaxies. We compare Holm 15A with N-body simulations of mergers between galaxies with black holes and find that the observed amount of tangential anisotropy and the shape of the light profile are consistent with a formation scenario where Holm 15A is the remnant of a merger between two ETGs with pre-existing depleted cores. We find that black hole masses in cored galaxies, including Holm 15A, scale inversely with the central stellar surface brightness and mass density, respectively. These correlation are independent of a specific parameterization of the light profile.

Using two years of data from the NOvA Near Detector at Fermilab, we report a seasonal variation of cosmic ray induced multiple-muon event rates which has an opposite phase to the seasonal variation in the atmospheric temperature. The strength of the seasonal multipl$ increase as a function of the muon multiplicity. However, no significant dependence of the strength of the seasonal variation of the multiple-muon variation is seen as a function of the muon zenith angle, or the spatial or angular separation between the correlated muons.

We report on a detailed optical spectroscopic follow-up of the black hole transient MAXI J1820+070 (ASASSN-18ey). The observations cover the main part of the X-ray binary outburst, when the source alternated between hard and soft states following the classical pattern widely seen in other systems. We focus the analysis on the He I emission lines at 5876 and 6678 Angs, as well as on Halpha. We detect clear accretion disk wind features (P-Cyg profiles and broad emission line wings) in the hard state, both during outburst rise and decay. These are not witnessed during the several months long soft state. However, our data suggest that the visibility of the outflow might be significantly affected by the ionisation state of the accretion disk. The terminal velocity of the wind is above ~ 1200 km/s, which is similar to outflow velocities derived from (hard-state) optical winds and (soft-state) X-ray winds in other systems. The wind signatures, in particular the P-Cyg profiles, are very shallow, and their detection has only been possible thanks to a combination of source brightness and intense monitoring at very high signal-to-noise. This study indicates that cold, optical winds are most likely a common feature of black hole accretion, and therefore, that wind-like outflows are a general mechanism of mass and angular momentum removal operating throughout the entire X-ray binary outburst.

Urea, NH2C(O)NH2, is a molecule of great importance in organic chemistry and biology. Two searches for urea in the interstellar medium were reported in the past, but neither were conclusive. We want to take advantage of the increased sensitivity and angular resolution provided by ALMA to search for urea toward the hot cores embedded in the high-mass star forming region Sgr B2(N). We use the new spectral line survey ReMoCA performed toward Sgr B2(N) with ALMA in its observing cycle 4. The spectra are analyzed under the LTE approximation. We construct a full synthetic spectrum that includes all the molecules identified so far. We use new spectroscopic predictions for urea in its vibrational ground state and first vibrationally excited state to search for this complex organic molecule in the ReMoCA data set. We employ the gas-grain chemical kinetics model MAGICKAL to interpret the astronomical observations. We report the secure detection of urea toward the hot core Sgr B2(N1) at a position called N1S slightly offset from the continuum peak, which avoids obscuration by the dust.. We derive a column density of 2.7x10^16 cm-2 for urea, two orders of magnitude lower than formamide, and one order of magnitude below methyl isocyanate, acetamide, and N-methylformamide. The latter molecule is reliably identified toward N1S with 60 clearly detected lines, confirming an earlier claim of its tentative interstellar detection. We report the first interstellar detections of NH2CH18O and 15NH2CHO. We also report the nondetection of urea toward the secondary hot core Sgr B2(N2) with an abundance relative to the other four species at least one order of magnitude lower than toward the main hot core. Our chemical model roughly reproduces the relative abundances of formamide, methyl isocyanate, acetamide, and N-methylformamide, but it overproduces urea by at least one order of magnitude.

The astrophysical $S$-factor for the radiative proton capture reaction on $^7$Be ($S_{17}$) at low energies is affected by the $s$-wave scattering lengths. We report the measurement of elastic and inelastic scattering cross sections for the $^7$Be+p system in the center-of-mass energy range 0.474 - 2.740 MeV and center-of-mass angular range of 70$^\circ$- 150$^\circ$. A radioactive $^7$Be beam produced at Oak Ridge National Laboratory's (ORNL) Holifield Radioactive Ion Beam Facility was accelerated and bombarded a thin polypropylene (CH$_{2}$)$_\text n$ target. Scattered ions were detected in the segmented Silicon Detector Array. Using an $\textit{R}$-matrix analysis of ORNL and Louvain-la-Neuve cross section data, the $s$-wave scattering lengths for channel spins 1 and 2 were determined to be 17.34$^{+1.11}_{-1.33}$ and -3.18$^{+0.55}_{-0.50}$ fm, respectively. The uncertainty in the $s$-wave scattering lengths reported in this work is smaller by a factor of 5-8 compared to the previous measurement, which may reduce the overall uncertainty in $S_{17}$ at zero energy. The level structure of $^8$B is discussed based upon the results from this work. Evidence for the existence of 0$^+$ and 2$^+$ levels in $^8$B at 1.9 and 2.21 MeV, respectively, is observed.

We use observations from the GEMINI-N/GMOS integral-field spectrograph (IFS) to obtain spatially resolved stellar kinematics of the central $\sim 1$ kpc of 20 early-type galaxies (ETGs) with stellar masses greater than $10^{11.7} M_\odot$ in the MASSIVE survey. Together with observations from the wide-field Mitchell IFS at McDonald Observatory in our earlier work, we obtain unprecedentedly detailed kinematic maps of local massive ETGs, covering a scale of $\sim 0.1-30$ kpc. The high ($\sim 120$) signal-to-noise of the GMOS spectra enable us to obtain two-dimensional maps of the line-of-sight velocity, velocity dispersion $\sigma$, as well as the skewness $h_3$ and kurtosis $h_4$ of the stellar velocity distributions. All but one galaxy in the sample have $\sigma(R)$ profiles that increase towards the center, whereas the slope of $\sigma(R)$ at one effective radius ($R_e$) can be of either sign. The $h_4$ is generally positive, with 14 of the 20 galaxies having positive $h_4$ within the GMOS aperture and 18 having positive $h_4$ within $1 R_e$. The positive $h_4$ and rising $\sigma(R)$ towards small radii are indicative of a central black hole and velocity anisotropy. We demonstrate the constraining power of the data on the mass distributions in ETGs by applying Jeans anisotropic modeling (JAM) to NGC~1453, the most regular fast rotator in the sample. Despite the limitations of JAM, we obtain a clear $\chi^2$ minimum in black hole mass, stellar mass-to-light ratio, velocity anisotropy parameters, and the circular velocity of the dark matter halo.

We study the impact of merging supermassive black holes (SMBHs) on the central regions of massive early-type galaxies (ETGs) using a series of merger simulations with varying mass ratios. The ETG models include realistic stellar and dark matter components and are evolved with the GADGET based regularized tree code KETJU. We show that observed key properties of the nuclear stellar populations of massive ETGs, namely flat stellar density distributions (cores), tangentially biased velocity distributions and kinematically decoupled (counter-)rotation can naturally result from a single process $-$ the scouring by SMBHs. Major mergers with mass ratios of $q>1/3$ produce flat, tangentially biased cores with kinematically distinct components. These features originate from spin reversals during the decay of the SMBH orbits caused by gravitational torques from tidally expelled nuclear material. Minor mergers ($q\lesssim1/3$) on the other hand, form non-rotating cores and the tidal reversal becomes less important. Low-density stellar cores scoured in (multiple) minor mergers are less tangentially biased. This implies that the nuclear stellar properties of massive ETGs can be solely explained by stellar dynamical processes during their final assembly without any need for `feedback' from accreting black holes. We predict a strong correlation between decoupled cores, central anisotropy and merger history: decoupled cores form in binary mergers and we predict them to occur in elliptical galaxies with the strongest central anisotropy. Measurements of the central orbital structure are the key to understanding the number of mergers a given galaxy has experienced.

We measure the stellar populations as a function of radius for 90 early-type galaxies (ETGs) in the MASSIVE survey, a volume-limited integral-field spectroscopic (IFS) galaxy survey targeting all northern-sky ETGs with absolute K-band magnitude M_K < -25.3 mag, or stellar mass M* 4x10^11 M_sun, within 108 Mpc. We are able to measure reliable stellar population parameters for individual galaxies out to 10-20 kpc (1-3 R_e) depending on the galaxy. Focusing on ~R_e (~10 kpc), we find significant correlations between the abundance ratios, sigma, and M* at large radius, but we also find that the abundance ratios saturate in the highest-mass bin. We see a strong correlation between the kurtosis of the line of sight velocity distribution (h4) and the stellar population parameters beyond R_e. Galaxies with higher radial anisotropy appear to be older, with metal-poorer stars and enhanced [alpha/Fe]. We suggest that the higher radial anisotropy may derive from more accretion of small satellites. Finally, we see some evidence for correlations between environmental metrics (measured locally and on >5 Mpc scales) and the stellar populations, as expected if satellites are quenched earlier in denser environments.

Given its velocity dispersion, the early-type galaxy NGC 1600 has an unusually massive ($M_\bullet = 1.7 \times 10^{10} M_\odot$) central supermassive black hole (SMBH), surrounded by a large core ($r_\mathrm{b} = 0.7$ kpc) with a tangentially biased stellar distribution. We present high-resolution equal-mass merger simulations including SMBHs to study the formation of such systems. The structural parameters of the progenitor ellipticals were chosen to produce merger remnants resembling NGC 1600. We test initial stellar density slopes of $\rho \propto r^{-1}$ and $\rho \propto r^{-3/2}$ and vary the initial SMBH masses from $8.5 \times 10^8$ to $8.5 \times 10^9$ $M_\odot$. With increasing SMBH mass the merger remnants show a systematic decrease in central surface brightness, an increasing core size, and an increasingly tangentially biased central velocity anisotropy. Two-dimensional kinematic maps reveal decoupled, rotating core regions for the most massive SMBHs. The stellar cores form rapidly as the SMBHs become bound, while the velocity anisotropy develops more slowly after the SMBH binaries become hard. The simulated merger remnants follow distinct relations between the core radius and the sphere-of-influence, and the SMBH mass, similar to observed systems. We find a systematic change in the relations as a function of the progenitor density slope, and present a simple scouring model reproducing this behavior. Finally, we find the best agreement with NGC 1600 using SMBH masses totaling the observed value of $M_\bullet = 1.7 \times 10^{10} M_\odot$. In general, density slopes of $\rho \propto r^{-3/2}$ for the progenitor galaxies are strongly favored for the equal-mass merger scenario.

Fossil groups (FGs) are galaxy aggregates with an extended and luminous X-ray halo, which are dominated by a very massive early-type galaxy and lack of L* objects. FGs are indeed characterized by a large magnitude gap between their central and surrounding galaxies. This is explained by either speculating that FGs are failed groups that formed without bright satellite galaxies and did not suffer any major merger, or by suggesting that FGs are very old systems that had enough time to exhaust their bright satellite galaxies through multiple major mergers. Since major mergers leave signatures in the stellar populations of the resulting galaxy, we study the stellar population parameters of the brightest central galaxies (BCGs) of FGs as a benchmark against which the formation and evolution scenarios of FGs can be compared. We present long-slit spectroscopic observations along different axes of NGC 6482 and NGC 7556, which are the BCGs of two nearby FGs. The measurements include spatially resolved stellar kinematics and radial profiles of line-strength indices, which we converted into stellar population parameters using single stellar-population models. NGC 6482 and NGC 7556 are very massive and large galaxies and host a centrally concentrated stellar population, which is significantly younger and more metal rich than the rest of the galaxy. The age gradients of both galaxies are somewhat larger than those of the other FG BCGs studied so far, whereas their metallicity gradients are similarly negative and shallow. They have negligible gradients of alpha-element abundance ratio. The measured metallicity gradients are less steep than those predicted for massive galaxies that formed monolithically and evolved without experiencing any major merger. We conclude that the observed FGs formed through major mergers rather than being failed groups that lacked bright satellite galaxies from the beginning.

MICADO will enable the ELT to perform diffraction limited near-infrared observations at first light. The instrument's capabilities focus on imaging (including astrometric and high contrast) as well as single object spectroscopy. This contribution looks at how requirements from the observing modes have driven the instrument design and functionality. Using examples from specific science cases, and making use of the data simulation tool, an outline is presented of what we can expect the instrument to achieve.

We use spatially resolved two-dimensional stellar velocity maps over a $107"\times 107"$ field of view to investigate the kinematic features of 90 early-type galaxies above stellar mass $10^{11.5}M_\odot$ in the MASSIVE survey. We measure the misalignment angle $\Psi$ between the kinematic and photometric axes and identify local features such as velocity twists and kinematically distinct components. We find 46% of the sample to be well aligned ($\Psi < 15^{\circ}$), 33% misaligned, and 21% without detectable rotation (non-rotators). Only 24% of the sample are fast rotators, the majority of which (91%) are aligned, whereas 57% of the slow rotators are misaligned with a nearly flat distribution of $\Psi$ from $15^{\circ}$ to $90^{\circ}$. 11 galaxies have $\Psi \gtrsim 60^{\circ}$ and thus exhibit minor-axis ("prolate") rotation in which the rotation is preferentially around the photometric major axis. Kinematic misalignments occur more frequently for lower galaxy spin or denser galaxy environments. Using the observed misalignment and ellipticity distributions, we infer the intrinsic shape distribution of our sample and find that MASSIVE slow rotators are consistent with being mildly triaxial, with mean axis ratios of $b/a=0.88$ and $c/a=0.65$. In terms of local kinematic features, 51% of the sample exhibit kinematic twists of larger than $20^{\circ}$, and 2 galaxies have kinematically distinct components. The frequency of misalignment and the broad distribution of $\Psi$ reported here suggest that the most massive early-type galaxies are mildly triaxial, and that formation processes resulting in kinematically misaligned slow rotators such as gas-poor mergers occur frequently in this mass range.

We present near-infrared observations of 35 of the most massive early-type galaxies in the local universe. The observations were made using the infrared channel of the Hubble Space Telescope Wide Field Camera 3 in the F110W (1.1 $\mu$m) filter. We measured surface brightness profiles and elliptical isophotal fit parameters from the nuclear regions out to a radius of ~10 kpc in most cases. We find that 37% (13) of the galaxies in our sample have isophotal position angle rotations greater than 20 degrees over the radial range imaged by WFC3/IR, which is often due to the presence of neighbors or multiple nuclei. Most galaxies in our sample are significantly rounder near the center than in the outer regions. This sample contains six fast rotators and 28 slow rotators. We find that all fast rotators are either disky or show no measurable deviation from purely elliptical isophotes. Among slow rotators, significantly disky and boxy galaxies occur with nearly equal frequency. The galaxies in our sample often exhibit changing isophotal shapes, sometimes showing both significantly disky and boxy isophotes at different radii. The fact that parameters vary widely between galaxies and within individual galaxies is evidence that these massive galaxies have complicated formation histories, and some of them have experienced recent mergers and have not fully relaxed. These data demonstrate the value of high spatial resolution IR imaging of galaxies and provide measurements necessary for determining stellar masses, dynamics, and black hole masses in high mass galaxies.

We measure the radial profiles of the stellar velocity dispersions, $\sigma(R)$, for 90 early-type galaxies (ETGs) in the MASSIVE survey, a volume-limited integral-field spectroscopic (IFS) galaxy survey targeting all northern-sky ETGs with absolute $K$-band magnitude $M_K < -25.3$ mag, or stellar mass $M_* > 4 \times 10^{11} M_\odot$, within 108 Mpc. Our wide-field 107" $\times$ 107" IFS data cover radii as large as 40 kpc, for which we quantify separately the inner (2 kpc) and outer (20 kpc) logarithmic slopes $\gamma_{\rm inner}$ and $\gamma_{\rm outer}$ of $\sigma(R)$. While $\gamma_{\rm inner}$ is mostly negative, of the 56 galaxies with sufficient radial coverage to determine $\gamma_{\rm outer}$ we find 36% to have rising outer dispersion profiles, 30% to be flat within the uncertainties, and 34% to be falling. The fraction of galaxies with rising outer profiles increases with $M_*$ and in denser galaxy environment, with 10 of the 11 most massive galaxies in our sample having flat or rising dispersion profiles. The strongest environmental correlations are with local density and halo mass, but a weaker correlation with large-scale density also exists. The average $\gamma_{\rm outer}$ is similar for brightest group galaxies, satellites, and isolated galaxies in our sample. We find a clear positive correlation between the gradients of the outer dispersion profile and the gradients of the velocity kurtosis $h_4$. Altogether, our kinematic results suggest that the increasing fraction of rising dispersion profiles in the most massive ETGs are caused (at least in part) by variations in the total mass profiles rather than in the velocity anisotropy alone.

Detection of ultra-high energy (UHE, $\gtrsim 1$ PeV) neutrinos from astrophysical sources will be a major advancement in identifying and understanding the sources of UHE cosmic rays (CRs) in nature. Long duration gamma-ray burst (GRB) blast waves have been considered as potential acceleration sites of UHECRs. These CRs are expected to interact with GRB afterglow photons, which is synchrotron radiation from relativistic electrons co-accelerated with CRs in the blast wave, and naturally produce UHE neutrinos. Fluxes of these neutrinos are uncertain, however, and crucially depend on the observed afterglow modeling. We have selected a sample of 23 long duration GRBs within redshift 0.5 for which adequate electromagnetic afterglow data are available and which could produce high flux of UHE afterglow neutrinos, being nearby. We fit optical, X-ray and $\gamma$-ray afterglow data with an adiabatic blast wave model in a constant density interstellar medium and in a wind environment where the density of the wind decreases as inverse square of the radius from the center of the GRB. The blast wave model parameters extracted from these fits are then used for calculating UHECR acceleration and $p\gamma$ interactions to produce UHE neutrino fluxes from these GRBs. We have also explored the detectability of these neutrinos by currently running and upcoming large area neutrino detectors, such as the Pierre Auger Observatory, IceCube Gen-2 and KM3NeT observatories. We find that our realistic flux models from nearby GRBs will be unconstrained in foreseeable future.

We present stellar-dynamical measurements of the central supermassive black hole (SMBH) in the S0 galaxy NGC 307, using adaptive-optics IFU data from VLT-SINFONI. We investigate the effects of including dark-matter haloes as well as multiple stellar components with different mass-to-light (M/L) ratios in the dynamical modeling. Models with no halo and a single stellar component yield a relatively poor fit with a low value for the SMBH mass ($7.0 \pm 1.0 \times 10^{7} M_{\odot}$) and a high stellar M/L ratio (K-band M/L = $1.3 \pm 0.1$). Adding a halo produces a much better fit, with a significantly larger SMBH mass ($2.0 \pm 0.5 \times 10^{8} M_{\odot}$) and a lower M/L ratio ($1.1 \pm 0.1$). A model with no halo but with separate bulge and disc components produces a similarly good fit, with a slightly larger SMBH mass ($3.0 \pm 0.5 \times 10^{8} M_{\odot}$) and an identical M/L ratio for the bulge component, though the disc M/L ratio is biased high (disc M/L $ = 1.9 \pm 0.1$). Adding a halo to the two-stellar-component model results in a much more plausible disc M/L ratio of $1.0 \pm 0.1$, but has only a modest effect on the SMBH mass ($2.2 \pm 0.6 \times 10^{8} M_{\odot}$) and leaves the bulge M/L ratio unchanged. This suggests that measuring SMBH masses in disc galaxies using just a single stellar component and no halo has the same drawbacks as it does for elliptical galaxies, but also that reasonably accurate SMBH masses and bulge M/L ratios can be recovered (without the added computational expense of modeling haloes) by using separate bulge and disc components.

We analyse the environmental properties of 370 local early-type galaxies (ETGs) in the MASSIVE and ATLAS3D surveys, two complementary volume-limited integral-field spectroscopic (IFS) galaxy surveys spanning absolute $K$-band magnitude $-21.5 > M_K > -26.6$, or stellar mass $8 \times 10^{9} < M_* < 2 \times 10^{12} M_\odot$. We find these galaxies to reside in a diverse range of environments measured by four methods: group membership (whether a galaxy is a brightest group/cluster galaxy, satellite, or isolated), halo mass, large-scale mass density (measured over a few Mpc), and local mass density (measured within the $N$th neighbour). The spatially resolved IFS stellar kinematics provide robust measurements of the spin parameter $\lambda_e$ and enable us to examine the relationship among $\lambda_e$, $M_*$, and galaxy environment. We find a strong correlation between $\lambda_e$ and $M_*$, where the average $\lambda_e$ decreases from $\sim 0.4$ to below 0.1 with increasing mass, and the fraction of slow rotators $f_{\rm slow}$ increases from $\sim 10$% to 90%. We show for the first time that at fixed $M_*$, there are almost no trends between galaxy spin and environment; the apparent kinematic morphology-density relation for ETGs is therefore primarily driven by $M_*$ and is accounted for by the joint correlations between $M_*$ and spin, and between $M_*$ and environment. A possible exception is that the increased $f_{\rm slow}$ at high local density is slightly more than expected based only on these joint correlations. Our results suggest that the physical processes responsible for building up the present-day stellar masses of massive galaxies are also very efficient at reducing their spin, in any environment.

We present the first systematic investigation of the existence, spatial distribution, and kinematics of warm ionized gas as traced by the [O II] 3727AA emission line in 74 of the most massive galaxies in the local Universe. All of our galaxies have deep integral field spectroscopy from the volume- and magnitude-limited MASSIVE survey of early-type galaxies with stellar mass log(M_*/M_sun) > 11.5 (M_K < -25.3 mag) and distance D < 108 Mpc. Of the 74 galaxies in our sample, we detect warm ionized gas in 28, which yields a global detection fraction of 38\pm6% down to a typical [O II] equivalent width limit of 2AA. MASSIVE fast rotators are more likely to have gas than MASSIVE slow rotators with detection fractions of 80\pm10% and 28\pm6%, respectively. The spatial extents span a wide range of radii (0.6 - 18.2 kpc; 0.1 - 4R_e), and the gas morphologies are diverse, with 17/28 = 61\pm9% being centrally concentrated, 8/28 = 29\pm9% exhibiting clear rotation out to several kpc, and 3/28 = 11\pm6% being extended but patchy. Three out of four fast rotators show kinematic alignment between the stars and gas, whereas the two slow rotators with robust kinematic measurements available exhibit kinematic misalignment. Our inferred warm ionized gas masses are roughly ~10^5M_sun. The emission line ratios and radial equivalent width profiles are generally consistent with excitation of the gas by the old underlying stellar population. We explore different gas origin scenarios for MASSIVE galaxies and find that a variety of physical processes are likely at play, including internal gas recycling, cooling out of the hot gaseous halo, and gas acquired via mergers.

New photometric and long-slit spectroscopic observations are presented for NGC 7113, PGC 1852, and PGC 67207 which are three bright galaxies residing in low-density environments. The surface-brightness distribution is analysed from the K_S-band images taken with adaptive optics at the Gemini North Telescope and the ugriz-band images from the Sloan Digital Sky Survey while the line-of-sight stellar velocity distribution and line-strength Lick indices inside the effective radius are measured along several position angles. The age, metallicity, and alpha-element abundance of the galaxies are estimated from single stellar-population models. In spite of the available morphological classification, images show that PGC 1852 is a barred spiral which we do not further consider for mass modelling. The structural parameters of the two early-type galaxies NGC 7113 and PGC 67207 are obtained from a two-dimensional photometric decomposition and the mass-to-light ratio of all the (luminous and dark) mass that follows the light is derived from orbit-based axisymmetric dynamical modelling together with the mass density of the dark matter halo. The dynamically derived mass that follows the light is about a factor of 2 larger than the stellar mass derived using stellar-population models with Kroupa initial mass function. Both galaxies have a lower content of halo dark matter with respect to early-type galaxies in high-density environments and in agreement with the predictions of semi-analytical models of galaxy formation.

We present spatially-resolved two-dimensional stellar kinematics for the 41 most massive early-type galaxies (MK <~ -25.7 mag, stellar mass M* >~ 10^11.8 Msun) of the volume-limited (D < 108 Mpc) MASSIVE survey. For each galaxy, we obtain high-quality spectra in the wavelength range of 3650 to 5850 A from the 246-fiber Mitchell integral-field spectrograph (IFS) at McDonald Observatory, covering a 107" x 107" field of view (often reaching 2 to 3 effective radii). We measure the 2-D spatial distribution of each galaxy's angular momentum (lambda and fast or slow rotator status), velocity dispersion (sigma), and higher-order non-Gaussian velocity features (Gauss-Hermite moments h3 to h6). Our sample contains a high fraction (~80% ) of slow and non-rotators with lambda <~ 0.2. When combined with the lower-mass ETGs in the ATLAS3D survey, we find the fraction of slow-rotators to increase dramatically with galaxy mass, reaching ~50% at MK ~ -25.5 mag and ~90% at MK <~ -26 mag. All of our fast rotators show a clear anti-correlation between h3 and V/sigma, and the slope of the anti-correlation is steeper in more round galaxies. The radial profiles of sigma show a clear luminosity and environmental dependence: the 12 most luminous galaxies in our sample (MK <~ -26 mag) are all brightest cluster/group galaxies (except NGC 4874) and all have rising or nearly flat sigma profiles, whereas five of the seven "isolated" galaxies are all fainter than MK = -25.8 mag and have falling sigma. All of our galaxies have positive average h4; the most luminous galaxies have average h4 ~ 0.05 while less luminous galaxies have a range of values between 0 and 0.05. Most of our galaxies show positive radial gradients in h4, and those galaxies also tend to have rising sigma profiles. We discuss the implications for the relationship among dynamical mass, sigma, h4, and velocity anisotropy for these massive galaxies.

Studies of the physical properties of local elliptical galaxies (e.g., gas temperatures, halo masses, stellar kinematics) are shedding new light on galaxy formation. Here we present the hot X-ray gas properties of 33 early-type systems within the MASSIVE galaxy survey sample that have archival Chandra X-ray observations. Through careful X-ray spectral modeling, we derive X-ray luminosities (L_X) and plasma temperatures (T_gas) for the diffuse gas components in these galaxies. We combine the MASSIVE sample with 41 galaxies from the ATLAS^3D survey to investigate the X-ray and optical properties of a statistically significant sample of nearby early-type galaxies across a wide-range of environments. We deduce that all early-type galaxies (independent of galaxy mass and rotational support) follow a universal scaling law such that L_X~T_gas^4.5. When X-ray measurements are performed consistently in apertures set by the galaxy stellar content, the wide-scale environment does not contribute to the intrinsic scatter (~0.5dex) within the scaling relation. We further demonstrate that the scatter in L_X around both K-band luminosity (L_K) and the galaxy stellar velocity dispersion is primarily driven by T_gas, with no clear trends with halo mass, radio power, or angular momentum of the stars. It is not trivial to tie the origin of the gas directly to either the stellar mass or the galaxy potential. Indeed, our data require a steeper relation between L_X, L_K, and sigma_e than predicted by standard mass-loss models. Finally, we find a statistically significant correlation between sigma_e and T_gas, suggesting that T_gas is set by the galaxy potential inside the optical effective radius. We conclude that within the inner-most 10-30kpc region, early-types maintain pressure-supported hot gas, with a minimum T_gas set by the virial temperature, but the majority show evidence for some additional heating.

We obtained adaptive-optics assisted SINFONI observations of the central regions of the giant elliptical galaxy NGC5419 with a spatial resolution of 0.2 arcsec ($\approx 55$ pc). NGC5419 has a large depleted stellar core with a radius of 1.58 arcsec (430 pc). HST and SINFONI images show a point source located at the galaxy's photocentre, which is likely associated with the low-luminosity AGN previously detected in NGC5419. Both the HST and SINFONI images also show a second nucleus, off-centred by 0.25 arcsec ($\approx 70$ pc). Outside of the central double nucleus, we measure an almost constant velocity dispersion of $\sigma \sim 350$ km/s. In the region where the double nucleus is located, the dispersion rises steeply to a peak value of $\sim 420$ km/s. In addition to the SINFONI data, we also obtained stellar kinematics at larger radii from the South African Large Telescope. While NGC5419 shows low rotation ($v < 50$ km/s), the central regions (inside $\sim 4 \, r_b$) clearly rotate in the opposite direction to the galaxy's outer parts. We use orbit-based dynamical models to measure the black hole mass of NGC5419 from the kinematical data outside of the double nuclear structure. The models imply M$_{\rm BH}=7.2^{+2.7}_{-1.9} \times 10^9$ M$_{\odot}$. The enhanced velocity dispersion in the region of the double nucleus suggests that NGC5419 possibly hosts two supermassive black holes at its centre, separated by only $\approx 70$ pc. Yet our measured M$_{\rm BH}$ is consistent with the black hole mass expected from the size of the galaxy's depleted stellar core. This suggests, that systematic uncertainties in M$_{\rm BH}$ related to the secondary nucleus are small.

MICADO will equip the E-ELT with a first light capability for diffraction limited imaging at near-infrared wavelengths. The instrument's observing modes focus on various flavours of imaging, including astrometric, high contrast, and time resolved. There is also a single object spectroscopic mode optimised for wavelength coverage at moderately high resolution. This contribution provides an overview of the key functionality of the instrument, outlining the scientific rationale for its observing modes. The interface between MICADO and the adaptive optics system MAORY that feeds it is summarised. The design of the instrument is discussed, focussing on the optics and mechanisms inside the cryostat, together with a brief overview of the other key sub-systems.

Quasars are associated with and powered by the accretion of material onto massive black holes; the detection of highly luminous quasars with redshifts greater than z = 6 suggests that black holes of up to ten billion solar masses already existed 13 billion years ago. Two possible present-day dormant descendants of this population of active black holes have been found in the galaxies NGC 3842 and NGC 4889 at the centres of the Leo and Coma galaxy clusters, which together form the central region of the Great Wall - the largest local structure of galaxies. The most luminous quasars, however, are not confined to such high-density regions of the early Universe; yet dormant black holes of this high mass have not yet been found outside of modern-day rich clusters. Here we report observations of the stellar velocity distribution in the galaxy NGC 1600 - a relatively isolated elliptical galaxy near the centre of a galaxy group at a distance of 64 Mpc from Earth. We use orbit superposition models to determine that the black hole at the centre of NGC 1600 has a mass of 17 billion solar masses. The spatial distribution of stars near the centre of NGC 1600 is rather diffuse. We find that the region of depleted stellar density in the cores of massive elliptical galaxies extends over the same radius as the gravitational sphere of influence of the central black holes, and interpret this as the dynamical imprint of the black holes.

We investigate the correlations between the black hole mass $M_{BH}$, the velocity dispersion $\sigma$, the bulge mass $M_{Bu}$, the bulge average spherical density $\rho_h$ and its spherical half mass radius $r_h$, constructing a database of 97 galaxies (31 core ellipticals, 17 power-law ellipticals, 30 classical bulges, 19 pseudo bulges) by joining 72 galaxies from the literature to 25 galaxies observed during our recent SINFONI black hole survey. For the first time we discuss the full error covariance matrix. We analyse the well known $M_{BH}-\sigma$ and $M_{BH}-M_{Bu}$ relations and establish the existence of statistically significant correlations between $M_{Bu}$ and $r_h$ and anti-correlations between $M_{Bu}$ and $\rho_h$. We establish five significant bivariate correlations ($M_{BH}-\sigma-\rho_h$, $M_{BH}-\sigma-r_h$, $M_{BH}-M_{Bu}-\sigma$, $M_{BH}-M_{Bu}-\rho_h$, $M_{BH}-M_{Bu}-r_h$) that predict $M_{BH}$ of 77 core and power-law ellipticals and classical bulges with measured and intrinsic scatter as small as $\approx 0.36$ dex and $\approx 0.33$ dex respectively, or 0.26 dex when the subsample of 45 galaxies defined by Kormendy and Ho (2013) is considered. In contrast, pseudo bulges have systematically lower $M_{BH}$, but approach the predictions of all the above relations at spherical densities $\rho_h\ge 10^{10} M_\odot/kpc^3$ or scale lengths $r_h\le 1$ kpc. These findings fit in a scenario of co-evolution of BH and classical-bulge masses, where core ellipticals are the product of dry mergers of power-law bulges and power-law Es and bulges the result of (early) gas-rich mergers and of disk galaxies. In contrast, the (secular) growth of BHs is decoupled from the growth of their pseudo bulge hosts, except when (gas) densities are high enough to trigger the feedback mechanism responsible for the existence of the correlations between $M_{BH}$ and galaxy structural parameters.

Isolating ground-based interferometric gravitational wave observatories from environmental disturbances is one of the great challenges of the advanced detector era. In order to directly observe gravitational waves, the detector components and test masses must be highly inertially decoupled from the ground motion not only to sense the faint strain of space-time induced by gravitational waves, but also to maintain the resonance of the very sensitive 4 km interferometers. This article presents the seismic isolation instrumentation and strategy developed for Advanced LIGO interferometers. It reviews over a decade of research on active isolation in the context of gravitational wave detection, and presents the performance recently achieved with the Advanced LIGO observatory. Lastly, it discusses prospects for future developments in active seismic isolation and the anticipated benefits to astrophysical gravitational wave searches. Beyond gravitational wave research, the goal of this article is to provide detailed isolation strategy guidelines for sensitive ground-based physics experiments that may benefit from similar levels of inertial isolation.

In this work we present CO(1-0) and CO(2-1) observations of a pilot sample of 15 early-type galaxies (ETGs) drawn from the MASSIVE galaxy survey, a volume-limited integral-field spectroscopic study of the most massive ETGs ($M_* >10^{11.5}M_\odot$) within 108 Mpc. These objects were selected because they showed signs of an interstellar medium and/or star formation. A large amount of gas ($>$2$\times$10$^8$ M$_{\odot}$) is present in 10 out of 15 objects, and these galaxies have gas fractions higher than expected based on extrapolation from lower mass samples. We tentatively interpret this as evidence that stellar mass loss and hot halo cooling may be starting to play a role in fuelling the most massive galaxies. These MASSIVE ETGs seem to have lower star-formation efficiencies (SFE=SFR/M$_{\rm H2}$) than spiral galaxies, but the SFEs derived are consistent with being drawn from the same distribution found in other lower mass ETG samples. This suggests that the SFE is not simply a function of stellar mass, but that local, internal processes are more important for regulating star formation. Finally we used the CO line profiles to investigate the high-mass end of the Tully-Fisher relation (TFR). We find that there is a break in the slope of the TFR for ETGs at high masses (consistent with previous studies). The strength of this break correlates with the stellar velocity dispersion of the host galaxies, suggesting it is caused by additional baryonic mass being present in the centre of massive ETGs. We speculate on the root cause of this change and its implications for galaxy formation theories.

MWC 314 is a bright candidate luminous blue variable that resides in a fairly close binary system, with an orbital period of 60.753$\pm$0.003 d. We observed MWC 314 with a combination of optical spectroscopy, broad-band ground- and space-based photometry, as well as with long baseline, near-infrared interferometry. We have revised the single-lined spectroscopic orbit and explored the photometric variability. The orbital light curve displays two minima each orbit that can be partially explained in terms of the tidal distortion of the primary that occurs around the time of periastron. The emission lines in the system are often double-peaked and stationary in their kinematics, indicative of a circumbinary disc. We find that the stellar wind or circumbinary disc is partially resolved in the K\prime-band with the longest baselines of the CHARA Array. From this analysis, we provide a simple, qualitative model in an attempt to explain the observations. From the assumption of Roche Lobe overflow and tidal synchronisation at periastron, we estimate the component masses to be M1 $\approx 5$ M$_\odot$ and M2$\approx 15$ M$_\odot$, which indicates a mass of the LBV that is extremely low. In addition to the orbital modulation, we discovered two pulsational modes with the MOST satellite. These modes are easily supported by a low-mass hydrogen-poor star, but cannot be easily supported by a star with the parameters of an LBV. The combination of these results provides evidence that the primary star was likely never a normal LBV, but rather is the product of binary interactions. As such, this system presents opportunities for studying mass-transfer and binary evolution with many observational techniques.

We examine stellar population gradients in ~100 massive early type galaxies spanning 180 < sigma* < 370 km/s and M_K of -22.5 to -26.5 mag, observed as part of the MASSIVE survey (Ma et al. 2014). Using integral-field spectroscopy from the Mitchell Spectrograph on the 2.7m telescope at McDonald Observatory, we create stacked spectra as a function of radius for galaxies binned by their stellar velocity dispersion, stellar mass, and group richness. With excellent sampling at the highest stellar mass, we examine radial trends in stellar population properties extending to beyond twice the effective radius (~2.5 R_e). Specifically, we examine trends in age, metallicity, and abundance ratios of Mg, C, N, and Ca, and discuss the implications for star formation histories and elemental yields. At a fixed physical radius of 3-6 kpc (the likely size of the galaxy cores formed at high redshift) stellar age and [alpha/Fe] increase with increasing sigma* and depend only weakly on stellar mass, as we might expect if denser galaxies form their central cores earlier and faster. If we instead focus on 1-1.5 R_e, the trends in abundance and abundance ratio are washed out, as might be expected if the stars at large radius were accreted by smaller galaxies. Finally, we show that when controlling for \sigmastar, there are only very subtle differences in stellar population properties or gradients as a function of group richness; even at large radius internal properties matter more than environment in determining star formation history.

We compare the performance of mass estimators for elliptical galaxies that rely on the directly observable surface brightness and velocity dispersion profiles, without invoking computationally expensive detailed modeling. These methods recover the mass at a specific radius where the mass estimate is expected to be least sensitive to the anisotropy of stellar orbits. One method (Wolf et al. 2010) uses the total luminosity-weighted velocity dispersion and evaluates the mass at a 3D half-light radius $r_{1/2}$, i.e., it depends on the GLOBAL galaxy properties. Another approach (Churazov et al. 2010) estimates the mass from the velocity dispersion at a radius $R_2$ where the surface brightness declines as $R^{-2}$, i.e., it depends on the LOCAL properties. We evaluate the accuracy of the two methods for analytical models, simulated galaxies and real elliptical galaxies that have already been modeled by the Schwarzschild's orbit-superposition technique. Both estimators recover an almost unbiased circular speed estimate with a modest RMS scatter ($\lesssim 10 \%$). Tests on analytical models and simulated galaxies indicate that the local estimator has a smaller RMS scatter than the global one. We show by examination of simulated galaxies that the projected velocity dispersion at $R_2$ could serve as a good proxy for the virial galaxy mass. For simulated galaxies the total halo mass scales with $\sigma_p(R_2)$ as $M_{vir} \left[M_{\odot}h^{-1}\right] \approx 6\cdot 10^{12} \left( \frac{\sigma_p(R_2)}{200\, \rm km\, s^{-1}} \right)^{4}$ with RMS scatter $\approx 40 \%$.

Under cosmic irradiation, the interstellar water ice mantles evolve towards a compact amorphous state. Crystalline ice amorphisation was previously monitored mainly in the keV to hundreds of keV ion energies. We experimentally investigate heavy ion irradiation amorphisation of crystalline ice, at high energies closer to true cosmic rays, and explore the water-ice sputtering yield. We irradiated thin crystalline ice films with MeV to GeV swift ion beams, produced at the GANIL accelerator. The ice infrared spectral evolution as a function of fluence is monitored with in-situ infrared spectroscopy (induced amorphisation of the initial crystalline state into a compact amorphous phase). The crystalline ice amorphisation cross-section is measured in the high electronic stopping-power range for different temperatures. At large fluence, the ice sputtering is measured on the infrared spectra, and the fitted sputtering-yield dependence, combined with previous measurements, is quadratic over three decades of electronic stopping power. The final state of cosmic ray irradiation for porous amorphous and crystalline ice, as monitored by infrared spectroscopy, is the same, but with a large difference in cross-section, hence in time scale in an astrophysical context. The cosmic ray water-ice sputtering rates compete with the UV photodesorption yields reported in the literature. The prevalence of direct cosmic ray sputtering over cosmic-ray induced photons photodesorption may be particularly true for ices strongly bonded to the ice mantles surfaces, such as hydrogen-bonded ice structures or more generally the so-called polar ices.

We study nine S0-Sb galaxies with (photometric) bulges consisting of two distinct components. The outer component is a flattened, kinematically cool, disklike structure: a "disky pseudobulge". Embedded inside is a rounder, kinematically hot spheroid: a "classical bulge". This indicates that pseudobulges and classical bulges are not mutually exclusive: some galaxies have both. The disky pseudobulges almost always have an exponential disk (scale lengths = 125-870 pc, mean $\sim 440$ pc) with disk-related subcomponents: nuclear rings, bars, and/or spiral arms. They constitute 11-59% of the galaxy stellar mass (mean PB/T = 0.33), with stellar masses $\sim 7 \times 10^{9}$-$9 \times 10^{10} M_{\odot}$. Classical-bulge components have Sersic indices of 0.9-2.2, effective radii of 25-430 pc and stellar masses of $5 \times 10^{8}$-$3 \times 10^{10} M_{\odot}$ (usually < 10% of the galaxy's stellar mass; mean B/T = 0.06). The classical bulges show rotation, but are kinematically hotter than the disky pseudobulges. Dynamical modeling of three systems indicates that velocity dispersions are isotropic in the classical bulges and equatorially biased in the disky pseudobulges. In the mass--radius and mass--stellar mass density planes, classical-bulge components follow sequences defined by ellipticals and (larger) classical bulges. Disky pseudobulges also fall on this sequence; they are more compact than similar-mass large-scale disks. Although some classical bulges are quite compact, they are distinct from nuclear star clusters in both size and mass, and coexist with nuclear clusters in at least two galaxies. Since almost all the galaxies in this study are barred, they probably also host boxy/peanut-shaped bulges (vertically thickened inner parts of bars). NGC 3368 shows evidence for such a zone outside its disky pseudobulge, making it a galaxy with all three types of "bulge".

Massive early-type galaxies represent the modern-day remnants of the earliest major star formation episodes in the history of the universe. These galaxies are central to our understanding of the evolution of cosmic structure, stellar populations, and supermassive black holes, but the details of their complex formation histories remain uncertain. To address this situation, we have initiated the MASSIVE Survey, a volume-limited, multi-wavelength, integral-field spectroscopic (IFS) and photometric survey of the structure and dynamics of the ~100 most massive early-type galaxies within a distance of 108 Mpc. This survey probes a stellar mass range M* > 10^{11.5} Msun and diverse galaxy environments that have not been systematically studied to date. Our wide-field IFS data cover about two effective radii of individual galaxies, and for a subset of them, we are acquiring additional IFS observations on sub-arcsecond scales with adaptive optics. We are also acquiring deep K-band imaging to trace the extended halos of the galaxies and measure accurate total magnitudes. Dynamical orbit modeling of the combined data will allow us to simultaneously determine the stellar, black hole, and dark matter halo masses. The primary goals of the project are to constrain the black hole scaling relations at high masses, investigate systematically the stellar initial mass function and dark matter distribution in massive galaxies, and probe the late-time assembly of ellipticals through stellar population and kinematical gradients. In this paper, we describe the MASSIVE sample selection, discuss the distinct demographics and structural and environmental properties of the selected galaxies, and provide an overview of our basic observational program, science goals and early survey results.

Observations of late-type main-sequence stars have revealed empirical scalings of coronal activity versus rotation period or Rossby number $Ro$ (a ratio of rotation period to convective turnover time) which has hitherto lacked explanation. For $Ro >> 1$, the activity observed as X-ray to bolometric flux varies as $Ro^{-q}$ with $2\le q \le 3$, whilst $|q| < 0.12$ for $Ro << 1$. Here we explain the transition between these two regimes and the power law in the $Ro >> 1$ regime by constructing an expression for the coronal luminosity based on dynamo magnetic field generation and magnetic buoyancy. We explain the $Ro<<1$ behavior from the inference that observed rotation is correlated with internal differential rotation and argue that once the shear time scale is shorter than the convective turnover time, eddies will be shredded on the shear time scale and so the eddy correlation time actually becomes the shear time and the convection time drops out of the equations. We explain the $Ro >> 1$ behavior using a dynamo saturation theory based on magnetic helicity buildup and buoyant loss.

Detailed imaging and spectroscopic analysis of the centers of nearby S0 and spiral galaxies shows the existence of "composite bulges", where both classical bulges and disky pseudobulges coexist in the same galaxy. As part of a search for supermassive black holes in nearby galaxy nuclei, we obtained VLT-SINFONI observations in adaptive-optics mode of several of these galaxies. Schwarzschild dynamical modeling enables us to disentangle the stellar orbital structure of the different central components, and to distinguish the differing contributions of kinematically hot (classical bulge) and kinematically cool (pseudobulge) components in the same galaxy.

This paper presents the results of the past seven years of experimental investigation and testing done on the two-stage twelve-axis vibration isolation platform for Advanced LIGO gravity waves observatories. This five-ton two-and-half-meter wide system supports more than a 1000 kg of very sensitive equipment. It provides positioning capability and seismic isolation in all directions of translation and rotation. To meet the very stringent requirements of Advanced LIGO, the system must provide more than three orders of magnitude of isolation over a very large bandwidth. It must bring the motion below 10^(-11) m/(Hz)^0.5 at 1 Hz and 10^(-12) m/(Hz)^0.5 at 10 Hz. A prototype of this system has been built in 2006. It has been extensively tested and analyzed during the following two years. This paper shows how the experimental results obtained with the prototype were used to engineer the final design. It highlights how the engineering solutions implemented not only improved the isolation performance but also greatly simplified the assembly, testing, and commissioning process. During the past two years, five units have been constructed, tested, installed and commissioned at each of the two LIGO observatories. Five other units are being built for an upcoming third observatory. The test results presented show that the system meets the motion requirements, and reach the sensor noise in the control bandwidth.

A sample of 1.53$\times$10$^{9}$ cosmic-ray-induced single muon events has been recorded at 225 meters-water-equivalent using the MINOS Near Detector. The underground muon rate is observed to be highly correlated with the effective atmospheric temperature. The coefficient $\alpha_{T}$, relating the change in the muon rate to the change in the vertical effective temperature, is determined to be 0.428$\pm$0.003(stat.)$\pm$0.059(syst.). An alternative description is provided by the weighted effective temperature, introduced to account for the differences in the temperature profile and muon flux as a function of zenith angle. Using the latter estimation of temperature, the coefficient is determined to be 0.352$\pm$0.003(stat.)$\pm$0.046(syst.).

We present new constraints on the density profiles of dark matter (DM) halos in seven nearby dwarf galaxies from measurements of their integrated stellar light and gas kinematics. The gas kinematics of low mass galaxies frequently suggest that they contain constant density DM cores, while N-body simulations instead predict a cuspy profile. We present a data set of high resolution integral field spectroscopy on seven galaxies and measure the stellar and gas kinematics simultaneously. Using Jeans modeling on our full sample, we examine whether gas kinematics in general produce shallower density profiles than are derived from the stars. Although 2/7 galaxies show some localized differences in their rotation curves between the two tracers, estimates of the central logarithmic slope of the DM density profile, gamma, are generally robust. The mean and standard deviation of the logarithmic slope for the population are gamma=0.67+/-0.10 when measured in the stars and gamma=0.58+/-0.24 when measured in the gas. We also find that the halos are not under concentrated at the radii of half their maximum velocities. Finally, we search for correlations of the DM density profile with stellar velocity anisotropy and other baryonic properties. Two popular mechanisms to explain cored DM halos are an exotic DM component or feedback models that strongly couple the energy of supernovae into repeatedly driving out gas and dynamically heating the DM halos. We investigate correlations that may eventually be used to test models. We do not find a secondary parameter that strongly correlates with the central DM density slope, but we do find some weak correlations. Determining the importance of these correlations will require further model developments and larger observational samples. (Abridged)

Most Milky Way globular clusters (GCs) exhibit measurable flattening, even if on a very low level. Both cluster rotation and tidal fields are thought to cause this flattening. Nevertheless, rotation has only been confirmed in a handful of GCs, based mostly on individual radial velocities at large radii. We are conducting a survey of the central kinematics of Galactic GCs using the new Integral Field Unit instrument VIRUS-W. We detect rotation in all 11 GCs that we have observed so far, rendering it likely that a large majority of the Milky Way GCs rotate. We use published catalogs of the ACS survey of GCs to derive central ellipticities and position angles. We show that in all cases where the central ellipticity permits an accurate measurement of the position angle, those angles are in excellent agreement with the kinematic position angles that we derive from the VIRUS-W velocity fields. We find an unexpected tight correlation between central rotation and outer ellipticity, indicating that rotation drives flattening for the objects in our sample. We also find a tight correlation between central rotation and published values for the central velocity dispersion, most likely due to rotation impacting the old dispersion measurements.

Previous studies have reported the existence of two counter-rotating stellar disks in the early-type spiral galaxy NGC7217. We have obtained high-resolution optical spectroscopic data (R ~ 9000) with the new fiber-based Integral Field Unit instrument VIRUS-W at the 2.7m telescope of the McDonald Observatory in Texas. Our analysis confirms the existence of two components. However, we find them to be co-rotating. The first component is the more luminous (~ 77% of the total light), has the higher velocity dispersion (~ 170 km/s) and rotates relatively slowly (projected $v_{max}$ = 50 km/s). The lower luminosity second component, (~ 23% of the total light), has a low velocity dispersion (~ 20 km/s) and rotates quickly (projected $v_{max}$ = 150 km/s). The difference in the kinematics of the two stellar components allows us to perform a kinematic decomposition and to measure the strengths of their Mg and Fe Lick indices separately. The rotational velocities and dispersions of the less luminous and faster component are very similar to those of the interstellar gas as measured from the [OIII] emission. Morphological evidence of active star formation in this component further suggests that NGC7217 may be in the process of (re)growing a disk inside a more massive and higher dispersion stellar halo. The kinematically cold and regular structure of the gas disk in combination with the central almost dust-free morphology allows us to compare the dynamical mass inside of the central 500pc with predictions from a stellar population analysis. We find agreement between the two if a Kroupa stellar initial mass function is assumed.

The most massive elliptical galaxies have low-density centers or cores that differ dramatically from the high-density centers of less massive ellipticals and bulges of disk galaxies. These cores have been interpreted as the result of mergers of supermassive black hole binaries, which depopulate galaxy centers by gravitationally slingshotting central stars toward large radii. Such binaries naturally form in mergers of luminous galaxies. Here, we analyze the population of central stellar orbits in 11 massive elliptical galaxies that we observed with the integral field spectrograph SINFONI at the European Southern Observatory Very Large Telescope. Our dynamical analysis is orbit-based and includes the effects of a central black hole, the mass distribution of the stars, and a dark matter halo. We show that the use of integral field kinematics and the inclusion of dark matter is important to conclude upon the distribution of stellar orbits in galaxy centers. Six of our galaxies are core galaxies. In these six galaxies, but not in the galaxies without cores, we detect a coherent lack of stars on radial orbits in the core region and a uniform excess of radial orbits outside of it: when scaled by the core radius, the radial profiles of the classical anisotropy parameter beta are nearly identical in core galaxies. Moreover, they match quantitatively the predictions of black hole binary simulations, providing the first convincing dynamical evidence for core scouring in the most massive elliptical galaxies.

We present an analysis of the H2 emission-line gas kinematics in the inner < 4 arcsec radius of six nearby spiral galaxies, based on AO-assisted integral-field observations obtained in the K-band with SINFONI/VLT. Four of the six galaxies in our sample display ordered H2 velocity fields, consistent with gas moving in the plane of the galaxy and rotating in the same direction as the stars. However, the gas kinematics is typically far from simple circular motion. We can classify the observed velocity fields into four different types of flows, ordered by increasing complexity: (1) circular motion in a disc (NGC3351); (2) oval motion in the galaxy plane (NGC3627 and NGC4536); (3) streaming motion superimposed on circular rotation (NGC4501); and (4) disordered streaming motions (NGC4569 and NGC4579). The H2 velocity dispersion in the galaxies is usually higher than 50 km/s in the inner 1-2 arcsec radii. The four galaxies with ordered kinematics have v/sigma < 1 at radii less than 40-80 pc. The radius at which v/sigma = 1 is independent of the type of nuclear activity. While the low values of v/sigma could be taken as an indication of a thick disc in the innermost regions of the galaxies, other lines of evidence (e.g. H2 morphologies and velocity fields) argue for a thin disc interpretation in the case of NGC3351 and NGC4536. We discuss the implications of the high values of velocity dispersion for the dynamics of the gaseous disc and suggest caution when interpreting the velocity dispersion of ionized and warm tracers as being entirely dynamical. Understanding the nature and role of the velocity dispersion in the gas dynamics, together with the full 2D information of the gas, is essential for obtaining accurate black hole masses from gas kinematics.

Shallow cores in bright, massive galaxies are commonly thought to be the result of scouring of stars by mergers of binary supermassive black holes. Past investigations have suggested correlations between the central black hole mass and the stellar light or mass deficit in the core, using proxy measurements of $M_{\rm BH}$ or stellar mass-to-light ratios ($\Upsilon$). Drawing on a wealth of dynamical models which provide both $M_{\rm BH}$ and $\Upsilon$, we identify cores in 23 galaxies, of which 20 have direct, reliable measurements of $M_{\rm BH}$ and dynamical stellar mass-to-light ratios ($\Upsilon_{\star,{\rm dyn}}$). These cores are identified and measured using Core-S\'ersic model fits to surface brightness profiles which extend out to large radii (typically more than the effective radius of the galaxy); for approximately one fourth of the galaxies, the best fit includes an outer (\sersic) envelope component. We find that the core radius is most strongly correlated with the black hole mass and that it correlates better with total galaxy luminosity than it does with velocity dispersion. The strong core-size-- $M_{\rm BH}$ correlation enables estimation of black hole masses (in core galaxies) with an accuracy comparable to the $M_{\rm BH}$--$\sigma$ relation (rms scatter of 0.30 dex in $\log M_{\rm BH}$), without the need for spectroscopy. The light and mass deficits correlate more strongly with galaxy velocity dispersion than they do with black hole mass. Stellar mass deficits span a range of 0.2--39 \mbh, with almost all (87%) being $< 10 \, M_{\rm BH}$; the median value is 2.2 $M_{\rm BH}$.

Adaptive-Optics assisted SINFONI observations of the central regions of ten early-type galaxies are presented. Based primarily on the SINFONI kinematics, ten black hole masses occupying the high-mass regime of the M_BH-sigma relation are derived using three-integral Schwarzschild models. The effect of dark matter inclusion on the black hole mass is explored. The omission of a dark matter halo in the model results in a higher stellar mass-to-light ratio, especially when extensive kinematic data are used in the model. However, when the diameter of the sphere of influence -- computed using the black hole mass derived without a dark halo -- is at least 10 times the PSF FWHM during the observations, it is safe to exclude a dark matter component in the dynamical modeling, i.e. the change in black hole mass is negligible. When the spatial resolution is marginal, restricting the mass-to-light ratio to the right value returns the correct M_BH although dark halo is not present in the model. Compared to the M_BH-sigma and M_BH-L relations of McConnell et al. (2011), the ten black holes are all more massive than expected from the luminosities and seven black hole masses are higher than expected from the stellar velocity dispersions of the host bulges. Using new fitted relations which include the ten galaxies, we find that the space density of the most massive black holes (M_BH >~ 10^9 Msun) estimated from the M_BH-L relation is higher than the estimate based on the M_BH-sigma relation and the latter is higher than model predictions based on quasar counts, each by about an order of magnitude.

Dynamical modeling and strong lensing data indicate that the total density profiles of early-type galaxies are close to isothermal, i.e., rho_tot ~ r^gamma with gamma approx -2. To understand the origin of this universal slope we study a set of simulated spheroids formed in isolated binary mergers as well as the formation within the cosmological framework. The total stellar plus dark matter density profiles can always be described by a power law with an index of gamma approx -2.1 with a tendency toward steeper slopes for more compact, lower-mass ellipticals. In the binary mergers the amount of gas involved in the merger determines the precise steepness of the slope. This agrees with results from the cosmological simulations where ellipticals with steeper slopes have a higher fraction of stars formed in situ. Each gas-poor merger event evolves the slope toward gamma ~ -2, once this slope is reached further merger events do not change it anymore. All our ellipticals have flat intrinsic combined stellar and dark matter velocity dispersion profiles. We conclude that flat velocity dispersion profiles and total density distributions with a slope of gamma ~ -2 for the combined system of stars and dark matter act as a natural attractor. The variety of complex formation histories as present in cosmological simulations, including major as well as minor merger events, is essential to generate the full range of observed density slopes seen for present-day elliptical galaxies.

In this paper we analyse a set of 33 optical spectra, which were acquired with the ARCES echelle spectrograph (R = 38,000) on the 3.5-m telescope at the Apache Point Observatory. We examine the H{\alpha} profile in each of these observations in order to determine the geometry and velocity structure of the previously discovered bipolar jet, which originates from the secondary star of HD 44179 located at the centre of the Red Rectangle nebula. Using a 3D geometric model we are able to determine the orbital coverage during which the jet occults the primary star. During the occultation, part of the H{\alpha} line profile appears in absorption. The velocity structure of the jet was determined by modelling the absorption line profile using the Sobolev approximation for each orbital phase during which we have observations. The results indicate the presence of a wide angle jet, likely responsible for observed biconical structure of the outer nebula. Furthermore, we were able to determine a likely velocity structure and rule out several others. We find that the jet is comprised of low-density, high-velocity, central region and a higher-density, lower-velocity, conical shell.

We studied the stellar populations, distribution of dark matter, and dynamical structure of a sample of 25 early-type galaxies in the Coma and Abell 262 clusters. We derived dynamical mass-to-light ratios and dark matter densities from orbit-based dynamical models, complemented by the ages, metallicities, and \alpha-elements abundances of the galaxies from single stellar population models. Most of the galaxies have a significant detection of dark matter and their halos are about 10 times denser than in spirals of the same stellar mass. Calibrating dark matter densities to cosmological simulations we find assembly redshifts z_{DM} \approx 1-3. The dynamical mass that follows the light is larger than expected for a Kroupa stellar initial mass function, especially in galaxies with high velocity dispersion \sigma_{eff} inside the effective radius r_{eff}. We now have 5 of 25 galaxies where mass follows light to 1-3 r_{eff}, the dynamical mass-to-light ratio of all the mass that follows the light is large (\approx 8-10 in the Kron-Cousins R band), the dark matter fraction is negligible to 1-3 r_{eff}. This could indicate a "massive" initial mass function in massive early-type galaxies. Alternatively, some of the dark matter in massive galaxies could follow the light very closely suggesting a significant degeneracy between luminous and dark matter.

We present the first results of an analysis of the properties of the molecular gas in the nuclear regions (r < 300 pc) of a sample of six nearby galaxies, based on new high spatial resolution observations obtained in the K-band with the near-infrared integral field spectrograph SINFONI at the Very Large Telescope. We derive two-dimensional distributions of the warm molecular and ionized gas from the H2, Br_gamma and HeI emission lines present in the spectra of the galaxies. We find a range of morphologies, including bar- and ring-like distributions and either centrally peaked or off-centre emission. The morphologies of the molecular and the ionized gas are not necessarily coincident. The observed emission-line ratios point towards thermal processes as the principal mechanism responsible for the H2 excitation in the nuclear and circumnuclear regions of the galaxies, independently of the presence of an active nucleus. We find that a rescaling of the H2 2.12 microns emission-line luminosity by a factor beta~1200 gives a good estimate (within a factor of 2) of the total (cold) molecular gas mass. The galaxies of the sample contain large quantities of molecular gas in their centres, with total masses in the ~ 105 - 108 Msol range. Never the less, these masses correspond to less than 3 per cent of the stellar masses derived for the galaxies in these regions, indicating that the presence of gas should not affect black hole mass estimates based on the dynamical modelling of the stars. The high-spatial resolution provided by the SINFONI data allowed us to resolve a circumnuclear ring (with a radius of ~270 pc) in the galaxy NGC 4536. The measured values of the Br_gamma equivalent width and the HeI/Br_gamma emission-line ratio suggests that bursts of star formation occurred throughout this ring as recently as 6.5 Myr ago.

We present radially resolved spectroscopy of 8 early-type galaxies in Abell~262, measuring rotation, velocity dispersion, $H_3$ and $H_4$ coefficients along three axes, and line-strength index profiles of Mg, Fe and H$\beta$. Ionized-gas velocity and velocity dispersion is included for 6 galaxies. We derive dynamical mass-to-light ratios and dark matter densities from orbit-based dynamical models, complemented by the galaxies' ages, metallicities, and $\alpha$-elements abundances. Four galaxies have significant dark matter with halos about 10 times denser than in spirals of the same stellar mass. Using dark matter densities and cosmological simulations, assembly redshifts $\zdm\approx 1-3$, which we found earlier for Coma. The dynamical mass following the light is larger than expected for a Kroupa stellar IMF, especially in galaxies with high velocity dispersion $\sigeff$ inside the effective radius $\reff$. This could indicate a `massive' IMF in massive galaxies. Alternatively, some dark matter in massive galaxies could follow the light closely. Combining with our comparison sample of Coma early-types, we now have 5 of 24 galaxies where (1) mass follows light to $1-3\,\reff$, (2) the dynamical mass-to-light ratio {of all the mass that follows the light is large ($\approx\,8-10$ in the Kron-Cousins $R$ band), (3) the dark matter fraction is negligible to $1-3\,\reff$. Unless the IMF in these galaxies is particularly `massive' and somehow coupled to the dark matter content, there seems a significant degeneracy between luminous and dark matter in some early-type galaxies. The role of violent relaxation is briefly discussed.

We revisit collisionless major and minor mergers of spheroidal galaxies in the context of the size evolution of elliptical galaxies. The simulations are performed as a series of mergers with mass-ratios of 1:1 and 1:10 for models representing pure bulges as well as bulges embedded in dark matter halos. For major and minor mergers, respectively, we identify and analyze two different processes, violent relaxation and stripping, leading to size evolution and a change of the dark matter fraction within the observable effective radius. Violent relaxation - which is the dominant mixing process for major mergers but less important for minor mergers - scatters relatively more dark matter particles than bulge particles to small. Stripping in minor mergers assembles stellar satellite particles at large radii in halo dominated regions of the massive host. This strongly increases the size of the bulge into regions with higher dark matter fractions leaving the inner host structure almost unchanged. A factor of two mass increase by minor mergers increases the dark matter fraction by 20 per cent. We present analytic corrections to simple one-component virial estimates for the evolution of the gravitational radii. If such a two-component system grows by minor mergers alone its size growth, $r_{\mathrm{e}} \propto M^\alpha$, reaches values of $\alpha \approx 2.4$, significantly exceeding the simple theoretical limit of $\alpha = 2$. For major mergers the sizes grow with $\alpha \lesssim 1$. Our results indicate that minor mergers of galaxies embedded in massive dark matter halos provide a potential mechanism for explaining the rapid size growth and the build-up of massive elliptical systems predicting significant dark matter fractions and radially biased velocity dispersions at large radii (abbreviated)

This work aims to study the distribution of luminous and dark matter in Coma early-type galaxies. Dynamical masses obtained under the assumption that mass follows light do not match with the masses of strong gravitational lens systems of similar velocity dispersions. Instead, dynamical fits with dark matter halos are in good agreement with lensing results. We derive mass-to-light ratios of the stellar populations from Lick absorption line indices, reproducing well the observed galaxy colours. Even in dynamical models with dark matter halos the amount of mass that follows the light increases more rapidly with galaxy velocity dispersion than expected for a constant stellar initial mass function (IMF). While galaxies around sigma ~ 200 km/s are consistent with a Kroupa IMF, the same IMF underpredicts luminous dynamical masses of galaxies with sigma ~ 300 km/s by a factor of two and more. A systematic variation of the stellar IMF with galaxy velocity dispersion could explain this trend with a Salpeter IMF for the most massive galaxies. If the IMF is instead constant, then some of the dark matter in high velocity dispersion galaxies must follow a spatial distribution very similar to that of the light. A combination of both, a varying IMF and a component of dark matter that follows the light is possible as well. For a subsample of galaxies with old stellar populations we show that the tilt in the fundamental plane can be explained by systematic variations of the total (stellar + dark) mass inside the effective radius. We tested commonly used mass estimator formulae, finding them accurate at the 20-30% level.

In this paper we examine the profiles of the complex Na I D-lines in the Red Rectangle. The spectra were acquired with the ARCES \'{e}chelle spectrograph R = 38,000 on the 3.5-m telescope at the Apache Point Observatory. Additional spectra taken with STIS were acquired from the Hubble Legacy Archive (HLA) and were used to independently confirm the spatial origin of the spectral features. The profile of a single D-line consists of double-peaked emission, red-shifted absorption and blue-shifted absorption. We find that the double-peaked emission originates from the bipolar outflow, the red-shifted absorption feature is due to the photospheric line, and the blue-shifted absorption arises from the bipolar outflow as seen against the photosphere of the luminous post-AGB component in HD 44179. In order to better understand the Na I D-line profile, we examined the periodically variable asymmetric photospheric absorption lines. The asymmetric lines are interpreted as a signature of slow self-accretion following enhanced mass-loss around periastron. An empirical model was constructed to remove the photospheric component from the Na I D-line profile in order to study the nebular emission and absorption of sodium along the line-of-sight to the primary. This paper also discusses the different origins of the single-peaked emission, the double-peaked emission and the blue-shifted and red-shifted absorption components.

Using a slitless spectroscopy method with (a) the 8.2 m Subaru telescope and its FOCAS Cassegrain spectrograph, and (b) the ESO Very Large Telescope (VLT) unit 1 (Antu) and its FORS2 Cassegrain spectrograph, we have detected 326 planetary nebulae (PNs) in the giant Virgo elliptical galaxy NGC 4649 (M 60), and we have measured their radial velocities. After rejecting some PNs more likely to belong to the companion galaxy NGC 4647, we have built a catalog with kinematic information for 298 PNs in M 60. Using these radial velocities we have concluded that they support the presence of a dark matter halo around M 60. Based on an isotropic, two-component Hernquist model, we estimate the dark matter halo mass within 3$R_{\rm e}$ to be 4$\times10^{11} M_{\odot}$, which is almost one half of the total mass of about $10^{12} M_{\odot}$ within 3$R_{\rm e}$. This total mass is similar to that estimated from globular cluster, XMM-Newton and Chandra observations. The dark matter becomes dominant outside. More detailed dynamical modeling of the PN data is being published in a companion paper. We have also measured the $m$(5007) magnitudes of many of these PNs, and built a statistically complete sample of 218 PNs. The resulting PN luminosity function (PNLF) was used to estimate a distance modulus of 30.7$\pm$0.2 mag, equivalent to 14$\pm$1 Mpc. This confirms an earlier PNLF distance measurement, based on a much smaller sample. The PNLF distance modulus remains smaller than the surface brightness fluctuation (SBF) distance modulus by 0.4 mag.

The masses of the most massive supermassive black holes (SMBHs) predicted by the M_BH-sigma and M_BH-luminosity relations appear to be in conflict. Which of the two relations is the more fundamental one remains an open question. NGC 1332 is an excellent example that represents the regime of conflict. It is a massive lenticular galaxy which has a bulge with a high velocity dispersion sigma of ~320 km/s; bulge--disc decomposition suggests that only 44% of the total light comes from the bulge. The M_BH-sigma and the M_BH-luminosity predictions for the central black hole mass of NGC 1332 differ by almost an order of magnitude. We present a stellar dynamical measurement of the SMBH mass using an axisymmetric orbit superposition method. Our SINFONI integral-field unit (IFU) observations of NGC 1332 resolve the SMBH's sphere of influence which has a diameter of ~0.76 arcsec. The sigma inside 0.2 arcsec reaches ~400 km/s. The IFU data allow us to increase the statistical significance of our results by modelling each of the four quadrants separately. We measure a SMBH mass of (1.45 \pm 0.20) x 10^9 M_sun with a bulge mass-to-light ratio of 7.08 \pm 0.39 in the R-band. With this mass, the SMBH of NGC 1332 is offset from the M_BH-luminosity relation by a full order of magnitude but is consistent with the M_BH-sigma relation.

This article describes the Schwarzschild orbit superposition method. It is the state-of-the-art dynamical modelling tool for early-type galaxies. Tests with analytic models show that masses and orbital anisotropies of not too face-on galaxies can be recovered with about 15 percent accuracy from typical observational data. Applying Schwarzschild models to a sample of Coma galaxies their dark matter halos were found to be 13 times denser than those of spirals with the same stellar mass. Since denser halos assembled earlier, this result indicates that the formation redshift 1+z of ellipticals is about two times higher than of spirals. Roughly half of the sample galaxies have halo assembly redshifts in agreement with their stellar-population ages. Galaxies where stars appear younger than the halos show strong phase-space density gradients in their orbital structure, indicative for dissipational evolution and possibly connected with secondary star-formation after the main halo assembly epoch. The importance of considering dark-matter in dynamical models aimed to measure black-hole masses is briefly discussed.

We present high quality long slit spectra along the major and minor axes out to 1.5-2 Re (14-22 kpc) of three bright elliptical galaxies (NGC1600, NGC4125, NGC7619) obtained at the Hobby-Eberly Telescope (HET). We derive stellar kinematic profiles and Lick/IDS indices (Hbeta, Mgb, Fe5015, Fe5270, Fe5335, Fe5406). Moreover, for NGC4125 we derive gas kinematics and emission line strengths. We model the absorption line strengths using Simple Stellar Populations models that take into account the variation of [\alpha/Fe] and derive ages, total metallicity and element abundances. Overall, we find that the three galaxies have old and [\alpha/Fe] overabundant stellar populations with no significant gradients. The metallicity is supersolar at the center with a strong negative radial gradient. For NGC4125, several pieces of evidence point to a recent dissipational merger event. We calculate the broad band color profiles with the help of SSP models. All of the colors show sharp peaks at the center of the galaxies, mainly caused by the metallicity gradients, and agree well with the measured colors. Using the Schwarzschild's axisymmetric orbit superposition technique, we model the stellar kinematics to constrain the dark halos of the galaxies. We use the tight correlation between the Mgb strength and local escape velocity to set limits on the extent of the halos by testing different halo sizes. Logarithmic halos - cut at 60 kpc -minimize the overall scatter of the Mgb-Vesc relation. Larger cutoff radii are found if the dark matter density profile is decreasing more steeply at large radii.

We report studies of ultra-high energy cosmic ray composition via analysis of depth of airshower maximum (Xmax), for airshower events collected by the High Resolution Fly's Eye (HiRes) observatory. The HiRes data are consistent with a constant elongation rate d<Xmax>/d(log(E)) of 47.9 +- 6.0 (stat.) +- 3.2 (syst.)g/cm^2/decade for energies between 1.6 EeV and 63 EeV, and are consistent with a predominantly protonic composition of cosmic rays when interpreted via the QGSJET01 and QGSJET-II high-energy hadronic interaction models. These measurements constrain models in which the galactic-to-extragalactic transition is the cause of the energy spectrum "ankle' at 4 EeV.

Stereo data collected by the HiRes experiment over a six year period are examined for large-scale anisotropy related to the inhomogeneous distribution of matter in the nearby Universe. We consider the generic case of small cosmic-ray deflections and a large number of sources tracing the matter distribution. In this matter tracer model the expected cosmic ray flux depends essentially on a single free parameter, the typical deflection angle theta. We find that the HiRes data with threshold energies of 40 EeV and 57 EeV are incompatible with the matter tracer model at a 95% confidence level unless theta is larger than 10 degrees and are compatible with an isotropic flux. The data set above 10 EeV is compatible with both the matter tracer model and an isotropic flux.

It is now well established that all galaxies with a massive bulge component harbour a central supermassive black hole (SMBH). The mass of the SMBH correlates with bulge properties such as the bulge mass and the velocity dispersion, which implies that the bulge and the SMBH of a galaxy have grown together during the formation process. The spiral galaxy NGC3368 and the S0 galaxy NGC3489 both host a pseudobulge and a much smaller classical bulge component at the centre. We present high resolution, near-infrared IFU data of these two galaxies, taken with SINFONI at the VLT, and use axisymmetric orbit models to determine the masses of the SMBHs. The SMBH mass of NGC3368 is M_BH=7.5x10^6 M_sun with an error of 1.5x10^6 M_sun, which mostly comes from the non-axisymmetry in the data. For NGC3489, a solution without black hole cannot be excluded when modelling the SINFONI data alone, but can be clearly ruled out when modelling a combination of SINFONI, OASIS and SAURON data, for which we obtain M_BH=6.00^{+0.56}_{-0.54} (stat) +/- 0.64 (sys) x 10^6 M_sun. Although both galaxies seem to be consistent with the M_BH-sigma relation, at face value they do not agree with the relation between bulge magnitude and black hole mass when the total bulge magnitude (i.e., including both classical bulge and pseudobulge) is considered; the agreement is better when only the small classical bulge components are considered. However, taking into account the ageing of the stellar population could change this conclusion.

We model the dynamical structure of M87 (NGC4486) using high spatial resolution long-slit observations of stellar light in the central regions, two-dimensional stellar light kinematics out to half of the effective radius, and globular cluster velocities out to 8 effective radii. We simultaneously fit for four parameters, black hole mass, dark halo core radius, dark halo circular velocity, and stellar mass-to-light ratio. We find a black hole mass of 6.4(+-0.5)x10^9 Msun(the uncertainty is 68% confidence marginalized over the other parameters). The stellar M/L_V=6.3+-0.8. The best-fitted dark halo core radius is 14+-2 kpc, assuming a cored logarithmic potential. The best-fitted dark halo circular velocity is 715+-15 km/s. Our black hole mass is over a factor of 2 larger than previous stellar dynamical measures, and our derived stellar M/L ratio is 2 times lower than previous dynamical measures. When we do not include a dark halo, we measure a black hole mass and stellar M/L ratio that is consistent with previous measures, implying that the major difference is in the model assumptions. The stellar M/L ratio from our models is very similar to that derived from stellar population models of M87. The reason for the difference in the black hole mass is because we allow the M/L ratio to change with radius. The dark halo is degenerate with the stellar M/L ratio, which is subsequently degenerate with the black hole mass. We argue that dynamical models of galaxies that do not include the contribution from a dark halo may produce a biased result for the black hole mass. This bias is especially large for a galaxy with a shallow light profile such as M87, and may not be as severe in galaxies with steeper light profiles unless they have a large stellar population change with radius.

The High Resolution Fly's Eye experiment has measured the flux of ultrahigh energy cosmic rays using the stereoscopic air fluorescence technique. The HiRes experiment consists of two detectors that observe cosmic ray showers via the fluorescence light they emit. HiRes data can be analyzed in monocular mode, where each detector is treated separately, or in stereoscopic mode where they are considered together. Using the monocular mode the HiRes collaboration measured the cosmic ray spectrum and made the first observation of the Greisen-Zatsepin-Kuzmin cutoff. In this paper we present the cosmic ray spectrum measured by the stereoscopic technique. Good agreement is found with the monocular spectrum in all details.

Recent progress in theoretical modeling of a sunspot is reviewed. The observed properties of umbral dots are well reproduced by realistic simulations of magnetoconvection in a vertical, monolithic magnetic field. To understand the penumbra, it is useful to distinguish between the inner penumbra, dominated by bright filaments containing slender dark cores, and the outer penumbra, made up of dark and bright filaments of comparable width with corresponding magnetic fields differing in inclination by some 30 degrees and strong Evershed flows in the dark filaments along nearly horizontal or downward-plunging magnetic fields. The role of magnetic flux pumping in submerging magnetic flux in the outer penumbra is examined through numerical experiments, and different geometric models of the penumbral magnetic field are discussed in the light of high-resolution observations. Recent, realistic numerical MHD simulations of an entire sunspot have succeeded in reproducing the salient features of the convective pattern in the umbra and the inner penumbra. The siphon-flow mechanism still provides the best explanation of the Evershed flow, particularly in the outer penumbra where it often consists of cool, supersonic downflows.

We use oblate axisymmetric dynamical models including dark halos to determine the orbital structure of intermediate mass to massive Coma early-type galaxies. We find a large variety of orbital compositions. Averaged over all sample galaxies the unordered stellar kinetic energy in the azimuthal and the radial direction are of the same order, but they can differ by up to 40 percent in individual systems. In contrast, both for rotating and non-rotating galaxies the vertical kinetic energy is on average smaller than in the other two directions. This implies that even most of the rotating ellipticals are flattened by an anisotropy in the stellar velocity dispersions. Using three-integral axisymmetric toy models we show that flattening by stellar anisotropy maximises the entropy for a given density distribution. Collisionless disk merger remnants are radially anisotropic. The apparent lack of strong radial anisotropy in observed early-type galaxies implies that they may not have formed from mergers of disks unless the influence of dissipational processes was significant.

Axisymmetric, orbit-based dynamical models are used to derive dark matter scaling relations for Coma early-type galaxies. From faint to bright galaxies halo core-radii and asymptotic circular velocities increase. Compared to spirals of the same brightness, the majority of Coma early-types -- those with old stellar populations -- have similar halo core-radii but more than 2 times larger asymptotic halo velocities. The average dark matter density inside 2 reff decreases with increasing luminosity and is 6.8 times larger than in disk galaxies of the same B-band luminosity. Compared at the same stellar mass, dark matter densities in ellipticals are 13.5 times higher than in spirals. Different baryon concentrations in ellipticals and spirals cannot explain the higher dark matter density in ellipticals. Instead, the assembly redshift (1+z) of Coma early-type halos is likely about two times larger than of comparably bright spirals. Assuming that local spirals typically assemble at a redshift of one, the majority of bright Coma early-type galaxy halos must have formed around z = 2-3. For about half of our Coma galaxies the assembly redshifts match with constraints derived from stellar populations. We find dark matter densities and estimated assembly redshifts of our observed Coma galaxies in reasonable agreement with recent semi-analytic galaxy formation models.

The radio galaxy Fornax A (NGC 1316) is a prominent merger remnant in the outskirts of the Fornax cluster. Its giant radio lobes suggest the presence of a powerful AGN and thus a central supermassive black hole (SMBH). We present high-resolution adaptive optics assisted integral-field data of Fornax A, taken with SINFONI at the Very Large Telescope in the K band. We use axisymmetric orbit models to determine the mass of the SMBH in the centre of Fornax A. The three-dimensional nature of our data provides the possibility to directly test the consistency of the data with axisymmetry by modelling each of the four quadrants separately. According to our dynamical models, consistent SMBH masses and dynamical Ks band mass-to-light ratios are obtained for all quadrants, with <M_BH>=1.3x10^8 M_\odot (rms(M_BH)=0.4x10^8 Msun) and <M/L>=0.68 (rms(M/L)=0.03), confirming the assumption of axisymmetry. For the folded and averaged data we find M_BH=(1.5+0.75-0.8)x10^8 Msun and M/L=(0.65+0.075-0.05) (3-sigma errors). Thus the black-hole mass of Fornax A is consistent within the error with the Tremaine (2002) M-sigma relation, but is a factor ~4 smaller than expected from its bulge mass and the Marconi&Hunt (2003) relation.

We have searched for correlations between the pointing directions of ultrahigh energy cosmic rays observed by the High Resolution Fly's Eye experiment and Active Galactic Nuclei (AGN) visible from its northern hemisphere location. No correlations, other than random correlations, have been found. We report our results using search parameters prescribed by the Pierre Auger collaboration. Using these parameters, the Auger collaboration concludes that a positive correlation exists for sources visible to their southern hemisphere location. We also describe results using two methods for determining the chance probability of correlations: one in which a hypothesis is formed from scanning one half of the data and tested on the second half, and another which involves a scan over the entire data set. The most significant correlation found occurred with a chance probability of 24%.

Air-fluorescence detectors such as the High Resolution Fly's Eye (HiRes) detector are very sensitive to upward-going, Earth-skimming ultrahigh energy electron-neutrino-induced showers. This is due to the relatively large interaction cross sections of these high-energy neutrinos and to the Landau-Pomeranchuk-Migdal (LPM) effect. The LPM effect causes a significant decrease in the cross sections for bremsstrahlung and pair production, allowing charged-current electron-neutrino-induced showers occurring deep in the Earth's crust to be detectable as they exit the Earth into the atmosphere. A search for upward-going neutrino-induced showers in the HiRes-II monocular dataset has yielded a null result. From an LPM calculation of the energy spectrum of charged particles as a function of primary energy and depth for electron-induced showers in rock, we calculate the shape of the resulting profile of these showers in air. We describe a full detector Monte Carlo simulation to determine the detector response to upward-going electron-neutrino-induced cascades and present an upper limit on the flux of electron-neutrinos.

Extended Red Emission (ERE) was recently attributed to the photo-luminescence of either doubly ionized Polycyclic Aromatic Hydrocarbons (PAH$^{++}$), or charged PAH dimers. We analysed the visible and mid-infrared (mid-IR) dust emission in the North-West and South photo-dissociation regions of the reflection nebula NGC 7023.Using a blind signal separation method, we extracted the map of ERE from images obtained with the Hubble Space Telescope, and at the Canada France Hawaii Telescope. We compared the extracted ERE image to the distribution maps of the mid-IR emission of Very Small Grains (VSGs), neutral and ionized PAHs (PAH$^0$ and PAH$^+$) obtained with the Spitzer Space Telescope and the Infrared Space Observatory. ERE is dominant in transition regions where VSGs are being photo-evaporated to form free PAH molecules, and is not observed in regions dominated by PAH$^+$. Its carrier makes a minor contribution to the mid-IR emission spectrum. These results suggest that the ERE carrier is a transition species formed during the destruction of VSGs. Singly ionized PAH dimers appear as good candidates but PAH$^{++}$ molecules seem to be excluded.

The long-slit spectra obtained along the minor axis, offset major axis and diagonal axis are presented for 12 E and S0 galaxies of the Coma cluster drawn from a magnitude-limited sample studied before. The rotation curves, velocity dispersion profiles and the H_3 and H_4 coefficients of the Hermite decomposition of the line of sight velocity distribution are derived. The radial profiles of the Hbeta, Mg, and Fe line strength indices are measured too. In addition, the surface photometry of the central regions of a subsample of 4 galaxies recently obtained with Hubble Space Telescope is presented. The data will be used to construct dynamical models of the galaxies and study their stellar populations.

Dynamical models for 17 Coma early-type galaxies are presented. The galaxy sample consists of flattened, rotating as well as non-rotating early-types including cD and S0 galaxies with luminosities between M=-18.79 and M=-22.56. Kinematical long-slit observations cover at least the major and minor axis and extend to 1-4 effective radii. Axisymmetric Schwarzschild models are used to derive stellar mass-to-light ratios and dark halo parameters. In every galaxy models with a dark matter halo match the data better than models without. The statistical significance is over 95 percent for 8 galaxies, around 90 percent for 5 galaxies and for four galaxies it is not significant. For the highly significant cases systematic deviations between observed and modelled kinematics are clearly seen; for the remaining galaxies differences are more statistical in nature. Best-fit models contain 10-50 percent dark matter inside the half-light radius. The central dark matter density is at least one order of magnitude lower than the luminous mass density. The central phase-space density of dark matter is often orders of magnitude lower than in the luminous component, especially when the halo core radius is large. The orbital system of the stars along the major-axis is slightly dominated by radial motions. Some galaxies show tangential anisotropy along the minor-axis, which is correlated with the minor-axis Gauss-Hermite coefficient H4. Changing the balance between data-fit and regularisation constraints does not change the reconstructed mass structure significantly. Model anisotropies tend to strengthen if the weight on regularisation is reduced, but the general property of a galaxy to be radially or tangentially anisotropic, respectively, does not change. (abridged)

Measurements are reported of the yield and spectrum of fluorescence, excited by a 28.5 GeV electron beam, in air at a range of pressures of interest to ultra-high energy cosmic ray detectors. The wavelength range was 300 - 420 nm. System calibration has been performed using Rayleigh scattering of a nitrogen laser beam. In atmospheric pressure dry air at 304 K the yield is 20.8 +/- 1.6 photons per MeV.

We model mock observations of collisionless N-body disc-disc mergers with an axisymmetric orbit superposition program that has been used to model Coma ellipticals. The remnants sample representatively the shapes of disc-disc mergers including prolate, triaxial and oblate objects. The aim is to better understand how the assumption of axial symmetry affects reconstructed masses and stellar motions of systems which are intrinsically not axisymmetric, whether it leads to a bias and how such a potential bias can be recognised in models of real galaxies. The mass recovery at the half-light radius depends on viewing-angle and intrinsic shape: edge-on views allow to reconstruct total masses with an accuracy between 20% (triaxial/prolate remnants) and 3% (oblate remnant). Masses of highly flattened, face-on systems are underestimated by up to 50%. Deviations in local mass densities can be larger where remnants are strongly triaxial or prolate. Luminous M/L are sensitive to box orbits in the remnants. Box orbits cause the central kinematics to vary with viewing-angle. Reconstructed luminous M/L and central masses follow this variation. Luminous M/L are always underestimated (up to a factor of 2.5). Respective dark halos in the models can be overestimated by about the same amount, depending again on viewing angle. Reconstructed velocity anisotropies depend on viewing angle and on the orbital composition of the remnant. We construct N-body realisations of the Schwarzschild models to discuss chaotic orbits and the virial equilibrium in our models. Apparently flattened, rotating ellipticals of intermediate mass are likely close to both, axial symmetry and edge-on orientation. Our results imply that Schwarzschild models allow a reconstruction of their masses and stellar anisotropies with high accuracy. (abridged)

The near-infrared integral field spectrograph SINFONI at the ESO VLT opens a new window for the study of central supermassive black holes. With a near-IR spatial resolution similar to HST optical and the ability to penetrate dust it provides the possibility to explore the low-mass end of the M-sigma relation (sigma<120km/s) where so far very few black hole masses were measured with stellar dynamics. With SINFONI we observed the central region of the low-luminosity elliptical galaxy NGC4486a at a spatial resolution of ~0.1arcsec in the K band. The stellar kinematics was measured with a maximum penalised likelihood method considering the region around the CO absorption band heads. We determined a black hole mass of M_BH=1.25^{+0.75}_{-0.79} x 10^7 M_sun (90% C.L.) using the Schwarzschild orbit superposition method including the full 2-dimensional spatial information. This mass agrees with the predictions of the M-sigma relation, strengthening its validity at the lower sigma end.

This work was conducted as part of the SPECPDR program, dedicated to the study of very small particles and astrochemistry, in Photo-Dissociation Regions (PDRs). We present the analysis of the mid-IR spectro-imagery observations of Ced 201, NCG 7023 East and North-West and rho-Ophiuchi West filament. Using the data from all four modules of the InfraRed Spectrograph onboard the Spitzer Space Telescope, we produced a spectral cube ranging from 5 to 35 um for each one of the observed PDRs. The resulting cubes were analysed using Blind Signal Separation methods (NMF and FastICA). For Ced 201, rho-Ophiuchi West filament and NGC 7023 East, we find that two signals can be extracted from the original data cubes, which are 5 to 35 um spectra. The main features of the first spectrum are a strong continuum emission at long wavelengths, and a broad 7.8 um band. On the contrary, the second spectrum exhibits the classical Aromatic Infrared Bands (AIBs) and no continuum. The reconstructed spatial distribution maps show that the latter spectrum is mainly present at the cloud surface, close to the star whereas the first one is located slightly deeper inside the PDR. The study of the spectral energy distribution of Ced 201 up to 100 um suggests that, in cool PDRs, the 5-25 um continuum is carried by Very Small Grains (VSGs). The AIB spectra in the observed objects can be interpreted as the contribution of neutral and positively-charged Polycyclic Aromatic Hydrocarbons (PAHs). We extracted the 5 to 25 um emission spectrum of VSGs in cool PDRs, these grains being most likely carbonaceous. We show that the variations of the mid-IR (5-35 um) spectra of PDRs can be explained by the photo-chemical processing of these VSGs and PAHs, VSGs being the progenitors of free PAH.

We report R~4300 VLT SINFONI adaptive optics integral field K-band spectroscopy of the nucleus of the Seyfert 1 galaxy NGC3227 at a spatial resolution of 0.085" (7pc). We present the morphologies and kinematics of emission lines and absorption features, and give the first derivation of a black hole mass in a Seyfert 1 nucleus from spatially resolved stellar dynamics. We show that the gas in the nucleus has a mean column density of order 10^{24}-10^{25}cm^{-2} and that it is geometrically thick, in agreement with the standard `molecular torus' scenario. We discuss which heating processes may be responsible for maintaining the vertical height of the torus. We have also resolved the nuclear stellar distribution, and find that within a few parsecs of the AGN there has been an intense starburst, the most recent episode of which began ~40Myr ago but has now ceased. The current luminosity of stars within 30pc of the AGN, ~3x10^9L_sun, is comparable to that of the AGN. Based on a comparison of the respective size scales, we argue that the star formation has been occuring in the obscuring torus. Finally, we present the first derivation of a black hole mass in a Seyfert 1 nucleus from stellar dynamics which marginally spatially resolve the black hole's sphere of influence. We apply Schwarzschild orbit superposition models to our full 2-dimensional data and derive the mass of the black hole, paying careful attention to the input parameters which are often uncertain: the contribution of the large scale bulge and its mass-to-light ratio; the recent star formation in the nucleus and its mass-to-light ratio; the contribution of the gas mass to the potential; and the inclination. Our models yield a 1sigma range for the black hole mass of M_{BH} = 7x10^6-2x10^7M_sun.

We present some preliminary (half-way) results on our adaptive optics spectroscopic survey of AGN at spatial scales down to 0.085arcsec. Most of the data were obtained with SINFONI which provides integral field capability at a spectral resolution of R~4000. The themes on which we focus in this contribution are: star formation around the AGN, the properties of the molecular gas and its relation to the torus, and the mass of the black hole.

Measurements are reported on the fluorescence of air as a function of depth in electromagnetic showers initiated by bunches of 28.5 GeV electrons. The light yield is compared with the expected and observed depth profiles of ionization in the showers. It validates the use of atmospheric fluorescence profiles in measuring ultra high energy cosmic rays.

We present HST spectroscopy of the nucleus of M31 obtained with STIS. Spectra taken around the CaT lines at 8500 see only the red giants in the double bright- ness peaks P1 and P2. In contrast, spectra taken at 3600-5100 A are sensitive to the tiny blue nucleus embedded in P2, the lower surface brightness red nucleus. P2 has a K-type spectrum, but the embedded blue nucleus has an A-type spectrum with strong Balmer absorption lines. Given the small likelihood for stellar collisions, a 200 Myr old starburst appears to be the most plausible origin of the blue nucleus. In stellar population, size, and velocity dispersion, the blue nucleus is so different from P1 and P2 that we call it P3. The line-of-sight velocity distributions of the red stars in P1+P2 strengthen the support for Tremaine s eccentric disk model. The kinematics of P3 is consistent with a circular stellar disk in Keplerian rotation around a super-massive black hole with M_bh = 1.4 x 10^8 M_sun. The P3 and the P1+P2 disks rotate in the same sense and are almost coplanar. The observed velocity dispersion of P3 is due to blurred rotation and has a maximum value of sigma = 1183+-201 km/s. The observed peak rotation velocity of P3 is V = 618+-81 km/s at radius 0.05" = 0.19 pc corresponding to a circular rotation velocity at this radius of ~1700 km/s. Any dark star cluster alternative to a black hole must have a half-mass radius <= 0.03" = 0.11 pc. We show that this excludes clusters of brown dwarfs or dead stars on astrophysical grounds.

Observational evidence for anisotropy in supernovae (SN) may signal the importance of angular momentum and differential rotation in the progenitors. Free energy in differential rotation and rotation can be extracted magnetically or via turbulent dissipation. The importance that magnetohydrodyamic jets and coronae may play in driving SN motivates understanding large scale dynamos in SN progenitors. We develop a dynamical large scale interface dynamo model in which the differential rotation and rotation deplete both through Poynting flux and turbulent diffusion. We apply the model to a differentially rotating core surrounded by a convection zone of a SN progenitor from a initial 15$M_\odot$ star. Unlike the Sun, the dynamo is transient because the differential rotation is primarily due to the initial collapse. Up to $\sim 10^{51}$erg can be drained into time-integrated Poynting flux and heat, the relative fraction of which depends on the relative amount of turbulence in the shear layer vs. convection zone and the fraction of the shear layer into which the magnetic field penetrates. Both sinks can help facilitate explosions and could lead to different levels of anisotropy and pulsar kicks. In all cases, the poloidal magnetic field is much weaker than the toroidal field, and the Poynting flux is lower than previous estimates which invoke the magnitude of the total magnetic energy. A signature of a large scale dynamo is that the oscillation of the associated Poynting flux on $\sim 1$ sec time scales, implying the same for the energy delivery to a SN.

This is the second in a series of papers dedicated to unveil the mass structure and orbital content of a sample of flattened early-type galaxies in the Coma cluster. The ability of our orbit libraries to reconstruct internal stellar motions and the mass composition of a typical elliptical in the sample is investigated by means of Monte-Carlo simulations of isotropic rotator models. The simulations allow a determination of the optimal amount of regularization needed in the orbit superpositions. It is shown that under realistic observational conditions and with the appropriate regularization internal velocity moments can be reconstructed to an accuracy of about 15 per cent; the same accuracy can be achieved for the circular velocity and dark matter fraction. In contrast, the flattening of the halo remains unconstrained. Regularized orbit superpositions are applied to a first galaxy in our sample, NGC 4807, for which stellar kinematical observations extend to 3 Reff. The galaxy seems dark matter dominated outside 2 Reff. Logarithmic dark matter potentials are consistent with the data, as well as NFW-profiles, mimicking logarithmic potentials over the observationally sampled radial range. In both cases, the derived stellar mass-to-light ratio agrees well with independently obtained mass-to-light ratios from stellar population analysis. Kinematically, NGC 4807 is characterized by mild radial anisotropy outside r>0.5 Reff, becoming isotropic towards the center. Our orbit models hint at either a distinct stellar component or weak triaxiality in the outer parts of the galaxy.

This is the first of a series of papers dedicated to unveil the mass composition and dynamical structure of a sample of flattened early type galaxies in the Coma cluster. We describe our modifications to the Schwarzschild code of Richstone et al. (in preparation). Applying a Voronoi tessellation in the surface of section we are able to assign accurate phase-space volumes to individual orbits and to reconstruct the full three-integral phase-space distribution function (DF) of any axisymmetric orbit library. Two types of tests have been performed to check the accuracy with which DFs can be represented by appropriate orbit libraries. First, by mapping DFs of spherical gamma-models and flattened Plummer models on to the library we show that the resulting line-of-sight velocity distributions and internal velocity moments of the library match those directly derived from the DF to a precision better than that of present day observational errors. Second, by fitting libraries to the projected kinematics of the same DFs we show that the distribution function reconstructed from the fitted library matches the input DF to a rms of about 15 per cent over a region in phase-space covering 90 per cent of the mass of the library. The achieved accuracy allows us to implement effective entropy-based regularisation to fit real, noisy and spatially incomplete data.

We describe a new mechanism for pulsations in evolved stars: relaxation oscillations driven by a coupling between the luminosity-dependent mass-loss rate and the H fuel abundance in a nuclear-burning shell. When mass loss is included, the outward flow of matter can modulate the flow of fuel into the shell when the stellar luminosity is close to the Eddington luminosity $L_{\rm Edd}$. When the luminosity drops below $L_{\rm Edd}$, the mass outflow declines and the shell is re-supplied with fuel. This process can be repetitive. We demonstrate the existence of such oscillations and discuss the dependence of the results on the stellar parameters. In particular, we show that the oscillation period scales specifically with the mass of the H-burning relaxation shell (HBRS), defined as the part of the H-burning shell above the minimum radius at which the luminosity from below first exceeds the Eddington threshold at the onset of the mass loss phase. For a stellar mass $M_*\sim 0.7\Msun$, luminosity $L_*\sim 10^4\Lsun$, and mass loss rate $|\dot M|\sim 10^{-5}\Msun$ yr$^{-1}$, the oscillations have a recurrence time $\sim 1400$ years $\sim 57\tau_{\rm fsm}$, where $\tau_{\rm fsm}$ is the timescale for modulation of the fuel supply in the HBRS by the varying mass-loss rate. This period agrees with the $\sim$ 1400-year period inferred for the spacings between the shells surrounding some planetary nebulae, and the the predictied shell thickness, of order 0.4 times the spacing, also agrees reasonably well.

Recently Perlmutter et al. suggested a positive value of Einstein's cosmological constant Lambda on the basis of luminosity distances from type-Ia supernovae. However, Lambda world models had earlier been proposed by Hoell & Priester and Liebscher et al. on the basis of quasar absorption-line data. Employing more general repulsive fluids ("dark energy") encompassing the Lambda component we quantitatively compare both approaches with each other. Fitting the SN-data by a minimum-component model consisting of dark energy + dust yields a closed universe with a large amount of dust exceeding the baryonic content constrained by big-bang nucleosynthesis. The nature of the dark energy is hardly constrained. Only when enforcing a flat universe there is a clear tendency to a dark-energy Lambda fluid and the `canonical' value Omega_M = 0.3 for dust. Conversely, fitting the quasar-data by a minimum-component model yields a sharply defined, slightly closed model with a low dust density ruling out significant pressureless dark matter. The dark-energy component obtains an equation-of-state P = -0.96 epsilon close to that of a Lambda-fluid. Omega_M = 0.3 or a precisely flat spatial geometry are inconsistent with minimum-component models. It is found that quasar and supernova data sets cannot be reconciled with each other via (repulsive ideal fluid+dust+radiation)-world models. Compatibility could be reached by drastic expansion of the parameter space with at least two exotic fluids added to dust and radiation as world constituents. If considering such solutions as far-fetched one has to conclude that the quasar absorption line and the SN-Ia constraints are incompatible.

Planetary nebulae are thought to be formed when a slow wind from the progenitor giant star is overtaken by a subsequent fast wind generated as the star enters its white dwarf stage$^{1}$. A shock forms near the boundary between the winds, which creates a relatively dense shell that provides the characteristic appearance of a planetary nebula. A spherically symmetric wind will produce a spherically symmetric shell, yet over half of known planetary nebulae are not spherical; rather, they are elliptical or bipolar in shape$^{2}$. While a magnetic field could launch and collimate a bipolar outflow, the origin of such a field has hitherto been unclear, as previous work suggested that a field could not be generated${^3}$. Here we show that an asymptotic-giant-branch (AGB) star can indeed generate a strong magnetic field, in a dynamo at the interface between a rapidly rotating core and the more slowly rotating envelope of the star. The field is strong enough to shape the bipolar outflows that produce the observed bipolar planetary nebulae. Magnetic braking of the stellar core during this process may also explain the puzzlingly$^{4}$ slow rotation of most white dwarf stars.

Various classification schemes exist for homogeneous and isotropic (CP) world models, which include pressureless matter (so-called dust) and Einstein's cosmological constant Lambda. We here classify the solutions of more general world models consisting of up to four non-interacting fluids, each with pressure P, energy density epsilon and an equation of state P = (gamma - 1) epsilon with 0 <= gamma <= 2. In addition to repulsive fluids with negative pressure and positive energy density, which generalize the classical repulsive (positive) Lambda component, we consider fluids with negative energy density as well. The latter generalize a negative Lambda component. This renders possible new types of models that do not occur among the classical classifications of world models. Singularity-free periodic solutions as well as further `hill-type', `hollow-type' and `shifting-type' models are feasible. However, if one only allows for three components (dust, radiation and one repulsive component) in a spatially flat universe the repulsive classical Lambda fluid (with Lambda > 0) tends to yield the smoothest fits of the Supernova Ia data from Perlmutter et al. (1999). Adopting the SN Ia constraints, exotic negative energy density components can be fittingly included only if the universe consists of four or more fluids.

Extensive results from the commissioning phase of PUEO, the adaptive optics instrument adaptor for the Canada-France-Hawaii telescope, are presented and discussed. Analyses of more than 750 images recorded with a CCD and a near-IR camera on 16 nights in wavelengths from B to H are used to derive the properties of the compensated wavefront and images in a variety of conditions. The performance characteristics of the system are analyzed and presented in several ways, in terms of delivered Strehl ratios, full-width-half-maxima, and quantities describing the improvements of both. A qualitative description is given of how the properties of the corrected images result from the structure function of the compensated phase. Under median seeing conditions, PUEO delivers essentially diffraction-limited images at H and K, images with FWHM ~0.1" at J and I, and provides significant gains down to B, with guide stars as faint as R = 14. During good conditions, substantial gains were realized with guide stars as faint as R = 17. A simple user-interface and software which automatically and continuously optimizes the mode gains during observations makes the operational efficiency extremely high. A few astronomical examples are briefly discussed.