Counterrotating Nuclear Disks in Arp 220

The ultraluminous infrared galaxy Arp 220 has been observed at 0.5" resolution in CO(2-1) and 1 mm continuum using the newly expanded Owens Valley Millimeter Array. The CO and continuum peaks at the double nuclei and the surrounding molecular gas disk are clearly resolved. We find steep velocity gradients across each nucleus (dV ~ 500 km/s within r= 0.3") whose directions are not aligned with each other and with that of the outer gas disk. We conclude that the double nuclei have their own gas disks (r ~ 100 pc). They are counterrotating with respect to each other and embedded in the outer gas disk (r ~ 1 kpc) rotating around the dynamical center of the system. The masses of each nucleus are M_dyn > 2* 10^9 M_sun based on the CO kinematics. Although there is no evidence of an old stellar population in the optical or near infrared spectroscopy of the nuclei (probably due to the much brighter young population), it seems likely that these nuclei were 'seeded' from the pre-merger nuclei in view of their counterrotating gas kinematics. The gas disks probably constitute a significant fraction (~ 50 %) of the mass in each nucleus. The CO and continuum brightness temperatures imply that the nuclear gas disks have high area filling factors (~ 0.5-1) and have extremely high visual extinctions (Av ~ 1000 mag). The molecular gas must be hot (>= 40 K) and dense (>= 10^4-5 cm^-3), given the large mass and small scale-height of the nuclear disks. The continuum data suggest that the large luminosity (be it starburst or AGN) must originate within 100 pc of the two nuclear gas disks which were presumably formed through concentration of gas from the progenitor outer galaxy disks.Comment: 20 pages, 5 figures. Accepted for publication in The Astrophysical Journa

Sharing and Preserving Computational Analyses for Posterity with encapsulator

Open data and open-source software may be part of the solution to science's "reproducibility crisis", but they are insufficient to guarantee reproducibility. Requiring minimal end-user expertise, encapsulator creates a "time capsule" with reproducible code in a self-contained computational environment. encapsulator provides end-users with a fully-featured desktop environment for reproducible research.Comment: 11 pages, 6 figure

Modelling CO emission from Mira's wind

We have modelled the circumstellar envelope of {\it o} Ceti (Mira) using new observational constraints. These are obtained from photospheric light scattered in near-IR vibrational-rotational lines of circumstellar CO molecules at 4.6 micron: absolute fluxes, the radial dependence of the scattered intensity, and two line ratios. Further observational constraints are provided by ISO observations of far-IR emission lines from highly excited rotational states of the ground vibrational state of CO, and radio observations of lines from rotational levels of low excitation of CO. A code based on the Monte-Carlo technique is used to model the circumstellar line emission. We find that it is possible to model the radio and ISO fluxes, as well as the highly asymmetric radio-line profiles, reasonably well with a spherically symmetric and smooth stellar wind model. However, it is not possible to reproduce the observed NIR line fluxes consistently with a `standard model' of the stellar wind. This is probably due to incorrectly specified conditions of the inner regions of the wind model, since the stellar flux needs to be larger than what is obtained from the standard model at the point of scattering, i.e., the intermediate regions at approximately 100-400 stellar radii (2"-7") away from the star. Thus, the optical depth in the vibrational-rotational lines from the star to the point of scattering has to be decreased. This can be accomplished in several ways. For instance, the gas close to the star (within approximately 2") could be in such a form that light is able to pass through, either due to the medium being clumpy or by the matter being in radial structures (which, further out, developes into more smooth or shell-like structures).Comment: 18 pages, 3 figures, accepted for publication in Ap

Mira's wind explored in scattering infrared CO lines

We have observed the intermediate regions of the circumstellar envelope of Mira (o Ceti) in photospheric light scattered by three vibration-rotation transitions of the fundamental band of CO, from low-excited rotational levels of the ground vibrational state, at an angular distance of beta = 2"-7" away from the star. The data were obtained with the Phoenix spectrometer mounted on the 4 m Mayall telescope at Kitt Peak. The spatial resolution is approximately 0.5" and seeing limited. Our observations provide absolute fluxes, leading to an independent new estimate of the mass-loss rate of approximately 3e-7 Msun/yr, as derived from a simple analytic wind model. We find that the scattered intensity from the wind of Mira for 2" < beta < 7" decreases as beta^-3, which suggests a time constant mass-loss rate, when averaged over 100 years, over the past 1200 years.Comment: accepted for publication in the Astrophysical Journa

A C2HC zinc finger is essential for the RING-E2 interaction of the ubiquitin ligase RNF125

The activity of RING ubiquitin ligases (E3s) depends on an interaction between the RING domain and ubiquitin conjugating enzymes (E2), but posttranslational events or additional structural elements, yet largely undefined, are frequently required to enhance or regulate activity. Here, we show for the ubiquitin ligase RNF125 that, in addition to the RING domain, a C2HC Zn finger (ZnF) is crucial for activity, and a short linker sequence (Li2(120-128)) enhances activity. The contribution of these regions was first shown with truncated proteins, and the essential role of the ZnF was confirmed with mutations at the Zn chelating Cys residues. Using NMR, we established that the C2HC ZnF/Li2(120-128) region is crucial for binding of the RING domain to the E2 UbcH5a. The partial X-ray structure of RNF125 revealed the presence of extensive intramolecular interactions between the RING and C2HC ZnF. A mutation at one of the contact residues in the C2HC ZnF, a highly conserved M112, resulted in the loss of ubiquitin ligase activity. Thus, we identified the structural basis for an essential role of the C2HC ZnF and conclude that this domain stabilizes the RING domain, and is therefore required for binding of RNF125 to an E2

The Complexity of Computing the Optimal Composition of Differential Privacy

In the study of differential privacy, composition theorems (starting with the original paper of Dwork, McSherry, Nissim, and Smith (TCC’06)) bound the degradation of privacy when composing several differentially private algorithms. Kairouz, Oh, and Viswanath (ICML’15) showed how to compute the optimal bound for composing k arbitrary ( , δ)- differentially private algorithms. We characterize the optimal composition for the more general case of k arbitrary ( 1, δ1), . . . ,( k, δk)-differentially private algorithms where the privacy parameters may differ for each algorithm in the composition. We show that computing the optimal composition in general is #P-complete. Since computing optimal composition exactly is infeasible (unless FP=#P), we give an approximation algorithm that computes the composition to arbitrary accuracy in polynomial time. The algorithm is a modification of Dyer’s dynamic programming approach to approximately counting solutions to knapsack problems (STOC’03).Engineering and Applied Science