The Micro-X sounding rocket flew for the first time on July 22, 2018, becoming the first program to fly Transition-Edge Sensors and multiplexing SQUID readout electronics in space. While a rocket pointing failure led to no time on-target, the success of the flight systems was demonstrated. The successful flight operation of the instrument puts the program in a position to modify the payload for indirect galactic dark matter searches. The payload modifications are motivated by the science requirements of this observation. Micro-X can achieve world-leading sensitivity in the keV regime with a single flight. Dark matter sensitivity projections have been updated to include recent observations and the expected sensitivity of Micro-X to these observed fluxes. If a signal is seen (as seen in the X-ray satellites), Micro-X can differentiate an atomic line from a dark matter signature.

Micro-X is a sounding rocket-borne instrument that uses a microcalorimeter array to perform high-resolution X-ray spectroscopy. This instrument flew for the first time on July 22nd, 2018 from the White Sands Missile Range, USA. This flight marks the first successful operation of a Transition-Edge Sensor array and its time division multiplexing read-out system in space. This launch was dedicated to the observation of the supernova remnant Cassiopeia A. A failure in the attitude control system prevented the rocket from pointing and led to no time on target. The on-board calibration source provided X-rays in flight, and it is used to compare detector performance during pre-flight integration, flight, and after the successful post-flight recovery. This calibration data demonstrates the capabilities of the detector in a space environment as well as its potential for future flights.

Observations of galaxies and primordial radiation suggest that the Universe is made mostly of non-luminous dark matter. Several types of new fundamental particles have been proposed as candidates for dark matter such as weakly interacting massive particles (WIMPs) but no definitive signal has been seen despite concerted efforts by many collaborations. One exception is the much-debated claim by the DAMA collaboration of a statistically significant annual modulation in the event rate of their experiment with a period and phase consistent with that expected from WIMP dark matter. Several groups have been working to develop experiments with the aim of reproducing DAMA's results using the same target medium. Here we report results from the initial operation of the COSINE-100 experiment. COSINE-100 uses sodium iodide as the target medium-the same medium as DAMA-and is designed to carry out a model-independent test of DAMA's claim. Initial data based on the first 59.5 days indicate that there is no excess of events over the expected background, confirming that DAMA's annual modulation signal is in severe tension with results from other experiments under the assumption of dark matter having spin independent interactions and the Standard Halo Model. COSINE-100 is now taking data to study the presence of dark matter-induced annual modulation in the event rate of the sodium iodide detectors.

COSINE-100 is a dark matter direct detection experiment designed to test the annual modulation signal observed by the DAMA/LIBRA experiment. COSINE-100 consists of 8 NaI(Tl) crystals with a total mass of 106 kg, a 2200 L liquid scintillator veto, and 37 muon detector panels. We present details of the data acquisition system of COSINE-100, including waveform storage using flash analog-to-digital converters for crystal events and integrated charge storage using charge-sensitive analog-to-digital converters for liquid scintillator and plastic scintillator muon veto events. We also discuss several trigger conditions developed in order to distinguish signal events from photomultiplier noise events. The total trigger rate observed for the crystal/liquid scintillator (plastic scintillator) detector is 15 Hz (24 Hz).

The COSINE-100 dark matter search experiment is an array of NaI(Tl) crystal detectors located in the Yangyang Underground Laboratory (Y2L). To understand measured backgrounds in the NaI(Tl) crystals we have performed Monte Carlo simulations using the Geant4 toolkit and developed background models for each crystal that consider contributions from both internal and external sources, including cosmogenic nuclides. The background models are based on comparisons of measurement data with Monte Carlo simulations that are guided by a campaign of material assays and are used to evaluate backgrounds and identify their sources. The average background level for the six crystals (70 kg total mass) that are studied is 3.5 counts/day/keV/kg in the (2-6) keV energy interval. The dominant contributors in this energy region are found to be $^{210}$Pb and $^3$H.

Meteoritical and astrophysical models of planet formation make contradictory predictions for dust concentration factors in chondrule forming regions of the solar nebula. Meteoritical and cosmochemical models strongly suggest that chondrules, a key component of the meteoritical record, formed in regions with solids-to-gas mass ratios orders of magnitude above background. However, models of dust grain dynamics in protoplanetary disks struggle to surpass factors of a few outside of very brief windows in the lifetime of the dust grains. Worse, those models do not predict significant concentration factors for dust grains the size of chondrule precursors. We briefly develop the difficulty in concentrating dust particles in the context of nebular chondrule formation and show that the disagreement is sufficiently stark that cosmochemists should explore ideas that might revise the concentration factor requirements downwards.

COSINE is a dark matter search experiment based on an array of low background NaI(Tl) crystals located at the Yangyang underground laboratory. The assembly of COSINE-100 was completed in the summer of 2016 and the detector is currently collecting physics quality data aimed at reproducing the DAMA/LIBRA experiment that reported an annual modulation signal. Stable operation has been achieved and will continue for at least two years. Here, we describe the design of COSINE-100, including the shielding arrangement, the configuration of the NaI(Tl) crystal detection elements, the veto systems, and the associated operational systems, and we show the current performance of the experiment.

