Aims: Investigate how the strength of episodic accretion bursts depends on eccentricity. Methods: We investigate the binary trigger hypothesis in longer-period (>20yr) binaries by carrying out three-dimensional magnetohydrodynamical (MHD) simulations of the formation of low-mass binary stars down to final separations of ~10AU, including the effects of gas turbulence and magnetic fields. We run two simulations with an initial turbulent gas core of one solar mass each and two different initial turbulent Mach numbers, M = 0.1 and 0.2, for 6500yr after protostar formation. Results: We observe bursts of accretion at periastron during the early stages when the eccentricity of the binary system is still high. We find that this correlation between bursts of accretion and passing periastron breaks down at later stages, because of the gradual circularisation of the orbits. For eccentricities greater than e=0.2, we observe episodic accretion triggered near periastron. However, we do not find any strong correlation between the strength (the ratio of the burst accretion ratio to the quiescent accretion rate) of episodic accretion and eccentricity. We determine that accretion events are likely triggered by torques between the rotation of the circumstellar disc and the approaching binary stars. We compare our results with observational data of episodic accretion in short-period binaries and find good agreement between our simulations and the observations. Conclusions: We conclude that episodic accretion is a universal mechanism operating in eccentric young binary-star systems, independent of separation, and should be observable in long-period binaries as well as in short-period binaries, but that the strength will depend on the torques, and hence the separation at periastron.

We address the problem of the origin of massive stars, namely the origin, path and timescale of the mass flows that create them. Based on extensive numerical simulations, we propose a scenario where massive stars are assembled by large-scale, converging, inertial flows that naturally occur in supersonic turbulence. We refer to this scenario of massive-star formation as the {\it Inertial-Inflow Model}. This model stems directly from the idea that the mass distribution of stars is primarily the result of {\it turbulent fragmentation}. Under this hypothesis, the statistical properties of the turbulence determine the formation timescale and mass of prestellar cores, posing definite constraints on the formation mechanism of massive stars. We quantify such constraints by the analysis of a simulation of supernova-driven turbulence in a 250-pc region of the interstellar medium, describing the formation of hundreds of massive stars over a time of approximately 30 Myr. Due to the large size of our statistical sample, we can say with full confidence that {\it massive stars in general do not form from the collapse of massive cores, nor from competitive accretion}, as both models are incompatible with the numerical results. We also compute synthetic continuum observables in Herschel and ALMA bands. We find that, depending on the distance of the observed regions, estimates of core mass based on commonly-used methods may exceed the actual core masses by up to two orders of magnitude, and that there is essentially no correlation between estimated and real core masses.

Jupiter and Saturn formed early, before the gas disk dispersed. The presence of gap-opening planets affects the dynamics of the gas and embedded solids and halts the inward drift of grains above a certain size. A drift barrier can explain the absence of calcium aluminium rich inclusions (CAIs) in chondrites originating from parent bodies that accreted in the inner solar system. Employing an interdisciplinary approach, we use a $\mu$-X-Ray-fluorescence scanner to search for large CAIs and a scanning electron microscope to search for small CAIs in the ordinary chondrite NWA 5697. We carry out long-term, two-dimensional simulations including gas, dust, and planets to characterize the transport of grains within the viscous $\alpha$-disk framework exploring the scenarios of a stand-alone Jupiter, Jupiter and Saturn \textit{in situ}, or Jupiter and Saturn in a 3:2 resonance. In each case, we find a critical grain size above which drift is halted as a function of the physical conditions in the disk. From the laboratory search we find four CAIs with a largest size of $\approx$200$\,\mu$m. \Combining models and data, we provide an estimate for the upper limit of the $\alpha$-viscosity and the surface density at the location of Jupiter, using reasonable assumptions about the stellar accretion rate during inward transport of CAIs, and assuming angular momentum transport to happen exclusively through viscous effects. Moreover, we find that the compound gap structure in the presence of Saturn in a 3:2 resonance favors inward transport of grains larger than CAIs currently detected in ordinary chondrites.

