We present high-resolution (30 mas or 130 au at 4.2 kpc) Atacama Large Millimeter/submillimeter Array observations at 1.2 mm of the disc around the forming O-type star AFGL 4176 mm1. The disc (AFGL 4176 mm1-main) has a radius of ~1000 au and contains significant structure, most notably a spiral arm on its redshifted side. We fitted the observed spiral with logarithmic and Archimedean spiral models. We find that both models can describe its structure, but the Archimedean spiral with a varying pitch angle fits its morphology marginally better. As well as signatures of rotation across the disc, we observe gas arcs in CH$_3$CN that connect to other millimetre continuum sources in the field, supporting the picture of interactions within a small cluster around AFGL 4176 mm1-main. Using local thermodynamic equilibrium modelling of the CH$_3$CN K-ladder, we determine the temperature and velocity field across the disc, and thus produce a map of the Toomre stability parameter. Our results indicate that the outer disc is gravitationally unstable and has already fragmented or is likely to fragment in the future, possibly producing further companions. These observations provide evidence that disc fragmentation is one possible pathway towards explaining the high fraction of multiple systems around high-mass stars.

In the disk-mediated accretion scenario for the formation of the most massive stars, gravitational instabilities in the disk can force it to fragment. We investigate the effects of inclination and spatial resolution on observable kinematics and stability of disks in high-mass star formation. We study a high-resolution 3D radiation-hydrodynamic simulation that leads to the fragmentation of a massive disk. Using RADMC-3D we produce 1.3 mm continuum and CH3CN line cubes at different inclinations. The model is set to different distances and synthetic observations are created for ALMA at ~80 mas resolution and NOEMA at ~0.3''. The synthetic ALMA observations resolve all fragments and their kinematics well. The synthetic NOEMA observations at 800 pc (~300 au resolution) are able to resolve the fragments, while at 2000 pc (~800 au resolution) only a single slightly elongated structure is observed. The position-velocity (PV) plots show the differential rotation of material best in the edge-on views. As the observations become less resolved, the inner high-velocity components of the disk become blended with the envelope and the PV plots resemble rigid-body-like rotation. Protostellar mass estimates from PV plots of poorly resolved observations are therefore overestimated. We fit the emission of CH3CN lines and produce maps of gas temperature with values in the range of 100-300 K. Studying the Toomre stability of the disks in the resolved observations, we find Q values below the critical value for stability against gravitational collapse at the positions of the fragments and the arms connecting the fragments. For the poorly resolved observations we find low Q values in the outskirts of the disk. Therefore we are able to predict that the disk is unstable and fragmenting even in poorly resolved observations. This conclusion is true regardless of knowledge about the inclination of the disk.

We present a detailed observational and modeling study of the hot core VLA 3 in the high-mass star-forming region AFGL 2591, which is a target region of the NOrthern Extended Millimeter Array (NOEMA) large program CORE. Using NOEMA observations at 1.37 mm with an angular resolution of ~0."42 (1 400 au at 3.33 kpc), we derived the physical and chemical structure of the source. We modeled the observed molecular abundances with the chemical evolution code MUSCLE (MUlti Stage ChemicaL codE). Results. With the kinetic temperature tracers CH3CN and H2CO we observe a temperature distribution with a power-law index of q = 0.41+-0.08. Using the visibilities of the continuum emission we derive a density structure with a power-law index of p = 1.7+-0.1. The hot core spectra reveal high molecular abundances and a rich diversity in complex molecules. The majority of the molecules have an asymmetric spatial distribution around the forming protostar(s), which indicates a complex physical structure on scales < 1 400 au. Using MUSCLE, we are able to explain the observed molecular abundance of 10 out of 14 modeled species at an estimated hot core chemical age of ~21 100 years. In contrast to the observational analysis, our chemical modeling predicts a lower density power-law index of p < 1.4. Reasons for this discrepancy are discussed. Conclusions. Combining high spatial resolution observations with detailed chemical modeling allows us to derive a concise picture of the physical and chemical structure of the famous AFGL 2591 hot core. The next steps are to conduct a similar analysis for the whole CORE sample, and then use this analysis to constrain the chemical diversity in high-mass star formation to a much greater depth.

