Magnetic fields play an important role in the dynamics of present-day molecular clouds. Recent work has shown that magnetic fields are equally important for primordial clouds, which form the first stars in the Universe. While the primordial magnetic field strength on cosmic scales is largely unconstrained, theoretical models strongly suggest that a weak seed field existed in the early Universe. We study how the amplification of such a weak field can influence the evolution of accretion discs around the first stars, and thus affect the primordial initial mass function (IMF). We perform a suite of 3D magneto-hydrodynamic (MHD) simulations with different initial field strengths and numerical resolutions. We find that, in simulations with sufficient spatial resolution to resolve the Jeans scale during the collapse, even initially weak magnetic fields grow exponentially to become dynamically important due to both the so-called \textit{small-scale turbulent dynamo} and the \textit{large-scale mean-field dynamo}. Capturing the small-scale dynamo action depends primarily on how well we resolve the Jeans length, while capturing the large-scale dynamo depends on the Jeans resolution as well as the maximum absolute resolution. Provided enough resolution, we find that fragmentation does not depend strongly on the initial field strength, because even weak fields grow to become strong. However, fragmentation in runs with magnetic fields differs significantly from those without magnetic fields. We conclude that the development of dynamically strong magnetic fields during the formation of the first stars is likely inevitable, and that these fields had a significant impact on the primordial IMF.

The subtype of hot algol semidetached binaries dubbed Double Periodic Variables (DPVs) are characterized by a photometric cycle longer than the orbital one, whose nature has been related to a magnetic dynamo in the donor component controlling the mass transfer rate. We aim to understand the morphologic changes observed in the light curve of OGLE-BLG-ECL-157529 that are linked to the long cycle. In particular, we want to explain the changes in relative depth of primary and secondary eclipses. We analyze $I$ and $V$-band OGLE photometric times series spanning 18.5 years and model the orbital light curve. We find that OGLE-BLG-ECL-157529 is a new eclipsing Galactic DPV of orbital period 24\fd8, and that its long cycle length decreases in amplitude and length during the time baseline. We show that the changes of the orbital light curve can be reproduced considering an accretion disk of variable thickness and radius, surrounding the hottest stellar component. Our models indicate changes in the temperatures of hot spot and bright spot during the long cycle, and also in the position of the bright spot. This, along with the changes in disk radius might indicate a variable mass transfer in this system.

Runaway stellar collisions in dense star clusters are invoked to explain the presence of very massive stars or blue stragglers in the center of those systems. This process has also been explored for the first star clusters in the Universe and shown to yield stars that may collapse at some points into an intermediate mass black hole. Although the early evolution of star clusters requires the explicit modeling of the gas out of which the stars form, these calculations would be extremely time-consuming and often the effects of the gas can be accurately treated by including a background potential to account for the extra gravitational force. We apply this approximation to model the early evolution of the first dense star clusters formed in the Universe by performing $N$-body simulations, our goal is to understand how the additional gravitational force affects the growth of a very massive star through stellar mergers in the central parts of the star cluster. Our results show that the background potential increases the velocities of the stars, causing an overall delay in the evolution of the clusters and in the runaway growth of a massive star at the center. The population of binary stars is lower due to the increased kinetic energy of the stars, initially reducing the number of stellar collisions, and we show that relaxation processes are also affected. Despite these effects, the external potential enhances the mass of the merger product by a factor $\sim$2 if the collisions are maintained for long times.

While the stellar Initial Mass Function (IMF) appears to be close to universal within the Milky Way galaxy, it is strongly suspected to be different in the primordial Universe, where molecular hydrogen cooling is less efficient and the gas temperature can be higher by a factor of 30. In between these extreme cases, the gas temperature varies depending on the environment, metallicity and radiation background. In this paper we explore if changes of the gas temperature affect the IMF of the stars considering fragmentation and accretion. The fragmentation behavior depends mostly on the Jeans mass at the turning point in the equation of state where a transition occurs from an approximately isothermal to an adiabatic regime due to dust opacities. The Jeans mass at this transition in the equation of state is always very similar, independent of the initial temperature, and therefore the initial mass of the fragments is very similar. Accretion on the other hand is strongly temperature dependent. We argue that the latter becomes the dominant process for star formation efficiencies above 5 - 7 %, increasing the average mass of the stars.

Cosmic rays are a global source of ionisation, and the ionisation fraction represents a fundamental parameter in the interstellar medium. Ions couple to magnetic fields, affect the chemistry, and the dynamics of star-forming regions as well as planetary atmospheres. However, the cosmic-ray ionisation rate represents one of the bottlenecks for astrochemical models, and its determination is one of the most puzzling problems in astrophysics. While for diffuse clouds reasonable values have been provided from H$_3^+$ observations, for dense clouds, due to the lack of rotational transitions, this is not possible, and estimates are strongly biased by the employed model. We present here an analytical expression, obtained from first principles, to estimate the cosmic-ray ionisation rate from observational quantities. The theoretical predictions are validated with high-resolution three-dimensional numerical simulations and applied to the well known core L1544; we obtained an estimate of $\zeta_2 \sim 2-3 \times 10^{-17}$ s$^{-1}$. Our results and the analytical formulae provided represent the first model-independent, robust tool to probe the cosmic-ray ionisation rate in the densest part of star-forming regions (on spatial scales of $R \leq 0.05$ pc). An error analysis is presented to give statistical relevance to our study.

Runaway collisions in dense clusters may lead to the formation of supermassive black hole (SMBH) seeds, and this process can be further enhanced by accretion, as recent models of SMBH seed formation in Population III star clusters have shown. This may explain the presence of supermassive black holes already at high redshift, $z>6$. However, in this context, mass loss during collisions was not considered and could play an important role for the formation of the SMBH seed. Here, we study the effect of mass loss, due to collisions of protostars, in the formation and evolution of a massive object in a dense primordial cluster. We consider both constant mass loss fractions as well as analytic models based on the stellar structure of the collision components. Our calculations indicate that mass loss can significantly affect the final mass of the possible SMBH seed. Considering a constant mass loss of 5% for every collision, we can lose between 60-80% of the total mass that is obtained if mass loss were not considered. Using instead analytical prescriptions for mass loss, the mass of the final object is reduced by 15-40%, depending on the accretion model for the cluster we study. Altogether, we obtain masses of the order of $10^4M_{\odot}$, which are still massive enough to be SMBH seeds.

The Event Horizon Telescope (EHT), now with its first ever image of the photon ring around the supermassive black hole of M87, provides a unique opportunity to probe the physics of supermassive black holes through Very Long Baseline Interferometry (VLBI), such as the existence of the event horizon, the accretion processes as well as jet formation in Low Luminosity AGN (LLAGN). We build a theoretical model which includes an Advection Dominated Accretion Flow (ADAF) with emission from thermal and non-thermal electrons in the flow and a simple radio jet outflow. The predicted spectral energy distribution (SED) of this model is compared to sub-arcsec resolution observations to get the best estimates of the model parameters. The model-predicted radial emission profiles at different frequency bands are used to predict whether the inflow can be resolved by the EHT or with telescopes such as the Global 3-mm VLBI array (GMVA). In this work the model is initially tested with high resolution SED data of M87 and then applied to our sample of 5 galaxies (Cen A, M84, NGC 4594, NGC 3998 and NGC 4278). The model then allows us to predict if one can detect and resolve the inflow for any of these galaxies using the EHT or GMVA within an 8 hour integration time.

Disentangling the different stages of the star-formation process, in particular in the high-mass regime, is a challenge in astrophysics. Chemical clocks could help alleviating this problem, but their evolution strongly depends on many parameters, leading to degeneracy in the interpretation of the observational data. One of these uncertainties is the degree of CO depletion. We present here the first self-consistent magneto-hydrodynamic simulations of high-mass star-forming regions at different scales, fully coupled with a non-equilibrium chemical network, which includes C-N-O bearing molecules. Depletion and desorption processes are treated time-dependently. The results show that full CO-depletion (i.e. all gas-phase CO frozen-out on the surface of dust grains), can be reached very quickly, in one third or even smaller fractions of the free-fall time, whether the collapse proceeds on slow or fast timescales. This leads to a high level of deuteration in a short time both for typical tracers like N$_2$H$^+$, as well as for the main ion H$_3^+$, the latter being in general larger and more extended. N$_2$ depletion is slightly less efficient, and no direct effects on N-bearing molecules and deuterium fractionation are observed. We show that CO depletion is not the only driver of deuteration, and that there is a strong impact on $D_{frac}$ when changing the grain-size. We finally apply a two-dimensional gaussian Point Spread Function to our results to mimic observations with single-dish and interferometers. Our findings suggest that the low-values observed in high-mass star-forming clumps are in reality masking a full-depletion stage in the inner 0.1 pc region.

Eclipsing time variations have been observed for a wide range of binary systems, including post-common-envelope binaries. A frequently proposed explanation, apart from the possibility of having a third body, is the effect of magnetic activity, which may alter the internal structure of the secondary star, particularly its quadrupole moment, and thereby cause quasi-periodic oscillations. Here we present two compressible non-ideal magneto-hydrodynamical (MHD) simulations of the magnetic dynamo in a solar mass star, one of them with three times the solar rotation rate ("slow rotator"), the other one with twenty times the solar rotation rate ("rapid rotator"), to account for the high rotational velocities in close binary systems. For the slow rotator, we find that both the magnetic field and the stellar quadrupole moment change in a quasi-periodic manner, leading to O-C (observed - corrected times of the eclipse) variations of ~0.025 s. For the rapid rotator, the behavior of the magnetic field as well as the quadrupole moment changes become considerably more complex, due to the less coherent dynamo solution. The resulting O-C variations are of the order 0.13 s. The observed system V471~Tau shows two modes of eclipsing time variations, with amplitudes of 151 s and 20 s, respectively. However, the current simulations may not capture all relevant effects due to the neglect of the centrifugal force and self-gravity. Considering the model limitations and that the rotation of V471 Tau is still a factor of 2.5 faster than our rapid rotator, it may be conceivable to reach the observed magnitudes.

We have performed cosmological hydrodynamic simulations to study the effect of fragmentation on the SMBH seed mass in the direct collapse formation scenario. We considered different background UV intensities, host halo spin, and halo merger histories. Our simulations in low-spin halos, in the presence of a strong UV background are consistent with the Direct Collapse Black Hole model, in which a single massive object $\sim10^5$ M$_{\odot}$ is formed in the center of a proto-galaxy. While in our simulations under the presence of a low UV background, we find fragmentation and the formation of various minor seeds. These fragments have masses of 10$^3$ - 10$^4$ M$_{\odot}$. These values are significant if we consider the potential mergers between them, and the fact that these minor objects are formed earlier in cosmic time compared to the massive single seeds. In one of our simulations, we observed gas fragmentation even in the presence of a strong UV intensity. Said structure arose in a dark matter halo that formed after various merger episodes, and the one with the highest spin value. The final mass obtained was $\sim 10^5$ M$_{\odot}$ in this run. From these results, we conclude that fragmentation in fact produces less massive objects, however, they are still prone to merge. In simulations that formed many fragments, they all approach the most massive one with time. We see no uniqueness in the strength of the UV intensity value required to achieve a DCBH, since it depends in other factors like the system dynamics in our cases.

We report the serendipitous discovery of a nitrogen-rich, mildly metal-poor ([Fe/H]=-1.08) giant star in a single-lined spectroscopic binary system found in the SDSS-IV Apache Point Observatory Galactic Evolution Experiment (APOGEE-2) survey, Data Release 14 (DR14). Previous work has assumed that the two percent of halo giants with unusual elemental abundances have been evaporated from globular clusters, but other origins for their abundance signatures, including binary mass transfer, must also be explored. We present the results of an abundance re-analysis of the APOGEE-2 high-resolution near-infrared spectrum of 2M12451043+1217401 with the Brussels Automatic Stellar Parameter (BACCHUS) automated spectral analysis code, and re-derive manually the main element families, namely the light elements (C, N), elements (O, Mg, Si), iron-peak element (Fe), \textit{s}-process element (Ce), and the light odd-Z element (Al). Our analysis confirm the N-rich nature of 2M12451043+1217401, which has a [N/Fe] ratio of $+0.69$, and shows that the abundances of C and Al are slightly discrepant from that of a typical mildly metal-poor RGB star, but exhibit Mg, Si, O and \textit{s}-process abundances (Ce) of typical field stars. We also detect a particularly large variability in its radial velocity over the period of the APOGEE-2 observations, and the most likely orbit fit to the radial velocity data has a period of $730.89\pm106.86$ days, a velocity semi-amplitude of $9.92 \pm 0.14$ km s$^{-1}$, and an eccentricity of $\sim 0.1276 \pm0.1174$, which support the hypothesis of a binary companion, and that has probably been polluted by a now-extinct AGB star.

We present elemental abundance analysis of high-resolution spectra for five giant stars, deriving Fe, Mg, Al, C, N, O, Si and Ce abundances, and spatially located within the innermost regions of the bulge globular cluster NGC 6522, based on H-band spectra taken with the multi-object APOGEE-north spectrograph from the SDSS-IV Apache Point Observatory Galactic Evolution Experiment (APOGEE) survey. Of the five cluster candidates, two previously unremarked stars are confirmed to have second-generation (SG) abundance patterns, with the basic pattern of depletion in C and Mg simultaneous with enrichment in N and Al as seen in other SG globular cluster populations at similar metallicity. } In agreement with the most recent optical studies, the NGC 6522 stars analyzed exhibit (when available) only mild overabundances of the s-process element Ce, contradicting the idea of the NGC 6522 stars being formed from gas enriched by spinstars and indicating that other stellar sources such as massive AGB stars could be the primary intra-cluster medium polluters. The peculiar abundance signature of SG stars have been observed in our data, confirming the presence of multiple generations of stars in NGC 6522.

Observations reveal the presence of supermassive black holes (SMBH) as early as ~700 million years after the Big Bang. Their formation path is still subject to current debate. We explore the influence of magnetic fields, which are strongly amplified via the turbulent small-scale dynamo, on the formation of SMBH seeds within the direct collapse scenario. In this study, we perform for the first time cosmological magnetohydrodynamic large eddy simulations that employ a model for unresolved, compressible MHD turbulence. In total we perform 36 simulations for 9 haloes each with two different initial magnetic field strengths, and with and without employing the unresolved turbulence model. We make use of the adaptive mesh refinement approach to achieve an effective spatial resolution of less than one proper astronomical unit. We consider a regime where cooling is regulated by atomic hydrogen and the molecular hydrogen gets dissociated by a strong radiation field. Our main finding is that the majority of the gas properties in the haloes at the final output are predominantly determined by the run-away gravitational collapse. Turbulence is supersonic and super-Alfv\'enic in all cases, and magnetic fields are amplified to an approximately dynamically relevant regime. Finally, fragmentation during the collapse is intermittent and mass accretion rates range from 0.2-3 Msun/yr. This suggests that the presence of strongly amplified magnetic fields and turbulence provides additional pressure support on small scales and make the direct collapse a viable scenario for the formation of massive objects under the required ambient conditions.

Young massive stars are usually found embedded in dense and massive molecular clumps and are known for being highly obscured and distant. During their formation process, deuteration is regarded as a potentially good indicator of the formation stage. Therefore, proper observations of such deuterated molecules are crucial, but still, hard to perform. In this work, we test the observability of the transition o-H$_2$D$^+(1_{10}$-$1_{11})$, using a synthetic source, to understand how the physical characteristics are reflected in observations through interferometers and single-dish telescopes. In order to perform such tests, we post-processed a magneto-hydrodynamic simulation of a collapsing magnetized core using the radiative transfer code POLARIS. Using the resulting intensity distributions as input, we performed single-dish (APEX) and interferometric (ALMA) synthetic observations at different evolutionary times, always mimicking realistic configurations. Finally, column densities were derived to compare our simulations with real observations previously performed. Our derivations for o-H$_2$D$^+$ are in agreement with values reported in the literature, in the range of 10$^{\!10-11}$cm$^{\!-2}$ and 10$^{\!12-13}$cm$^{\!-2}$ for single-dish and interferometric measurements, respectively.

We present results from our 47-night imaging campaign of Comet 41P/Tuttle-Giacobini-Kresak conducted from Lowell Observatory between 2017 February 16 and July 2. Coma morphology revealed gas jets, whose appearance and motion as a function of time yielded the rotation period and other properties. All narrowband CN images exhibited either one or two jets; one jet appeared as a partial face-on spiral with clockwise rotation while the second jet evolved from a side-on corkscrew, through face-on, and finally corkscrew again, with only a slow evolution throughout the apparition due to progressive viewing geometry changes. A total of 78 period determinations were made over a 7-week interval, yielding a smooth and accelerating rotation period starting at 24 hr (March 21&22) and passing 48 hr on April 28. While this is by far the fastest rate of change ever measured for a comet nucleus, the torque required is readily within what can exist given likely properties of the nucleus. If the torque remained constant, we estimate that the nucleus could have stopped rotating and/or began to tumble as soon as only two months following perihelion, and will certainly reach this stage by early in the next apparition. Working backwards in time, Tuttle-Giacobini-Kresak would have been rotating near its rotational break-up velocity 3-4 orbits earlier, suggesting that its extreme 7-magnitude outburst observed in 2001 might have been caused by a partial fragmentation at that time, as might the pair of 1973 8-magnitude outbursts if there had been an earlier spin-down and spin-up cycle.

Context: In cold and dense gas prior to the formation of young stellar objects, heavy molecular species (including CO) are accreted onto dust grains. Under these conditions H$_3^+$ and its deuterated isotopologues become more abundant, enhancing the deuterium fraction of molecules such as N$_2$H$^+$ that are formed via ion-neutral reactions. Because this process is extremely temperature sensitive, the abundance of these species is likely linked to the evolutionary stage of the source. Aims: We investigate how the abundances of o-H$_2$D$^+$ and N$_2$D$^+$ vary with evolution in high-mass clumps. Methods: We observed with APEX the ground-state transitions of o-H$_2$D$^+$ near 372 GHz, and N$_2$D$^+$(3-2) near 231 GHz for three massive clumps in different evolutionary stages. The sources were selected within the G351.77-0.51 complex to minimise the variation of initial chemical conditions, and to remove distance effects. We modelled their dust continuum emission to estimate their physical properties, and also modelled their spectra under the assumption of local thermodynamic equilibrium to calculate beam-averaged abundances. Results: We find an anticorrelation between the abundance of o-H$_2$D$^+$ and that of N$_2$D$^+$, with the former decreasing and the latter increasing with evolution. With the new observations we are also able to provide a qualitative upper limit to the age of the youngest clump of about 10$^5$ yr, comparable to its current free-fall time. Conclusions: We can explain the evolution of the two tracers with simple considerations on the chemical formation paths, depletion of heavy elements, and evaporation from the grains. We therefore propose that the joint observation and the relative abundance of o-H$_2$D$^+$ and N$_2$D$^+$ can act as an efficient tracer of the evolutionary stages of the star-formation process.

We present a simple model for the Double Periodic Variable (DPV) V495 Cen, which evolves as a binary system of intermediate mass, where the gainer cannot accrete at high rate, limited by the Eddington accretion rate, leading to the formation of an accretion disc. The theoretical model begins at the zero age main sequence considering the rotation for both stars. For this purpose we used the stellar evolution code MESA, developed to calculate the evolution of stars in a wide range of parameters. We started the model adjusting fundamental parameters published for this system through a chi-square optimization algorithm, and adopting an initial orbital period of 3.9 days and initial masses for the primary component $M_{i,d} = 3.40$ $M_{\odot}$ and $M_{i,g}= 3.18$ $M_{\odot}$ for the gainer, with a metallicity associated to this type of DPV of $Z = 0.02$. The method converged successfully for eight free degrees and 5% of confidence with a chi-square of $\Delta \chi^{2}_{0.95,8}= 0.212$. We describe each evolutionary stage of both components until that the donor reaches 20% core helium depletion as stop criterion. We offer a complementary analysis for understanding the mechanism of the magnetic dynamos inside the donor star using the Tayler-Spruit formalism. Currently, the theoretical model is consistent with the fundamental parameters published for V495 Cen and we discuss how our predictions can help to develop efficients theoretical models for DPV stars.

The Event Horizon Telescope (EHT) provides a unique opportunity to probe the physics of supermassive black holes through Very Large Baseline Interferometry (VLBI), such as the existence of the event horizon, the accretion processes as well as jet formation. We build a theoretical model which includes an Advection Dominated Accretion Flow (ADAF) and a simple radio jet outflow. The predicted spectral energy distribution (SED) of this model is compared to observations to get the best estimates of the model parameters. Also the model-predicted radial emission profiles at different frequency bands can be used to predict whether the inflow can be resolved by the EHT.

Black holes with masses of $\rm 10^6-10^9~M_{\odot}$ dwell in the centers of most galaxies, but their formation mechanisms are not well known. A subdominant dissipative component of dark matter with similar properties to the ordinary baryons, known as mirror dark matter, may collapse to form massive black holes during the epoch of first galaxies formation. In this study, we explore the possibility of massive black hole formation via this alternative scenario. We perform three-dimensional cosmological simulations for four distinct halos and compare their thermal, chemical and dynamical evolution in both the ordinary and the mirror sectors. We find that the collapse of halos is significantly delayed in the mirror sector due to the lack of $\rm H_2$ cooling and only halos with masses above $ \rm \geq 10^7~ M_{\odot}$ are formed. Overall, the mass inflow rates are $\rm \geq 10^{-2}~M_{\odot}/yr$ and there is less fragmentation. This suggests that the conditions for the formation of massive objects, including black holes, are more favorable in the mirror sector.