The COSINE-100 dark matter search experiment has started taking physics data with the goal of performing an independent measurement of the annual modulation signal observed by DAMA/LIBRA. A muon detector was constructed by using plastic scintillator panels in the outermost layer of the shield surrounding the COSINE-100 detector. It is used to detect cosmic ray muons in order to understand the impact of the muon annual modulation on dark matter analysis. Assembly and initial performance test of each module have been performed at a ground laboratory. The installation of the detector in Yangyang Underground Laboratory (Y2L) was completed in the summer of 2016. Using three months of data, the muon underground flux was measured to be 328 $\pm$ 1(stat.)$\pm$ 10(syst.) muons/m$^2$/day. In this report, the assembly of the muon detector and the results from the analysis are presented.

Current protoplanetary dust coagulation theory does not predict dry silicate planetesimals, in tension with the Earth. While remedies to this predicament have been proposed, they have generally failed numerical studies, or are in tension with the Earth's (low, volatility dependent) volatile and moderately volatile elemental abundances. Expanding on the work of Boley et al. (2014), we examine the implications of molten grain collisions and find that they may provide a solution to the dry silicate planetesimal problem. Further, the source of the heating, be it the hot inner disk or an FU Orionis scale accretion event, would dictate the location of the resulting planetesimals, potentially controlling subsequent planetary system architectures. We hypothesize that systems which did undergo FU Orionis scale accretion events host planetary systems similar to our own, while ones that did not instead host very close in, tightly packed planets such as seen by Kepler.

We present the first search for a dark matter annual modulation signal in the Southern Hemisphere conducted with NaI(Tl) detectors, performed by the DM-Ice17 experiment. Nuclear recoils from dark matter interactions are expected to yield an annually modulated signal independent of location within the Earth's hemispheres. DM-Ice17, the first step in the DM-Ice experimental program, consists of 17 kg of NaI(Tl) located at the South Pole under 2200 m.w.e. overburden of Antarctic glacial ice. Taken over 3.6 years for a total exposure of 60.8 kg yr, DM-Ice17 data are consistent with no modulation in the energy range of 4-20 keV, providing the strongest limits on weakly interacting massive particle dark matter from a direct detection experiment located in the Southern Hemisphere. The successful deployment and stable long-term operation of DM-Ice17 establishes the South Pole ice as a viable location for future dark matter searches and in particular for a high-sensitivity NaI(Tl) dark matter experiment to directly test the DAMA/LIBRA claim of the observation of dark matter.

Ices, including water ice, prefer to recondense onto pre-existing nuclei rather than spontaneously forming grains from a cloud of vapor. Interestingly, different potential recondensation nuclei have very different propensities to actually nucleate water ice at the temperatures associated with freeze-out in protoplanetary discs. Therefore, if a region in a disc is warmed and then recooled, water vapor should not be expected to refreeze evenly onto all available grains. Instead it will preferentially recondense onto the most favorable grains. When the recooling is slow enough, only the most favorable grains will nucleate ice, allowing them to recondense thick ice mantles. We quantify the conditions for preferential recondensation to rapidly create pebble-sized grains in protoplanetary discs and show that FU Orionis type outbursts have the appropriate cooling rates to drive pebble creation in a band about 5 astronomical units wide outside of the quiescent frost line from approximately Jupiter's orbit to Saturn's (about 4 to 10 au). Those pebbles could be of the appropriate size to proceed to planetesimal formation via the Streaming Instability, or to contribute to the growth of planetesimals through pebble accretion. We suggest that this phenomenon contributed to the formation of the gas giants in our own Solar System.

An intriguing aspect of chondritic meteorites is that they are complementary: while their separate components have wildly varying abundances, bulk chondrites have nearly solar composition. This implies that the nearly-solar reservoirs in which chondrites were born were in turn assembled from sub-reservoirs of differing compositions that birthed the different components. We focus on explaining the potassium abundance variations between chondrules even within a single chondrite, while maintaining the observed CI $^{41}$K to $^{39}$K ratios. This requires physically separating potassium and chondrules while the temperature is high enough for K to be in the gas phase. We examine several mechanisms which could drive the dust through gas and show that to do so locally would have required long (sub-orbital to many orbits) time scales; with shortest potassium depletion time scales occurring in a scenario where chondrules formed high above the midplane and settled out of the evaporated potassium. While orbital time scales are at odds with laboratory chondrule cooling rate estimates, any other model for the origin for the potassium abundance variation has to wrestle with the severe logistical difficulty of generating a plethora of correlated reservoirs which varied strongly in their potassium abundances, but not in their potassium isotope ratios.