The formation process of high-mass stars (>8M$_\odot$) is poorly constrained, particularly, the effects of clump fragmentation creating multiple systems and the mechanism of mass accretion onto the cores. We study the fragmentation of dense gas clumps, and trace the circumstellar rotation and outflows by analyzing observations of the high-mass (~500M$_\odot$) star-forming region IRAS 23033+5951. Using the Northern Extended Millimeter Array (NOEMA) in three configurations and the IRAM 30-m single-dish telescope at 220GHz, we probe the gas and dust emission at an angular resolution of ~0.45arcsec, corresponding to 1900au. In the mm continuum emission, we identify a protostellar cluster with at least four mm-sources, where three of them show a significantly higher peak intensity well above a signal-to-noise ratio of 100. Hierarchical fragmentation from large to small spatial scales is discussed. Two fragments are embedded in rotating structures and drive molecular outflows, traced by $^{13}$CO (2-1) emission. The velocity profiles across two of the cores are similar to Keplerian but are missing the highest velocity components close to the center of rotation, which is a common phenomena from observations like these, and other rotation scenarios are not excluded entirely. Position-velocity diagrams suggest protostellar masses of ~6 and 19M$_\sun$. Rotational temperatures from fitting CH$_3$CN ($12_K-11_K$) spectra are used for estimating the gas temperature and by that the disk stability against gravitational fragmentation, utilizing Toomre's $Q$ parameter. [We] identify only one candidate disk to be unstable against gravitational instability caused by axisymmetric perturbations. The dominant sources cover different evolutionary stages within the same maternal gas clump. The appearance of rotation and outflows of the cores are similar to those found in low-mass star-forming regions.

We present the highest angular resolution (20x15mas - 44x33au) Atacama Large Millimeter/sub-millimeter Array (ALMA) observations currently possible of the proto-O-star G17.64+0.16 in Band 6. The Cycle 5 observations with baselines out to 16km probes scales <50au and reveal the rotating disc around G17.64+0.16, a massive forming O-type star. The disc has a ring-like enhancement in the dust emission, especially visible as arc structures to the north and south. The Keplerian kinematics are most prominently seen in the vibrationally excited water line, H2O (Eu=3461.9K). The mass of the central source found by modelling the Keplerian rotation is consistent with 45+/-10Mo. The H30alpha (231.9GHz) radio-recombination line and the SiO (5-4) molecular line were detected at up to the 10 sigma$ level. The estimated disc mass is 0.6-2.6Mo under the optically thin assumption. Analysis of the Toomre Q parameter, in the optically thin regime, indicates that the disc stability is highly dependent on temperature. The disc currently appears stable for temperatures >150K, this does not preclude that the substructures formed earlier through disc fragmentation.

This study is part of the project ``CORE'', an IRAM/NOEMA large program consisting of observations of the millimeter continuum and molecular line emission towards 20 selected high-mass star forming regions. We focus on IRAS23385+6053, which is believed to be the least evolved source of the CORE sample. The observations were performed at ~1.4 mm and employed three configurations of NOEMA and additional single-dish maps, merged with the interferometric data to recover the extended emission. Our correlator setup covered a number of lines from well-known hot core tracers and a few outflow tracers. The angular (~0.45"$-$0.9") and spectral (0.5 km/s) resolutions were sufficient to resolve the clump in IRAS23385+6053 and investigate the existence of large-scale motions due to rotation, infall, or expansion. We find that the clump splits into six distinct cores when observed at sub-arcsecond resolution. These are identified through their 1.4 mm continuum and molecular line emission. We produce maps of the velocity, line width, and rotational temperature from the methanol and methyl cyanide lines, which allow us to investigate the cores and reveal a velocity and temperature gradient in the most massive core. We also find evidence of a bipolar outflow, possibly powered by a low-mass star. We present the tentative detection of a circumstellar self-gravitating disk lying in the most massive core and powering a large-scale outflow previously known in the literature. In our scenario, the star powering the flow is responsible for most of the luminosity of IRAS23385+6053 (~$3000~L_\odot$). The other cores, albeit with masses below the corresponding virial masses, appear to be accreting material from their molecular surroundings and are possibly collapsing or on the verge of collapse. We conclude that we are observing a sample of star-forming cores that is bound to turn into a cluster of massive stars.