We present Lampray: a multi-group long characteristics ray tracing method for adaptive mesh radiation hydrodynamics in the Ramses code. It avoids diffusion, captures shadows, and treats colliding beams correctly, and therefore complements existing moment-based ray tracing in Ramses. Lampray includes different options for interpolation between ray and cell domain, and use either integral, Fourier, or an implicit method for hydrogen ionization to solve the radiative transfer. The opacity can either be tabulated or computed through a coupling to the general non-equilibrium astro-chemistry framework Krome. We use an H-He-C-O network with 36 species and 240 reactions to track the photo-chemistry in the interstellar medium across 6 and 10 orders of magnitude in temperature and density. Self-shielding prescriptions for H$_2$ and CO are used together with a new model for the diffuse interstellar UV-field. We also track the dust temperature, formation of H$_2$ on grains, and H$_2$O and CO ices in detail. Lampray is tested against standard benchmarks for molecular cloud and star formation simulations, including the formation of a Str\"omgren sphere, the expansion of an ionization front, the photo-evaporation of a dense clump, and the H-He-C-O chemistry in a static photo-dissociation front. Efficient parallelisation is achieved with a separate domain decomposition for rays where points along a ray reside in the same memory space, and data movement from cell- to ray-domain is done with a direct hash-table lookup algorithm. Point sources are treated without splitting rays, and therefore the method currently only scales to a few point sources, while diffuse radiation has excellent scaling.

(abridged) Using the particularly long gravitational microlensing event OGLE-2014-BLG-1186 with a time-scale $t_\mathrm{E}$ ~ 300 d, we present a methodology for identifying the nature of localised deviations from single-lens point-source light curves, which ensures that 1) the claimed signal is substantially above the noise floor, 2) the inferred properties are robustly determined and their estimation not subject to confusion with systematic noise in the photometry, 3) there are no alternative viable solutions within the model framework that might have been missed. Annual parallax and binarity could be separated and robustly measured from the wing and the peak data, respectively. We find matching model light curves that involve either a binary lens or a binary source. Our binary-lens models indicate a planet of mass $M_2$ = (45 $\pm$ 9) $M_\oplus$, orbiting a star of mass $M_1$ = (0.35 $\pm$ 0.06) $M_\odot$, located at a distance $D_\mathrm{L}$ = (1.7 $\pm$ 0.3) kpc from Earth, whereas our binary-source models suggest a brown-dwarf lens of $M$ = (0.046 $\pm$ 0.007) $M_\odot$, located at a distance $D_\mathrm{L}$ = (5.7 $\pm$ 0.9) kpc, with the source potentially being a (partially) eclipsing binary involving stars predicted to be of similar colour given the ratios between the luminosities and radii. The ambiguity in the interpretation would be resolved in favour of a lens binary by observing the luminous lens star separating from the source at the predicted proper motion of $\mu$ = (1.6 $\pm$ 0.3) mas yr$^{-1}$, whereas it would be resolved in favour of a source binary if the source could be shown to be a (partially) eclipsing binary matching the obtained model parameters. We experienced that close binary source stars pose a challenge for claiming the detection of planets by microlensing in events where the source passes very close to the lens star hosting the planet.

The current study was developed to provide a database of relatively simple numerical simulations of protostellar collapse, as a template library for observations of cores and very young protostars, and for researchers who wish to test their chemical modeling under dynamic astrophysical conditions. It was also designed to identify statistical trends that may appear when running many models of the formation of low-mass stars by varying the initial conditions. A large set of 143 calculations of the gravitational collapse of an isolated sphere of gas with uniform temperature and a Bonnor-Ebert like density profile was undertaken using a 1D fully implicit Lagrangian radiation hydrodynamics code. The parameter space covered initial masses from 0.2 to 8 Msun, temperatures of 5-30 K and radii between 3000 and 30,000 AU. A spread in the thermal evolutionary tracks of the runs was found, due to differing initial conditions and optical depths. Within less than an order of magnitude, all first and second Larson cores had masses and radii independent of the initial conditions. The time elapsed between the formation of the first and second cores was found to strongly depend on the first core mass accretion rate, and no first core in our grid of models lived for longer than 2000 years, before the onset of the second collapse. The end product of a protostellar cloud collapse, the second Larson core, is, at birth, a canonical object with a mass and radius of about 3 Mjup and 8 Rjup, independent of its initial conditions. The evolution sequence which brings the gas to stellar densities can however proceed in a variety of scenarios, on different timescales, along different isentropes, but each story line can largely be predicted by the initial conditions. All the data from the simulations are publicly available at this address: http://starformation.hpc.ku.dk/grid-of-protostars.