Recent studies of chondrites have found a tungsten isotopic anomaly between chondrules and matrix. Given the refractory nature of tungsten, this implies that W was carried into the solar nebula by at least two distinct families of pre-solar grains. The observed chondrule/matrix split requires that the distinct families were kept separate during the dust coagulation process, and that the two families of grain interacted with the chondrule formation mechanism differently. We take the co-existence of different families of solids in the same general orbital region at the chondrule-precursor size as given, and explore the requirements for them to have interacted with the chondrule formation process at significantly different rates. We show that this sorting of families of solids into chondrule and matrix destined dust had to have been at least as powerful a sorting mechanism as the relative settling of aerodynamically distinct grains at at least two scale heights above the midplane. The requirement that the chondrule formation mechanism was correlated in some fashion with a dust grain sorting mechanism argues strongly for spatially localized chondrule formation mechanisms such as turbulent dissipation in non-thermally ionized disk surface layers, and argues against volume filling mechanisms such as planetesimal bow shocks.

The meteoritical record shows both iron partitioning and tungsten isotopic partitioning between matrix and chondrules. Tungsten is not abundant enough to have driven its own isotopic partitioning, but if tungsten were correlated with iron, then ferromagnetic interactions grains could help explain both observations. We derive a practical parameterization for the increase in particle-particle collision rates caused by mutually attracting particle magnetic dipole moments. While the appropriate magnetic parameters remain uncertain, we show that ambient magnetic fields in protoplanetary disks are expected to be strong enough to magnetize iron metal bearing dust grains sufficiently to drive large increases in their collision rates. Such increased collision rates between iron metal rich grains could help preserve primordial iron and W isotopic inhomogeneities; and would help explain why the meteoritical record shows their partitioning in the solar nebula. The importance of magnetic interactions for larger grains whose growth is balanced by fragmentation is less clear, and will require future laboratory or numerical studies.

We report the measurement of muons and muon-induced phosphorescence in DM-Ice17, a NaI(Tl) direct detection dark matter experiment at the South Pole. Muon interactions in the crystal are identified by their observed pulse shape and large energy depositions. The measured muon rate in DM-Ice17 is 2.93 +/- 0.04 muons/crystal/day with a modulation amplitude of 12.3 +/- 1.7%, consistent with expectation. Following muon interactions, we observe long-lived phosphorescence in the NaI(Tl) crystals with a decay time of 5.5 +/- 0.5 s. The prompt energy deposited by a muon is correlated to the amount of delayed phosphorescence, the brightest of which consist of tens of millions of photons. These photons are distributed over tens of seconds with a rate and arrival timing that do not mimic a scintillation signal above 2 keVee. While the properties of phosphorescence vary among individual crystals, the annually-modulating signal observed by DAMA cannot be accounted for by phosphorescence with the characteristics observed in DM-Ice17.

Turbulence acting on mixes of gas and particles generally evenly diffuses the latter through the former. However, in the presence of background gas temperature gradients a phenomenon known as turbulent thermal diffusion appears as a particle drift velocity (rather than a diffusive term). This process moves particles from hot regions to cold ones. We rederive turbulent thermal diffusion using astrophysical language and demonstrate that it could play a major role in protoplanetary discs by concentrating particles by factors of tens. Such a concentration would set the stage for collective behavior such as the streaming instability and hence planetesimal formation.

We derive expressions for the photophoretic force on opaque spherical particles in a dilute gas in the optically thick regime where the radiation field is in local thermal equilibrium. Under those conditions, the radiation field has a simple form, leading to well defined analytical approximations for the photophoretic force that also consider both the internal thermal conduction within the particle, and the effects of heat conduction and radiation to the surrounding gas. We derive these results for homogeneous spherical particles; and for the double layered spheres appropriate for modeling solid grains with porous aggregate mantles. Then, as a specific astrophysical application of these general physical results, we explore the parameter space relevant to the photophoresis driven drift of dust in protoplanetary disks. We show that highly porous silicate grains have sufficiently low thermal conductivities that photophoretic effects, such as significant relative velocities between particles with differing porosity or levitation above the midplane, are expected to occur.

About 4-5% of chondrules are compound: two separate chondrules stuck together. This is commonly believed to be the result of the two component chondrules having collided shortly after forming, while still molten. This allows high velocity impacts to result in sticking. However, at T ~ 1100K, the temperature below which chondrules collide as solids (and hence usually bounce), coalescence times for droplets of appropriate composition are measured in tens of seconds. Even at 1025K, at which temperature theory predicts that the chondrules must have collided extremely slowly to have stuck together, the coalescence time scale is still less than an hour. These coalescence time scales are too short for the collision of molten chondrules to explain the observed frequency of compound chondrules. We suggest instead a scenario where chondrules stuck together in slow collisions while fully solid; and the resulting chondrule pair was subsequently briefly heated to a temperature in the range of 900-1025K. In that temperature window the coalescence time is finite but long, covering a span of hours to a decade. This is particularly interesting because those temperatures are precisely the critical window for thermally ionized MRI activity, so compound chondrules provide a possible probe into that vital regime.