To explain the observed population of supermassive black holes at z~7, very massive seed black holes or, alternatively, super-Eddington scenarios are needed to reach final masses of the order of 10^9 solar masses. A popular explanation for massive seeds has been the direct collapse model, which predicts the formation of a single massive object due to the direct collapse of a massive gas cloud. Simulations over the last years have however shown that such a scenario is very difficult to achieve. A realistic model of black hole formation should therefore take fragmentation into account, and consider the interaction between stellar-dynamical and gas-dynamical processes. We present here numerical simulations pursued with the AMUSE code, employing an approximate treatment of the gas. Based on these simulations, we show that very massive black holes of 10^4-10^5 solar masses may form depending on the gas supply and the accretion onto the protostars.

The presence of supermassive black holes (SMBHs) of a few billion solar masses at very high redshift has motivated us to study how these massive objects formed during the first billion years after the Big Bang. The most promising model that has been proposed to explain this is the direct collapse of protogalactic gas clouds. In this scenario, very high accretion rates are needed to form massive objects early on and the suppression of $\mathrm H_2$ cooling is important in regulating the fragmentation. Recent studies have shown that if we use a strong radiation background, the hydrogen molecules are destroyed, favoring the high accretion rates and therefore producing objects of very high mass. In this work we study the impact of UV radiation fields in a primordial gas cloud using the recently coupled code GRADSPH-KROME for the modeling of gravitational collapse including primordial chemistry to explore the fragmentation in AU scales and hence the formation of first SMBHs. We found that to suppress the formation of $\mathrm H_2$ a very high value of $J_{21}$ is required, because of that we conclude that the direct collapse black holes (DCBHs) are very unlikely to be an explanation for the formation of the first SMBHs.

DQ Velorum is a galactic double periodic variable (DPV), this system is a semi-detached binary comprised of a B-type gainer and an A-type donor star plus an extended accretion disc around the gainer. The system also presents an orbital period of $6.08337$~days and a long period of $189$~days whose origin is still under debate. Here we studied the possibility that this period may be driven by a magnetic dynamo investigating the entire evolution of the system. The model matches in a very good way the current state of the system and it can potentially be used to describe the evolution of DQ Velorum. It also predicts an increase of the dynamo number of the donor during epochs of high mass transfer in this system, and a theoretical long/orbital period ratio very close to the observed one at the present system age.

We present here a three-dimesional hydrodynamical simulation for star formation. Our aim is to explore the effect of the metal-line cooling on the thermodynamics of the star-formation process. We explore the effect of changing the metallicty of the gas from $Z/Z_{\odot}=10^{-4}$ to $Z/Z_{\odot}=10^{-2}$. Furthermore, we explore the implications of using the observational abundance pattern of a CEMP-no star, which have been considered to be the missing second-generation stars, the so-called Pop. III.2 stars. In order to pursue our aim, we modelled the microphysics by employing the public astrochemistry package KROME, using a chemical network which includes sixteen chemical species (H, H$^{+}$, H$^{-}$, He, He$^{+}$, He$^{++}$, e$^{-}$, H$_{2}$, H$_{2}^{+}$, C, C$^{+}$, O, O$^{+}$, Si, Si$^{+}$, and Si$^{++}$). We couple KROME with the fully three-dimensional Smoothed-particle hydrodynamics (SPH) code GRADSPH. With this framework we investigate the collapse of a metal-enhanced cloud, exploring the fragmentation process and the formation of stars. We found that the metallicity has a clear impact on the thermodynamics of the collapse, allowing the cloud to reach the CMB temperature floor for a metallicity $Z/Z_{\odot}=10^{-2}$, which is in agreement with previous work. Moreover, we found that adopting the abundance pattern given by the star SMSS J031300.36-670839.3 the thermodynamics behavior is very similar to simulations with a metallicity of $Z/Z_{\odot}=10^{-2}$, due to the high carbon abundance. As long as only metal line cooling is considered, our results support the metallicity threshold proposed by previous works, which will very likely regulate the first episode of fragmentation and potentially determine the masses of the resulting star clusters.

Since its proposal in 1992, the Applegate mechanism has been discussed as a potential intrinsical mechanism to explain transit timing variations in various kinds of close binary systems. Most analytical arguments presented so far focused on the energetic feasibility of the mechanism, while applying rather crude one- or two-zone prescriptions to describe the exchange of angular momentum within the star. In this paper, we present the most detailed approach to date to describe the physics giving rise to the modulation period from kinetic and magnetic fluctuations. Assuming moderate levels of stellar parameter fluctuations, we find that the resulting binary period variations are one or two orders of magnitude lower than the observed values in RS-CVn like systems, supporting the conclusion of existing theoretical work that the Applegate mechanism may not suffice to produce the observed variations in these systems. The most promising Applegate candidates are low-mass post-common-envelope binaries (PCEBs) with binary separations $\lesssim 1~\mathrm{R}_\odot$ and secondary masses in the range of $0.30~\mathrm{M}_\odot$ and $0.36~\mathrm{M}_\odot$.

In this chapter, we review the current observational status of the first supermassive black holes. It is clear that such a review can hardly be complete, due to the wealth of surveys that has been pursued, including different wavelengths and different observational techniques. This chapter will focus on the main results that have been obtained, considering the detections of z~6 supermassive black holes in large surveys such as SDSS, CFHQS and Pan-STARRS. In addition, we will discuss upper limits and constraints on the population of the first black holes that can be derived from observational data, in particular in the X-ray regime, as these provide additional relevant information for the comparison with formation scenarios.

We report on narrowband photometry and extensive imaging observations of comet C/Lulin (2007 N3) obtained at Lowell Observatory during 2008 and 2009. Enhanced CN images revealed a double corkscrew morphology with two near-polar jets oriented approximately east-west, and both CN and dust images showed nightly rotational variability and seasonal changes in bulk morphology. We determined a rotational pole direction of RA/Dec = 81/+29 with an obliquity of 97 deg, and a sidereal rotation period of 41.45 +/- 0.05 hr. Monte Carlo numerical modeling best replicated the observed CN features with an eastern source area at lat/long -80/125 and a 10 deg radius and a western source area at lat/long +77/245 and a 20 deg radius, 4 times larger than the eastern source. An additional small, near-equatorial source was necessary to reproduce some dust features. Water morphology, based on OH, was quite different than that of the carbon-bearing species, implying a different driver for the polar jets such as CO or CO2. Ion tails were detected in decontaminated images from both the dust and NH filters, likely being H2O+ and OH+, respectively. We measured water production both before and after perihelion, and extrapolate peak water production at perihelion to be about 1.0*10^29 molecules/s. We estimate an active fraction of only 4-5% and a nucleus radius of up to ~8 km. Our data suggest that Lulin, defined as dynamically new in a statistical sense, behaves more like a long-period comet due to its nearly asteroidal early appearance, isolated source regions, and dust properties.

V393 Scorpii is a member of the subclass of Algols dubbed Double Periodic Variables (DPVs). These are semidetached binaries with B-type primaries showing a long-photometric cycle lasting in average 33 times the orbital period. We describe the behavior of unreported metallic emission lines in the cool stellar component of this system. The emissions can be single or double for a same line and sometimes show velocity shifts regarding the velocity of the center of mass of the star. In addition, these lines are stronger during the high state. This behavior suggests the presence of active regions in the surface of the rapidly rotating A7 donor covering a fraction of the visible hemisphere, which have larger emissivity during the high state. Our finding supports the recently proposed dynamo model for the long cycle of DPVs proposed by Schleicher & Mennickent. The model predicts an increase of the dynamo number of the donor during epochs of mass transfer in this system, and a theoretical long/orbital period ratio very close to the observed one at the present system age.

Eclipsing time variations (ETVs) are observed in many close binary systems. In particular, for several post-common-envelope binaries (PCEBs) that consist of a white dwarf and a main sequence star, the O-C diagram suggests that real or apparent orbital period variations are driven by Jupiter-mass planets or as a result of magnetic activity, the so-called Applegate mechanism. The latter explains orbital period variations as a result of changes in the stellar quadrupole moment due to magnetic activity. We explore the feasibility of driving ETVs via the Applegate mechanism for a sample of PCEB systems, including a range of different rotations. Using the MESA code we evolve 12 stars with different masses and rotation rates. We apply a simple dynamo model to their radial profiles to investigate on which scale the predicted activity cycle matches the observed modulation period, and quantify the uncertainty, and further calculate the required energies to drive que Applegate mechanism. We show that the Applegate mechanism is energetically feasible in 5 PCEB systems, and note that these are the systems with the highest rotation rate compared to the critical rotation rate of the main-sequence star. The results suggest that the ratio of physical to critical rotation in the main sequence star is an important indicator for the feasibility of Applegate's mechanism, but exploring larger samples will be necessary to probe this hypothesis.

The first stars are known to form in primordial gas, either in minihalos with about $10^6$~M$_\odot$ or so-called atomic cooling halos of about $10^8$~M$_\odot$. Simulations have shown that gravitational collapse and disk formation in primordial gas yield dense stellar clusters. In this paper, we focus particularly on the formation of protostellar binary systems, and aim to quantify their properties during the early stage of their evolution. For this purpose, we combine the smoothed particle hydrodynamics code GRADSPH with the astrochemistry package KROME. The GRADSPH-KROME framework is employed to investigate the collapse of primordial clouds in the high-density regime, exploring the fragmentation process and the formation of binary systems. We observe a strong dependence of fragmentation on the strength of the turbulent Mach number $\mathcal{M}$ and the rotational support parameter $\beta{}$. Rotating clouds show significant fragmentation, and have produced several Pop.~III proto-binary systems. We report maximum and minimum mass accretion rates of $2.31 \times 10^{-1}$~M$_{\odot}$ yr$^{-1}$ and $2.18\times 10^{-4}$~M$_{\odot}$ yr$^{-1}$. The mass spectrum of the individual Pop III proto-binary components ranges from $0.88$~M$_{\odot}$ to $31.96$~M$_{\odot}$ and has a sensitive dependence on the Mach number $\mathcal{M}$ as well as on the rotational parameter $\beta{}$. We also report a range from $\sim0.01$ to $\sim1$ for the mass ratio of our proto-binary systems.

Although timing variations in close binary systems have been studied for a long time, their underlying causes are still unclear. A possible explanation is the so-called Applegate mechanism, where a strong, variable magnetic field can periodically change the gravitational quadrupole moment of a stellar component, thus causing observable period changes. One of the systems exhibiting such strong orbital variations is the RS CVn binary HR 1099, whose activity cycle has been studied by various authors via photospheric and chromospheric activity indicators, resulting in contradicting periods. We aim at independently determining the magnetic activity cycle of HR 1099 using archival X-ray data to allow for a comparison to orbital period variations. Archival X-ray data from 80 different observations of HR 1099 acquired with 12 different X-ray facilities and covering almost four decades were used to determine X-ray fluxes in the energy range of 2-10 keV via spectral fitting and flux conversion. Via the Lomb-Scargle periodogram we analyze the resulting long-term X-ray light curve to search for periodicities. We do not detect any statistically significant periodicities within the X-ray data. An analysis of optical data of HR 1099 shows that the derivation of such periods is strongly dependent on the time coverage of available data, since the observed optical variations strongly deviate from a pure sine wave. We argue that this offers an explanation as to why other authors derive such a wide range of activity cycle periods based on optical data. We conclude that our analysis constitutes the longest stellar X-ray activity light curve acquired to date, yet the still rather sparse sampling of the X-ray data, along with stochastic flaring activity, does not allow for the independent determination of an X-ray activity cycle.

The origin of very low-mass stars (VLMS) and brown dwarfs (BDs) is still an unresolved topic of star formation. We here present numerical simulations of the formation of VLMS, BDs, and planet mass objects (planemos) resulting from the gravitational collapse and fragmentation of solar mass molecular cores with varying rotation rates and initial density perturbations. Our simulations yield various types of binary systems including the combinations VLMS-VLMS, BD-BD, planemo-planemo, VLMS-BD, VLMS-planemos, BD-planemo. Our scheme successfully addresses the formation of wide VLMS and BD binaries with semi-major axis up to 441 AU and produces a spectrum of mass ratios closer to the observed mass ratio distribution (q > 0.5). Molecular cores with moderate values of the ratio of kinetic to gravitational potential energy (0.16 <= beta <= 0.21) produce planemos. Solar mass cores with rotational parameters beta outside of this range yield either VLMS/BDs or a combination of both. With regard to the mass ratios we find that for both types of binary systems the mass ratio distribution varies in the range 0.31 <= q <= 0.74. We note that in the presence of radiative feedback, the length scale of fragmentation would increase by approximately two orders of magnitude, implying that the formation of binaries may be efficient for wide orbits, while being suppressed for short-orbit systems.

We use the POLARIS radiative transport code to generate predictions of the two main observables directly sensitive to the magnetic field morphology and strength in filaments: dust polarization and gas Zeeman line splitting. We simulate generic gas filaments with power-law density profiles assuming two density-field strength dependencies, six different filament inclinations, and nine distinct magnetic field morphologies, including helical, toroidal, and warped magnetic field geometries. We present idealized spatially resolved dust polarization and Zeeman-derived field strengths and directions maps. Under the assumption that dust grains are aligned by radiative torques (RATs), dust polarization traces the projected plane-of-the-sky magnetic field morphology. Zeeman line splitting delivers simultaneously the intensity-weighted line-of-sight field strength and direction. We show that linear dust polarization alone is unable to uniquely constrain the 3D field morphology. We demonstrate that these ambiguities are ameliorated or resolved with the addition of the Zeeman directional information. Thus, observations of both the dust polarization and Zeeman splitting together provide the most promising means for obtaining constraints of the 3D magnetic field configuration. We find that the Zeeman-derived field strengths are at least a factor of a few below the input field strengths due to line-of-sight averaging through the filament density gradient. Future observations of both dust polarization and Zeeman splitting are essential for gaining insights into the role of magnetic fields in star and cluster forming filaments.

We have observed 104 gravitationally-lensed quasars at $z\sim1-4$ with Herschel/SPIRE, the largest such sample ever studied. By targeting gravitational lenses, we probe intrinsic far-infrared (FIR) luminosities and star formation rates (SFRs) more typical of the population than the extremely luminous sources that are otherwise accessible. We detect 72 objects with Herschel/SPIRE and find 66 percent (69 sources) of the sample have spectral energy distributions (SEDs) characteristic of dust emission. For 53 objects with sufficiently constrained SEDs, we find a median effective dust temperature of $38^{+12}_{-5}$ K. By applying the radio-infrared correlation, we find no evidence for an FIR excess which is consistent with star-formation-heated dust. We derive a median magnification-corrected FIR luminosity of $3.6^{+4.8}_{-2.4}~\times 10^{11}~{\rm L_{\odot}}$ and median SFR of $120^{+160}_{-80}~{\rm M_{\odot}~yr^{-1}}$ for 94 quasars with redshifts. We find $\sim10$ percent of our sample have FIR properties similar to typical dusty star-forming galaxies at $z\sim2-3$ and a range of SFRs $<20-10000~{\rm M_{\odot}~yr^{-1}}$ for our sample as a whole. These results are in line with current models of quasar evolution and suggests a coexistence of dust-obscured star formation and AGN activity is typical of most quasars. We do not find a statistically-significant difference in the FIR luminosities of quasars in our sample with a radio excess relative to the radio-infrared correlation. Synchrotron emission is found to dominate at FIR wavelengths for $<15$ percent of those sources classified as powerful radio galaxies.

Dark matter annihilation can have a strong impact on many astrophysical processes in the Universe. In the case of Sommerfeld-enhanced annihilation cross sections, the annihilation rates are enhanced at late times, thus enhancing the potential annihilation signatures. We here calculate the Sommerfeld-enhanced annihilation signatures during the epoch of helium reionization, the epoch where helium becomes fully ionized due to energetic photons. When considering the upper limits on the energy injection from the CMB, we find that the resulting abundance of He$^{++}$ becomes independent of the dark matter particle mass. The resulting enhancement compared to a standard scenario is thus 1-2 orders of magnitude higher. For realistic scenarios compatible with CMB constraints, there is no significant shift in the epoch of helium reionization, which is completed between redshifts $3$ and $4$. While it is thus difficult to disentangle dark matter annihilation from astrophysical contributions (active galactic nuclei), a potential detection of dark matter particles and its interactions using the Large Hadron Collider (LHC) would allow one to quantify the dark matter contribution.

Collisions were suggested to potentially play a role in the formation of massive stars in present day clusters, and have likely been relevant during the formation of massive stars and intermediate mass black holes within the first star clusters. In the early Universe, the first stellar clusters were particularly dense, as fragmentation typically only occurred at densities above $10^9$cm$^{-3}$, and the radii of the protostars were enhanced due to the larger accretion rates, suggesting a potentially more relevant role of stellar collisions. We present here a detailed parameter study to assess how the number of collisions as well as the mass growth of the most massive object depends on the properties of the cluster, and we characterize the time evolution with three effective parameters, the time when most collisions occur, the duration of the collisions period, as well as the normalization required to obtain the total number of collisions. We apply our results to typical Population III (Pop.III) clusters of about $1000$M$_\odot$, finding that a moderate enhancement of the mass of the most massive star by a factor of a few can be expected. For more massive Pop.III clusters as expected in the first atomic cooling halos, we expect a more significant enhancement by a factor of $15-32$. We therefore conclude that collisions in massive Pop.III clusters were likely relevant to form the first intermediate mass black holes.

Models aiming to explain the formation of massive black hole seeds, and in particular the direct collapse scenario, face substantial difficulties. These are rooted in rather ad hoc and fine-tuned initial conditions, such as the simultaneous requirements of extremely low metallicities and strong radiation backgrounds. Here we explore a modification of such scenarios where a massive primordial star cluster is initially produced. Subsequent stellar collisions give rise to the formation of massive (10^4 - 10^5 solar mass) objects. Our calculations demonstrate that the interplay between stellar dynamics, gas accretion and protostellar evolution is particularly relevant. Gas accretion onto the protostars enhances their radii, resulting in an enhanced collisional cross section. We show that the fraction of collisions can increase from 0.1-1% of the initial population to about 10% when compared to gas-free models or models of protostellar clusters in the local Universe. We conclude that very massive objects can form in spite of initial fragmentation, making the first massive protostellar clusters viable candidate birth places for observed supermassive black holes.

By using a novel interface between the modern smoothed particle hydrodynamics code GASOLINE2 and the chemistry package KROME, we follow the hydrodynamical and chemical evolution of an isolated galaxy. In order to assess the relevance of different physical parameters and prescriptions, we constructed a suite of ten simulations, in which we vary the chemical network (primordial and metal species), how metal cooling is modelled (non-equilibrium versus equilibrium; optically thin versus thick approximation), the initial gas metallicity (from ten to hundred per cent solar), and how molecular hydrogen forms on dust. This is the first work in which metal injection from supernovae, turbulent metal diffusion, and a metal network with non-equilibrium metal cooling are self-consistently included in a galaxy simulation. We find that properly modelling the chemical evolution of several metal species and the corresponding non-equilibrium metal cooling has important effects on the thermodynamics of the gas, the chemical abundances, and the appearance of the galaxy: the gas is typically warmer, has a larger molecular gas mass fraction, and has a smoother disc. We also conclude that, at relatively high metallicity, the choice of molecular-hydrogen formation rates on dust is not crucial. Moreover, we confirm that a higher initial metallicity produces a colder gas and a larger fraction of molecular gas, with the low-metallicity simulation best matching the observed molecular Kennicutt-Schmidt relation. Finally, our simulations agree quite well with observations which link star formation rate to metal emission lines.

Observations show a large spread in the luminosities of young protostars, which are frequently explained in the context of episodic accretion. We here test this scenario using numerical simulations following the collapse of a solar mass molecular cloud using the GRADSPH code, varying the strength of the initial perturbations and the temperature of the cores. A specific emphasis of this paper is to investigate the role of binaries and multiple systems in the context of episodic accretion, and to compare their evolution to the evolution in isolated fragments. Our models form a variety of low mass protostellar objects including single, binary and triple systems with binaries more active in exhibiting episodic accretion than isolated protostars. We also find a general decreasing trend for the average mass accretion rate over time, suggesting that the majority of the protostellar mass is accreted within the first 10^5 years. This result can potentially help to explain the surprisingly low average luminosities in the majority of the protostellar population.

The Orion nebula is a prime example of a massive star-forming region in our galaxy. Observations have shown that gravitational and magnetic energy are comparable in its integral shaped filament (ISF) on a scale of ~1 pc, and that the population of pre-main sequence stars appears dynamically heated compared to the protostars. These results have been attributed to a slingshot mechanism resulting from the oscillation of the filament (Stutz & Gould 2016). In this paper, we show that radially contracting filaments naturally evolve toward a state where gravitational, magnetic, and rotational energy are comparable. While the contraction of the filament will preferentially amplify the axial component of the magnetic field, the presence of rotation leads to a helical field structure. We show how magnetic tension can give rise to a filament oscillation, and estimate a typical timescale of 0.7 million years for the motion of the filament to the position of maximum displacement, consistent with the characteristic timescale of the ejected stars. Furthermore, the presence of helical magnetic fields is expected to give rise to magneto-hydrodynamical instabilities. We show here that the presence of a magnetic field significantly enhances the overall instability, which operates on a characteristic scale of about 1 pc. We expect the physics discussed here to be generally relevant in massive star forming regions, and encourage further investigations in the future.