We describe a simple and effective algorithm for solving Poisson's equation in the context of self-gravity within the DISPATCH astrophysical fluid framework. The algorithm leverages the fact that DISPATCH stores multiple time slices and uses asynchronous time-stepping to produce a scheme that does not require any explicit global communication or sub-cycling, only the normal, local communication between patches and the iterative solution to Poisson's equation. We demonstrate that the implementation is suitable for both collections of patches of a single resolution and for hierarchies of adaptively resolved patches. Benchmarks are presented that demonstrate the accuracy, effectiveness and efficiency of the scheme.

Using zoom-simulations carried out with the adaptive mesh-refinement code RAMSES with a dynamic range of up to $2^{27} \approx 1.34 \times 10^8$ we investigate the accretion profiles around six stars embedded in different environments inside a (40 pc)$^3$ giant molecular cloud, the role of mass infall and disc instabilities on the accretion profile, and thus on the luminosity of the forming protostar. Our results show that the environment in which the protostar is embedded determines the overall accretion profile of the protostar. Infall on to the circumstellar disc may trigger gravitational disc instabilities in the disc at distances of around ~10 to ~50 au leading to rapid transport of angular momentum and strong accretion bursts. These bursts typically last for about ~10 to a ~100 yr, consistent with typical orbital times at the location of the instability, and enhance the luminosity of the protostar. Calculations with the stellar evolution code mesa show that the accretion bursts induce significant changes in the protostellar proper- ties, such as the stellar temperature and radius. We apply the obtained protostellar properties to produce synthetic observables with RADMC3D and predict that accretion bursts lead to ob- servable enhancements around 20 to 200 $\mu$m in the spectral energy distribution of Class 0 type young stellar objects.

We address the turbulent fragmentation scenario for the origin of the stellar initial mass function (IMF), using a large set of numerical simulations of randomly driven supersonic MHD turbulence. The turbulent fragmentation model successfully predicts the main features of the observed stellar IMF assuming an isothermal equation of state without any stellar feedback. As a test of the model, we focus on the case of a magnetized isothermal gas, neglecting stellar feedback, while pursuing a large dynamic range in both space and timescales covering the full spectrum of stellar masses from brown dwarfs to massive stars. Our simulations represent a generic 4 pc region within a typical Galactic molecular cloud, with a mass of 3000 Msun and an rms velocity 10 times the isothermal sound speed and 5 times the average Alfven velocity, in agreement with observations. We achieve a maximum resolution of 50 au and a maximum duration of star formation of 4.0 Myr, forming up to a thousand sink particles whose mass distribution closely matches the observed stellar IMF. A large set of medium-size simulations is used to test the sink particle algorithm, while larger simulations are used to test the numerical convergence of the IMF and the dependence of the IMF turnover on physical parameters predicted by the turbulent fragmentation model. We find a clear trend toward numerical convergence and strong support for the model predictions, including the initial time evolution of the IMF. We conclude that the physics of isothermal MHD turbulence is sufficient to explain the origin of the IMF.

Hertzsprung-Russell diagrams of star forming regions show a large luminosity spread. This is incompatible with well-defined isochrones based on classic non-accreting protostellar evo- lution models. Protostars do not evolve in isolation of their environment, but grow through accretion of gas. In addition, while an age can be defined for a star forming region, the ages of individual stars in the region will vary. We show how the combined effect of a protostellar age spread, a consequence of sustained star formation in the molecular cloud, and time-varying protostellar accretion for individual protostars can explain the observed luminosity spread. We use a global MHD simulation including a sub-scale sink particle model of a star forming region to follow the accretion process of each star. The accretion profiles are used to compute stellar evolution models for each star, incorporating a model of how the accretion energy is distributed to the disk, radiated away at the accretion shock, or incorporated into the outer layers of the protostar. Using a modelled cluster age of 5 Myr we naturally reproduce the lumi- nosity spread and find good agreement with observations of the Collinder 69 cluster, and the Orion Nebular Cluster. It is shown how stars in binary and multiple systems can be externally forced creating recurrent episodic accretion events. We find that in a realistic global molecular cloud model massive stars build up mass over relatively long time-scales. This leads to an important conceptual change compared to the classic picture of non-accreting stellar evolution segmented in to low-mass Hayashi tracks and high-mass Henyey tracks.