Early in the study of viscous accretion disks it was realized that energy transfers from distant sources must be important, not least because the flow at the disk midplane in the bulk of the disk is likely outwards, out of the gravitational potential well. If the source of the viscosity is powered by accretion, such as in the case of the magneto-rotational instability, such distant energy sources must lie in the innermost regions of the disk, where accretion occurs even at the midplane. We argue here that modulations in this energy supply can alter the accretion rate on dynamical, rather than far longer viscous, time scales. This means that both the steady state value of and fluctuations in the inner disk's accretion rate, depending on the details of the inner boundary condition and occurring on the inner disk's rapid evolution time, can affect the outer disk. This is particularly interesting because observations have shown that disk accretion is not steady (e.g.~EX Lupi type objects). We also note that the power supplied to shearing boxes is set by the boxes themselves rather than the physical energy fluxes in a global disk. That is, their saturated magnetic field is not subject to the full set of energy constraints present in an actual disk. Our analysis suggests that large scale radial transport of energy has a critical impact on the evolution and variability of accretion disks.

In light of new results, the one-dimensional mean-field dynamo model of Brandenburg & Kapyla (2007) with dynamical quenching and a nonlocal Babcock-Leighton alpha effect is re-examined for the solar dynamo. We extend the one-dimensional model to include the effects of turbulent downward pumping (Kitchatinov & Olemskoy 2011), and to combine dynamical quenching with shear. We use both the conventional dynamical quenching model of Kleeorin & Ruzmaikin (1982) and the alternate one of Hubbard & Brandenburg (2011), and confirm that with varying levels of non-locality in the alpha effect, and possibly shear as well, the saturation field strength can be independent of the magnetic Reynolds number.

We report the first analysis of background data from DM-Ice17, a direct-detection dark matter experiment consisting of 17 kg of NaI(Tl) target material. It was codeployed with IceCube 2457 m deep in the South Pole glacial ice in December 2010 and is the first such detector operating in the Southern Hemisphere. The background rate in the 6.5 - 8.0 keVee region is measured to be 7.9 +/- 0.4 counts/day/keV/kg. This is consistent with the expected background from the detector assemblies with negligible contributions from the surrounding ice. The successful deployment and operation of DM-Ice17 establishes the South Pole ice as a viable location for future underground, low-background experiments in the Southern Hemisphere. The detector assembly and deployment are described here, as well as the analysis of the DM-Ice17 backgrounds based on data from the first two years of operation after commissioning, July 2011 - June 2013.

We consider the evidence presented by the LL3.0 chondrite Semarkona, including its chondrule fraction, chondrule size distribution and matrix thermal history. We show that no more than a modest fraction of the ambient matrix material in the Solar Nebula could have been melted into chondrules; and that much of the unprocessed matrix material must have been filtered out at some stage of Semarkona's parent body formation process. We conclude that agglomerations of many chondrules must have formed in the Solar Nebula, which implies that chondrules and matrix grains had quite different collisional sticking parameters. Further, we note that the absence of large melted objects in Semarkona means that chondrules must have exited the melting zone rapidly, before the chondrule agglomerations could form. The simplest explanation for this rapid exit is that chondrule melting occurred in surface layers of the disk. The newly formed, compact, chondrules then settled out of those layers on short time scales.

The magnetorotational instability (MRI) drives magnetized turbulence in sufficiently ionized regions of protoplanetary disks, leading to mass accretion. The dissipation of the potential energy associated with this accretion determines the thermal structure of accreting regions. Until recently, the heating from the turbulence has only been treated in an azimuthally averaged sense, neglecting local fluctuations. However, magnetized turbulence dissipates its energy intermittently in current sheet structures. We study this intermittent energy dissipation using high resolution numerical models including a treatment of radiative thermal diffusion in an optically thick regime. Our models predict that these turbulent current sheets drive order unity temperature variations even where the MRI is damped strongly by Ohmic resistivity. This implies that the current sheet structures where energy dissipation occurs must be well resolved to correctly capture the flow structure in numerical models. Higher resolutions are required to resolve energy dissipation than to resolve the magnetic field strength or accretion stresses. The temperature variations are large enough to have major consequences for mineral formation in disks, including melting chondrules, remelting calcium-aluminum rich inclusions, and annealing silicates; and may drive hysteresis: current sheets in MRI active regions could be significantly more conductive than the remainder of the disk.

In protoplanetary disks, dust grains rich in metallic iron can attract each other magnetically. If they are magnetized to values near saturation, the magnetically induced collision speeds are high enough to knock off the non-magnetized, loosely bound silicates. This process enriches the surviving portions of the dust grains in metallic iron, which further enhances the magnetically mediated collisions. The magnetic enhancement to the collisional cross-section between the iron rich dust results in rapid grain growth leading to planetesimal formation. While this process of knocking off silicates, which we term "magnetic erosion", occurs only in a very limited portion of a protoplanetary disk, it is a possible explanation for Mercury's disproportionately large iron core.