Deuterium fractionation, i.e. the enhancement of deuterated species with respect to the non-deuterated ones, is considered to be a reliable chemical clock of star-forming regions. This process is strongly affected by the ortho-to-para (o-p) H$_2$ ratio. In this letter we explore the effect of the o-p H$_2$ conversion on grains on the deuteration timescale in fully depleted dense cores, including the most relevant uncertainties that affect this complex process. We show that (i) the o-p H$_2$ conversion on grains is not strongly influenced by the uncertainties on the conversion time and the sticking coefficient and (ii) that the process is controlled by the temperature and the residence time of ortho-H$_2$ on the surface, i.e. by the binding energy. We find that for binding energies in between 330-550 K, depending on the temperature, the o-p H$_2$ conversion on grains can shorten the deuterium fractionation timescale by orders of magnitude, opening a new route to explain the large observed deuteration fraction $D_\mathrm{frac}$ in dense molecular cloud cores. Our results suggest that the star formation timescale, when estimated through the timescale to reach the observed deuteration fractions, might be shorter than previously proposed. However, more accurate measurements of the binding energy are needed to better assess the overall role of this process.

We present observations of comet 67P/Churyumov-Gerasimenko acquired in support of the $Rosetta$ mission. We obtained usable data on 68 nights from 2014 September until 2016 May, with data acquired regularly whenever the comet was observable. We collected an extensive set of near-IR $J$, $H$, and $Ks$ data throughout the apparition plus visible-light images in $g'$, $r'$, $i'$, and $z'$ when the comet was fainter. We also obtained broadband $R$ and narrowband $CN$ filter observations when the comet was brightest using telescopes at Lowell Observatory. The appearance was dominated by a central condensation and the tail until 2015 June. From 2015 August onwards there were clear asymmetries in the coma, which enhancements revealed to be due to the presence of up to three features (i.e., jets). The features were similar in all broadband filters; $CN$ images did not show these features but were instead broadly enhanced in the southeastern hemisphere. Modeling using the parameters from Vincent et al. (2013) replicated the dust morphology reasonably well, indicating that the pole orientation and locations of active areas have been relatively unchanged over at least the last three apparitions. The dust production, as measured by $A(0^{\circ})f{\rho}$ peaked $\sim$30 days after perihelion and was consistent with predictions from previous apparitions. $A(0^{\circ})f{\rho}$ as a function of heliocentric distance was well fit by a power-law with slope $-$4.2 from 35-120 days post-perihelion. We detected photometric evidence of apparent outbursts on 2015 August 22 and 2015 September 19, although neither was discernible morphologically in this dataset.

We analyzed images of comet 49P/Arend-Rigaux on 33 nights between 2012 January and May and obtained R-band lightcurves of the nucleus. Through usual phasing of the data we found a double-peaked lightcurve having a synodic rotation period of 13.450 +/- 0.005 hr. Similarly, phase dispersion minimization and the Lomb-Scargle method both revealed rotation periods of 13.452 hr. Throughout the 2011/12 apparition, the rotation period was found to increase by a small amount, consistent with a retrograde rotation of the nucleus. We also reanalyzed the publicly available data from the 1984/85 apparition by applying the same techniques, finding a rotation period of 13.45 +/- 0.01 hr. Based on these findings we show that the change in rotation period is less than 14 seconds per apparition. Furthermore, the amplitudes of the light curves from the two apparitions are comparable, to within reasonable errors, even though the viewing geometries differ, implying that we are seeing the comet at a similar sub-Earth latitude. We detected the presence of a short term jet-like feature in 2012 March which appears to have been created by a short duration burst of activity on March 15. Production rates obtained in 2004/05, along with reanalysis of previous results from 1984/85 imply a strong seasonal effect and a very steep fall-off after perihelion. This, in turn, implies that a single source region dominates activity, rather than leakage from the entire nucleus.

Deuterium fractionation is considered as an important process to infer the chemical ages of prestellar cores in filaments. We present here the first magneto-hydrodynamical simulations including a chemical network to study deuterium fractionation in magnetized and turbulent filaments and their substructures. The filaments typically show widespread deuterium fractionation with average values $\gtrsim0.01$. For individual cores of similar age, we observe the deuteration fraction to increase with time, but also to be independent of their average properties such as density, virial or mass-to-magnetic flux ratio. We further find a correlation of the deuteration fraction with core mass, average H$_2$ density and virial parameter only at late evolutionary stages of the filament and attribute this to the lifetime of the individual cores. Specifically, chemically old cores reveal higher deuteration fractions. Within the radial profiles of selected cores, we notice differences in the structure of the deuteration fraction or surface density, which we can attribute to their different turbulent properties. High deuteration fractions of the order $0.01-0.1$ may be reached within approximately $200$~kyrs, corresponding to two free-fall times, as defined for cylindrical systems, of the filaments

Astrochemistry plays a central role during the process of star formation, both in the primordial regime as well as in the present-day Universe. We revisit here the chemistry in both regimes, focusing first on the chemistry under close to primordial conditions, as observed in the so-called Caffau star SDSS J102915+172927, and subsequently discuss deuteration processes in present-day star-forming cores. In models of the high-redshift Universe, the chemistry is particularly relevant to determine the cooling, while it also serves as an important diagnostic in the case of present-day star formation.

Observations of the Orion-A integral shaped filament (ISF) have shown indications of an oscillatory motion of the gas filament. This evidence is based on both the wave-like morphology of the filament as well as the kinematics of the gas and stars, where the characteristic velocities of the stars require a dynamical heating mechanism. As proposed by Stutz and Gould (2016), such a heating mechanism (the "Slingshot") may be the result of an oscillating gas filament in a gas-dominated (as opposed to stellar-mass dominated) system. Here we test this hypothesis with the first stellar-dynamical simulations in which the stars are subjected to the influence of an oscillating cylindrical potential. The accelerating, cylindrical background potential is populated with a narrow distribution of stars. By coupling the potential to N-body dynamics, we are able to measure the influence of the potential on the stellar distribution. The simulations provide evidence that the slingshot mechanism can successfully reproduce several stringent observational constraints. These include the stellar spread (both in projected position and in velocity) around the filament, the symmetry in these distributions, and a bulk motion of the stars with respect to the filament. Using simple considerations we show that star-star interactions are incapable of reproducing these spreads on their own when properly accounting for the gas potential. Thus, properly accounting for the gas potential is essential for understanding the dynamical evolution of star forming filamentary systems in the era of Gaia.

We report the peculiar chemical abundance patterns of eleven atypical Milky Way (MW) field red giant stars observed by the Apache Point Observatory Galactic Evolution Experiment (APOGEE). These atypical giants exhibit strong Al and N enhancements accompanied by C and Mg depletions, strikingly similar to those observed in the so-called second-generation (SG) stars of globular clusters (GCs). Remarkably, we find low-Mg abundances ([Mg/Fe]$<$0.0) together with strong Al and N overabundances in the majority (5/7) of the metal-rich ([Fe/H]$\gtrsim - 1.0$) sample stars, which is at odds with actual observations of SG stars in Galactic CGs of similar metallicities. This chemical pattern is unique and unprecedented among MW stars, posing urgent questions about its origin. These atypical stars could be former SG stars of dissolved GCs formed with intrinsically lower abundances of Mg and enriched Al (subsequently self-polluted by massive AGB stars) or the result of exotic binary systems. We speculate that the stars Mg-deficiency as well as the orbital properties suggest that they could have an extragalactic origin. This discovery should guide future dedicated spectroscopic searches of atypical stellar chemical patterns in our Galaxy; a fundamental step forward to understand the Galactic formation and evolution.