In this paper the non-linear effect of massive neutrinos on cosmological structures is studied in a conceptually new way. We have solved the non-linear continuity and Euler equations for the neutrinos on a grid in real space in $N$-body simulations, and closed the Boltzmann hierarchy at the non-linear Euler equation using the stress and pressure perturbations from linear theory. By comparing with state-of-the art cosmological neutrino simulations, we are able to simulate the non-linear neutrino power spectrum very accurately. This translates into a negligible error in the matter power spectrum, and so our CONCEPT code is ideally suited for extracting the neutrino mass from future high precision non-linear observational probes such as Euclid.

Stars and their corresponding protoplanetary disks form in diverse environments. To account for these natural variations, we investigate the formation process around nine solar mass stars with a maximum resolution of 2 AU in a Giant Molecular Cloud of (40 pc)$^3$ in volume by using the adaptive mesh refinement code \ramses. The magnetohydrodynamic simulations reveal that the accretion process is heterogeneous in time, in space, and among protostars of otherwise similar mass. During the first roughly 100 kyr of a protostar evolving to about a solar mass, the accretion rates peak around $10^{-5}$ to $10^{-4}$ M$_{\odot}$ yr$^{-1}$ shortly after its birth, declining with time after that. The different environments also affect the spatial accretion, and infall of material to the star-disk system is mostly through filaments and sheets. Furthermore, the formation and evolution of disks varies significantly from star to star. We interpret the variety in disk formation as a consequence of the differences in the combined effects of magnetic fields and turbulence that may cause differences in the efficiency of magnetic braking, as well as differences in the strength and distribution of specific angular momentum.

We compute the star formation rate (SFR) in molecular clouds (MCs) that originate {\it ab initio} in a new, higher-resolution simulation of supernova-driven turbulence. Because of the large number of well-resolved clouds with self-consistent boundary and initial conditions, we obtain a large range of cloud physical parameters with realistic statistical distributions, an unprecedented sample of star-forming regions to test SFR models and to interpret observational surveys. We confirm the dependence of the SFR per free-fall time, $SFR_{\rm ff}$, on the virial parameter, $\alpha_{\rm vir}$, found in previous simulations, and compare a revised version of our turbulent fragmentation model with the numerical results. The dependences on Mach number, ${\cal M}$, gas to magnetic pressure ratio, $\beta$, and compressive to solenoidal power ratio, $\chi$ at fixed $\alpha_{\rm vir}$ are not well constrained, because of random scatter due to time and cloud-to-cloud variations in $SFR_{\rm ff}$. We find that $SFR_{\rm ff}$ in MCs can take any value in the range $0 \le SFR_{\rm ff} \lesssim 0.2$, and its probability distribution peaks at a value $SFR_{\rm ff}\approx 0.025$, consistent with observations. The values of $SFR_{\rm ff}$ and the scatter in the $SFR_{\rm ff}$--$\alpha_{\rm vir}$ relation are consistent with recent measurements in nearby MCs and in clouds near the Galactic center. Although not explicitly modeled by the theory, the scatter is consistent with the physical assumptions of our revised model and may also result in part from a lack of statistical equilibrium of the turbulence, due to the transient nature of MCs.

We present a combined analysis of the observations of the gravitational microlensing event OGLE-2015-BLG-0479 taken both from the ground and by the {\it Spitzer Space Telescope}. The light curves seen from the ground and from space exhibit a time offset of $\sim 13$ days between the caustic spikes, indicating that the relative lens-source positions seen from the two places are displaced by parallax effects. From modeling the light curves, we measure the space-based microlens parallax. Combined with the angular Einstein radius measured by analyzing the caustic crossings, we determine the mass and distance of the lens. We find that the lens is a binary composed of two G-type stars with masses $\sim 1.0\ M_\odot$ and $\sim 0.9\ M_\odot$ located at a distance $\sim 3$ kpc. In addition, we are able to constrain the complete orbital parameters of the lens thanks to the precise measurement of the microlens parallax derived from the joint analysis. In contrast to the binary event OGLE-2014-BLG-1050, which was also observed by {\it Spitzer}, we find that the interpretation of OGLE-2015-BLG-0479 does not suffer from the degeneracy between $(\pm,\pm)$ and $(\pm,\mp)$ solutions, confirming that the four-fold parallax degeneracy in single-lens events collapses into the two-fold degeneracy for the general case of binary-lens events. The location of the blend in the color-magnitude diagram is consistent with the lens properties, suggesting that the blend is the lens itself. The blend is bright enough for spectroscopy and thus this possibility can be checked from future follow-up observations.