Blackman & Brandenburg argued that magnetic helicity conservation in dynamo theory can in principle be captured by diagrams of mean field dynamos when the magnetic fields are represented by ribbons or tubes, but not by lines. Here we present such a schematic ribbon diagram for the $\alpha^2$ dynamo that tracks magnetic helicity and provides distinct scales of large scale magnetic helicity, small scale magnetic helicity, and kinetic helicity involved in the process. This also motivates our construction of a new "2.5 scale" minimalist generalization of the helicity-evolving equations for the \alpha^2 dynamo that separately allows for these three distinct length scales while keeping only two dynamical equations. We solve these equations and, as in previous studies, find that the large scale field first grows at a rate independent of the magnetic Reynolds number R_M before quenching to an R_M dependent regime. But we also show that the larger the ratio of the wavenumber where the small scale current helicity resides to that of the forcing scale, the earlier the non-linear dynamo quenching occurs, and the weaker the large scale field is at the turnoff from linear growth. The harmony between the theory and the schematic diagram exemplifies a general lesson that magnetic fields in MHD are better visualized as two-dimensional ribbons (or pairs of lines) rather than single lines.

The Earth is known to be depleted in volatile lithophile elements in a fashion that defies easy explanation. We resolve this anomaly with a model that combines the porosity of collisionally grown dust grains in protoplanetary disks with heating from FU Orionis events that dramatically raise protoplanetary disk temperatures. The heating from an FU Orionis event alters the aerodynamical properties of the dust while evaporating the volatiles. This causes the dust to settle, abandoning those volatiles. The success of this model in explaining the elemental composition of the Earth is a strong argument in favor of highly porous collisionally grown dust grains in protoplanetary disks outside our Solar System. Further, it demonstrates how thermal (or condensation based) alterations of dust porosity, and hence aerodynamics, can be a strong factor in planet formation, leading to the onset of rapid gravitational instabilities in the dust disk and the subsequent collapse that forms planetesimals.

This letter expands the stability criterion for radially stratified, vertically {unstratified} accretion disks incorporating thermal relaxation. We find a linear amplification of epicyclic oscillations in these disks that depends on the effective cooling time, i.e. an overstability. The growth rates of the overstability vanish for both extreme cases, e.g. infinite cooling time and instantaneous cooling, i.e. the adiabatic and fully isothermal cases. However, for thermal relaxation times $\tau$ on the order of the orbital frequency, $\tau\Omega \sim 1$, modes grow at a rate proportional to the square of the Brunt-V\"ais\"al\"a frequency. The overstability is based on epicyclic motions, with the thermal relaxation causing gas to heat while radially displaced inwards, and cool while radially displaced outwards. This causes the gas to have a lower density when moving outwards compared to when it moves inwards, so it feels the outwards directed pressure force more strongly on that leg of the journey. We suggest the term ``Convective Overstability" for the phenomenon that has already been numerically studied in the non-linear regime in the context of amplifying vortices in disks, under the name ``Subcritical Baroclinic Instability". The point of the present paper is to make clear that vortex formation in three-dimensional disks is neither subcritical, i.e. does not need a finite perturbation, nor is it baroclinic in the sense of geophysical fluid dynamics, which requires on vertical shear. We find that Convective Overstability is a linear instability that will operate under a wide range of physical conditions for circumstellar disks.

Meteoritic chondrules were formed in the early solar system by brief heating of silicate dust to melting temperatures. Some highly refractory grains (Type B calcium-aluminum-rich inclusions, CAIs) also show signs of transient heating. A similar process may occur in other protoplanetary disks, as evidenced by observations of spectra characteristic of crystalline silicates. One possible environment for this process is the turbulent magnetohydrodynamic flow thought to drive accretion in these disks. Such flows generally form thin current sheets, which are sites of magnetic reconnection, and dissipate the magnetic fields amplified by a disk dynamo. We suggest that it is possible to heat precursor grains for chondrules and other high-temperature minerals in current sheets that have been concentrated by our recently described short-circuit instability. We extend our work on this process by including the effects of radiative cooling, taking into account the temperature dependence of the opacity; and by examining current sheet geometry in three-dimensional, global models of magnetorotational instability. We find that temperatures above 1600 K can be reached for favorable parameters that match the ideal global models. This mechanism could provide an efficient means of tapping the gravitational potential energy of the protoplanetary disk to heat grains strongly enough to form high-temperature minerals. The volume-filling nature of turbulent magnetic reconnection is compatible with constraints from chondrule-matrix complementarity, chondrule-chondrule complementarity, the occurrence of igneous rims, and compound chondrules. The same short-circuit mechanism may perform other high-temperature mineral processing in protoplanetary disks such as the production of crystalline silicates and CAIs.

