Why haven't loose globular clusters collapsed yet?

We report on the discovery of a surprising observed correlation between the slope of the low-mass stellar global mass function (GMF) of globular clusters (GCs) and their central concentration parameter c=log(r_t/r_c), i.e. the logarithmic ratio of tidal and core radii. This result is based on the analysis of a sample of twenty Galactic GCs with solid GMF measurements from deep HST or VLT data. All the high-concentration clusters in the sample have a steep GMF, most likely reflecting their initial mass function. Conversely, low-concentration clusters tend to have a flatter GMF implying that they have lost many stars via evaporation or tidal stripping. No GCs are found with a flat GMF and high central concentration. This finding appears counter-intuitive, since the same two-body relaxation mechanism that causes stars to evaporate and the cluster to eventually dissolve should also lead to higher central density and possibly core-collapse. Therefore, more concentrated clusters should have lost proportionately more stars and have a shallower GMF than low concentration clusters, contrary to what is observed. It is possible that severely depleted GCs have also undergone core collapse and have already recovered a normal radial density profile. It is, however, more likely that GCs with a flat GMF have a much denser and smaller core than suggested by their surface brightness profile and may well be undergoing collapse at present. In either case, we may have so far seriously underestimated the number of post core-collapse clusters and many may be lurking in the Milky Way.Comment: Four pages, one figure, accepted for publication in ApJ Letter

Greedy kernel methods for accelerating implicit integrators for parametric ODEs

We present a novel acceleration method for the solution of parametric ODEs by single-step implicit solvers by means of greedy kernel-based surrogate models. In an offline phase, a set of trajectories is precomputed with a high-accuracy ODE solver for a selected set of parameter samples, and used to train a kernel model which predicts the next point in the trajectory as a function of the last one. This model is cheap to evaluate, and it is used in an online phase for new parameter samples to provide a good initialization point for the nonlinear solver of the implicit integrator. The accuracy of the surrogate reflects into a reduction of the number of iterations until convergence of the solver, thus providing an overall speedup of the full simulation. Interestingly, in addition to providing an acceleration, the accuracy of the solution is maintained, since the ODE solver is still used to guarantee the required precision. Although the method can be applied to a large variety of solvers and different ODEs, we will present in details its use with the Implicit Euler method for the solution of the Burgers equation, which results to be a meaningful test case to demonstrate the method's features

Star Formation in the Starburst Cluster in NGC 3603

We have used new, deep, visible and near infrared observations of the compact starburst cluster in the giant HII region NGC 3603 and its surroundings with the WFC3 on HST and HAWK-I on the VLT to study in detail the physical properties of its intermediate mass (~ 1 - 3 M_sun) stellar population. We show that after correction for differential extinction and actively accreting stars, and the study of field star contamination, strong evidence remains for a continuous spread in the ages of pre-main sequence stars in the range ~ 2 to ~ 30 Myr within the temporal resolution available. Existing differences among presently available theoretical models account for the largest possible variation in shape of the measured age histograms within these limits. We also find that this isochronal age spread in the near infrared and visible Colour-Magnitude Diagrams cannot be reproduced by any other presently known source of astrophysical or instrumental scatter that could mimic the luminosity spread seen in our observations except, possibly, episodic accretion. The measured age spread and the stellar spatial distribution in the cluster are consistent with the hypothesis that star formation started at least 20-30 Myrs ago progressing slowly but continuously up to at least a few million years ago. All the stars in the considered mass range are distributed in a flattened oblate spheroidal pattern with the major axis oriented in an approximate South-East - North-West direction, and with the length of the equatorial axis decreasing with increasing age. This asymmetry is most likely due to the fact that star formation occurred along a filament of gas and dust in the natal molecular cloud oriented locally in this direction.Comment: 21 pages, 19 figures, accepted for publication in Astrophysics & Space Scienc

WIYN Open Cluster Study 1: Deep Photometry of NGC 188

We have employed precise V and I photometry of NGC 188 at WIYN to explore the cluster luminosity function (LF) and study the cluster white dwarfs (WDs). Our photometry is offset by V = 0.052 (fainter) from Sandage (1962) and Eggen & Sandage (1969). All published photometry for the past three decades have been based on these two calibrations, which are in error by 0.05 +- 0.01. We employ the Pinsonneault etal (1998) fiducial main sequence to derive a cluster distance modulus of 11.43 +- 0.08. We report observations that are >= 50% complete to V = 24.6 and find that the cluster central-field LF peaks at M_I ~ 3 to 4. This is unlike the solar neighborhood LF and unlike the LFs of dynamically unevolved portions of open and globular clusters, which rise continuously until M_I ~ 9.5. Although we find that >= 50% of the unresolved cluster objects are multiple systems, their presence cannot account for the shape of the NGC 188 LF. For theoretical reasons (Terlevich 1987; Vesperini & Heggie 1997) having to do with the survivability of NGC 188 we believe the cluster is highly dynamically evolved and that the missing low luminosity stars are either in the cluster outskirts or have left the cluster altogether. We identify nine candidate WDs, of which we expect three to six are bona fide cluster WDs. The luminosities of the faintest likely WD indicates an age (Bergeron, Wesemael, & Beauchamp 1995) of 1.14 +- 0.09 Gyrs. This is a lower limit to the cluster age and observations probing to V = 27 or 28 will be necessary to find the faintest cluster WDs and independently determine the cluster age. While our age limit is not surprising for this ~6 Gyr old cluster, our result demonstrates the value of the WD age technique with its very low internal errors. (abridged)Comment: 26 pages, uuencoded gunzip'ed latex + 16 postscrip figures, to be published in A