We show the benefits of using Electron-Multiplying CCDs and the shift-and-add technique as a tool to minimise the effects of the atmospheric turbulence such as blending between stars in crowded fields and to avoid saturated stars in the fields observed. We intend to complete, or improve, the census of the variable star population in globular cluster NGC~6715. Our aim is to obtain high-precision time-series photometry of the very crowded central region of this stellar system via the collection of better angular resolution images than has been previously achieved with conventional CCDs on ground-based telescopes. Observations were carried out using the Danish 1.54-m Telescope at the ESO La Silla observatory in Chile. The telescope is equipped with an Electron-Multiplying CCD that allowed to obtain short-exposure-time images (ten images per second) that were stacked using the shift-and-add technique to produce the normal-exposure-time images (minutes). The high precision photometry was performed via difference image analysis employing the DanDIA pipeline. We attempted automatic detection of variable stars in the field. We statistically analysed the light curves of 1405 stars in the crowded central region of NGC~6715 to automatically identify the variable stars present in this cluster. We found light curves for 17 previously known variable stars near the edges of our reference image (16 RR Lyrae and 1 semi-regular) and we discovered 67 new variables (30 RR Lyrae, 21 long-period irregular, 3 semi-regular, 1 W Virginis, 1 eclipsing binary, and 11 unclassified). Photometric measurements for these stars are available in electronic form through the Strasbourg Astronomical Data Centre.

We present a comparison of molecular clouds (MCs) from a simulation of supernova-driven interstellar medium (ISM) turbulence with real MCs from the Outer Galaxy Survey. The radiative transfer calculations to compute synthetic CO spectra are carried out assuming the CO relative abundance depends only on gas density, according to four different models. Synthetic MCs are selected above a threshold brightness temperature value, $T_{\rm B,min}=1.4$ K, of the $J=1-0$ $^{12}$CO line, generating 16 synthetic catalogs (four different spatial resolutions and four CO abundance models), each containing up to several thousands MCs. The comparison with the observations focuses on the mass and size distributions and on the velocity-size and mass-size Larson relations. The mass and size distributions are found to be consistent with the observations, with no significant variations with spatial resolution or chemical model, except in the case of the unrealistic model with constant CO abundance. The velocity-size relation is slightly too steep for some of the models, while the mass-size relation is a bit too shallow for all models only at a spatial resolution $dx\approx 1$ pc. The normalizations of the Larson relations show a clear dependence on spatial resolution, for both the synthetic and the real MCs. The comparison of the velocity-size normalization suggests that the SN rate in the Perseus arm is approximately 70\% or less of the rate adopted in the simulation. Overall, the realistic properties of the synthetic clouds confirm that supernova-driven turbulence can explain the origin and dynamics of MCs.

(abridged) Context. Wide binaries are a potential pathway for the formation of hot Jupiters. The binary fraction among host stars is an important discriminator between competing formation theories, but has not been well characterised. Additionally, contaminating light from unresolved stars can significantly affect the accuracy of photometric and spectroscopic measurements in studies of transiting exoplanets. Aims. We observed 101 transiting exoplanet host systems in the Southern hemisphere in order to create a homogeneous catalogue of both bound companion stars and contaminating background stars. We investigate the binary fraction among the host stars in order to test theories for the formation of hot Jupiters, in an area of the sky where transiting exoplanetary systems have not been systematically searched for stellar companions. Methods. Lucky imaging observations from the Two Colour Instrument on the Danish 1.54m telescope at La Silla were used to search for previously unresolved stars at small angular separations. The separations and relative magnitudes of all detected stars were measured. For 12 candidate companions to 10 host stars, previous astrometric measurements were used to evaluate how likely the companions are to be physically associated. Results. We provide measurements of 499 candidate companions within 20 arcseconds of our sample of 101 planet host stars. 51 candidates are located within 5 arcseconds of a host star, and we provide the first published measurements for 27 of these. Calibrations for the plate scale and colour performance of the Two Colour Instrument are presented. Conclusions. We find that the overall multiplicity rate of the host stars is 38 +17 -13%, consistent with the rate among solar-type stars in our sensitivity range, suggesting that planet formation does not preferentially occur in long period binaries compared to a random sample of field stars.