In this paper, we present a novel method for computing the relative entropy as well as the expected relative entropy using an MCMC chain. The relative entropy from information theory can be used to quantify differences in posterior distributions of a pair of experiments. In cosmology, the relative entropy has been proposed as an interesting tool for model selection, experiment design, forecasting and measuring information gain from subsequent experiments. In contrast to Gaussian distributions, these quantities are not generally available analytically and one needs to use numerical methods to estimate them which are certainly computationally expensive. We propose a method and provide its python package to estimate the relative entropy as well as expected relative entropy from a posterior sample. We consider the linear Gaussian model to check the accuracy of our code. Our results indicate that the relative error is below $0.2\%$ for sample size larger than $10^5$ in the linear Gaussian model. In addition, we study the robustness of our code in estimating the expected relative entropy in this model.

We present high angular resolution 0.2 arcsec continuum and molecular emission line Atacama Large Millimeter/sub-millimeter Array (ALMA) observations of G17.64+0.16 in Band 6 (220GHz) taken as part of a campaign in search of circumstellar discs around (proto)-O-stars. At a resolution of 400au the main continuum core is essentially unresolved and isolated from other strong and compact emission peaks. At a resolution of 400au the main continuum core is essentially unresolved and isolated from other strong and compact emission peaks. We detect SiO (5-4) emission that is marginally resolved and elongated in a direction perpendicular to the large-scale outflow seen in the 13CO (2-1) line using the main ALMA array in conjunction with the Atacama Compact Array (ACA). Morphologically, the SiO appears to represent a disc-like structure. Using parametric models we show that the position-velocity profile of the SiO is consistent with the Keplerian rotation of a disc around an object between 10-30Mo in mass, only if there is also radial expansion from a separate structure. The radial motion component can be interpreted as a disc wind from the disc surface. Models with a central stellar object mass between 20 and 30Mo are the most consistent with the stellar luminosity (100000 Lo) and indicative of an O-type star. The H30a millimetre recombination line (231.9GHz) is also detected, but spatially unresolved, and is indicative of a very compact, hot, ionised region co-spatial with the dust continuum core. Accounting for all observables, we suggest that G17.64 is consistent with a O-type young stellar object in the final stages of protostellar assembly, driving a wind, but that has not yet developed into a compact HII region. The existance and detection of the disc in G17.64 is likely related to its isolated and possibly more evolved nature, traits which may underpin discs in similar sources.

The fragmentation mode of high-mass molecular clumps and the properties of the central rotating structures surrounding the most luminous objects have yet to be comprehensively characterised. Using the IRAM NOrthern Extended Millimeter Array (NOEMA) and the IRAM 30-m telescope, the CORE survey has obtained high-resolution observations of 20 well-known highly luminous star-forming regions in the 1.37 mm wavelength regime in both line and dust continuum emission. We present the spectral line setup of the CORE survey and a case study for W3(H2O). At ~0.35" (700 AU at 2 kpc) resolution, the W3(H2O) clump fragments into two cores (West and East), separated by ~2300 AU. Velocity shifts of a few km/s are observed in the dense-gas tracer, CH3CN, across both cores, consistent with rotation and perpendicular to the directions of two bipolar outflows, one emanating from each core. The kinematics of the rotating structure about W3(H2O) W shows signs of differential rotation of material, possibly in a disk-like object. The observed rotational signature around W3(H2O) E may be due to a disk-like object, an unresolved binary (or multiple) system, or a combination of both. We fit the emission of CH3CN (12-11) K = 4-6 and derive a gas temperature map with a median temperature of ~165 K across W3(H2O). We create a Toomre Q map to study the stability of the rotating structures against gravitational instability. The rotating structures appear to be Toomre unstable close to their outer boundaries, with a possibility of further fragmentation in the differentially-rotating core W3(H2O) W. Rapid cooling in the Toomre-unstable regions supports the fragmentation scenario. Combining millimeter dust continuum and spectral line data toward the famous high-mass star-forming region W3(H2O), we identify core fragmentation on large scales, and indications for possible disk fragmentation on smaller spatial scales.