The Evolution of Globular Clusters in the Galaxy

We investigate the evolution of globular clusters using N-body calculations and anisotropic Fokker-Planck (FP) calculations. The models include a mass spectrum, mass loss due to stellar evolution, and the tidal field of the parent galaxy. Recent N-body calculations have revealed a serious discrepancy between the results of N-body calculations and isotropic FP calculations. The main reason for the discrepancy is an oversimplified treatment of the tidal field employed in the isotropic FP models. In this paper we perform a series of calculations with anisotropic FP models with a better treatment of the tidal boundary and compare these with N-body calculations. The new tidal boundary condition in our FP model includes one free parameter. We find that a single value of this parameter gives satisfactory agreement between the N-body and FP models over a wide range of initial conditions. Using the improved FP model, we carry out an extensive survey of the evolution of globular clusters over a wide range of initial conditions varying the slope of the mass function, the central concentration, and the relaxation time. The evolution of clusters is followed up to the moment of core collapse or the disruption of the clusters in the tidal field of the parent galaxy. In general, our model clusters, calculated with the anisotropic FP model with the improved treatment for the tidal boundary, live longer than isotropic models. The difference in the lifetime between the isotropic and anisotropic models is particularly large when the effect of mass loss via stellar evolution is rather significant. On the other hand the difference is small for relaxation- dominated clusters which initially have steep mass functions and high central concentrations.Comment: 36 pages, 11 figures, LaTeX; added figures and tables; accepted by Ap

Function-Based Mapping of Industrial Assistance Systems to User Groups in Production

Abstract By looking at the last few decades, industrial production has undergone great changes. Industry 4.0, also called the fourth industrial revolution, describes the change in the entire value chain through the digital networking of systems, machines, and products. In addition, product variety and complexity in assembly increased due to customization. Big Data Analytics, Internet of Things, Horizontal and Vertical Data Integration and Cyber-Physical Production Systems are just some examples of technologies that find their way from research into industrial practice. However, the most important resource is often neglected, when talking about industry: the human. When we look at companies, we find different types of personnel in production, each with different requirements and capabilities. Assistance systems can be used to counteract these new challenges and offer adequate support to each individual worker. In the past, much research has been done to develop new worker assistance systems, while the analysis of specific needs of user groups in production has been ignored. This paper presents a function-based mapping of industrial worker assistance systems to different user groups and proposes a method for selecting the most appropriate assistance system to each user group

Small crater populations on Vesta

The NASA Dawn mission has extensively examined the surface of asteroid Vesta, the second most massive body in the main belt. The high quality of the gathered data provides us with an unique opportunity to determine the surface and internal properties of one of the most important and intriguing main belt asteroids (MBAs). In this paper, we focus on the size frequency distributions (SFDs) of sub-kilometer impact craters observed at high spatial resolution on several selected young terrains on Vesta. These small crater populations offer an excellent opportunity to determine the nature of their asteroidal precursors (namely MBAs) at sizes that are not directly observable from ground-based telescopes (i.e., below ~100 m diameter). Moreover, unlike many other MBA surfaces observed by spacecraft thus far, the young terrains examined had crater spatial densities that were far from empirical saturation. Overall, we find that the cumulative power-law index (slope) of small crater SFDs on Vesta is fairly consistent with predictions derived from current collisional and dynamical models down to a projectile size of ~10 m diameter (Bottke et al., 2005a,b). The shape of the impactor SFD for small projectile sizes does not appear to have changed over the last several billions of years, and an argument can be made that the absolute number of small MBAs has remained roughly constant (within a factor of 2) over the same time period. The apparent steady state nature of the main belt population potentially provides us with a set of intriguing constraints that can be used to glean insights into the physical evolution of individual MBAs as well as the main belt as an ensemble.Comment: Accepted by PSS, to appear on Vesta cratering special issu

Modulation induced frequency shifts in a CPT-based atomic clock

We investigate systematic errors associated with a common modulation technique used for phase sensitive detection of a coherent population trapping (CPT) resonance. In particular, we show that modification of the CPT resonance lineshape due to the presence of off-resonant fields leads to frequency shifts which may limit the stability of CPT-based atomic clocks. We also demonstrate that an alternative demodulation technique greatly reduces these effects.Comment: 14 pages, 7 figure