We study the collision rates and velocities for point-particles of different sizes in turbulent flows. We construct fits for the collision rates at specified velocities (effectively a collisional velocity probability distribution) for particle stopping time ratios up to four; already by that point the collisional partners are very poorly correlated and so the results should be robust for even larger stopping time ratios. Significantly, we find that while particles of very different masses have approximately Maxwellian collisional statistics, as the mass ratio shrinks the distribution changes dramatically. At small stopping time ratios, the collisional partners are highly correlated and we find a population of high number density (clustered), low relative-velocity particle pairs. Unlike in the case of identical stopping time collisional partners, this low relative-velocity clustered population is collisional, but the clustering is barely adequate to trigger bulk effects such as the streaming instability. We conclude our analysis by constructing a master fit to the collisional statistics as a function only of the stopping time ratio. Together with our previous work for identical stopping time particle pairs, this provides a recipe for including collisional velocity probability distributions in dust coagulation models for protoplanetary disks. We also include our recipe for determining particle collisional diagnostics from numerical simulations.

Jets from active galaxies propagate from the central black hole out to the radio lobes on scales of hundreds of kiloparsecs. The jets may encounter giant stars with strong stellar winds and produce observable signatures. For strong winds and weak jets, the interaction may truncate the jet flow during its transit via the mass loading. For weaker jets, the interaction can produce knots in the jet. We present recent 3DMHD numerical simulations to model the evolution of this jet-wind interaction and its observational consequences. We explore (i) the relative mechanical luminosity of the radio jets and the stellar winds (ii) the impact parameter between the jets' axis and the stellar orbital path (iii) the relative magnetic field strength of the jets and the stellar winds.

Many astrophysical systems of interest, including protoplanetary accretion disks, are made of turbu- lent magnetized gas with near solar metallicity. Thermal ionization of alkali metals in such gas exceeds non-thermal ionization when temperatures climb above roughly 1000 K. As a result, the conductiv- ity, proportional to the ionization fraction, gains a strong, positive dependence on temperature. In this paper, we demonstrate that this relation between the temperature and the conductivity triggers an exponential instability that acts similarly to an electrical short, where the increased conductivity concentrates the current and locally increases the Ohmic heating. This contrasts with the resistiv- ity increase expected in an ideal magnetic reconnection region. The instability acts to focus narrow current sheets into even narrower sheets with far higher currents and temparatures. We lay out the basic principles of this behavior in this paper using protoplanetary disks as our example host system, motivated by observations of chondritic meteorites and their ancestors, dust grains in protoplanetary disks, that reveal the existence of strong, frequent heating events that this instability could explain.

To understand the earliest stages of planet formation, it is crucial to be able to predict the rate and the outcome of dust grains collisions, be it sticking and growth, bouncing, or fragmentation. The outcome of such collisions depends on the collision speed, so we need a solid understanding of the rate and velocity distribution of turbulence-induced dust grain collisions. The rate of the collisions depends both on the speed of the collisions and the degree of clustering experienced by the dust grains, which is a known outcome of turbulence. We evolve the motion of dust grains in simulated turbulence, an approach that allows a large turbulent inertial range making it possible to investigate the effect of turbulence on meso-scale grains (millimeter and centimeter). We find three populations of dust grains: one highly clustered, cold and collisionless; one warm; and the third "hot". Our results can be fit by a simple formula, and predict both significantly slower typical collisional velocities for a given turbulent strength than previously considered, and modest effective clustering of the collisional populations, easing difficulties associated with bouncing and fragmentation barriers to dust grain growth. Nonetheless, the rate of high velocity collisions falls off merely exponentially with relative velocity so some mid- or high-velocity collisions will still occur, promising some fragmentation.

At large magnetic Reynolds numbers, magnetic helicity evolution plays an important role in astrophysical large-scale dynamos. The recognition of this fact led to the development of the dynamical alpha quenching formalism, which predicts catastrophically low mean fields in open systems. Here we show that in oscillatory alpha-Omega dynamos this formalism predicts an unphysical magnetic helicity transfer between scales. An alternative technique is proposed where this artifact is removed by using the evolution equation for the magnetic helicity of the total field in the shearing--advective gauge. In the traditional dynamical alpha quenching formalism, this can be described as an additional magnetic helicity flux of small-scale fields that does not appear in homogeneous alpha-squared dynamos. In alpha-Omega dynamos, the alternative formalism is shown to lead to larger saturation fields than previously obtained with the traditional formalism.