We present a summary of the campaign of remote observations that supported the European Space Agency's Rosetta mission. Telescopes across the globe (and in space) followed comet 67P/Churyumov-Gerasimenko from before Rosetta's arrival until nearly the end of mission in September 2016. These provided essential data for mission planning, large-scale context information for the coma and tails beyond the spacecraft, and a way to directly compare 67P with other comets. The observations revealed 67P to be a relatively `well behaved' comet, typical of Jupiter family comets and with activity patterns that repeat from orbit-to-orbit. Comparison between this large collection of telescopic observations and the in situ results from Rosetta will allow us to better understand comet coma chemistry and structure. This work is just beginning as the mission ends -- in this paper we present a summary of the ground-based observations and early results, and point to many questions that will be addressed in future studies.

While supermassive black holes are known to co-evolve with their host galaxy, the precise nature and origin of this co-evolution is not clear. We here explore the possible connection between star formation and black hole growth in the circumnuclear disk (CND) to probe this connection in the vicinity close to the black hole. We adopt here the circumnuclear disk model developed by Kawakatu & Wada (2008) and Wutschik et al. (2013), and explore both the dependence on the star formation recipe as well as the role of the gravitational field, which can be dominated by the central black hole, the CND itself or the host galaxy. A specific emphasis is put on the turbulence regulated star formation model by Krumholz et al. (2005) to explore the impact of a realistic star formation recipe. It is shown that this model helps to introduce realistic fluctuations in the black hole and star formation rate, without overestimating them. Consistent with previous works, we show that the final black hole masses are rather insensitive to the masses of the initial seeds, even for seed masses of up to 10^6 M_sol. In addition, we apply our model to the formation of high-redshift quasars, as well as to the nearby system NGC 6951, where a tentative comparison is made in spite of the presence of a bar in the galaxy. We show that our model can reproduce the high black hole masses of the high-redshift quasars within a sufficiently short time, provided a high mass supply rate from the host galaxy. In addition, it reproduces several of the properties observed in NGC 6951. With respect to the latter system, our analysis suggests that supernova feedback may be important to create the observed fluctuations in the star formation history as a result of negative feedback effects.

The class of Double Period Variables (DPVs) consists of close interacting binaries, with a characteristic long period that is an order of magnitude longer than the corresponding orbital period, many of them with a characteristic ratio of about 3.5x10^1. We consider here the possibility that the accretion flow is modulated as a result of a magnetic dynamo cycle. Due to the short binary separations, we expect the rotation of the donor star to be synchronized with the rotation of the binary due to tidal locking. We here present a model to estimate the dynamo number and the resulting relation between the activity cycle length and the orbital period, as well as an estimate for the modulation of the mass transfer rate. The latter is based on Applegate's scenario, implying cyclic changes in the radius of the donor star and thus in the mass transfer rate as a result of magnetic activity. Our model is applied to a sample of 17 systems with known physical parameters, 10 also with known orbital periods. In spite of the uncertainties of our simplified framework, the results show a reasonable agreement, indicating that a dynamo interpretation is potentially feasible. At the same time, we note that the orbital period variations resulting from Applegate's model are sufficiently small to be consistent with the data. We conclude that both larger samples with known physical parameters as well as potential direct probes of the magnetism of the donor star, including cold spots as well as polarization, will be valuable to further constrain the nature of these systems.

High-mass stars are expected to form from dense prestellar cores. Their precise formation conditions are widely discussed, including their virial condition, which results in slow collapse for super-virial cores with strong support by turbulence or magnetic fields, or fast collapse for sub-virial sources. To disentangle their formation processes, measurements of the deuterium fractions are frequently employed to approximately estimate the ages of these cores and to obtain constraints on their dynamical evolution. We here present 3D magneto-hydrodynamical simulations including for the first time an accurate non-equilibrium chemical network with 21 gas-phase species plus dust grains and 213 reactions. With this network we model the deuteration process in fully depleted prestellar cores in great detail and determine its response to variations in the initial conditions. We explore the dependence on the initial gas column density, the turbulent Mach number, the mass-to-magnetic flux ratio and the distribution of the magnetic field, as well as the initial ortho-to-para ratio of H2. We find excellent agreement with recent observations of deuterium fractions in quiescent sources. Our results show that deuteration is rather efficient, even when assuming a conservative ortho-to-para ratio of 3 and highly sub-virial initial conditions, leading to large deuterium fractions already within roughly a free-fall time. We discuss the implications of our results and give an outlook to relevant future investigations.

Large eddy simulations (LES) are a powerful tool in understanding processes that are inaccessible by direct simulations due to their complexity, for example, in the highly turbulent regime. However, their accuracy and success depends on a proper subgrid-scale (SGS) model that accounts for the unresolved scales in the simulation. We evaluate the applicability of two traditional SGS models, namely the eddy-viscosity (EV) and the scale-similarity (SS) model, and one recently proposed nonlinear (NL) SGS model in the realm of compressible MHD turbulence. Using 209 simulations of decaying, supersonic (initial sonic Mach number of ~3) MHD turbulence with a shock-capturing scheme and varying resolution, SGS model and filter, we analyze the ensemble statistics of kinetic and magnetic energy spectra and structure functions. Furthermore, we compare the temporal evolution of lower and higher order statistical moments of the spatial distributions of kinetic and magnetic energy, vorticity, current density, and dilatation magnitudes. We find no statistical influence on the evolution of the flow by any model if grid-scale quantities are used to calculate SGS contributions. In addition, the SS models, which employ an explicit filter, have no impact in general. On the contrary, both EV and NL models change the statistics if an explicit filter is used. For example, they slightly increase the dissipation on the smallest scales. We demonstrate that the nonlinear model improves higher order statistics already with a small explicit filter, i.e. a three-point stencil. The results of e.g. the structure functions or the skewness and kurtosis of the current density distribution are closer to the ones obtained from simulations at higher resolution. We conclude that the nonlinear model with a small explicit filter is suitable for application in more complex scenarios when higher order statistics are important.

A key for understanding the evolution of galaxies and in particular their star formation history will be future ultra-deep radio surveys. While star formation rates (SFRs) are regularly estimated with phenomenological formulas based on the local FIR-radio correlation, we present here a physically motivated model to relate star formation with radio fluxes. Such a relation holds only in frequency ranges where the flux is dominated by synchrotron emission, as this radiation originates from cosmic rays produced in supernova remnants, therefore reflecting recent star formation. At low frequencies synchrotron emission can be absorbed by the free-free mechanism. This suppression becomes stronger with increasing number density of the gas, more precisely of the free electrons. We estimate the critical observing frequency below which radio emission is not tracing the SFR, and use the three well-studied local galaxies M 51, M 82, and Arp 220 as test cases for our model. If the observed galaxy is at high redshift, this critical frequency moves along with other spectral features to lower values in the observing frame. In the absence of systematic evolutionary effects, one would therefore expect that the method can be applied at lower observing frequencies for high redshift observations. However, in case of a strong increase of the typical gas column densities towards high redshift, the increasing free-free absorption may erase the star formation signatures at low frequencies. At high radio frequencies both, free-free emission and the thermal bump, can dominate the spectrum, also limiting the applicability of this method.

Dust plays a key role in the evolution of the ISM and its correct modelling in numerical simulations is therefore fundamental. We present a new and self-consistent model that treats grain thermal coupling with the gas, radiation balance, and surface chemistry for molecular hydrogen. This method can be applied to any dust distribution with an arbitrary number of grain types without affecting the overall computational cost. In this paper we describe in detail the physics and the algorithm behind our approach, and in order to test the methodology, we present some examples of astrophysical interest, namely (i) a one-zone collapse with complete gas chemistry and thermochemical processes, (ii) a 3D model of a low-metallicity collapse of a minihalo starting from cosmological initial conditions, and (iii) a turbulent molecular cloud with H-C-O chemistry (277 reactions), together with self-consistent cooling and heating solved on the fly. Although these examples employ the publicly available code KROME, our approach can be easily integrated into any computational framework.

Understanding the formation of the extremely metal poor star SDSS-J102915+172927 is of fundamental importance to improve our knowledge on the transition between the first and second generation of stars in the Universe. In this paper, we perform three-dimensional cosmological hydrodynamical simulations of dust-enriched halos during the early stages of the collapse process including a detailed treatment of the dust physics. We employ the astrochemistry package \krome coupled with the hydrodynamical code \textsc{enzo} assuming grain size distributions produced by the explosion of core-collapse supernovae of 20 and 35 M$_\odot$ primordial stars which are suitable to reproduce the chemical pattern of the SDSS-J102915+172927 star. We find that the dust mass yield produced from Population III supernovae explosions is the most important factor which drives the thermal evolution and the dynamical properties of the halos. Hence, for the specific distributions relevant in this context, the composition, the dust optical properties, and the size-range have only minor effects on the results due to similar cooling functions. We also show that the critical dust mass to enable fragmentation provided by semi-analytical models should be revised, as we obtain values one order of magnitude larger. This determines the transition from disk fragmentation to a more filamentary fragmentation mode, and suggests that likely more than one single supernova event or efficient dust growth should be invoked to get such a high dust content.

Galaxy mergers are expected to play a central role for the evolution of galaxies, and may have a strong impact on their magnetic fields. We present the first grid-based 3D magneto-hydrodynamical simulations investigating the evolution of magnetic fields during merger events. For this purpose, we employ a simplified model considering the merger event of magnetized gaseous disks in the absence of stellar feedback and without a stellar or dark matter component. We show that our model naturally leads to the production of two peaks in the evolution of the average magnetic field strength within 5 kpc, within 25 kpc and on scales in between 5 and 25 kpc. The latter is consistent with the peak in the magnetic field strength reported by Drzazga et al. (2011) in a merger sequence of observed galaxies. We show that the peak on the galactic scale and in the outer regions is likely due to geometrical effects, as the core of one galaxy enters the outskirts of the other one. In addition, there is a physical enhancement of the magnetic field within the central ~5 kpc, reflecting the enhancement in density due to efficient angular momentum transport. We conclude that high-resolution observations of the central regions will be particularly relevant to probe the evolution of magnetic field structures during merger events.

The far-infrared - radio correlation connects star formation and magnetic fields in galaxies, and has been confirmed over a large range of far-infrared luminosities. Recent investigations indicate that it may even hold in the regime of local dwarf galaxies, and we explore here the expected behavior in the regime of star formation surface densities below 0.1 M_sun kpc^{-2} yr^{-1}. We derive two conditions that can be particularly relevant for inducing a change in the expected correlation: a critical star formation surface density to maintain the correlation between star formation rate and the magnetic field, and a critical star formation surface density below which cosmic ray diffusion losses dominate over their injection via supernova explosions. For rotation periods shorter than 1.5x10^7 (H/kpc)^2 yrs, with H the scale height of the disk, the first correlation will break down before diffusion losses are relevant, as higher star formation rates are required to maintain the correlation between star formation rate and magnetic field strength. For high star formation surface densities Sigma_SFR, we derive a characteristic scaling of the non-thermal radio to the far-infrared / infrared emission with Sigma_SFR^{1/3}, corresponding to a scaling of the non-thermal radio luminosity L_s with the infrared luminosity L_{th} as L_{th}^{4/3}. The latter is expected to change when the above processes are no longer steadily maintained. In the regime of long rotation periods, we expect a transition towards a steeper scaling with Sigma_SFR^{2/3}, implying L_s~L_th^{5/3}, while the regime of fast rotation is expected to show a considerably enhanced scatter. These scaling relations explain the increasing thermal fraction of the radio emission observed within local dwarfs, and can be tested with future observations by the SKA and its precursor radio telescopes.

Galactic magnetic fields in the local Universe are strong and omnipresent. There is mounting evidence that galaxies were magnetized already in the early Universe. Theoretical scenarios including the turbulent small-scale dynamo predict magnetic energy densities comparable to the one of turbulence. Based on the assumption of this energy equipartition, we determine the galactic synchrotron flux as a function of redshift z. The conditions in the early Universe are different from the present day, in particular the galaxies have more intense star formation. To cover a large range of conditions we consider models based on two different systems: one model galaxy comparable to the Milky Way and one typical high-z starburst galaxy. We include a model of the steady state cosmic ray spectrum and find that synchrotron emission can be detected up to cosmological redshifts with current and future radio telescopes. Turbulent dynamo theory is in agreement with the origin of the observed correlation between the far-infrared (FIR) luminosity L_FIR and the radio luminosity L_radio. Our model reproduces this correlation well at z=0. We extrapolate the FIR-radio correlation to higher redshift and predict a time evolution with a significant deviation from its present-day appearance already at z~2 for a gas density that increases strongly with z. In particular, we predict a decrease of the radio luminosity with redshift which is caused by the increase of cosmic ray energy losses at high z. The result is an increase of the ratio between L_FIR and L_radio. Simultaneously, we predict that the slope of the FIR-radio correlation becomes shallower with redshift. This behavior of the correlation could be observed in the near future with ultra-deep radio surveys.

Compressible magnetohydrodynamic (MHD) turbulence is ubiquitous in astrophysical phenomena ranging from the intergalactic to the stellar scales. In studying them, numerical simulations are nearly inescapable, due to the large degree of nonlinearity involved. However the dynamical ranges of these phenomena are much larger than what is computationally accessible. In large eddy simulations (LES), the resulting limited resolution effects are addressed explicitly by introducing to the equations of motion additional terms associated with the unresolved, subgrid-scale (SGS) dynamics. This renders the system unclosed. We derive a set of nonlinear structural closures for the ideal MHD LES equations with particular emphasis on the effects of compressibility. The closures are based on a gradient expansion of the finite-resolution operator (W.K. Yeo CUP 1993, ed. Galperin & Orszag) and require no assumptions about the nature of the flow or magnetic field. Thus the scope of their applicability ranges from the sub- to the hyper-sonic and -Alfvenic regimes. The closures support spectral energy cascades both up and down-scale, as well as direct transfer between kinetic and magnetic resolved and unresolved energy budgets. They implicitly take into account the local geometry, and in particular the anisotropy, of the flow. Their properties are $\textit{a priori}$ validated in an accompanying article (Grete et al. Phys. Plasmas, 2016) against alternative closures available in the literature with respect to a wide range of simulation data of homogeneous and isotropic turbulence.

Even though compressible plasma turbulence is encountered in many astrophysical phenomena, its effect is often not well understood. Furthermore, direct numerical simulations are typically not able to reach the extreme parameters of these processes. For this reason, large-eddy simulations (LES), which only simulate large and intermediate scales directly, are employed. The smallest, unresolved scales and the interactions between small and large scales are introduced by means of a subgrid-scale (SGS) model. We propose and verify a new set of nonlinear SGS closures for future application as an SGS model in LES of compressible magnetohydrodynamics (MHD). We use 15 simulations (without explicit SGS model) of forced, isotropic, homogeneous turbulence with varying sonic Mach number $\mathrm{M_s} = 0.2$ to $20$ as reference data for the most extensive \textit{a priori} tests performed so far in literature. In these tests we explicitly filter the reference data and compare the performance of the new closures against the most widely tested closures. These include eddy-viscosity and scale-similarity type closures with different normalizations. Performance indicators are correlations with the turbulent energy and cross-helicty flux, the average SGS dissipation, the topological structure and the ability to reproduce the correct magnitude and direction of the SGS vectors. We find that only the new nonlinear closures exhibit consistently high correlations (median value \textgreater$0.8$) with the data over the entire parameter space and outperform the other closures in all tests. Moreover, we show that these results are independent of resolution and chosen filter scale. Additionally, the new closures are effectively coefficient-free with a deviation of less than $20\%$.

We present results from our observing campaign of Comet 209P/LINEAR during its exceptionally close approach to Earth during May of 2014, the third smallest perigee of any comet in two centuries. These circumstances permitted us to pursue several studies of this intrinsically faint object, including measurements of gas and dust production rates, searching for coma morphology, and direct detection of the nucleus to measure its properties. Indeed, we successfully measured the lowest water production rates of an intact comet in over 35 years and a corresponding smallest active area, ~0.007 km^2. When combined with the nucleus size found from radar (Howell et al. 2014), this also yields the smallest active fraction for any comet, ~0.024%. In all, this strongly suggests that 209P/LINEAR is on its way to becoming an inert object. The nucleus was detected but could not easily be disentangled from the inner coma due to seeing variations and changing spatial scales. Even so, we were able to measure a double-peaked lightcurve consistent with the shorter of two viable rotational periods found by Hergenrother (2014). Radial profiles of the dust coma are quite steep, similar to that observed for some other very anemic comets and suggestive that vaporizing icy grains are present.

Observations of quasars at $ z > 6$ suggest the presence of black holes with a few times $\rm 10^9 ~M_{\odot}$. Numerous models have been proposed to explain their existence including the direct collapse which provides massive seeds of $\rm 10^5~M_{\odot}$. The isothermal direct collapse requires a strong Lyman-Werner flux to quench $\rm H_2$ formation in massive primordial halos. In this study, we explore the impact of trace amounts of metals and dust enrichment. We perform three dimensional cosmological simulations for two halos of $\rm > 10^7~M_{\odot}$ with $\rm Z/Z_{\odot}= 10^{-4}-10^{-6}$ illuminated by an intense Lyman Werner flux of $\rm J_{21}=10^5$. Our results show that initially the collapse proceeds isothermally with $\rm T \sim 8000$ K but dust cooling becomes effective at densities of $\rm 10^{8}-10^{12} ~cm^{-3}$ and brings the gas temperature down to a few 100-1000 K for $\rm Z/Z_{\odot} \geq 10^{-6}$. No gravitationally bound clumps are found in $\rm Z/Z_{\odot} \leq 10^{-5}$ cases by the end of our simulations in contrast to the case with $\rm Z/Z_{\odot} = 10^{-4}$. Large inflow rates of $\rm \geq 0.1~M_{\odot}/yr$ are observed for $\rm Z/Z_{\odot} \leq 10^{-5}$ similar to a zero-metallicity case while for $\rm Z/Z_{\odot} = 10^{-4}$ the inflow rate starts to decline earlier due to the dust cooling and fragmentation. For given large inflow rates a central star of $\rm \sim 10^4~M_{\odot}$ may form for $\rm Z/Z_{\odot} \leq 10^{-5}$.

We present and test chemical models for three-dimensional hydrodynamical simulations of galaxies. We explore the effect of changing key parameters such as metallicity, radiation and non-equilibrium versus equilibrium metal cooling approximations on the transition between the gas phases in the interstellar medium. The microphysics is modelled by employing the public chemistry package KROME and the chemical networks have been tested to work in a wide range of densities and temperatures. We describe a simple H/He network following the formation of H$_2$, and a more sophisticated network which includes metals. Photochemistry, thermal processes, and different prescriptions for the H$_2$ catalysis on dust are presented and tested within a one-zone framework. The resulting network is made publicly available on the KROME webpage. We find that employing an accurate treatment of the dust-related processes induces a faster HI--H$_2$ transition. In addition, we show when the equilibrium assumption for metal cooling holds, and how a non-equilibrium approach affects the thermal evolution of the gas and the HII--HI transition. These models can be employed in any hydrodynamical code via an interface to KROME and can be applied to different problems including isolated galaxies, cosmological simulations of galaxy formation and evolution, supernova explosions in molecular clouds, and the modelling of star-forming regions. The metal network can be used for a comparison with observational data of CII 158 $\mu$m emission both for high-redshift as well as for local galaxies.

The detection of $\rm z>6$ quasars reveals the existence of supermassive black holes of a few $\rm 10^9~M_{\odot}$. One of the potential pathways to explain their formation in the infant universe is the so-called direct collapse model which provides massive seeds of $\rm 10^5-10^6~M_{\odot}$. An isothermal direct collapse mandates that halos should be of a primordial composition and the formation of molecular hydrogen remains suppressed in the presence of a strong Lyman Werner flux. In this study, we perform high resolution cosmological simulations for two massive primordial halos employing a detailed chemical model which includes $\rm H^-$ cooling as well as realistic opacities for both the bound-free $\rm H^-$ emission and the Rayleigh scattering of hydrogen atoms. We are able to resolve the collapse up to unprecedentedly high densities of $\rm \sim 10^{-3}~g/cm^3$ and to scales of about $\rm 10^{-4}$ AU. Our results show that the gas cools down to $\rm \sim $ 5000 K in the presence of $\rm H^-$ cooling, and induces fragmentation at scales of about 8000 AU in one of the two simulated halos, which may lead to the formation of a binary. In addition, fragmentation also occurs on the AU scale in one of the halos but the clumps are expected to merge on short time scales. Our results confirm that $\rm H^-$ cooling does not prevent the formation of a supermassive star and the trapping of cooling radiation stabilises the collapse on small scales.

The observed eclipsing time variations in post-common-envelope binaries (PCEBs) can be interpreted as potential evidence for massive Jupiter-like planets, or as a result of magnetic activity, leading to quasi-periodic changes in the quadrupole moment of the secondary star. The latter is commonly referred to as the Applegate mechanism. Following Brinkworth et al. (2006), we employ here an improved version of Applegate's model including the angular momentum exchange between a finite shell and the core of the star. The framework is employed to derive the general conditions under which the Applegate mechanism can work, and is subsequently applied to a sample of 16 close binary systems with potential planets, including 11 PCEBs. Further, we present a detailed derivation and study of analytical models which allow for an straightforward extension to other systems. Using our full numerical framework, we show that the Applegate mechanism can clearly explain the observed eclipsing time variations in 4 of the systems, while the required energy to produce the quadrupole moment variations is too high in at least 8 systems. In the remaining 4 systems, the required energy is comparable to the available energy produced by the star, which we consider as borderline cases. Therefore, the Applegate mechanism cannot uniquely explain the observed period time variations for this entire population. Even in systems where the required energy is too high, the Applegate mechanism may provide an additional scatter, which needs to be considered in the derivation and analysis of planetary models.

Magnetic fields are considered as a vital ingredient of contemporary star formation, and may have been important during the formation of the first stars in the presence of an efficient amplification mechanism. Initial seed fields are provided via plasma fluctuations, and are subsequently amplified by the small-scale dynamo, leading to a strong tangled magnetic field. Here we explore how the magnetic field provided by the small-scale dynamo is further amplified via the $\alpha-\Omega$ dynamo in a protostellar disk and assess its implications. For this purpose, we consider two characteristic cases, a typical Pop.~III star with $10$~M$_\odot$ and an accretion rate of $10^{-3}$~M$_\odot$~yr$^{-1}$, and a supermassive star with $10^5$~M$_\odot$ and an accretion rate of $10^{-1}$~M$_\odot$~yr$^{-1}$. For the $10$~M$_\odot$ Pop.~III star, we find that coherent magnetic fields can be produced on scales of at least $100$~AU, which are sufficient to drive a jet with a luminosity of $100$~L$_\odot$ and a mass outflow rate of $10^{-3.7}$~M$_\odot$~yr$^{-1}$. For the supermassive star, the dynamical timescales in its environment are even shorter, implying smaller orbital timescales and an efficient magnetization out to at least $1000$~AU. The jet luminosity corresponds to $\sim10^{6.0}$~L$_\odot$, and a mass outflow rate of $10^{-2.1}$~M$_\odot$~yr$^{-1}$. We expect that the feedback from the supermassive star can have a relevant impact on its host galaxy.

Numerical simulations show the formation of self-gravitating primordial disks during the assembly of the first structures in the Universe, in particular during the formation of Pop.~III and supermassive stars. Their subsequent evolution is expected to be crucial to determine the mass scale of the first cosmological objects, which depends on the temperature of the gas and the dominant cooling mechanism. Here, we derive a one-zone framework to explore the chemical evolution of such disks and show that viscous heating leads to the collisional dissociation of an initially molecular gas. The effect is relevant on scales of 10 AU (1000 AU) for a central mass of 10 M_sun (10^4 M_sun) at an accretion rate of 10^{-1} M_sun yr^{-1}, and provides a substantial heat input to stabilize the disk. If the gas is initially atomic, it remains atomic during the further evolution, and the effect of viscous heating is less significant. The additional thermal support is particularly relevant for the formation of very massive objects, such as the progenitors of the first supermassive black holes. The stabilizing impact of viscous heating thus alleviates the need for a strong radiation background as a means of keeping the gas atomic.

The origin of strong magnetic fields in the Universe can be explained by amplifying weak seed fields via turbulent motions on small spatial scales and subsequently transporting the magnetic energy to larger scales. This process is known as the turbulent dynamo and depends on the properties of turbulence, i.e. on the hydrodynamical Reynolds number and the compressibility of the gas, and on the magnetic diffusivity. While we know the growth rate the magnetic energy in the linear regime, the saturation level, i.e. the ratio of magnetic energy to turbulent kinetic energy that can be reached, is not known from analytical calculations. In this paper we present the first scale-dependent saturation model based on an effective turbulent resistivity which is determined by the turnover timescale of turbulent eddies and the magnetic energy density. The magnetic resistivity increases compared to the Spitzer value and the effective scale on which the magnetic energy spectrum is at its maximum moves to larger spatial scales. This process ends when the peak reaches a characteristic wavenumber k* which is determined by the critical magnetic Reynolds number. The saturation level of the dynamo also depends on the type of turbulence and differs for the limits of large and small magnetic Prandtl numbers Pm. With our model we find saturation levels between 43.8% and 1.3% for Pm>>1 and between 2.43% and 0.135% for Pm<<1, where the higher values refer to incompressible turbulence and the lower ones to highly compressible turbulence.

This paper is the result of the International Cometary Workshop, held in Toulouse, France in April 2014, where the participants came together to assess our knowledge of comets prior to the ESA Rosetta Mission. In this paper, we look at the composition of the gas and dust from the comae of comets. With the gas, we cover the various taxonomic studies that have broken comets into groups and compare what is seen at all wavelengths. We also discuss what has been learned from mass spectrometers during flybys. A few caveats for our interpretation are discussed. With dust, much of our information comes from flybys. They include {\it in situ} analyses as well as samples returned to Earth for laboratory measurements. Remote sensing IR observations and polarimetry are also discussed. For both gas and dust, we discuss what instruments the Rosetta spacecraft and Philae lander will bring to bear to improve our understanding of comet 67P/Churyumov-Gerasimenko as "ground-truth" for our previous comprehensive studies. Finally, we summarize some of the initial Rosetta Mission findings.

The reionization of helium describes the transition from its singly ionized state to a doubly-ionized state in the intergalactic medium (IGM). This process is important for the thermal evolution of the IGM and influences the mean free path of photons with energies above $54.4$~eV. While it is well-known that helium reionization is mostly driven by the contribution of energetic quasars at $z<6$, we study here how helium reionization proceeds if there is an additional contribution due to the annihilation of dark matter. We explore the effects of different dark matter profiles for the dark matter clumping factor, which can significantly enhance the annihilation rate at late times. We find that the presence of dark matter annihilation enhances the He$^{++}$ abundance at early stages where it would be zero within the standard model, and it can further increase during structure formation, reflecting the increase of the dark matter clumping factor. The latter is, however, degenerate with the build-up of the quasar contribution, and we therefore expect no significant changes at late times. We expect that future studies of the He$^+$ Lyman $\alpha$ forest may help to assess whether the evolution is consistent with the contribution from quasars alone, or if an additional component may be required.

We present new analyses of the photometric lightcurve of Comet 1P/Halley during its 1985/86 apparition. As part of a world-wide campaign coordinated by the International Halley Watch (IHW), narrowband photometry was obtained with telescopes at 18 observatories. Following submissions to and basic reductions by the Photometry and Polarimetry Network of the IHW, we computed production rates and created composite lightcurves for each species. These were used to measure how the apparent rotational period (~7.35 day), along with its shape, evolved with time during the apparition. The lightcurve shape systematically varied from double-peaked to triple-peaked and back again every 8-9 weeks, due to Halley's non-principal axis (complex) rotation and the associated component periods. Unexpectedly, we found a phase shift of one-half cycle also took place during this interval, and therefore the actual beat frequency between the component periods is twice this interval or 16-18 weeks. Preliminary modeling suggests that a single source might produce the entire post-perihelion lightcurve variability and associated evolution. The detailed evolution of the apparent period varied in a non-smooth manner between 7.2 and 7.6 day, likely due to a combination of synodic effects and the interaction of solar illumination with isolated source regions on a body in complex rotation. The need to simultaneously reproduce each of these characteristics will provide very strong additional constraints on Halley's component periods associated with its complex rotation. To assist in these and future analyses, we created a synthetic lightcurve based directly on the measured data. We unexpectedly discovered a strong correlation of ion tail disconnection event start times with minima in the comet's gas production, implying that a decrease in outgassing is another cause of these events.

We present an analysis of Kitt Peak National Observatory and Lowell Observatory observations of comet 103P/Hartley 2 obtained from August through December 2010. The results are then compared with contemporaneous observations made by the EPOXI spacecraft. Each ground-based dataset has previously been investigated individually; the combined dataset has complementary coverage that reduces the time between observing runs and allows us to determine additional apparent periods at intermediate times. We compare CN coma morphology between ground-based datasets, making nine new measurements of apparent periods. The first five are consistent with the roughly linearly increasing apparent period during the apparition found by previous authors. The final four suggest that the change in apparent period slowed or stopped by late November. We also measure an inner coma lightcurve in both CN and R-band ground-based images, finding a single-peaked lightcurve which repeats in phase with the coma morphology. The apparent period from the lightcurve had significantly larger uncertainties than from the coma morphology, but varied over the apparition in a similar manner. Our ground-based lightcurve aligns with the published EPOXI lightcurve, indicating that the lightcurve represents changing activity rather than viewing geometry of structures in the coma. The EPOXI lightcurve can best be phased by a triple-peaked period near 54-55 hr that increases from October to November. This phasing reveals that the spacing between maxima is not constant, and that the overall lightcurve shape evolves from one triple-peaked cycle to the next. These behaviors suggest that much of the scatter in apparent periods derived from ground-based datasets acquired at similar epochs are likely due to limited sampling of the data.

The astonishing diversity in the observed planetary population requires theoretical efforts and advances in planet formation theories. Numerical approaches provide a method to tackle the weaknesses of current planet formation models and are an important tool to close gaps in poorly constrained areas. We present a global disk setup to model the first stages of giant planet formation via gravitational instabilities (GI) in 3D with the block-structured adaptive mesh refinement (AMR) hydrodynamics code ENZO. With this setup, we explore the impact of AMR techniques on the fragmentation and clumping due to large-scale instabilities using different AMR configurations. Additionally, we seek to derive general resolution criteria for global simulations of self-gravitating disks of variable extent. We run a grid of simulations with varying AMR settings, including runs with a static grid for comparison, and study the effects of varying the disk radius. Adopting a marginally stable disk profile (Q_init=1), we validate the numerical robustness of our model for different spatial extensions, from compact to larger, extended disks (R_disk = 10, 100 and 300 AU, M_disk ~ 0.05 M_Sun, M_star = 0.646 M_Sun). By combining our findings from the resolution and parameter studies we find a lower limit of the resolution to be able to resolve GI induced fragmentation features and distinct, turbulence inducing clumps. Irrespective of the physical extension of the disk, topologically disconnected clump features are only resolved if the fragmentation-active zone of the disk is resolved with at least 100 cells, which holds as a minimum requirement for all global disk setups. Our simulations illustrate the capabilities of AMR-based modeling techniques for planet formation simulations and underline the importance of balanced refinement settings to reproduce fragmenting structures.

Massive primordial halos exposed to moderate UV backgrounds are the potential birthplaces of supermassive black holes. In such a halo, an initially isothermal collapse will occur, leading to high accretion rates of $\sim0.1$~M$_\odot$~yr$^{-1}$. During the collapse, the gas in the interior will turn into a molecular state, and form an accretion disk due to the conservation of angular momentum. We consider here the structure of such an accretion disk and the role of viscous heating in the presence of high accretion rates for a central star of $10$, $100$ and $10^4$~M$_\odot$. Our results show that the temperature in the disk increases considerably due to viscous heating, leading to a transition from the molecular to the atomic cooling phase. We found that the atomic cooling regime may extend out to several $100$~AU for a $10^4$~M$_\odot$ central star and provides substantial support to stabilize the disk. It therefore favors the formation of a massive central object. The comparison of clump migration and contraction time scales shows that stellar feedback from these clumps may occur during the later stages of the evolution. Overall, viscous heating provides an important pathway to obtain an atomic gas phase within the center of the halo, and helps in the formation of very massive objects. The latter may collapse to form a massive black hole of about $\geq 10^4$~M$_\odot$.

Turbulence in compressible plasma plays a key role in many areas of astrophysics and engineering. The extreme plasma parameters in these environments, e.g. high Reynolds numbers, supersonic and super-Alfvenic flows, however, make direct numerical simulations computationally intractable even for the simplest treatment -- magnetohydrodynamics (MHD). To overcome this problem one can use subgrid-scale (SGS) closures -- models for the influence of unresolved, subgrid-scales on the resolved ones. In this work we propose and validate a set of constant coefficient closures for the resolved, compressible, ideal MHD equations. The subgrid-scale energies are modeled by Smagorinsky-like equilibrium closures. The turbulent stresses and the electromotive force (EMF) are described by expressions that are nonlinear in terms of large scale velocity and magnetic field gradients. To verify the closures we conduct a priori tests over 137 simulation snapshots from two different codes with varying ratios of thermal to magnetic pressure ($\beta_\mathrm{p} = 0.25, 1, 2.5, 5, 25$) and sonic Mach numbers ($M_s = 2, 2.5, 4$). Furthermore, we make a comparison to traditional, phenomenological eddy-viscosity and $\alpha-\beta-\gamma$ closures. We find only mediocre performance of the kinetic eddy-viscosity and $\alpha-\beta-\gamma$ closures, and that the magnetic eddy-viscosity closure is poorly correlated with the simulation data. Moreover, three of five coefficients of the traditional closures exhibit a significant spread in values. In contrast, our new closures demonstrate consistently high correlation and constant coefficient values over time and and over the wide range of parameters tested. Important aspects in compressible MHD turbulence such as the bi-directional energy cascade, turbulent magnetic pressure and proper alignment of the EMF are well described by our new closures.

Deep surveys with the SKA1-MID array offer for the first time the opportunity to systematically explore the polarization properties of the microJy source population. Our knowledge of the polarized sky approaching these levels is still very limited. In total intensity the population will be dominated by star-forming and normal galaxies to intermediate redshifts ($z \sim1-2$), and low-luminosity AGN to high redshift. The polarized emission from these objects is a powerful probe of their intrinsic magnetic fields and of their magnetic environments. For redshift of order 1 and above the broad bandwidth of the mid-bands span the Faraday thick and thin regimes allowing study of the intrinsic polarization properties of these objects as well as depolarization from embedded and foreground plasmas. The deep field polarization images will provide Rotation Measures data with very high solid angle density allowing a sensitive statistical analysis of the angular variation of RM on critical arc-minute scales from a magnetic component of Large Scale Structure of the Universe.

To understand the evolution of planetary systems, it is important to investigate planets in highly evolved stellar systems, and to explore the implications of their observed properties with respect to potential formation scenarios. Observations suggest the presence of giant planets in post-common-envelope binaries (PCEBs). A particularly well-studied system with planetary masses of 1.7 M_J and 7.0 M_J is NN Ser. We show here that a pure first-generation scenario where the planets form before the common envelope (CE) phase and the orbits evolve due to the changes in the gravitational potential is inconsistent with the current data. We propose a second-generation scenario where the planets are formed from the material that is ejected during the CE, which may naturally explain the observed planetary masses. In addition, hybrid scenarios where the planets form before the CE and evolve due to the accretion of the ejected gas appear as a realistic possibility.

Relativistic jets in active galactic nuclei (AGN) are among the most powerful astrophysical objects discovered to date. Indeed, jetted AGN studies have been considered a prominent science case for SKA, and were included in several different chapters of the previous SKA Science Book (Carilli & Rawlings 2004). Most of the fundamental questions about the physics of relativistic jets still remain unanswered, and await high-sensitivity radio instruments such as SKA to solve them. These questions will be addressed specially through analysis of the massive data sets arising from the deep, all-sky surveys (both total and polarimetric flux) from SKA1. Wide-field very-long-baseline-interferometric survey observations involving SKA1 will serve as a unique tool for distinguishing between extragalactic relativistic jets and star forming galaxies via brightness temperature measurements. Subsequent SKA1 studies of relativistic jets at different resolutions will allow for unprecedented cosmological studies of AGN jets up to the epoch of re-ionization, enabling detailed characterization of the jet composition, magnetic field, particle populations, and plasma properties on all scales. SKA will enable us to study the dependence of jet power and star formation on other properties of the AGN system. SKA1 will enable such studies for large samples of jets, while VLBI observations involving SKA1 will provide the sensitivity for pc-scale imaging, and SKA2 (with its extraordinary sensitivity and dynamic range) will allow us for the first time to resolve and model the weakest radio structures in the most powerful radio-loud AGN.

Our understanding of population III star formation is still in its infancy. They are formed in dark matter minihalos of $\rm 10^5-10^6 M_{\odot}$ at $z=20-30$. Recent high resolution cosmological simulations show that a protostellar disk forms as a consequence of gravitational collapse and fragments into multiple clumps. However, it is not entirely clear if these clumps will be able to survive to form multiple stars as simulations are unable to follow the disk evolution for longer times. In this study, we employ a simple analytical model to derive the properties of marginally stable steady-state disks. Our results show that the stability of the disk depends on the critical value of the viscous parameter $\alpha$. For $\alpha_{crit} = 1$, the disk is stable for an accretion rate of $\rm \leq 10^{-3} M_{\odot}/yr$ and becomes unstable at radii about $\rm \geq 100 AU$ in the presence of an accretion rate of $\rm 10^{-2} M_{\odot}/yr$. For $0.06 < \alpha_{crit} < 1$, the disk can be unstable for both accretion rates. The comparison of the migration and the Kelvin-Helmholtz time scales shows that clumps are expected to migrate inward before reaching the main sequence. Furthermore, in the presence of a massive central star the clumps within the central 1 AU will be tidally disrupted. We also find that UV feedback from the central star is unable to disrupt the disk, and that photo-evaporation becomes important only once the accretion rate has dropped to $\rm 2 \times 10^{-4} M_{\odot}/yr$. As a result, the central star may reach a mass of 100 $\rm M_{\odot}$ or even higher.

The turbulent dynamo may explain the origin of cosmic magnetism. While the exponential amplification of magnetic fields has been studied for incompressible gases, little is known about dynamo action in highly-compressible, supersonic plasmas, such as the interstellar medium of galaxies and the early Universe. Here we perform the first quantitative comparison of theoretical models of the dynamo growth rate and saturation level with three-dimensional magnetohydrodynamical simulations of supersonic turbulence with grid resolutions of up to 1024^3 cells. We obtain numerical convergence and find that dynamo action occurs for both low and high magnetic Prandtl numbers Pm = nu/eta = 0.1-10 (the ratio of viscous to magnetic dissipation), which had so far only been seen for Pm >= 1 in supersonic turbulence. We measure the critical magnetic Reynolds number, Rm_crit = 129 (+43, -31), showing that the compressible dynamo is almost as efficient as in incompressible gas. Considering the physical conditions of the present and early Universe, we conclude that magnetic fields need to be taken into account during structure formation from the early to the present cosmic ages, because they suppress gas fragmentation and drive powerful jets and outflows, both greatly affecting the initial mass function of stars.

Observations of high redshift quasars at $z>6$ indicate that they harbor supermassive black holes (SMBHs) of a billion solar masses. The direct collapse scenario has emerged as the most plausible way to assemble SMBHs. The nurseries for the direct collapse black holes are massive primordial halos illuminated with an intense UV flux emitted by population II (Pop II) stars. In this study, we compute the critical value of such a flux ($J_{21}^{\rm crit}$) for realistic spectra of Pop II stars through three-dimensional cosmological simulations. We derive the dependence of $J_{21}^{\rm crit}$ on the radiation spectra, on variations from halo to halo, and on the impact of X-ray ionization. Our findings show that the value of $J_{21}^{\rm crit}$ is a few times $\rm 10^4$ and only weakly depends on the adopted radiation spectra in the range between $T_{\rm rad}=2 \times 10^4-10^5$ K. For three simulated halos of a few times $\rm 10^{7}$~M$_{\odot}$, $J_{21}^{\rm crit}$ varies from $\rm 2 \times 10^4 - 5 \times 10^4$. The impact of X-ray ionization is almost negligible and within the expected scatter of $J_{21}^{\rm crit}$ for background fluxes of $J_{\rm X,21} \leq 0.1$. The computed estimates of $J_{21}^{\rm crit}$ have profound implications for the quasar abundance at $z=10$ as it lowers the number density of black holes forming through an isothermal direct collapse by a few orders of magnitude below the observed black holes density. However, the sites with moderate amounts of $\rm H_2$ cooling may still form massive objects sufficient to be compatible with observations.

Context. The seeds of the first supermassive black holes may have resulted from the direct collapse of hot primordial gas in $\gtrsim 10^4$ K haloes, forming a supermassive or quasistar as an intermediate stage. Aims. We explore the formation of a protostar resulting from the collapse of primordial gas in the presence of a strong Lyman-Werner radiation background. Particularly, we investigate the impact of turbulence and rotation on the fragmentation behaviour of the gas cloud. We accomplish this goal by varying the initial turbulent and rotational velocities. Methods. We performed 3D adaptive mesh refinement simulations with a resolution of 64 cells per Jeans length using the ENZO code, simulating the formation of a protostar up to unprecedentedly high central densities of $10^{21}$ cm$^{-3}$, and spatial scales of a few solar radii. To achieve this goal, we employed the KROME package to improve modelling of the chemical and thermal processes. Results. We find that the physical properties of the simulated gas clouds become similar on small scales, irrespective of the initial amount of turbulence and rotation. After the highest level of refinement was reached, the simulations have been evolved for an additional ~5 freefall times. A single bound clump with a radius of $2 \times 10^{-2}$ AU and a mass of ~$7 \times 10^{-2}$ M$_{\odot}$ is formed at the end of each simulation, marking the onset of protostar formation. No strong fragmentation is observed by the end of the simulations, regardless of the initial amount of turbulence or rotation, and high accretion rates of a few solar masses per year are found. Conclusions. Given such high accretion rates, a quasistar of $10^5$ M$_{\odot}$ is expected to form within $10^5$ years.

We observed dynamically new sungrazing comet ISON (C/2012 S1) extensively at Lowell Observatory throughout 2013 in order to characterize its behavior prior to perihelion. ISON had "typical" abundances for an Oort Cloud comet. Its dust production, as measured by Afrho, remained nearly constant during the apparition but its CN gas production increased by ~50x. The minimum active area necessary to support observed water production rates exceeded the likely surface area of the nucleus and suggests a population of icy grains in the coma. Together with the flattening of the dust radial profile over time, this is consistent with ejection of a large quantity of slow moving dust and icy grains in the coma at large heliocentric distance. The dust morphology was dominated by the tail, but a faint sunward dust fan was detected in March, April, May, and September. We imaged multiple gas species in September, October, and November. Excess CN signal was observed in the sunward hemisphere in September and early October. In November the excess CN signal was in the tailward hemisphere and two faint CN features appeared approximately orthogonal to the tail with position angles varying by about +/-20 degrees from night to night. Using numerical modeling, we best reproduced the orientation and shape of these features as well as the bulk brightness with a pole oriented approximately towards the Sun and a single source located within ~35 degrees of the equator. The production rates and coma morphology suggest a nucleus that was active over nearly its entire sunward facing hemisphere in September and October but which underwent a significant mass loss event, potentially including fragmentation, shortly before November 1. Significant mass loss likely continued at the same site over subsequent days/weeks and may have catastrophically weakened the nucleus prior to perihelion.

We study the Initial Mass Function (IMF) and host halo properties of Intermediate Mass Black Holes (IMBH, 10^{4-6} Msun) formed inside metal-free, UV illuminated atomic cooling haloes (virial temperature T_vir > 10^4 K) either via the direct collapse of the gas or via an intermediate Super Massive Star (SMS) stage. We achieve this goal in three steps: (a) we derive the gas accretion rate for a proto-SMS to undergo General Relativity instability and produce a direct collapse black hole (DCBH) or to enter the ZAMS and later collapse into a IMBH; (b) we use merger-tree simulations to select atomic cooling halos in which either a DCBH or SMS can form and grow, accounting for metal enrichment and major mergers that halt the growth of the proto-SMS by gas fragmentation. We derive the properties of the host halos and the mass distribution of black holes at this stage, and dub it the "Birth Mass Function"; (c) we follow the further growth of the DCBH due to accretion of leftover gas in the parent halo and compute the final IMBH mass.We consider two extreme cases in which minihalos (T_vir < 10^4 K) can (fertile) or cannot (sterile) form stars and pollute their gas leading to a different IMBH IMF. In the (fiducial) fertile case the IMF is bimodal extending over a broad range of masses, M= (0.5-20)x10^5 Msun, and the DCBH accretion phase lasts from 10 to 100 Myr. If minihalos are sterile, the IMF spans the narrower mass range M= (1-2.8)x10^6 Msun, and the DCBH accretion phase is more extended (70-120 Myr). We conclude that a good seeding prescription is to populate halos (a) of mass 7.5 < log (M_h/Msun) < 8, (b) in the redshift range 8 < z < 17, (c) with IMBH in the mass range 4.75 < log (M_BH/Msun) < 6.25.

The ability of metal free gas to cool by molecular hydrogen in primordial halos is strongly associated with the strength of ultraviolet (UV) flux produced by the stellar populations in the first galaxies. Depending on the stellar spectrum, these UV photons can either dissociate $\rm H_{2}$ molecules directly or indirectly by photo-detachment of $\rm H^{-}$ as the latter provides the main pathway for $\rm H_{2}$ formation in the early universe. In this study, we aim to determine the critical strength of the UV flux above which the formation of molecular hydrogen remains suppressed for a sample of five distinct halos at $z>10$ by employing a higher order chemical solver and a Jeans resolution of 32 cells. We presume that such flux is emitted by PopII stars implying atmospheric temperatures of $\rm 10^{4}$~K. We performed three-dimensional cosmological simulations and varied the strength of the UV flux below the Lyman limit in units of $\rm J_{21}$. Our findings show that the value of $\rm J_{21}^{crit}$ varies from halo to halo and is sensitive to the local thermal conditions of the gas. For the simulated halos it varies from 400-700 with the exception of one halo where $\rm J_{21}^{crit} \geq 1500$. This has important implications for the formation of direct collapse black holes and their estimated population at z > 6. It reduces the number density of direct collapse black holes by almost three orders of magnitude compared to the previous estimates.

Recent discoveries of carbon-enhanced metal-poor stars like SMSS J031300.36-670839.3 provide increasing observational insights into the formation conditions of the first second-generation stars in the Universe, reflecting the chemical conditions after the first supernova explosion. Here, we present the first cosmological simulations with a detailed chemical network including primordial species as well as C, C$^+$, O, O$^+$, Si, Si$^+$, and Si$^{2+}$ following the formation of carbon-enhanced metal poor stars. The presence of background UV flux delays the collapse from $z=21$ to $z=15$ and cool the gas down to the CMB temperature for a metallicity of Z/Z$_\odot$=10$^{-3}$. This can potentially lead to the formation of lower mass stars. Overall, we find that the metals have a stronger effect on the collapse than the radiation, yielding a comparable thermal structure for large variations in the radiative background. We further find that radiative backgrounds are not able to delay the collapse for Z/Z$_\odot$=10$^{-2}$ or a carbon abundance as in SMSS J031300.36-670839.3.

Radiative feedback from populations II stars played a vital role in early structure formation. Particularly, photons below the Lyman limit can escape the star forming regions and produce a background ultraviolet (UV) flux which consequently may influence the pristine halos far away from the radiation sources. These photons can quench the formation of molecular hydrogen by photo-detachment of $\rm H^{-}$. In this study, we explore the impact of such UV radiation on fragmentation in massive primordial halos of a few times $\rm 10^{7}$~M${_\odot}$. To accomplish this goal, we perform high resolution cosmological simulations for two distinct halos and vary the strength of the impinging background UV field in units of $\rm J_{21}$. We further make use of sink particles to follow the evolution for 10,000 years after reaching the maximum refinement level. No vigorous fragmentation is observed in UV illuminated halos while the accretion rate changes according to the thermal properties. Our findings show that a few 100-10, 000 solar mass protostars are formed when halos are irradiated by $\rm J_{21}=10-500$ at $\rm z>10$ and suggest a strong relation between the strength of UV flux and mass of a protostar. This mode of star formation is quite different from minihalos, as higher accretion rates of about $\rm 0.01-0.1$ M$_{\odot}$/yr are observed by the end of our simulations. The resulting massive stars are the potential cradles for the formation of intermediate mass black holes at earlier cosmic times and contribute to the formation of a global X-ray background.

Low-metallicity star formation poses a central problem of cosmology, as it determines the characteristic mass scale and distribution for the first and second generations of stars forming in our Universe. Here, we present a comprehensive investigation assessing the relative impact of metals and magnetic fields, which may both be present during low-metallicity star formation. We show that the presence of magnetic fields generated via the small-scale dynamo stabilises the protostellar disc and provides some degree of support against fragmentation. In the absence of magnetic fields, the fragmentation timescale in our model decreases by a factor of ~10 at the transition from Z=0 to Z>0, with subsequently only a weak dependence on metallicity. Similarly, the accretion timescale of the cluster is set by the large-scale dynamics rather than the local thermodynamics. In the presence of magnetic fields, the primordial disc can become completely stable, therefore forming only one central fragment. At Z>0, the number of fragments is somewhat reduced in the presence of magnetic fields, though the shape of the mass spectrum is not strongly affected in the limits of the statistical uncertainties. The fragmentation timescale, however, increases by roughly a factor of 3 in the presence of magnetic fields. Indeed, our results indicate comparable fragmentation timescales in primordial runs without magnetic fields and Z>0 runs with magnetic fields.

In this paper we show that the Universe is already strongly magnetized at very early epochs during cosmic evolution. Our calculations are based on the efficient amplification of weak magnetic seed fields, which are unavoidably present in the early Universe, by the turbulent small-scale dynamo. We identify two mechanisms for the generation of turbulence in the radiation dominated epoch where velocity fluctuations are produced by the primordial density perturbation and by possible first-order phase transitions at the electroweak or QCD scales. We show that all the necessities for the small-scale dynamo to work are fulfilled. Hence, this mechanism, operating due to primordial density perturbations, guarantees fields with comoving field strength $B_0\sim10^{-6}\varepsilon^{1/2}$ nG on scales up to $\lambda_c\sim0.1$ pc, where $\varepsilon$ is the saturation efficiency. The amplification of magnetic seed fields could be even larger if there are first-order phase transitions in the early Universe. Where, on scales up to $\lambda_c\sim100$ pc, the comoving field strength due to this mechanism will be $B_0\sim10^{-3}\varepsilon^{1/2}$ nG at the present time. Such fields, albeit on small scales, can play an important role in structure formation and could provide an explanation to the apparently observed magnetic fields in the voids of the large-scale structure.

The large-scale structure of the Universe formed from initially small perturbations in the cosmic density field, leading to galaxy clusters with up to 10^15 Msun at the present day. Here, we review the formation of structures in the Universe, considering the first primordial galaxies and the most massive galaxy clusters as extreme cases of structure formation where fundamental processes such as gravity, turbulence, cooling and feedback are particularly relevant. The first non-linear objects in the Universe formed in dark matter halos with 10^5-10^8 Msun at redshifts 10-30, leading to the first stars and massive black holes. At later stages, larger scales became non-linear, leading to the formation of galaxy clusters, the most massive objects in the Universe. We describe here their formation via gravitational processes, including the self-similar scaling relations, as well as the observed deviations from such self-similarity and the related non-gravitational physics (cooling, stellar feedback, AGN). While on intermediate cluster scales the self-similar model is in good agreement with the observations, deviations from such self-similarity are apparent in the core regions, where numerical simulations do not reproduce the current observational results. The latter indicates that the interaction of different feedback processes may not be correctly accounted for in current simulations. Both in the most massive clusters of galaxies as well as during the formation of the first objects in the Universe, turbulent structures and shock waves appear to be common, suggesting them to be ubiquitous in the non-linear regime.

The evolution of magnetic fields in galaxies is still an open problem in astrophysics. In nearby galaxies the far-infrared-radio correlation indicates the coupling between magnetic fields and star formation. The correlation arises from the synchrotron emission of cosmic ray electrons traveling through the interstellar magnetic fields. However, with an increase of the interstellar radiation field (ISRF), inverse Compton scattering becomes the dominant energy loss mechanism of cosmic ray electrons with a typical emission frequency in the X-ray regime. The ISRF depends on the one hand on the star formation rate and becomes stronger in starburst galaxies, and on the other hand increases with redshift due to the evolution of the cosmic microwave background. With a model for the star formation rate of galaxies, the ISRF, and the cosmic ray spectrum, we can calculate the expected X-ray luminosity resulting from the inverse Compton emission. Except for galaxies with an active galactic nucleus the main additional contribution to the X-ray luminosity comes from X-ray binaries. We estimate this contribution with an analytical model as well as with an observational relation, and compare it to the pure inverse Compton luminosity. Using data from the Chandra Deep Field Survey and far-infrared observations from ALMA we then determine upper limits for the cosmic ray energy. Assuming that the magnetic energy in a galaxy is in equipartition with the energy density of the cosmic rays, we obtain upper limits for the magnetic field strength. Our results suggest that the mean magnetic energy of young galaxies is similar to the one in local galaxies. This points toward an early generation of galactic magnetic fields, which is in agreement with current dynamo evolution models.

While Population III stars are typically thought to be massive, pathways towards lower-mass Pop III stars may exist when the cooling of the gas is particularly enhanced. A possible route is enhanced HD cooling during the merging of dark-matter halos. The mergers can lead to a high ionization degree catalysing the formation of HD molecules and may cool the gas down to the cosmic microwave background (CMB) temperature. In this paper, we investigate the merging of mini-halos with masses of a few 10$^5$ M$_\odot$ and explore the feasibility of this scenario. We have performed three-dimensional cosmological hydrodynamics calculations with the ENZO code, solving the thermal and chemical evolution of the gas by employing the astrochemistry package KROME. Our results show that the HD abundance is increased by two orders of magnitude compared to the no-merging case and the halo cools down to $\sim$60 K triggering fragmentation. Based on Jeans estimates the expected stellar masses are about 10 M$_\odot$. Our findings show that the merging scenario is a potential pathway for the formation of low-mass stars.

Baryonic streaming motions produced prior to the epoch of recombination became supersonic during the cosmic dark ages. Various studies suggest that such streaming velocities change the halo statistics and also influence the formation of Population III stars. In this study, we aim to explore the impact of streaming velocities on the formation of supermassive black holes at $z>10$ via the direct collapse scenario. To accomplish this goal, we perform cosmological large eddy simulations for two halos of a few times $\rm 10^{7} M_{\odot}$ with initial streaming velocities of 3, 6 and 9 $\rm km/s$. These massive primordial halos illuminated by the strong Lyman Werner flux are the potential cradles for the formation of direct collapse seed black holes. To study the evolution for longer times, we employ sink particles and track the accretion for 10,000 years. Our findings show that higher streaming velocities increase the circular velocities from about 14 $\rm km/s$ to 16 $\rm km/s$. They also delay the collapse of halos for a few million years, but do not have any significant impact on the halo properties such as turbulent energy, radial velocity, density and accretion rates. Sink particles of about $\rm \sim 10^5 M_{\odot}$ are formed at the end of our simulations and no clear distribution of sink masses is observed in the presence of streaming motions. It is further found that the impact of streaming velocities is less severe in massive halos compared to the minihalos as reported in the previous studies.

Observations of quasars at $\rm z> 6$ report the existence of a billion solar mass black holes. Comprehending their formation in such a short time scale is a matter of ongoing research. One of the most promising scenarios to assemble supermassive black holes is a monolithic collapse of protogalactic gas clouds in atomic cooling halos with $\rm T_{vir} \geq 10^{4} K$. In this article, we study the amplification and impact of magnetic fields during the formation of seed black holes in massive primordial halos. We perform high resolution cosmological magnetohydrodynamics simulations for four distinct halos and follow their collapse for a few free-fall times until the simulations reach a peak density of $\rm 7 \times 10^{-10} g/cm^{3}$. Our findings show that irrespective of the initial seed field, the magnetic field strength reaches a saturated state in the presence of strong accretion shocks. Under such conditions, the growth time becomes very short and amplification occurs rapidly within a small fraction of the free-fall time. We find that the presence of such strong magnetic fields provides additional support against gravity and helps in suppressing fragmentation. Massive clumps of a few hundred solar masses are formed at the end of our simulations and high accretion rates of $\rm 1 M_{\odot}/yr$ are observed. We expect that in the presence of such accretion rates, the clumps will grow to form supermassive stars of $\rm \sim 10^{5} M_{\odot}$. Overall, the role of the magnetic fields seems supportive for the formation of massive black holes.

We used the UltraViolet-Optical Telescope on board Swift to observe the dynamically young comet C/2009 P1 (Garradd) from a heliocentric distance of 3.5 AU pre-perihelion until 4.0 AU outbound. At 3.5 AU pre-perihelion, comet Garradd had one of the highest dust-to-gas ratios ever observed, matched only by comet Hale-Bopp. The evolving morphology of the dust in its coma suggests an outburst that ended around 2.2 AU pre-perihelion. Comparing slit-based measurements and observations acquired with larger fields of view indicated that between 3 AU and 2 AU pre-perihelion a significant extended source started producing water in the coma. We demonstrate that this source, which could be due to icy grains, disappeared quickly around perihelion. Water production by the nucleus may be attributed to a constantly active source of at least 75 km$^2$, estimated to be more than 20 percent of the surface. Based on our measurements, the comet lost $4x10^{11}$ kg of ice and dust during this apparition, corresponding to at most a few meters of its surface.Even though this was likely not Garradd's first passage through the inner solar system, the activity of the comet was complex and changed significantly during the time it was observed.

The first supernova explosions are potentially relevant sources for the production of the first large-scale magnetic fields. For this reason we present a set of high resolution simulations studying the effect of supernova explosions on magnetized, primordial halos. We focus on the evolution of an initially small-scale magnetic field formed during the collapse of the halo. We vary the degree of magnetization, the halo mass, and the amount of explosion energy in order to account for expected variations as well as to infer systematical dependencies of the results on initial conditions. Our simulations suggest that core collapse supernovae with an explosion energy of 10^51 erg and more violent pair instability supernovae with 10^53 erg are able to disrupt halos with masses up to a about 10^6 and 10^7 M_sun, respectively. The peak of the magnetic field spectra shows a continuous shift towards smaller k-values, i.e. larger length scales, over time reaching values as low as k = 4. On small scales the magnetic energy decreases at the cost of the energy on large scales resulting in a well-ordered magnetic field with a strength up to ~ 10^-8 G depending on the initial conditions. The coherence length of the magnetic field inferred from the spectra reaches values up to 250 pc in agreement with those obtained from autocorrelation functions. We find the coherence length to be as large as $50\%$ of the radius of the supernova bubble. Extrapolating this relation to later stages we suggest that significantly strong magnetic fields with coherence lengths as large as 1.5 kpc could be created. We discuss possible implications of our results on processes like recollapse of the halo, first galaxy formation, and the magnetization of the intergalactic medium.

The close binary system NN Serpentis must have gone through a common envelope phase before the formation of its white dwarf. During this phase, a substantial amount of mass was lost from the envelope. The recently detected orbits of circumbinary planets are likely inconsistent with planet formation before the mass loss.We explore whether new planets may have formed from the ejecta of the common envelope and derive the expected planetary mass as a function of radius.We employed the Kashi & Soker model to estimate the amount of mass that is retained during the ejection event and inferred the properties of the resulting disk from the conservation of mass and angular momentum. The resulting planetary masses were estimated from models with and without radiative feedback. We show that the observed planetary masses can be reproduced for appropriate model parameters. Photoheating can stabilize the disks in the interior, potentially explaining the observed planetary orbits on scales of a few AU. We compare the expected mass scale of planets for 11 additional systems with observational results and find hints of two populations, one consistent with planet formation from the ejecta of common envelopes and the other a separate population that may have formed earlier. The formation of the observed planets from the ejecta of common envelopes seems feasible. The model proposed here can be tested through refined observations of additional post-common envelope systems. While it appears observationally challenging to distinguish between the accretion on pre-existing planets and their growth from new fragments, it may be possible to further constrain the properties of the protoplanetary disk through additional observations of current planetary candidates and post-common envelope binary systems.

Chemistry plays a key role in many astrophysical situations regulating the cooling and the thermal properties of the gas, which are relevant during gravitational collapse, the evolution of disks and the fragmentation process. In order to simplify the usage of chemical networks in large numerical simulations, we present the chemistry package KROME, consisting of a Python pre-processor which generates a subroutine for the solution of chemical networks which can be embedded in any numerical code. For the solution of the rate equations, we make use of the high-order solver DLSODES, which was shown to be both accurate and efficient for sparse networks, which are typical in astrophysical applications. KROME also provides a large set of physical processes connected to chemistry, including photochemistry, cooling, heating, dust treatment, and reverse kinetics. The package presented here already contains a network for primordial chemistry, a small metal network appropriate for the modelling of low metallicities environments, a detailed network for the modelling of molecular clouds, a network for planetary atmospheres, as well as a framework for the modelling of the dust grain population. In this paper, we present an extended test suite ranging from one-zone and 1D-models to first applications including cosmological simulations with ENZO and RAMSES and 3D collapse simulations with the FLASH code. The package presented here is publicly available at http://kromepackage.org/ and https://bitbucket.org/krome/krome_stable

Population III stars are the first stars in the Universe to form at z=20-30 out of a pure hydrogen and helium gas in minihalos of 10^5-10^6 M$_\odot$ . Cooling and fragmentation is thus regulated via molecular hydrogen. At densities above 10^8 cm$^{-3}$, the three-body H2 formation rates are particularly important for making the gas fully molecular. These rates were considered to be uncertain by at least a few orders of magnitude. We explore the impact of new accurate three-body H2 formation rates derived by Forrey (2013) for three different minihalos, and compare to the results obtained with three-body rates employed in previous studies. The calculations are performed with the cosmological hydrodynamics code ENZO (release 2.2) coupled with the chemistry package KROME (including a network for primordial chemistry), which was previously shown to be accurate in high resolution simulations. While the new rates can shift the point where the gas becomes fully molecular, leading to a different thermal evolution, there is no trivial trend in how this occurs. While one might naively expect the results to be inbetween the calculations based on Palla et al. (1983) and Abel et al. (2002), the behavior can be close to the former or the latter depending on the dark matter halo that is explored. We conclude that employing the correct three-body rates is about as equally important as the use of appropriate initial conditions, and that the resulting thermal evolution needs to be calculated for every halo individually.

The High Resolution Instrument Infrared Spectrometer (HRI-IR) onboard the Deep Impact Flyby spacecraft detected H2O, CO2, and CO in the coma of the dynamically young Oort cloud comet C/2009 P1 (Garradd) post-perihelion at a heliocentric distance of 2 AU. Production rates were derived for the parent volatiles, QH2O = 4.6e28, QCO2 = 3.9e27, and QCO = 2.9e28 molecules s-1, and are consistent with the trends seen by other observers and within the error bars of measurements acquired during a similar time period. When compiled with other observations of the dominant volatiles of Garradd, unexpected behavior was seen in the release of CO. The H2O outgassing of Garradd, increasing and peaking pre-perihelion and then steadily decreasing, is more typical than that of CO, which monotonically increased throughout the entire apparition. Due to the temporal asymmetry in volatile release, Garradd exhibited the highest CO to H2O abundance ratio ever observed for any comet inside the water snow line at 60 percent during the HRI-IR observations. Also, the HRI-IR made the only direct measurement of CO2, giving a typical cometary abundance ratio of CO2 to H2O of 8 percent but, with only one measurement, no sense of how it varied with orbital position.

The Universe at present is highly magnetized, with fields of the order of a few 10^-5 G and coherence lengths larger than 10 kpc in typical galaxies like the Milky Way. We propose that the magnetic field was amplified to this values already during the formation and the early evolution of the galaxies. Turbulence in young galaxies is driven by accretion as well as by supernova (SN) explosions of the first generation of stars. The small-scale dynamo can convert the turbulent kinetic energy into magnetic energy and amplify very weak primordial magnetic seed fields on short timescales. The amplification takes place in two phases: in the kinematic phase the magnetic field grows exponentially, with the largest growth on the smallest non-resistive scale. In the following non-linear phase the magnetic energy is shifted towards larger scales until the dynamo saturates on the turbulent forcing scale. To describe the amplification of the magnetic field quantitatively we model the microphysics in the interstellar medium (ISM) of young galaxies and determine the growth rate of the small-scale dynamo. We estimate the resulting saturation field strengths and dynamo timescales for two turbulent forcing mechanisms: accretion-driven turbulence and SN-driven turbulence. We compare them to the field strength that is reached, when only stellar magnetic fields are distributed by SN explosions. We find that the small-scale dynamo is much more efficient in magnetizing the ISM of young galaxies. In the case of accretion-driven turbulence a magnetic field strength of the order of 10^-6 G is reached after a time of 24-270 Myr, while in SN-driven turbulence the dynamo saturates at field strengths of typically 10^-5 G after only 4-15 Myr. This is considerably shorter than the Hubble time. Our work can help to understand why present-day galaxies are highly magnetized.

Black holes of a billion solar masses are observed in the infant universe a few hundred million years after the Big Bang. The direct collapse of protogalactic gas clouds in primordial halos with $\rm T_{vir} \geq 10^{4} K$ provides the most promising way to assemble massive black holes. In this study, we aim to determine the characteristic mass scale of seed black holes and the time evolution of the accretion rates resulting from the direct collapse model. We explore the formation of supermassive black holes via cosmological large eddy simulations (LES) by employing sink particles and following their evolution for twenty thousand years after the formation of the first sink. As the resulting protostars were shown to have cool atmospheres in the presence of strong accretion, we assume here that UV feedback is negligible during this calculation. We confirm this result in a comparison run without sinks. Our findings show that black hole seeds with characteristic mass of $\rm 10^{5} M_{\odot}$ are formed in the presence of strong Lyman Werner flux which leads to an isothermal collapse. The characteristic mass is a about two times higher in LES compared to the implicit large eddy simulations (ILES). The accretion rates increase with time and reach a maximum value of 10 $\rm M_{\odot}/yr$ after $\rm 10^{4}$ years. Our results show that the direct collapse model is clearly feasible as it provides the expected mass of the seed black holes.

We explore the evolution of supermassive black holes (SMBH) centered in a circumnuclear disk (CND) as a function of the mass supply from the host galaxy and considering different star formation laws, which may give rise to a self-regulation via the injection of supernova-driven turbulence. A system of equations describing star formation, black hole accretion and angular momentum transport was solved for an axisymmetric disk in which the gravitational potential includes contributions from the black hole, the disk and the hosting galaxy. Our model extends the framework provided by Kawakatu et al. (2008) by separately considering the inner and outer part of the disk, and by introducing a potentially non-linear dependence of the star formation rate on the gas surface density and the turbulent velocity. The star formation recipes are calibrated using observational data for NGC 1097, while the accretion model is based on turbulent viscosity as a source of angular momentum transport in a thin viscous accretion disk. We find that current data provide no strong constraint on the star formation recipe, and can in particular not distinguish between models entirely regulated by the surface density, and models including a dependence on the turbulent velocity. The evolution of the black hole mass, on the other hand, strongly depends on the applied star formation law, as well as the mass supply from the host galaxy. We suggest to explore the star formation process in local AGN with high-resolution ALMA observations to break the degeneracy between different star formation models.

We present new lightcurve measurements of Comet 10P/Tempel 2 carried out with Lowell Observatory's Discovery Channel Telescope in early 2013 when the comet was at aphelion. These data represent some of the first science obtained with this new 4.3-m facility. With Tempel 2 having been observed to exhibit a small but ongoing spin-down in its rotation period for over two decades, our primary goals at this time were two-fold. First, to determine its current rotation period and compare it to that measured shortly after its most recent perihelion passage in 2010, and second, to disentangle the spin-down from synodic effects due to the solar day and the Earth's orbital motion and to determine the sense of rotation, i.e. prograde or retrograde. At our midpoint of 2013 Feb 24, the observed synodic period is 8.948+/-0.001 hr, exactly matching the predicted prograde rotation solution based on 2010 results, and yields a sidereal period of the identical value due to the solar and Earth synodic components just canceling out during the interval of the 2013 observations. The retrograde solution is ruled out because the associated sidereal periods in 2010 and 2013 are quite different even though we know that extremely little outgassing, needed to produce torques, occurred in this interval. With a definitive sense of rotation, the specific amounts of spin-down to the sidereal period could be assessed. The nominal values imply that the rate of spin-down has decreased over time, consistent with the secular drop in water production since 1988. Our data also exhibited an unexpectedly small lightcurve amplitude which appears to be associated with viewing from a large, negative sub-Earth latitude, and a lightcurve shape deviating from a simple sinusoid implying a highly irregularly shaped nucleus.

(Abrigded) Observations of galaxies up to z 2 show a tight correlation between far-infrared and radio continuum emission. We explain the far-infrared - radio continuum correlation by relating star formation and magnetic field strength in terms of turbulent magnetic field amplification, where turbulence is injected by supernova explosions from massive stars. We calculate the expected amount of turbulence in galaxies based on their star formation rates, and infer the expected magnetic field strength due to turbulent dynamo amplification. We estimate the timescales for cosmic ray energy losses via synchrotron emission, inverse Compton scattering, ionization and bremsstrahlung emission, probing up to which redshift strong synchrotron emission can be maintained. We find that the correlation between star formation rate and magnetic field strength in the local Universe can be understood as a result of turbulent magnetic field amplification. If the typical gas density in the interstellar medium increases at high z, we expect an increase of the magnetic field strength and the radio emission, as indicated by current observations. Such an increase would imply a modification of the far-infrared - radio correlation. We expect a breakdown when inverse Compton losses start dominating over synchrotron emission. For a given star formation surface density, we calculate the redshift where the breakdown occurs, yielding z (Sigma_SFR/0.0045 M_solar kpc^{-2} yr^{-1})^{1/(6-alpha/2)}. In this relation, the parameter \alpha describes the evolution of the characteristic ISM density in galaxies as (1+z)^\alpha. Both the possible raise of the radio emission at high redshift and the final breakdown of the far-infrared -- radio correlation at a critical redshift will be probed by the Square Kilometre Array (SKA) and its pathfinders, while the typical ISM density in galaxies will be probed with ALMA.

Population III stars forming in the infant universe at z=30 heralded the end of the cosmic dark ages. They are presumed to be assembled in so-called minihaloes with virial temperatures of a few thousand K where collapse is triggered by molecular hydrogen cooling. A central question concerns their final masses, and whether fragmentation occurs during their formation. While studies employing Lagrangian codes suggest fragmentation via a self-gravitating disk, recent high resolution simulations indicated that disk formation is suppressed. Here we report the first high-resolution large-eddy simulations performed with the Eulerian grid-based code Enzo following the evolution beyond the formation of the first peak, to investigate the accretion of the central massive clump and potential fragmentation. For a total of 3 halos, we see that a disk forms around the first clump. The central clump reaches $\sim10$ solar masses after 40 years, while subsequent accretion is expected at a rate of $10^{-2}$ solar masses per year. In one of these halos, additional clumps form as a result of fragmentation which proceeds at larger scales. We note that subgrid-scale turbulence yields relevant contributions to the stability of the protostellar disks. We conclude that the first protostar may reach masses up to $\rm 40-100 M_{\odot}$, which are only limited by the effect of radiative feedback.

Supermassive stars and quasi-stars (massive stars with a central black hole) are both considered as potential progenitors for the formation of supermassive black holes. They are expected to form from rapidly accreting protostars in massive primordial halos. We explore how long rapidly accreting protostars remain on the Hayashi track, implying large protostellar radii and weak accretion luminosity feedback. We assess the potential role of energy production in the nuclear core, and determine what regulates the evolution of such protostars into quasi-stars or supermassive stars. We follow the contraction of characteristic mass scales in rapidly accreting protostars, and infer the timescales for them to reach nuclear densities. We compare the characteristic timescales for nuclear burning with those for which the extended protostellar envelope can be maintained. We find that the extended envelope can be maintained up to protostellar masses of 3.6x10^8 \dot{m}^3 solar, where \dot{m} denotes the accretion rate in solar masses per year. We expect the nuclear core to exhaust its hydrogen content in 7x10^6 yrs. If accretion rates \dot{m}>>0.14 can still be maintained at this point, a black hole may form within the accreting envelope, leading to a quasi-star. Alternatively, the accreting object will gravitationally contract to become a main-sequence supermassive star. Due to the limited gas reservoir in dark matter halos with 10^7 solar masses, the accretion rate onto the central object may drop at late times, implying the formation of supermassive stars as the typical outcome of direct collapse. However, if high accretion rates are maintained, a quasi-star with an interior black hole may form.

The formation of the first stars in the Universe is regulated by a sensitive interplay of chemistry and cooling with the dynamics of a self-gravitating system. As the outcome of the collapse and the final stellar masses depend sensitively on the thermal evolution, it is necessary to accurately model the thermal evolution in high resolution simulations. As previous investigations raised doubts regarding the convergence of the temperature at high resolution, we investigate the role of the numerical method employed to model the chemistry and the thermodynamics. Here we compare the standard implementation in the adaptive-mesh refinement code \verb|ENZO|, employing a first order backward differentiation formula (BDF), with the 5th order accurate BDF solver \verb|DLSODES|. While the standard implementation in \verb|ENZO| shows a strong dependence on the employed resolution, the results obtained with \verb|DLSODES| are considerably more robust, both with respect to the chemistry and thermodynamics, but also for dynamical quantities such as density, total energy or the accretion rate. We conclude that an accurate modeling of the chemistry and thermodynamics is central for primordial star formation.

Supermassive black holes with up to a $\rm 10^{9}~M_{\odot}$ dwell in the centers of present-day galaxies, and their presence has been confirmed at z $\geq$ 6. Their formation at such early epochs is still an enigma. Different pathways have been suggested to assemble supermassive black holes in the first billion years after the Big Bang. Direct collapse has emerged as a highly plausible scenario to form black holes as it provides seed masses of $\rm 10^{5}-10^{6}~M_{\odot}$. Gravitational collapse in atomic cooling haloes with virial temperatures T$_{vir} \geq 10^{4}$~K may lead to the formation of massive seed black holes in the presence of an intense background UV flux. Turbulence plays a central role in regulating accretion and transporting angular momentum. We present here the highest resolution cosmological large-eddy simulations to date which track the evolution of high-density regions on scales of $0.25$~AU beyond the formation of the first peak, and study the impact of subgrid-scale turbulence. The peak density reached in these simulations is $\rm 1.2 \times 10^{-8}~g~cm^{-3}$. Our findings show that while fragmentation occasionally occurs, it does not prevent the growth of a central massive object resulting from turbulent accretion and occasional mergers. The central object reaches $\rm \sim 1000~M_{\odot}$ within $4$ free-fall times, and we expect further growth up to $\rm 10^{6}~M_{\odot}$ through accretion in about 1 million years. The direct collapse model thus provides a viable pathway of forming high-mass black holes at early cosmic times.

While present standard model of cosmology yields no clear prediction for the initial magnetic field strength, efficient dynamo action may compensate for initially weak seed fields via rapid amplification. In particular, the small-scale dynamo is expected to exponentially amplify any weak magnetic field in the presence of turbulence. We explore whether this scenario is viable using cosmological magneto-hydrodynamics simulations modeling the formation of the first galaxies, which are expected to form in so-called atomic cooling halos with virial temperatures $\rm T_{vir} \geq 10^{4}$ K. As previous calculations have shown that a high Jeans resolution is needed to resolve turbulent structures and dynamo effects, our calculations employ resolutions of up to 128 cells per Jeans length. The presence of the dynamo can be clearly confirmed for resolutions of at least 64 cells per Jeans length, while saturation occurs at approximate equipartition with turbulent energy. As a result of the large Reynolds numbers in primordial galaxies, we expect saturation to occur at early stages, implying magnetic field strengths of \sim0.1 $\mu$G at densities of 10^4 cm^{-3}.

Chemistry has a key role in the evolution of the interstellar medium (ISM), so it is highly desirable to follow its evolution in numerical simulations. However, it may easily dominate the computational cost when applied to large systems. In this paper we discuss two approaches to reduce these costs: (i) based on computational strategies, and (ii) based on the properties and on the topology of the chemical network. The first methods are more robust, while the second are meant to be giving important information on the structure of large, complex networks. To this aim we first discuss the numerical solvers for integrating the system of ordinary differential equations (ODE) associated with the chemical network. We then propose a buffer method that decreases the computational time spent in solving the ODE system. We further discuss a flux-based method that allows one to determine and then cut on the fly the less active reactions. In addition we also present a topological approach for selecting the most probable species that will be active during the chemical evolution, thus gaining information on the chemical network that otherwise would be difficult to retrieve. This topological technique can also be used as an a priori reduction method for any size network. We implemented these methods into a 1D Lagrangian hydrodynamical code to test their effects: both classes lead to large computational speed-ups, ranging from x2 to x5. We have also tested some hybrid approaches finding that coupling the flux method with a buffer strategy gives the best trade-off between robustness and speed-up of calculations.

(Abridged) The small-scale dynamo may play a substantial role in magnetizing the Universe under a large range of conditions, including subsonic turbulence at low Mach numbers, highly supersonic turbulence at high Mach numbers and a large range of magnetic Prandtl numbers Pm, i.e. the ratio of kinetic viscosity to magnetic resistivity. Low Mach numbers may in particular lead to the well-known, incompressible Kolmogorov turbulence, while for high Mach numbers, we are in the highly compressible regime, thus close to Burgers turbulence. In this study, we explore whether in this large range of conditions, a universal behavior can be expected. Our starting point are previous investigations in the kinematic regime. Here, analytic studies based on the Kazantsev model have shown that the behavior of the dynamo depends significantly on Pm and the type of turbulence, and numerical simulations indicate a strong dependence of the growth rate on the Mach number of the flow. Once the magnetic field saturates on the current amplification scale, backreactions occur and the growth is shifted to the next-larger scale. We employ a Fokker-Planck model to calculate the magnetic field amplification during the non-linear regime, and find a resulting power-law growth that depends on the type of turbulence invoked. For Kolmogorov turbulence, we confirm previous results suggesting a linear growth of magnetic energy. For more general turbulent spectra, where the turbulent velocity v_t scales with the characteristic length scale as u_\ell\propto \ell^{\vartheta}, we find that the magnetic energy grows as (t/T_{ed})^{2\vartheta/(1-\vartheta)}, with t the time-coordinate and T_{ed} the eddy-turnover time on the forcing scale of turbulence. For Burgers turbulence, \vartheta=1/2, a quadratic rather than linear growth may thus be expected, and a larger timescale until saturation is reached.

The Square Kilometre Array (SKA) is the most ambitious radio telescope ever planned. With a collecting area of about a square kilometre, the SKA will be far superior in sensitivity and observing speed to all current radio facilities. The scientific capability promised by the SKA and its technological challenges provide an ideal base for interdisciplinary research, technology transfer, and collaboration between universities, research centres and industry. The SKA in the radio regime and the European Extreme Large Telescope (E-ELT) in the optical band are on the roadmap of the European Strategy Forum for Research Infrastructures (ESFRI) and have been recognised as the essential facilities for European research in astronomy. This "White Paper" outlines the German science and R&D interests in the SKA project and will provide the basis for future funding applications to secure German involvement in the Square Kilometre Array.

The present-day Universe is highly magnetized, even though the first magnetic seed fields were most probably extremely weak. To explain the growth of the magnetic field strength over many orders of magnitude fast amplification processes need to operate. The most efficient mechanism known today is the small-scale dynamo, which converts turbulent kinetic energy into magnetic energy leading to an exponential growth of the magnetic field. The efficiency of the dynamo depends on the type of turbulence indicated by the slope of the turbulence spectrum v(l) \propto l^{theta}, where v(l) is the eddy velocity at a scale l. We explore turbulent spectra ranging from incompressible Kolmogorov turbulence with theta = 1/3 to highly compressible Burgers turbulence with theta = 1/2. In this work we analyze the properties of the small-scale dynamo for low magnetic Prandtl numbers Pm, which denotes the ratio of the magnetic Reynolds number, Rm, to the hydrodynamical one, Re. We solve the Kazantsev equation, which describes the evolution of the small-scale magnetic field, using the WKB approximation. In the limit of low magnetic Prandtl numbers the growth rate is proportional to Rm^{(1-theta)/(1+theta)}. We furthermore discuss the critical magnetic Reynolds number Rm_crit, which is required for small-scale dynamo action. The value of Rm_crit is roughly 100 for Kolmogorov turbulence and 2700 for Burgers. Furthermore, we discuss that Rm_crit provides a stronger constraint in the limit of low Pm than it does for large Pm. We conclude that the small-scale dynamo can operate in the regime of low magnetic Prandtl numbers, if the magnetic Reynolds number is large enough. Thus, the magnetic field amplification on small scales can take place in a broad range of physical environments and amplify week magnetic seed fields on short timescales.

Atomic cooling haloes with virial temperatures $\rm T_{vir} \geq 10^{4}$ K are the most plausible sites for the formation of the first galaxies and the first intermediate mass black holes. It is therefore important to assess whether one can obtain robust results concerning their main properties from numerical simulations. A major uncertainty is the presence of turbulence, which is barely resolved in cosmological simulations. We explore the latter both by pursuing high-resolution simulations with up to 64 cells per Jeans length and by incorporating a subgrid-scale turbulence model to account for turbulent pressure and viscosity on unresolved scales. We find that the main physical quantities in the halo, in particular the density, temperature and energy density profile, are approximately converged. However, the morphologies in the central 500 AU change significantly with increasing resolution and appear considerably more turbulent. In a systematic comparison of three different haloes, we further found that the turbulence subgrid-scale model gives rise to more compact central structures, and decreases the amount of vorticity. Such compact morphologies may in particular favor the accretion onto the central object.

We present the first 3D simulations to include the effects of dark matter annihilation feedback during the collapse of primordial mini-halos. We begin our simulations from cosmological initial conditions and account for dark matter annihilation in our treatment of the chemical and thermal evolution of the gas. The dark matter is modelled using an analytical density profile that responds to changes in the peak gas density. We find that the gas can collapse to high densities despite the additional energy input from the dark matter. No objects supported purely by dark matter annihilation heating are formed in our simulations. However, we find that the dark matter annihilation heating has a large effect on the evolution of the gas following the formation of the first protostar. Previous simulations without dark matter annihilation found that protostellar discs around Population III stars rapidly fragmented, forming multiple protostars that underwent mergers or ejections. When dark matter annihilation is included, however, these discs become stable to radii of 1000 AU or more. In the cases where fragmentation does occur, it is a wide binary that is formed.

The small-scale dynamo provides a highly efficient mechanism for the conversion of turbulent into magnetic energy. In astrophysical environments, such turbulence often occurs at high Mach numbers, implying steep slopes in the turbulent spectra. It is thus a central question whether the small-scale dynamo can amplify magnetic fields in the interstellar or intergalactic media, where such Mach numbers occur. To address this long-standing issue, we employ the Kazantsev model for turbulent magnetic field amplification, systematically exploring the effect of different turbulent slopes, as expected for Kolmogorov, Burgers, the Larson laws and results derived from numerical simulations. With the framework employed here, we give the first solution encompassing the complete range of magnetic Prandtl numbers, including Pm << 1, Pm ~ 1 and Pm >> 1. We derive scaling laws of the growth rate as a function of hydrodynamic and magnetic Reynolds number for Pm << 1 and Pm >> 1 for all types of turbulence. A central result concerns the regime of Pm ~ 1, where the magnetic field amplification rate increases rapidly as a function of Pm. This phenomenon occurs for all types of turbulence we explored. We further find that the dynamo growth rate can be decreased by a few orders of magnitude for turbulence spectra steeper than Kolmogorov. We calculate the critical magnetic Reynolds number Rm_c for magnetic field amplification, which is highest for the Burgers case. As expected, our calculation shows a linear behaviour of the amplification rate close to the threshold proportional to Rm-Rm_c. Based on the Kazantsev model, we therefore expect the existence of the small-scale dynamo for any given value of Pm, as long as the magnetic Reynolds number is above the critical threshold.

We explore the amplification of magnetic fields in the high-redshift Universe. For this purpose, we perform high-resolution cosmological simulations following the formation of primordial halos with \sim10^7 M_solar, revealing the presence of turbulent structures and complex morphologies at resolutions of at least 32 cells per Jeans length. Employing a turbulence subgrid-scale model, we quantify the amount of unresolved turbulence and show that the resulting turbulent viscosity has a significant impact on the gas morphology, suppressing the formation of low-mass clumps. We further demonstrate that such turbulence implies the efficient amplification of magnetic fields via the small-scale dynamo. We discuss the properties of the dynamo in the kinematic and non-linear regime, and explore the resulting magnetic field amplification during primordial star formation. We show that field strengths of \sim10^{-5} G can be expected at number densities of \sim5 cm^{-3}.

The first galaxies form due to gravitational collapse of primordial halos. During this collapse, weak magnetic seed fields get amplified exponentially by the small-scale dynamo - a process converting kinetic energy from turbulence into magnetic energy. We use the Kazantsev theory, which describes the small-scale dynamo analytically, to study magnetic field amplification for different turbulent velocity correlation functions. For incompressible turbulence (Kolmogorov turbulence), we find that the growth rate is proportional to the square root of the hydrodynamic Reynolds number, Re^(1/2). In the case of highly compressible turbulence (Burgers turbulence) the growth rate increases proportional to Re^(1/3). With a detailed chemical network we are able to follow the chemical evolution and determine the kinetic and magnetic viscosities (due to Ohmic and ambipolar diffusion) during the collapse of the halo. This way, we can calculate the growth rate of the small-scale dynamo quantitatively and predict the evolution of the small-scale magnetic field. As the magnetic energy is transported to larger scales on the local eddy-timescale, we obtain an estimate for the magnetic field on the Jeans scale. Even there, we find that equipartition with the kinetic energy is reached on small timescales. Dynamically relevant field structures can thus be expected already during the formation of the first objects in the Universe.

We report on photometry and imaging of Comet 10P/Tempel 2 obtained at Lowell Observatory from 1983 through 2011. We measured a nucleus rotation period of 8.950 +/- 0.002 hr from 2010 September to 2011 January. This rotation period is longer than the period we previously measured in 1999, which was itself longer than the period measured in 1988. A nearly linear jet was observed which varied little during a rotation cycle in both R and CN images acquired during the 1999 and 2010 apparitions. We measured the projected direction of this jet throughout the two apparitions and, under the assumption that the source region of the jet was near the comet's pole, determined a rotational pole direction of RA/Dec = 151deg/+59deg from CN measurements and RA/Dec = 173deg/+57deg from dust measurements (we estimate a circular uncertainty of 3deg for CN and 4deg for dust). Different combinations of effects likely bias both gas and dust solutions and we elected to average these solutions for a final pole of RA/Dec = 162 +/- 11deg/+58 +/- 1deg. Photoelectric photometry was acquired in 1983, 1988, 1999/2000, and 2010/2011. The activity exhibited a steep turn-on ~3 months prior to perihelion (the exact timing of which varies) and a relatively smooth decline after perihelion. The activity during the 1999 and 2010 apparitions was similar; limited data in 1983 and 1988 were systematically higher and the difference cannot be explained entirely by the smaller perihelion distance. We measured a "typical" composition, in agreement with previous investigators. Monte Carlo numerical modeling with our pole solution best replicated the observed coma morphology for a source region located near a comet latitude of +80deg and having a radius of ~10deg. Our model reproduced the seasonal changes in activity, suggesting that the majority of Tempel 2's activity originates from a small active region located near the pole.

Stars form by the gravitational collapse of interstellar gas. The thermodynamic response of the gas can be characterized by an effective equation of state. It determines how gas heats up or cools as it gets compressed, and hence plays a key role in regulating the process of stellar birth on virtually all scales, ranging from individual star clusters up to the galaxy as a whole. We present a systematic study of the impact of thermodynamics on gravitational collapse in the context of high-redshift star formation, but argue that our findings are also relevant for present-day star formation in molecular clouds. We consider a polytropic equation of state, P = k rho^Gamma, with both sub-isothermal exponents Gamma < 1 and super-isothermal exponents Gamma > 1. We find significant differences between these two cases. For Gamma > 1, pressure gradients slow down the contraction and lead to the formation of a virialized, turbulent core. Weak magnetic fields are strongly tangled and efficiently amplified via the small-scale turbulent dynamo on timescales corresponding to the eddy-turnover time at the viscous scale. For Gamma < 1, on the other hand, pressure support is not sufficient for the formation of such a core. Gravitational contraction proceeds much more rapidly and the flow develops very strong shocks, creating a network of intersecting sheets and extended filaments. The resulting magnetic field lines are very coherent and exhibit a considerable degree of order. Nevertheless, even under these conditions we still find exponential growth of the magnetic energy density in the kinematic regime.

Star formation in the centers of galaxies is thought to yield massive stars with a possibly top-heavy stellar mass distribution. It is likely that magnetic fields play a crucial role in the distribution of stellar masses inside star-forming molecular clouds. In this context, we explore the effects of magnetic fields, with a typical field strength of 38 microG, such as in RCW 38, and a field strength of 135 microG, similar to NGC 2024 and the infrared dark cloud G28.34+0.06, on the initial mass function (IMF) near (< 10 pc) a 10^7 solar mass black hole. Using these conditions, we perform a series of numerical simulations with the hydrodynamical code FLASH to elucidate the impact of magnetic fields on the IMF and the star-formation efficiency (SFE) emerging from an 800 solar mass cloud. We find that the collapse of a gravitationally unstable molecular cloud is slowed down with increasing magnetic field strength and that stars form along the field lines. The total number of stars formed during the simulations increases by a factor of 1.5-2 with magnetic fields. The main component of the IMF has a lognormal shape, with its peak shifted to sub-solar (< 0.3 M_sun) masses in the presence of magnetic fields, due to a decrease in the accretion rates from the gas reservoir. In addition, we see a top-heavy, nearly flat IMF above ~2 solar masses, from regions that were supported by magnetic pressure until high masses are reached. We also consider the effects of X-ray irradiation if the central black hole is active. X-ray feedback inhibits the formation of sub-solar masses and decreases the SFEs even further. Thus, the second contribution is no longer visible. We conclude that magnetic fields potentially change the SFE and the IMF both in active and inactive galaxies, and need to be taken into account in such calculations.

We report on photometry and imaging of Comet 103P/Hartley 2 obtained at Lowell Observatory from 1991 through 2011. We acquired photoelectric photometry on two nights in 1991, four nights in 1997/98, and 13 nights in 2010/11. We observed a strong secular decrease in water and all other observed species production in 2010/11 from the 1991 and 1997/98 levels. We see evidence for a strong asymmetry with respect to perihelion in the production rates of our usual bandpasses, with peak production occurring ~10 days post-perihelion and production rates considerably higher post-perihelion. The composition was "typical", in agreement with the findings of other investigators. We obtained imaging on 39 nights from 2010 July until 2011 January. We find that, after accounting for their varying parentage and lifetimes, the C2 and C3 coma morphology resemble the CN morphology we reported previously. These species exhibited an hourglass shape in October and November, and the morphology changed with rotation and evolved over time. The OH and NH coma morphology showed hints of an hourglass shape near the nucleus, but was also enhanced in the anti-sunward hemisphere. This tailward brightness enhancement did not vary significantly with rotation and evolved with the viewing geometry. We conclude that all five gas species likely originate from the same source regions on the nucleus, but that OH and NH were derived from small grains of water and ammonia ice that survived long enough to be affected by radiation pressure and driven in the anti-sunward direction. We detected the faint, sunward facing dust jet reported by other authors, and did not detect a corresponding gas feature. This jet varied little during a night but exhibited some variations from night to night, suggesting it is located near the total angular momentum vector.

We study the amplification of magnetic fields during the formation of primordial halos. The turbulence generated by gravitational infall motions during the formation of the first stars and galaxies can amplify magnetic fields very efficiently and on short timescales up to dynamically significant values. Using the Kazantsev theory, which describes the so-called small-scale dynamo - a magnetohydrodynamical process converting kinetic energy from turbulence into magnetic energy - we can then calculate the growth rate of the small-scale magnetic field. Our calculations are based on a detailed chemical network and we include non-ideal magnetohydrodynamical effects such as ambipolar diffusion and Ohmic dissipation. We follow the evolution of the magnetic field up to larger scales until saturation occurs on the Jeans scale. Assuming a weak magnetic seed field generated by the Biermann battery process, both Burgers and Kolmogorov turbulence lead to saturation within a rather small density range. Such fields are likely to become relevant after the formation of a protostellar disk and, thus, could influence the formation of the first stars and galaxies in the Universe.

Numerical simulations suggest that the first galaxies are formed in protogalactic halos with virial temperatures >= 10^4 K. It is likely that such halos are polluted with trace amounts of metals produced by the first generation of stars. The presence of dust can significantly change the chemistry and dynamics of early galaxies. In this article, we aim to assess the role of dust on the thermal and dynamical evolution of the first galaxies in the presence of a background UV flux, and its implications for the observability of Lyman alpha emitters and sub-mm sources. We have performed high resolution cosmological simulations using the adaptive mesh refinement code FLASH to accomplish this goal. We have developed a chemical network appropriate for these conditions and coupled it with the FLASH code. The main ingredients of our chemical model include the formation of molecules, a multi-level treatment of atomic hydrogen, line trapping of Lyman alpha photons and, photoionization and photodissociation processes in a UV background. We found that the formation of molecules (H_{2} and HD) is significantly enhanced in the presence of dust grains. The presence of a background UV flux strongly influences the formation of molecules by photodissociating them. We explore the evolution after a major merger, leading to the formation of a binary disk. These disks have gas masses of ~10^{7} M_sun at a redshift of 5.4. Each disk lies in a separate subhalo as a result of the merger event. The disks are supported by turbulent pressure due to the highly supersonic turbulence present in the halo. The presence of dust does not significantly reduce the Lyman alpha emission. The emission of Lyman alpha is extended and originates from the envelope of the halo due to line trapping effects. We also find that dust masses of 10^8 M_sun are required to observe the dust continuum emission from z 5 galaxies with ALMA.

Here we summarize our recent results of high-resolution computer simulations on the turbulent amplification of weak magnetic seed fields showing that such fields will be exponentially amplified also during the gravitational collapse reminiscent to the situation during primordial star formation. The exponential magnetic field amplification is driven by the turbulent small-scale dynamo that can be only observed in computer simulations if the turbulent motions in the central core are sufficiently resolved. We find that the Jeans length, which determines the central core region, has to be resolved by at least 30 grid cells to capture the dynamo activity. We conclude from our studies that strong magnetic fields will be unavoidably created already during the formation of the first stars in the Universe, potentially influencing their evolution and mass distribution.

The small-scale dynamo is a process by which turbulent kinetic energy is converted into magnetic energy, and thus is expected to depend crucially on the nature of turbulence. In this work, we present a model for the small-scale dynamo that takes into account the slope of the turbulent velocity spectrum v(l) ~ l^theta, where l and v(l) are the size of a turbulent fluctuation and the typical velocity on that scale. The time evolution of the fluctuation component of the magnetic field, i.e., the small-scale field, is described by the Kazantsev equation. We solve this linear differential equation for its eigenvalues with the quantum-mechanical WKB-approximation. The validity of this method is estimated as a function of the magnetic Prandtl number Pm. We calculate the minimal magnetic Reynolds number for dynamo action, Rm_crit, using our model of the turbulent velocity correlation function. For Kolmogorov turbulence (theta=1/3), we find that the critical magnetic Reynolds number is approximately 110 and for Burgers turbulence (theta=1/2) approximately 2700. Furthermore, we derive that the growth rate of the small-scale magnetic field for a general type of turbulence is Gamma ~ Re^((1-theta)/(1+theta)) in the limit of infinite magnetic Prandtl numbers. For decreasing magnetic Prandtl number (down to Pm approximately larger than 10), the growth rate of the small-scale dynamo decreases. The details of this drop depend on the WKB-approximation, which becomes invalid for a magnetic Prandtl number of about unity.

We study the influence of initial conditions on the magnetic field amplification during the collapse of a magnetised gas cloud. We focus on the dependence of the growth and saturation level of the dynamo generated field on the turbulent properties of the collapsing cloud. In particular, we explore the effect of varying the initial strength and injection scale of turbulence and the initial uniform rotation of the collapsing magnetised cloud. In order to follow the evolution of the magnetic field in both the kinematic and the nonlinear regime, we choose an initial field strength of $\simeq 1\,\mkG$ with the magnetic to kinetic energy ratio, $E_{\rm m}/E_{\rm k} \sim 10^{-4}$. Both gravitational compression and the small-scale dynamo initially amplify the magnetic field. Further into the evolution, the dynamo-generated magnetic field saturates but the total magnetic field continues to grow because of compression. The saturation of the small-scale dynamo is marked by a change in the slope of $B/\rho^{2/3}$ and by a shift in the peak of the magnetic energy spectrum from small scales to larger scales. For the range of initial Mach numbers explored in this study, the dynamo growth rate increases as the Mach number increases from $v_{\rm rms}/c_{\rm s}\sim 0.2$ to 0.4 and then starts decreasing from $v_{\rm rms}/c_{\rm s}\sim 1.0$. We obtain saturation values of $E_{\rm m}/E_{\rm k} = 0.2 - 0.3$ for these runs. Simulations with different initial injection scales of turbulence also show saturation at similar levels. For runs with different initial rotation of the cloud, the magnetic energy saturates at $E_{\rm m}/E_{\rm k}\sim 0.2 - 0.4$ of the equipartition value. (Abridged)

We discuss the amplification of magnetic fields by the small-scale dynamo, a process that could efficiently produce strong magnetic fields in the first galaxies. In addition, we derive constraints on the primordial field strength from the epoch of reionization.

Primordial magnetic fields generated in the early universe are subject of considerable investigation, and observational limits on their strength are required to constrain the theory. Due to their impact on the reionization process, the strength of primordial fields can be limited using the latest data on reionization and the observed UV-luminosity function of high-redshift galaxies. Given the steep faint-end slope of the luminosity function, faint galaxies contribute substantial ionizing photons, and the low-luminosity cutoff has an impact on the total budget thereof. Magnetic pressure from primordial fields affects such cutoff by preventing collapse in halos with mass below 10^{10} M_solar (B_0 / 3 nG)^3, with B_0 the co-moving field strength. In this letter, the implications of these effects are consistently incorporated in a simplified model for reionization, and the uncertainties due to the cosmological parameters, the reionization parameters and the observed UV luminosity function are addressed. We show that the observed ionization degree at z\sim7 leads to the strongest upper limit of B_0\lsim 2-3nG. Stronger limits could follow from measurements of high ionization degree at z>7.

We review current ideas on the origin of galactic and extragalactic magnetic fields. We begin by summarizing observations of magnetic fields at cosmological redshifts and on cosmological scales. These observations translate into constraints on the strength and scale magnetic fields must have during the early stages of galaxy formation in order to seed the galactic dynamo. We examine mechanisms for the generation of magnetic fields that operate prior during inflation and during subsequent phase transitions such as electroweak symmetry breaking and the quark-hadron phase transition. The implications of strong primordial magnetic fields for the reionization epoch as well as the first generation of stars is discussed in detail. The exotic, early-Universe mechanisms are contrasted with astrophysical processes that generate fields after recombination. For example, a Biermann-type battery can operate in a proto-galaxy during the early stages of structure formation. Moreover, magnetic fields in either an early generation of stars or active galactic nuclei can be dispersed into the intergalactic medium.

Magnetic fields appear to be ubiquitous in astrophysical environments. Their existence in the intracluster medium is established through observations of synchrotron emission and Faraday rotation. On the other hand, the nature of magnetic fields outside of clusters, where observations are scarce and controversial, remains largely unknown. In this chapter, we review recent developments in our understanding of the nature and origin of intergalactic magnetic fields, and in particular, intercluster fields. A plausible scenario for the origin of galactic and intergalactic magnetic fields is for seed fields, created in the early universe, to be amplified by turbulent flows induced during the formation of the large scale structure. We present several mechanisms for the generation of seed fields both before and after recombination. We then discuss the evolution and role of magnetic fields during the formation of the first starts. We describe the turbulent amplification of seed fields during the formation of large scale structure and the nature of the magnetic fields that arise. Finally, we discuss implications of intergalactic magnetic fields.

We study the growth rate and saturation level of the turbulent dynamo in magnetohydrodynamical simulations of turbulence, driven with solenoidal (divergence-free) or compressive (curl-free) forcing. For models with Mach numbers ranging from 0.02 to 20, we find significantly different magnetic field geometries, amplification rates, and saturation levels, decreasing strongly at the transition from subsonic to supersonic flows, due to the development of shocks. Both extreme types of turbulent forcing drive the dynamo, but solenoidal forcing is more efficient, because it produces more vorticity.

The Lyman alpha line is a robust tracer of high redshift galaxies. We present estimates of Lyman alpha emission from a protogalactic halo illuminated by UV background radiation fields with various intensities. For this purpose, we performed cosmological hydrodynamics simulations with the adaptive mesh refinement code FLASH, including a detailed network for primordial chemistry,comprising the formation of primordial molecules, a multi-level model for the hydrogen atom as well as the photo-ionization and photo-dissociation processes in a UV background. We find that the presence of a background radiation field J_21 excites the emission of Lyman alpha photons, increasing the Lyman alpha luminosity up to two orders of magnitude. For a halo of \sim 10^10 M_sun, we find that a maximum flux of 5 \times 10^-15 erg cm^-2 s^-1 is obtained for J21 \times f_esc = 0.1, where f_esc is the escape fraction of the ionizing radiation. Depending on the environmental conditions, the flux may vary by three orders of magnitude. For J_21 \times f_esc > 0.1 the Lyman alpha luminosity decreases as the atomic hydrogen abundance becomes rather small. The fluxes derived here can be probed using Subaru and the upcoming James Webb Space Telescope. The emission of Lyman alpha photons is extended and comes from the envelope of the halo rather than its core. In the center of the halo, line trapping becomes effective above columns of 10^22 cm^-2 and suppresses the emission of Lyman alpha. In addition, cooling by primordial molecules may decrease the gas temperature in the central region, which further reduces Lyman alpha emission. In the central core, H_2 is photo-dissociated for a background flux of J_21 \geq 1000. For weaker radiation fields, i.e. J_21 < 0.1, H_2 and HD cooling are particularly strong in the center of the halo, leading to gas temperatures as low as \sim 100 K.

We report on narrowband CN imaging of Comet 103P/Hartley 2 obtained at Lowell Observatory on 39 nights from 2010 July until 2011 January. We observed two features, one generally to the north and the other generally to the south. The CN morphology varied during the apparition: no morphology was seen in July; in August and September the northern feature dominated and appeared as a mostly face-on spiral; in October, November, and December the northern and southern features were roughly equal in brightness and looked like more side-on corkscrews; in January the southern feature was dominant but the morphology was indistinct due to very low signal. The morphology changed smoothly during each night and similar morphology was seen from night to night. However, the morphology did not exactly repeat each rotation cycle, suggesting that there is a small non-principal axis rotation. Based on the repetition of the morphology, we find evidence that the fundamental rotation period was increasing: 16.7 hr from August 13-17, 17.2 hr from September 10-13, 18.2 hr from October 12-19, and 18.7 hr from October 31-November 7. We conducted Monte Carlo jet modeling to constrain the pole orientation and locations of the active regions based on the observed morphology. Our preliminary, self-consistent pole solution has an obliquity of 10 deg relative to the comet's orbital plane (i.e., it is centered near RA = 257 deg and Dec=+67 deg with an uncertainty around this position of about 15 deg) and has two mid-latitude sources, one in each hemisphere.

Recent FERMI observations provide a lower limit of 10^{-15} G for the magnetic field strength in the intergalactic medium (IGM). This is consistent with theoretical expectations based on the Biermann battery effect, which predicts such IGM fields already at redshifts z~10. During gravitational collapse, such magnetic fields can be amplified by compression and by turbulence, giving rise to the small-scale dynamo. On scales below the Jeans length, the eddy turnover timescale is much shorter than the free-fall timescale, so that saturation can be reached during collapse. This scenario has been tested and confirmed with magneto-hydrodynamical simulations following the collapse of a turbulent, weakly magnetized cloud. Based on a spectral analysis, we confirm that turbulence is injected on the Jeans scale. For the power spectrum of the magnetic field, we obtain the Kazantsev slope which is characteristic for the small-scale dynamo. A calculation of the critical length scales for ambipolar diffusion and Ohmic dissipation shows that these scales are always small enough to allow significant amplification of the magnetic field by small-scale eddies. We discuss potential implications for the protostellar accretion disk, with particular focus on the magneto-rotational instability, which may change the morphology of the disk and reduce the accretion rate by a factor of a few.

Cosmic structure formation is characterized by the complex interplay between gravity, turbulence, and magnetic fields. The processes by which gravitational energy is converted into turbulent and magnetic energies, however, remain poorly understood. Here, we show with high-resolution, adaptive-mesh simulations that MHD turbulence is efficiently driven by extracting energy from the gravitational potential during the collapse of a dense gas cloud. Compressible motions generated during the contraction are converted into solenoidal, turbulent motions, leading to a natural energy ratio of E_sol/E_tot of approximately 2/3. We find that the energy injection scale of gravity-driven turbulence is close to the local Jeans scale. If small seeds of the magnetic field are present, they are amplified exponentially fast via the small-scale dynamo process. The magnetic field grows most efficiently on the smallest scales, for which the stretching, twisting, and folding of field lines, and the turbulent vortices are sufficiently resolved. We find that this scale corresponds to about 30 grid cells in the simulations. We thus suggest a new minimum resolution criterion of 30 cells per Jeans length in (magneto)hydrodynamical simulations of self-gravitating gas, in order to resolve turbulence on the Jeans scale, and to capture minimum dynamo amplification of the magnetic field. Due to numerical diffusion, however, any existing simulation today can at best provide lower limits on the physical growth rates. We conclude that a small, initial magnetic field can grow to dynamically important strength on time scales significantly shorter than the free-fall time of the cloud.

The formation of the first galaxies is accompanied by large accretion flows and virialization shocks, during which the gas is shock-heated to temperatures of $\sim10^4$ K, leading to potentially strong fluxes in the Lyman alpha line. Indeed, a number of Lyman alpha blobs has been detected at high redshift. In this letter, we explore the origin of such Lyman alpha emission using cosmological hydrodynamical simulations that include a detailed model of atomic hydrogen as a multi-level atom and the effects of line trapping with the adaptive mesh refinement code FLASH. We see that baryons fall into the center of a halo through cold streams of gas, giving rise to a Lyman alpha luminosity of at least $\rm 10^{44} erg s^{-1}$ at $\rm z=4.7$, similar to observed Lyman alpha blobs. We find that a Lyman alpha flux of $\rm 5.0\times 10^{-17} erg cm^{-2} s^{-1}$ emerges from the envelope of the halo rather than its center, where the photons are efficiently trapped. Such emission can be probed in detail with the upcoming James Webb Space Telescope (JWST) and will constitute an important probe of gas infall and accretion.

We imaged comet 10P/Tempel 2 on 32 nights from 1999 April through 2000 March. R-band lightcurves were obtained on 11 of these nights from 1999 April through 1999 June, prior to both the onset of significant coma activity and perihelion. Phasing of the data yields a double-peaked lightcurve and indicates a nucleus rotational period of 8.941 +/- 0.002 hr with a peak-to-peak amplitude of ~0.75 mag. Our data are sufficient to rule out all other possible double-peaked solutions as well as the single- and triple- peaked solutions. This rotation period agrees with one of five possible solutions found in post-perihelion data from 1994 by Mueller and Ferrin (1996, Icarus, 123, 463-477), and unambiguously eliminates their remaining four solutions. We applied our same techniques to published lightcurves from 1988 which were obtained at an equivalent orbital position and viewing geometry as in 1999. We found a rotation period of 8.932 +/- 0.001 hr in 1988, consistent with the findings of previous authors and incompatible with our 1999 solution. This reveals that Tempel 2 spun-down by ~32 s between 1988 and 1999 (two intervening perihelion passages). If the spin-down is due to a systematic torque, then the rotation period prior to perihelion during the 2010 apparition is expected to be an additional 32 s longer than in 1999.

We explore the amplification of magnetic seed fields during the formation of the first stars and galaxies. During gravitational collapse, turbulence is created from accretion shocks, which may act to amplify weak magnetic fields in the protostellar cloud. Numerical simulations showed that such turbulence is sub-sonic in the first star-forming minihalos, and highly supersonic in the first galaxies with virial temperatures larger than 10^4 K. We investigate the magnetic field amplification during the collapse both for Kolmogorov and Burgers-type turbulence with a semi-analytic model that incorporates the effects of gravitational compression and small-scale dynamo amplification. We find that the magnetic field may be substantially amplified before the formation of a disk. On scales of 1/10 of the Jeans length, saturation occurs after ~10^8 yr. Although the saturation behaviour of the small-scale dynamo is still somewhat uncertain, we expect a saturation field strength of the order ~10^{-7} n^{0.5} G in the first star-forming halos, with n the number density in cgs units. In the first galaxies with higher turbulent velocities, the magnetic field strength may be increased by an order of magnitude, and saturation may occur after 10^6 to 10^7 yr. In the Kolmogorov case, the magnetic field strength on the integral scale (i.e. the scale with most magnetic power) is higher due to the characteristic power-law indices, but the difference is less than a factor of 2 in the saturated phase. Our results thus indicate that the precise scaling of the turbulent velocity with length scale is of minor importance. They further imply that magnetic fields will be significantly enhanced before the formation of a protostellar disk, where they may change the fragmentation properties of the gas and the accretion rate.

Cosmological hydrodynamical simulations of primordial star formation suggest that the gas within the first star-forming halos is turbulent. This has strong implications on the subsequent evolution, in particular on the generation of magnetic fields. Using high-resolution numerical simulations, we show that in the presence of turbulence, weak seed magnetic fields are exponentially amplified by the small-scale dynamo during the formation of the first stars. We conclude that strong magnetic fields are generated during the birth of the first stars in the universe, potentially modifying the mass distribution of these stars and influencing the subsequent cosmic evolution. We find that the presence of the small-scale turbulent dynamo can only be identified in numerical simulations in which the turbulent motions in the central core are resolved with at least 32 grid cells.

We discuss the formation of supermassive black holes in the early universe, and how to probe their subsequent evolution with the upcoming mm/sub-mm telescope ALMA. We first focus on the chemical and radiative conditions for black hole formation, in particular considering radiation trapping and molecular dissociation effects. We then turn our attention towards the magnetic properties in the halos where the first black holes form, and show that the presence of turbulence may lead to a magnetic dynamo, which could support the black hole formation process by providing an efficient means of transporting the angular momentum. We finally focus on observable properties of high-redshift black holes with respect to ALMA, and discuss how to distinguish between chemistry driven by the starburst and chemistry driven by X-rays from the black hole.

In massive primordial galaxies, the gas may directly collapse and form a single central massive object if cooling is suppressed. Line cooling by molecular hydrogen can be suppressed in the presence of a strong soft-ultraviolet radiation field, but the role played by other cooling mechanisms is less clear. In optically thin gas, Lyman-Alpha cooling can be very effective, maintaining the gas temperature below 10^4 K over many orders of magnitude in density. However, the large neutral hydrogen column densities present in primordial galaxies render them highly optically thick to Lyman-Alpha photons. In this letter, we examine in detail the effects of the trapping of these Lyman-Alpha photons on the thermal and chemical evolution of the gas. We show that despite the high optical depth in the Lyman series lines, cooling is not strongly suppressed, and proceeds via other atomic hydrogen transitions, in particular the 2s-1s and the 3-2 transitions. At densities larger than 10^9 cm^{-3}, collisional dissociation of molecular hydrogen becomes the dominant cooling process and decreases the gas temperature to about 5000 K. The gas temperature evolves with density as $T \propto \rho^{\gamma_{\rm eff} - 1}$, with $\gamma_{\rm eff} = 0.97-0.98$. The evolution is thus very close to isothermal, and so fragmentation is possible, but unlikely to occur during the initial collapse. However, after the formation of a massive central object, we expect that later-infalling, higher angular momentum material will form an accretion disk that may be unstable to fragmentation, which may give rise to star formation with a top-heavy IMF.

(abridged) We explore how ALMA observations can probe high-redshift galaxies in unprecedented detail. We discuss the main observables that are excited by the large-scale starburst, and formulate expectations for the chemistry and the fluxes in the center of active galaxies, where chemistry may be driven by the absorption of X-ray photons. We show that such X-ray dominated regions (XDRs) should be large enough to be resolved with ALMA, and predict the expected amount of emission in CO and various fine-structure lines. We discuss how such XDRs can be distinguished from a strong starburst on the same spatial scales based on the CO line SED. Our models are compared to known sources like NGC 1068 and APM 08279. We also analyze the properties of the z=6.42 quasar SDSS J114816.64+525150.3, and find that the observed emission in CO, [CII] and [CI] requires a dense warm and a low-density cold gas component. We estimate the expected number of sources at redshifts higher than 6, finding that one could expect one black hole with $10^6$ solar masses per arcmin$^2$.

The formation of the first galaxies at redshifts z ~ 10-15 signaled the transition from the simple initial state of the universe to one of ever increasing complexity. We here review recent progress in understanding their assembly process with numerical simulations, starting with cosmological initial conditions and modelling the detailed physics of star formation. In this context we emphasize the importance and influence of selecting appropriate initial conditions for the star formation process. We revisit the notion of a critical metallicity resulting in the transition from primordial to present-day initial mass functions and highlight its dependence on additional cooling mechanisms and the exact initial conditions. We also review recent work on the ability of dust cooling to provide the transition to present-day low-mass star formation. In particular, we highlight the extreme conditions under which this transition mechanism occurs, with violent fragmentation in dense gas resulting in tightly packed clusters.

We explore the effects of magnetic energy dissipation on the formation of the first stars. For this purpose, we follow the evolution of primordial chemistry in the presence of magnetic fields in the post-recombination universe until the formation of the first virialized halos. From the point of virialization, we follow the protostellar collapse up to densities of $\sim10^{12}$ cm$^{-3}$ in a one-zone model. In the intergalactic medium (IGM), comoving field strengths of $\gtrsim0.1$ nG lead to Jeans masses of $10^8 M_\odot$ or more and thus delay gravitational collapse in the first halos until they are sufficiently massive. During protostellar collapse, we find that the temperature minimum at densities of $\sim10^3$ cm$^{-3}$ does not change significantly, such that the characteristic mass scale for fragmentation is not affected. However, we find a significant temperature increase at higher densities for comoving field strengths of $\gtrsim0.1$ nG. This may delay gravitational collapse, in particular at densities of $\sim10^9$ cm$^{-3}$, where the proton abundance drops rapidly and the main contribution to the ambipolar diffusion resistivity is due to collisions with Li$^+$. After the formation of the protostar, the increased gas temperatures may enhance the protostellar accretion rate. Our model confirms that initial weak magnetic fields may be amplified considerably during gravitational collapse and become dynamically relevant. For instance, a comoving field strength above $10^{-5}$ nG will be amplified above the critical value for the onset of jets which can magnetize the IGM.

Dark stars powered by dark matter annihilation have been proposed as the first luminous sources in the universe. These stars are believed to form in the central dark matter cusp of low-mass minihalos. Recent calculations indicate stellar masses up to \sim1000 solar masses and/or have very long lifetimes. The UV photons from these objects could therefore contribute significantly to cosmic reionization. Here we show that such dark star models would require a somewhat artificial reionization history, based on a double-reionization phase and a late star-burst near redshift $z\sim6$, in order to fulfill the WMAP constraint on the optical depth as well as the Gunn-Peterson constraint at $z\sim6$. This suggests that, if dark stars were common in the early universe, then models are preferred which predict a number of UV photons similar to conventional Pop. III stars. This excludes dark stars with 100 solar masses that enter a main-sequence phase and other models that lead to a strong increase in the number of UV photons. We also derive constraints for massive as well as light dark matter candidates from the observed X-ray, gamma-ray and neutrino background, considering dark matter profiles which have been steepened during the formation of dark stars. This increases the clumping factor at high redshift and gives rise to a higher dark matter annihilation rate in the early universe. We furthermore estimate the potential contribution from the annihilation products in the remnants of dark stars, which may provide a promising path to constrain such models further, but which is currently still uncertain.

Light dark matter annihilating into electron-positron pairs emits a significant amount of internal bremsstrahlung that may contribute to the cosmic gamma-ray background. The amount of emitted gamma-rays depends on the dark matter clumping factor. Recent calculations indicate that this value should be of order $10^6-10^7$. That allows us to calculate the expected gamma-ray background contribution from dark matter annihilation. We find that the light dark matter model can be ruled out if a constant thermally-averaged cross section is assumed (s-wave annihilation). For more massive dark matter candidates like neutralinos, however, cosmic constraints are weaker.

We study the prospects of finding the first quasars in the universe with ALMA and JWST. For this purpose, we derive a model for the high-redshift black hole population based on observed relations between the black hole mass and the host galaxy. We re-address previous constraints from the X-ray background with particular focus on black hole luminosities below the Eddington limit as observed in many local AGN. For such luminosities, up to 20% of high-redshift black holes can be active quasars. We then discuss the observables of high-redshift black holes for ALMA and JWST by adopting NGC 1068 as a reference system. We calculate the expected flux of different fine-structure lines for a similar system at higher redshift, and provide further predictions for high-J CO lines. We discuss the expected fluxes from stellar light, the AGN continuum and the Lyman $\alpha$ line for JWST. Line fluxes observed with ALMA can be used to derive detailed properties of high-redshift sources. We suggest two observational strategies to find potential AGN at high redshift and estimate the expected number of sources, which is between 1-10 for ALMA with a field of view of $\sim(1')^2$ searching for line emission and 100-1000 for JWST with a field of view of $(2.16')^2$ searching for continuum radiation. We find that both telescopes can probe high-redshift quasars down to redshift 10 and beyond, and therefore truely detect the first quasars in the universe.

We calculate the reionization history for different models of the stellar population and explore the effects of primordial magnetic fields, dark matter decay and dark matter annihilation on reionization. We find that stellar populations based on a Scalo-type initial mass function for Population II stars can be ruled out as sole sources for reionization, unless star formation efficiencies of more than 10% or very high photon escape fractions from the parental halo are adopted. When considering primordial magnetic fields, we find that the additional heat injection from ambipolar diffusion and decaying MHD turbulence has significant impact on the thermal evolution and the ionization history of the post-recombination universe and on structure formation. The magnetic Jeans mass changes the typical mass scale of the star forming halos, and depending on the adopted stellar model we derive upper limits to the magnetic field strength between 0.7 and $5 $nG (comoving). For dark matter annihilation, we find an upper limit to the thermally averaged mass-weighted cross section of $10^{-33} \mathrm{cm}^3\mathrm{/s/eV}$. For dark matter decay, our calculations yield a lower limit to the lifetime of dark matter particles of $3\times10^{23}$ s. These limits are in agreement with constraints from recombination and the X-ray background and provide an independent confirmation at a much later epoch.

Magnetic fields in the early universe can significantly alter the thermal evolution and the ionization history during the dark ages. This is reflected in the 21 cm line of atomic hydrogen, which is coupled to the gas temperature through collisions at high redshifts, and through the Wouthuysen-Field effect at low redshifts. We present a semi-analytic model for star formation and the build-up of a Lyman alpha background in the presence of magnetic fields, and calculate the evolution of the mean 21 cm brightness temperature and its frequency gradient as a function of redshift. We further discuss the evolution of linear fluctuations in temperature and ionization in the presence of magnetic fields and calculate the effect on the 21 cm power spectrum. At high redshifts, the signal is increased compared to the non-magnetic case due to the additional heat input into the IGM from ambipolar diffusion and the decay of MHD turbulence. At lower redshifts, the formation of luminous objects and the build-up of a Lyman alpha background can be delayed by a redshift interval of 10 due to the strong increase of the filtering mass scale in the presence of magnetic fields. This tends to decrease the 21 cm signal compared to the zero-field case. In summary, we find that 21 cm observations may become a promising tool to constrain primordial magnetic fields.

Previous works have demonstrated that the generation of secondary CMB anisotropies due to the molecular optical depth is likely too small to be observed. In this paper, we examine additional ways in which primordial chemistry and the dark ages might influence the CMB. We present a detailed and updated chemical network and give an overview of the interactions of molecules with the CMB. We consider the optical depth due to line absorption, photoionization, photodissociation and free-free processes, and estimate the resulting changes in the CMB temperature and its power spectrum. The most promising results are obtained for the negative hydrogen ion \HM and the \HeHII molecule. The free-free process of \HM yields a relative change in the CMB temperature of up to $2\times10^{-11}$, and leads to a frequency-dependent change in the power spectrum of the order $10^{-7}$ at 30 GHz. With a change of the order $10^{-10}$ in the power spectrum, our result for the bound-free process of \HM is significantly below a previous suggestion. \HeHII efficiently scatters CMB photons and smears out primordial fluctuations, leading to a change in the power spectrum of the order $10^{-8}$. Improvements in the accuracy of future CMB experiments may thus help to constrain and finally detect these interesting signals from the dark ages of the universe.