Obtain time-series photometry of the very crowded central regions of Galactic globular clusters with better angular resolution than previously achieved with conventional CCDs on ground-based telescopes to complete, or improve, the census of the variable star population in those stellar systems. Images were taken using the Danish 1.54-m Telescope at the ESO observatory at La Silla in Chile. The telescope was equipped with an electron-multiplying CCD and the short-exposure-time images obtained (10 images per second) were stacked using the shift-and-add technique to produce the normal-exposure-time images (minutes). Photometry was performed via difference image analysis. Automatic detection of variable stars in the field was attempted. The light curves of 12541 stars in the cores of 10 globular clusters were statistically analysed in order to automatically extract the variable stars. We obtained light curves for 31 previously known variable stars (3 L, 2 SR, 20 RR Lyrae, 1 SX Phe, 3 cataclysmic variables, 1 EW and 1 NC) and we discovered 30 new variables (16 L, 7 SR, 4 RR Lyrae, 1 SX Phe and 2 NC).

Context: Understanding how protostars accrete their mass is a central question of star formation. One aspect of this is trying to understand whether the time evolution of accretion rates in deeply embedded objects is best characterised by a smooth decline from early to late stages or by intermittent bursts of high accretion. Aims: We create synthetic observations of deeply embedded protostars in a large numerical simulation of a molecular cloud, which are compared directly to real observations. The goal is to compare episodic accretion events in the simulation to observations and to test the methodology used for analysing the observations. Methods: Simple freeze-out and sublimation chemistry is added to the simulation, and synthetic C$^{18}$O line cubes are created for a large number of simulated protostars. The spatial extent of C$^{18}$O is measured for the simulated protostars and compared directly to a sample of 16 deeply embedded protostars observed with the Submillimeter Array. If CO is distributed over a larger area than predicted based on the protostellar luminosity, it may indicate that the luminosity has been higher in the past and that CO is still in the process of refreezing. Results: Approximately 1% of the protostars in the simulation show extended C$^{18}$O emission, as opposed to approximately 50% in the observations, indicating that the magnitude and frequency of episodic accretion events in the simulation is too low relative to observations. The protostellar accretion rates in the simulation are primarily modulated by infall from the larger scales of the molecular cloud, and do not include any disk physics. The discrepancy between simulation and observations is taken as support for the necessity of disks, even in deeply embedded objects, to produce episodic accretion events of sufficient frequency and amplitude.

(abridged) Context: Both observations and simulations of embedded protostars have progressed rapidly in recent years. Bringing them together is an important step in advancing our knowledge about the earliest phases of star formation. Aims: To compare synthetic continuum images and SEDs, calculated from large-scale numerical simulations, to observational studies, thereby aiding in both the interpretation of the observations and in testing the fidelity of the simulations. Methods: The radiative transfer code RADMC-3D is used to create synthetic continuum images and SEDs of protostellar systems in a large numerical simulation of a molecular cloud. More than 13000 unique radiative transfer models are produced of a variety of different protostellar systems. Results: Over the course of 0.76 Myr the simulation forms more than 500 protostars, primarily within two sub-clusters. Synthetic SEDs are used to the calculate evolutionary tracers Tbol and Lsmm/Lbol. It is shown that, while the observed distributions of the tracers are well matched by the simulation, they generally do a poor job of tracking the protostellar ages. Disks form early in the simulation, with 40 % of the Class 0 protostars being encircled by one. The flux emission from the simulated disks is found to be, on average, a factor of 6 too low relative to real observations. The distribution of protostellar luminosities spans more than three order of magnitudes, similar to the observed distribution. Cores and protostars are found to be closely associated with one another, with the distance distribution between them being in excellent agreement with observations. Conclusions: The analysis and statistical comparison of synthetic observations to real ones is established as a powerful tool in the interpretation of observational results.