Context. Helically forced magneto-hydrodynamic shearing-sheet turbulence can support different large-scale dynamo modes, although the {\alpha} - {\Omega} mode is generally expected to dominate because it is the fastest growing. In an {\alpha} - {\Omega} dynamo, most of the field amplification is produced by the shear. As differential rotation is an ubiquitous source of shear in astrophysics, such dynamos are believed to be the source of most astrophysical large-scale magnetic fields. Aims. We study the stability of oscillatory migratory {\alpha} - {\Omega} type dynamos in turbulence simulations. Methods. We use shearing-sheet simulations of hydromagnetic turbulence that is helically forced at a wavenumber that is about three times larger than the lowest wavenumber in the domain so that both {\alpha} - {\Omega} and {\alpha}2 dynamo action is possible. Results. After initial dominance and saturation, the {\alpha} - {\Omega} mode is found to be destroyed by an orthogonal {\alpha}2 mode sustained by the helical turbulence alone. We show that there are at least two processes through which this transition can occur. Conclusions. The fratricide of {\alpha} - {\Omega} dynamos by its {\alpha}2 sibling is discussed in the context of grand minima of solar and stellar activity. However, the genesis of {\alpha} - {\Omega} dynamos from an {\alpha}2 dynamo has not yet been found.

Magnetic helicity fluxes are investigated in a family of gauges in which the contribution from ideal magnetohydrodynamics takes the form of a purely advective flux. Numerical simulations of magnetohydrodynamic turbulence in this advective gauge family exhibit instabilities triggered by the build-up of unphysical irrotational contributions to the magnetic vector potential. As a remedy, the vector potential is evolved in a numerically well behaved gauge, from which the advective vector potential is obtained by a gauge transformation. In the kinematic regime, the magnetic helicity density evolves similarly to a passive scalar when resistivity is small and turbulent mixing is mild, i.e. when the fluid Reynolds number is not too large. In the dynamical regime, resistive contributions to the magnetic helicity flux in the advective gauge are found to be significant owing to the development of small length scales in the irrotational part of the magnetic vector potential.

Magnetic helicity has risen to be a major player in dynamo theory, with the helicity of the small-scale field being linked to the dynamo saturation process for the large-scale field. It is a nearly conserved quantity, which allows its evolution equation to be written in terms of production and flux terms. The flux term can be decomposed in a variety of fashions. One particular contribution that has been expected to play a significant role in dynamos in the presence of mean shear was isolated by Vishniac & Cho (2001, ApJ 550, 752). Magnetic helicity fluxes are explicitly gauge dependent however, and the correlations that have come to be called the Vishniac-Cho flux were determined in the Coulomb gauge, which turns out to be fraught with complications in shearing systems. While the fluxes of small-scale helicity are explicitly gauge dependent, their divergences can be gauge independent. We use this property to investigate magnetic helicity fluxes of small-scale field through direct numerical simulations in a shearing-box system and find that in a numerically usable gauge the divergence of the small-scale helicity flux vanishes, while the divergence of the Vishniac-Cho flux remains finite. We attribute this seeming contradiction to the existence of horizontal fluxes of small-scale magnetic helicity with finite divergences even in our shearing-periodic domain.

The predictive power of mean-field theory is emphasized by comparing theory with simulations under controlled conditions. The recently developed test-field method is used to extract turbulent transport coefficients both in kinematic as well as nonlinear and quasi-kinematic cases. A striking example of the quasi-kinematic method is provided by magnetic buoyancy-driven flows that produce an alpha effect and turbulent diffusion.

We present the results of simulations of forced turbulence in a slab where the mean kinetic helicity has a maximum near the mid-plane, generating gradients of magnetic helicity of both large and small-scale fields. We also study systems that have poorly conducting buffer zones away from the midplane in order to assess the effects of boundaries. The dynamical alpha quenching phenomenology requires that the magnetic helicity in the small-scale fields approaches a nearly static, gauge independent state. To stress-test this steady state condition we choose a system with a uniform sign of kinetic helicity, so that the total magnetic helicity can reach a steady state value only through fluxes through the boundary, which are themselves suppressed by the velocity boundary conditions. Even with such a set up, the small-scale magnetic helicity is found to reach a steady state. In agreement with earlier work, the magnetic helicity fluxes of small-scale fields are found to be turbulently diffusive. By comparing results with and without halos, we show that artificial constraints on magnetic helicity at the boundary do not have a significant impact on the evolution of the magnetic helicity, except that "softer" (halo) boundary conditions give a lower energy of the saturated mean magnetic field.

Estimates for the nonlinear alpha effect in helical turbulence with an applied magnetic field are presented using two different approaches: the imposed-field method where the electromotive force owing to the applied field is used, and the test-field method where separate evolution equations are solved for a set of different test fields. Both approaches agree for stronger fields, but there are apparent discrepancies for weaker fields that can be explained by the influence of dynamo-generated magnetic fields on the scale of the domain that are referred to as meso-scale magnetic fields. Examples are discussed where these meso-scale fields can lead to both drastically overestimated and underestimated values of alpha compared with the kinematic case. It is demonstrated that the kinematic value can be recovered by resetting the fluctuating magnetic field to zero in regular time intervals. It is concluded that this is the preferred technique both for the imposed-field and the test-field methods.