As part of our effort to search for circumstellar disks around high-mass stellar objects, we observed the well-known core G31.41+0.31 with ALMA at 1.4 mm with an angular resolution of~0.22" (~1700 au). The dust continuum emission has been resolved into two cores namely Main and NE. The Main core, which has the stronger emission and is the more chemically rich, has a diameter of ~5300 au, and is associated with two free-free continuum sources. The Main core looks featureless and homogeneous in dust continuum emission and does not present any hint of fragmentation. Each transition of CH3CN and CH3OCHO, both ground and vibrationally excited, as well as those of CH3CN isotopologues, shows a clear velocity gradient along the NE-SW direction, with velocity linearly increasing with distance from the center, consistent with solid-body rotation. However, when comparing the velocity field of transitions with different upper level energies, the rotation velocity increases with increasing energy of the transition, which suggests that the rotation speeds up towards the center. Spectral lines towards the dust continuum peak show an inverse P-Cygni profile that supports the existence of infall in the core. The infall velocity increases with the energy of the transition suggesting that the infall is accelerating towards the center of the core, consistent with gravitational collapse. Despite the monolithic appearance of the Main core, the presence of red-shifted absorption, the existence of two embedded free-free sources at the center, and the rotational spin-up are consistent with an unstable core undergoing fragmentation with infall and differential rotation due to conservation of angular momentum. Therefore, the most likely explanation for the monolithic morphology is that the large opacity of the dust emission prevents the detection of any inhomogeneity in the core.

Aims: We aim to understand the fragmentation as well as the disk formation, outflow generation and chemical processes during high-mass star formation on spatial scales of individual cores. Methods: Using the IRAM Northern Extended Millimeter Array (NOEMA) in combination with the 30m telescope, we have observed in the IRAM large program CORE the 1.37mm continuum and spectral line emission at high angular resolution (~0.4'') for a sample of 20 well-known high-mass star-forming regions with distances below 5.5kpc and luminosities larger than 10^4Lsun. Results: We present the overall survey scope, the selected sample, the observational setup and the main goals of CORE. Scientifically, we concentrate on the mm continuum emission on scales on the order of 1000AU. We detect strong mm continuum emission from all regions, mostly due to the emission from cold dust. The fragmentation properties of the sample are diverse. We see extremes where some regions are dominated by a single high-mass core whereas others fragment into as many as 20 cores. A minimum-spanning-tree analysis finds fragmentation at scales on the order of the thermal Jeans length or smaller suggesting that turbulent fragmentation is less important than thermal gravitational fragmentation. The diversity of highly fragmented versus singular regions can be explained by varying initial density structures and/or different initial magnetic field strengths. Conclusions: The smallest observed separations between cores are found around the angular resolution limit which indicates that further fragmentation likely takes place on even smaller spatial scales. The CORE project with its numerous spectral line detections will address a diverse set of important physical and chemical questions in the field of high-mass star formation.

The nearby open cluster NGC 752 presents a rare opportunity to study stellar properties at ages >1 Gyr. However, constructing a membership catalog for it is challenging; most surveys have been limited to identifying its giants and dwarf members earlier than mid-K. We supplement past membership catalogs with candidates selected with updated photometric and proper-motion criteria, generating a list of 258 members, a >50% increase over previous catalogs. Using a Bayesian framework to fit MESA Isochrones & Stellar Tracks evolutionary models to literature photometry and the Tycho-Gaia Astrometric Solution data available for 59 cluster members, we infer the age of, and distance to, NGC 752: 1.34$\pm$0.06 Gyr and 438$_{-6}^{+8}$ pc. We also report the results of our optical monitoring of the cluster using the Palomar Transient Factory. We obtain rotation periods for 12 K and M cluster members, the first periods measured for such low-mass stars with a well-constrained age >1 Gyr. We compare these new periods to data from the younger clusters Praesepe and NGC 6811, and to a theoretical model for angular-momentum loss, to examine stellar spin down for low-mass stars over their first 1.3 Gyr. While on average NGC 752 stars are rotating more slowly than their younger counterparts, the difference is not significant. Finally, we use our spectroscopic observations to measure Halpha for cluster stars, finding that members earlier than $\approx$M2 are magnetically inactive, as expected at this age. Forthcoming Gaia data should solidify and extend the membership of NGC 752 to lower masses, thereby increasing its importance for studies of low-mass stars.