We have carefully studied how local measurements of the Hubble constant, $H_0$, can be influenced by a variety of different parameters related to survey geometry, depth, and size, as well as observer position in space. Our study is based on N-body simulations of structure in the standard $\Lambda$CDM model and our conclusion is that the expected variance in measurements of $H_0$ is far too small to explain the current discrepancy between the low value of $H_0$ inferred from measurements of the cosmic microwave background (CMB) by the Planck collaboration and the value measured directly in the local universe by use of Type Ia supernovae. This conclusion is very robust and does not change with different assumptions about survey geometry and effective sky coverage or observer position in space.

In this paper, we present the analysis of time-series observations from 2013 and 2014 of five metal rich ([Fe/H] $>$ -1) globular clusters: NGC~6388, NGC~6441, NGC~6528, NGC~6638, and NGC~6652. The data have been used to perform a census of the variable stars in the central parts of these clusters. The observations were made with the electron multiplying CCD (EMCCD) camera at the Danish 1.54m Telescope at La Silla, Chile, and they were analysed using difference image analysis (DIA) to obtain high-precision light curves of the variable stars. It was possible to identify and classify all of the previously known or suspected variable stars in the central regions of the five clusters. Furthermore, we were able to identify, and in most cases classify 48, 49, 7, 8, and 2 previously unknown variables in NGC~6388, NGC~6441, NGC~6528, NGC~6638, and NGC~6652, respectively. Especially interesting is the case of NGC~6441, for which the variable star population of about 150 stars has been thoroughly examined by previous studies, including a Hubble Space Telescope study. In this paper we are able to present 49 new variable stars for this cluster, of which one (possibly two) are RR Lyrae stars, two are W Virginis stars, and the rest are long period semi-regular/irregular variables on the red giant branch. We have also detected the first double mode RR Lyrae in the cluster.

We investigate the role of mass infall in the formation and evolution of protostars. To avoid ad hoc initial and boundary conditions, we consider the infall resulting self-consistently from modeling the formation of stellar clusters in turbulent molecular clouds. We show that infall rates in turbulent clouds are comparable to accretion rates inferred from protostellar luminosities or measured in pre-main-sequence stars. They should not be neglected in modeling the luminosity of protostars and the evolution of disks, even after the embedded protostellar phase. We find large variations of infall rates from protostar to protostar, and large fluctuations during the evolution of individuals protostars. In most cases, the infall rate is initially of order 10$^{-5}$\msun\ yr$^{-1}$, and may either decay rapidly in the formation of low-mass stars, or remain relatively large when more massive stars are formed. The simulation reproduces well the observed characteristic values and scatter of protostellar luminosities and matches the observed protostellar luminosity function. The luminosity problem is therefore solved once realistic protostellar infall histories are accounted for, with no need for extreme accretion episodes. These results are based on a simulation of randomly-driven magneto-hydrodynamic turbulence on a scale of 4pc, including self-gravity, adaptive-mesh refinement to a resolution of 50AU, and accreting sink particles. The simulation yields a low star formation rate, consistent with the observations, and a mass distribution of sink particles consistent with the observed stellar initial mass function during the whole duration of the simulation, forming nearly 1,300 sink particles over 3.2 Myr.

Context: The R CrA cloud hosts a handful of Class 0/I low-mass young stellar objects. The chemistry and physics at scales $>500$ AU in this cloud are dominated by the irradiation from the nearby Herbig Be star R CrA. The luminous large-scale emission makes it necessary to use high-resolution spectral imaging to study the chemistry and dynamics of the inner envelopes and discs of the protostars. Aims: We aim to better understand the structure of the inner regions of these protostars and, in particular, the interplay between the chemistry and the presence of discs. Methods: Using Atacama Large Millimeter/submillimeter Array (ALMA) high-resolution spectral imaging interferometry observations, we study the molecular line and dust continuum emission at submillimetre wavelengths. Results: We detect dust continuum emission from four circumstellar discs around Class 0/I objects within the R CrA cloud. Towards IRS7B we detect C$^{17}$O emission showing a rotation curve consistent with a Keplerian disc with a well-defined edge that gives a good estimate for the disc radius at 50 AU. We derive the central object mass to $2.3M_{\odot}$ and the disc mass to $0.024M_{\odot}$. The observations are also consistent with a model of material infalling under conservation of angular momentum; however, this model provides a worse fit to the data. We also report a likely detection of faint CH$_3$OH emission towards this point source, as well as more luminous CH$_3$OH emission in an outflow orthogonal to the major axis of the C$^{17}$O emission. Conclusions: The faint CH$_3$OH emission seen towards IRS7B can be explained by a flat density profile of the inner envelope caused by the disc with a radius $\lesssim50$ AU. We propose that the regions of the envelopes where complex organic molecules are present in Class 0/I young stellar objects can become quenched as the disc grows.