In the mean-field theory of magnetic fields, turbulent transport, i.e. the turbulent electromotive force, is described by a combination of the alpha effect and turbulent magnetic diffusion, which are usually assumed to be proportional respectively to the mean field and its spatial derivatives. For a passive scalar there is just turbulent diffusion, where the mean flux of concentration depends on the gradient of the mean concentration. However, these proportionalities are approximations that are valid only if the mean field or the mean concentration vary slowly in time. Examples are presented where turbulent transport possesses memory, i.e. where it depends crucially on the past history of the mean field. Such effects are captured by replacing turbulent transport coefficients with time integral kernels, resulting in transport coefficients that depend effectively on the frequency or the growth rate of the mean field itself. In this paper we perform numerical experiments to find the characteristic timescale (or memory length) of this effect as well as simple analytical models of the integral kernels in the case of passive scalar concentrations and kinematic dynamos. The integral kernels can then be used to find self-consistent growth or decay rates of the mean fields. In mean-field dynamos the growth rates and cycle periods based on steady state values of alpha effect and turbulent diffusivity can be quite different from the actual values.

Turbulent viscosity is frequently used in accretion disk theory to replace the microphysical viscosity in order to accomodate the observational need for in- stabilities in disks that lead to enhanced transport. However, simply replacing the microphysical transport coefficient by a single turbulent transport coeffi- cient hides the fact that the procedure should formally arise as part of a closure in which the hydrodynamic or magnetohydrodynamic equations are averaged, and correlations of turbulent fluctuations are replaced by transport coefficients. Here we show how a mean field approach leads quite naturally two transport coefficients, not one, that govern mass and angular momentum transport. In particular, we highlight that the conventional approach suffers from a seemingly inconsistent neglect of turbulent diffusion in the surface density equation. We constrain these new transport coefficients for specific cases of inward, outward, and zero net mass transport. In addition, we find that one of the new transport terms can lead to oscillations in the mean surface density which then requires a constant or small inverse Rossby number for disks to maintain a monotonic power-law surface density.

It is difficult to imagine a planet formation model that does not at some stage include a gravitationally unstable disc. Initially unstable gas-dust discs may form planets directly, but the high surface density required has motivated the alternative that gravitational instability occurs in a dust sub-layer only after grains have grown large enough by electrostatic sticking. Although such growth up to the instability stage is efficient for laminar discs, concern has mounted as to whether realistic disc turbulence catastrophically increases the settling time, thereby requiring additional processes to facilitate planet formation on the needed time scales. To evaluate this concern, we develop a model for grain growth that accounts for the influence of turbulence on the collisional velocity of grains and on the scale height of the dust layer. The relative effect on these quantities depends on the grain size. The model produces a disc-radius dependent time scale to reach a gravitationally unstable phase of planet formation. For a range of dust sticking and disc parameters, we find that for viscosity parameters $\alpha \le 10^{-3}$, this time scale is short enough over a significant range in radii $R$ that turbulence does not catastrophically slow the early phases of planet formation, even in the absence of agglomeration enhancement agents like vortices.

Active Galactic Nuclei can produce extremely powerful jets. While tightly collimated, the scale of these jets and the stellar density at galactic centers implies that there will be many jet/star interactions, which can mass-load the jet through stellar winds. Previous work employed modest wind mass outflow rates, but this does not apply when mass loading is provided by a small number of high mass-loss stars. We construct a framework for jet mass-loading by stellar winds for a broader spectrum of wind mass-loss rates than has been previously considered. Given the observed stellar mass distributions in galactic centers, we find that even highly efficient (0.1 Eddington luminosity) jets from supermassive black holes of masses $M_{BH} \la 10^4M_{\odot}$ are rapidly mass loaded and quenched by stellar winds. For $10^4 M_{\odot}<M_{BH}<10^8 M_{\odot}$, the quenching length of highly efficient jets is independent of the jet's mechanical luminosity. Stellar wind mass-loading is unable to quench efficient jets from more massive engines, but can account for the observed truncation of the inefficient M87 jet, and implies a baryon dominated composition on scales $\ga 2$ kpc therein even if the jet is initially pair plasma dominated.

We consider the possibility that an infalling bulge core or stellar cluster could form an eccentric disk, following tidal disruption by a massive black hole in the center of a galaxy. As a function of central black hole mass, we constrain the core radii and central densities of cluster progenitors capable of becoming nearly Keplerian disks which can support lopsided slow modes. We find that progenitor stellar clusters with core radii less than a pc and densities above a few times 10^5 Msol/pc^3 are likely eccentric disk progenitors near a massive black hole of mass 10^7 to 10^8 Msol. Lower density and larger progenitor cores are capable of causing eccentric stellar disks near more massive black holes. Our constraints on the progenitor cores are consistent with existing N-body simulations, which in one case has produced an eccentric disk. For M31 and NGC 4486B, the estimated progenitor cluster cores are dense and compact compared to Galactic globular clusters, however the cores of nearby galaxies such as M33, M32 and M31 itself are in the right regime. If galaxy mergers can create eccentric disks, then they would be a natural consequence of hierarchical galaxy formation.