Stellar rotation periods measured in open clusters have proved to be extremely useful for studying stars' angular momentum content and rotationally driven magnetic activity, which are both age- and mass-dependent processes. While period measurements have been obtained for hundreds of solar-mass members of the Pleiades, period measurements exist for only a few low-mass ($<$0.5 M$_{\odot}$) members of this key laboratory for stellar evolution theory. To fill this gap, we report rotation periods for 132 low-mass Pleiades members (including nearly 100 with M $\leq$ 0.45 M$_{\odot}$), measured from photometric monitoring of the cluster conducted by the Palomar Transient Factory in late 2011 and early 2012. These periods extend the portrait of stellar rotation at 125 Myr to the lowest-mass stars and re-establish the Pleiades as a key benchmark for models of the transport and evolution of stellar angular momentum. Combining our new rotation periods with precise BVIJHK photometry reported by Stauffer et al. and Kamai et al., we investigate known anomalies in the photometric properties of K and M Pleiades members. We confirm the correlation detected by Kamai et al. between a star's rotation period and position relative to the main sequence in the cluster's color-magnitude diagram. We find that rapid rotators have redder (V-K) colors than slower rotators at the same V, indicating that rapid and slow rotators have different binary frequencies and/or photospheric properties. We find no difference in the photometric amplitudes of rapid and slow rotators, indicating that asymmetries in the longitudinal distribution of starspots do not scale grossly with rotation rate.

Located within 10 degrees of the North Pole, northern Ellesmere Island offers continuous darkness in the winter months. This capability can greatly enhance the detection efficiency of planetary transit surveys and other time domain astronomy programs. We deployed two wide-field cameras at 80 degrees North, near Eureka, Nunavut, for a 152-hour observing campaign in February 2012. The 16-megapixel-camera systems were based on commercial f/1.2 lenses with 70mm and 42mm apertures, and they continuously imaged 504 and 1,295 square degrees respectively. In total, the cameras took over 44,000 images and produced better-than-1% precision light curves for approximately 10,000 stars. We describe a new high-speed astrometric and photometric data reduction pipeline designed for the systems, test several methods for the precision flat-fielding of images from very-wide-angle cameras, and evaluate the cameras' image qualities. We achieved a scintillation-limited photometric precision of 1-2% in each 10s exposure. Binning the short exposures into 10-min chunks provided a photometric stability of 2-3 millimagnitudes, sufficient for the detection of transiting exoplanets around the bright stars targeted by our survey. We estimate that the cameras, when operated over the full arctic winter, will be capable of discovering several transiting exoplanets around bright (V<9) stars.

Observations from near the Eureka station on Ellesmere Island, in the Canadian High Arctic at 80 degrees North, benefit from 24-hour darkness combined with dark skies and long cloud-free periods during the winter. Our first astronomical surveys conducted at the site are aimed at transiting exoplanets; compared to mid-latitude sites, the continuous darkness during the Arctic winter greatly improves the survey's detection efficiency for longer-period transiting planets. We detail the design, construction, and testing of the first two instruments: a robotic telescope, and a set of very wide-field imaging cameras. The 0.5m Dunlap Institute Arctic Telescope has a 0.8-square-degree field of view and is designed to search for potentially habitable exoplanets around low-mass stars. The very wide field cameras have several-hundred-square-degree fields of view pointed at Polaris, are designed to search for transiting planets around bright stars, and were tested at the site in February 2012. Finally, we present a conceptual design for the Compound Arctic Telescope Survey (CATS), a multiplexed transient and transit search system which can produce a 10,000-square-degree snapshot image every few minutes throughout the Arctic winter.