The primary focus of this paper is on the particle acceleration mechanism in solar coronal three-dimensional reconnection null-point regions. Starting from a potential field extrapolation of a Solar and Heliospheric Observatory (SOHO) magnetogram taken on 2002 November 16, we first performed magnetohydrodynamics (MHD) simulations with horizontal motions observed by SOHO applied to the photospheric boundary of the computational box. After a build-up of electric current in the fan-plane of the null-point, a sub-section of the evolved MHD data was used as initial and boundary conditions for a kinetic particle-in-cell model of the plasma. We find that sub-relativistic electron acceleration is mainly driven by a systematic electric field in the current sheet. A non-thermal population of electrons with a power-law distribution in energy forms in the simulated pre-flare phase, featuring a power-law index of about -1.78. This work provides a first step towards bridging the gap between macroscopic scales on the order of hundreds of Mm and kinetic scales on the order of cm in the solar corona, and explains how to achieve such a cross-scale coupling by utilizing either physical modifications or (equivalent) modifications of the constants of nature. With their exceptionally high resolution - up to 135 billion particles and 3.5 billion grid cells of size 17.5 km - these simulations offer a new opportunity to study particle acceleration in solar-like settings.

We present a novel method for implementing massive neutrinos in N-body simulations. Instead of sampling the neutrino velocity distribution by individual point particles we take neutrino free-streaming into account by treating it as an effective redshift dependent sound speed in a perfect isothermal fluid, and assume a relation between the sound speed and velocity dispersion of the neutrinos. Although the method fails to accurately model the true neutrino power spectrum, it is able to calculate the total matter power spectrum to the same accuracy as more complex hybrid neutrino methods, except on very small scales. We also present an easy way to update the publicly available Gadget-2 version with this neutrino approximation.

We present an analysis of peculiar velocities and their effect on supernova cosmology. In particular, we study (a) the corrections due to our own motion, (b) the effects of correlations in peculiar velocities induced by large-scale structure, and (c) uncertainties arising from a possible local under- or over-density. For all of these effects we present a case study of their impact on the cosmology derived by the Sloan Digital Sky Survey-II Supernova Survey (SDSS-II SN Survey). Correcting supernova redshifts for the CMB dipole slightly over-corrects nearby supernovae that share some of our local motion. We show that while neglecting the CMB dipole would cause a shift in the derived equation of state of Delta w ~ 0.04 (at fixed matter density) the additional local-motion correction is currently negligible (Delta w<0.01). We use a covariance-matrix approach to statistically account for correlated peculiar velocities. This down-weights nearby supernovae and effectively acts as a graduated version of the usual sharp low-redshift cut. Neglecting coherent velocities in the current sample causes a systematic shift of ~2% in the preferred value of w and will therefore have to be considered carefully when future surveys aim for percent-level accuracy. Finally, we perform n-body simulations to estimate the likely magnitude of any local density fluctuation (monopole) and estimate the impact as a function of the low-redshift cutoff. We see that for this aspect the low-z cutoff of z=0.02 is well-justified theoretically, but that living in a putative local density fluctuation leaves an indelible imprint on the magnitude-redshift relation.

If the diffuse extragalactic gamma ray emission traces the large scale structures of the universe, peculiar anisotropy patterns are expected in the gamma ray sky. In particular, because of the cutoff distance introduced by the absorption of 0.1-10 TeV photons on the infrared/optical background, prominent correlations with the local structures within a range of few hundreds Mpc should be present. We provide detailed predictions of the signal based on the PSCz map of the local universe. We also use mock N-body catalogues complemented with the halo model of structures to study some statistical features of the expected signatures. The results are largely independent from cosmological details, and depend mostly on the index of correlation (or bias) of the sources with respect to the large scale distribution of galaxies. For instance, the predicted signal in the case of a quadratic correlation (as it may happen for a dark matter annihilation contribution to the diffuse gamma flux) differs substantially from a linear correlation case, providing a complementary tool to unveil the nature of the sources of the diffuse gamma ray emission. The chances of the present and future space and ground based observatories to measure these features are discussed.