In simulations of viscously evolving accretion disks, the inner boundary condition is particularly important. If treated incorrectly, it induces incorrect behavior very quickly, because the viscous time is shortest near the inner boundary. Recent work has determined the correct inner boundary in Eulerian simulations. But in quasi-Lagrangian simulations (e.g., SPH, moving mesh, and mesh-less), where the inner boundary is modeled by removing mass within a finite zone, the inner density profile typically becomes anomalously depleted. Here we show how the boundary condition should be applied in such codes, via a simple modification of the usual approach: when one removes mass, one must speed up the remaining material so that the disk's angular momentum is unchanged. We show with both 1D and 2D moving-mesh (AREPO) simulations that this scheme works as desired in viscously evolving disks. It produces no spurious density depletions and is independent of the mass removal rate, provided that the disk is adequately resolved and that the mass removal rate is not so extreme as to trigger instabilities. This "torque-free" mass removal technique permits the use of quasi-Lagrangian codes to simulate viscously evolving disks, while including a variety of additional effects. As an example, we apply our scheme to a 2D simulation of an accretion disk perturbed by a very massive planet, in which the disk is evolved to viscous steady state.

Characterizing the surface material of an asteroid is important for understanding its geology and for informing mission decisions, such as the selection of a sample site. Diurnal surface temperature amplitudes are directly related to the thermal properties of the materials on the surface. We describe a numerical model for studying the thermal conductivity of particulate regolith in vacuum. Heat diffusion and surface-to-surface radiation calculations are performed using the finite element (FE) method in three-dimensional meshed geometries of randomly packed spherical particles. We validate the model for test cases where the total solid and radiative conductivity values of particulates with monodisperse particle size frequency distributions (SFDs) are determined at steady-state thermal conditions. Then, we use the model to study the bulk radiative thermal conductivity of particulates with polydisperse, cumulative power-law particle SFDs. We show that for each polydisperse particulate geometry tested, there is a corresponding monodisperse geometry with some effective particle diameter that has an identical radiative thermal conductivity. These effective diameters are found to correspond very well to the Sauter mean particle diameter, which is essentially the surface-area-weighted mean. Next, we show that the thermal conductivity of the particle material can have an important effect on the radiative component of the thermal conductivity of particulates, especially if the particle material conductivity is very low or the spheres are relatively large, owing to non-isothermality in each particle. We provide an empirical correlation to predict the effects of non-isothermality on radiative thermal conductivity in both monodisperse and polydisperse particulates.

We carry out 2D viscous hydrodynamics simulations of circumbinary disk (CBD) accretion using {\footnotesize AREPO}. We resolve the accretion flow from a large-scale CBD down to the streamers and disks around individual binary components. Extending our recent studies \citep{mun19}, we consider circular binaries with various mass ratios ($0.1\leq q_{\rm{b}}\leq1$) and study accretion from ``infinite'', steady-supply disks and from finite-sized, viscously spreading tori. For ``infinite'' disks, a global steady state can be reached, and the accretion variability has a dominant frequency ${\sim}0.2\Omega_{\rm{b}}$ for $q_{\rm{b}}>0.5$ and $\Omega_{\rm{b}}$ for $q_{\rm{b}}<0.5$, ($\Omega_{\rm{b}}$ is the binary angular frequency). We find that the accretion ``eigenvalue'' $l_0$ -- the net angular momentum transfer from the disk to the binary per unit accreted mass -- is always positive and falls in the range ($0.65$-$0.85)a_{\rm b}^2\Omega_{\rm{b}}$ (with $a_{\rm{b}}$ the binary separation), depending weakly on the mass ratio and viscosity. This leads to binary expansion when $q_{\rm{b}}\gtrsim0.3$. Accretion from a finite torus can be separated into two phases: an initial transient phase, corresponding to the filling of the binary cavity, followed by a viscous pseudo-stationary phase, during which the torus viscously spreads and accretes onto the binary. In the viscous phase, the net torque on the binary per unit accreted mass is close to $l_0$, the value derived for ``infinite'' disks. We conclude that similar-mass binaries accreting from CBDs gain angular momentum and expand over long time scales. This result significantly impacts the coalescence of supermassive binary black holes and newly formed binary stars. We offer a word of caution against conclusions drawn from simulations of transient accretion onto empty circumbinary cavities.

We carry out 2D viscous hydrodynamical simulations of circumbinary accretion using the AREPO code. We self-consistently compute the accretion flow over a wide range of spatial scales, from the circumbinary disk (CBD) far from the central binary, through accretion streamers, to the disks around individual binary components, resolving the flow down to 2% of the binary separation. We focus on equal-mass binaries with arbitrary eccentricities. We evolve the flow over long (viscous) timescales until a quasi-steady is reached, in which the mass supply rate at large distances $\dot{M}_0$ (assumed constant) equals the time-averaged mass transfer rate across the disk and the total mass accretion rate onto the binary components. This quasi-steady state allows us to compute the secular angular momentum transfer rate onto the binary, $\langle\dot{J}_b\rangle$, and the resulting orbital evolution. Through direct computation of the gravitational and accretion torques on the binary, we find that $\langle\dot{J}_b\rangle$ is consistently positive (i.e., the binary gains angular momentum), with $l_0\equiv\langle\dot{J}_b\rangle/\dot M_0$ in the range of $(0.4-0.8)a_b^2\Omega_b$, depending on the binary eccentricity (where $a_b,~\Omega_b$ are the binary semi-major axis and angular frequency); we also find that this $\langle\dot{J}_b\rangle$ is equal to the net angular momentum current across the CBD, indicating that global angular momentum balance is achieved in our simulations. We compute the time-averaged rate of change of the binary orbital energy for eccentric binaries, and thus obtain the secular rates $\langle\dot a_b\rangle$ and $\langle \dot{e}_b\rangle$. In all cases, $\langle\dot{a}_b\rangle$ is positive, i.e., the binary expands while accreting. We discuss the implications of our results for the merger of supermassive binary black holes and for the formation of close stellar binaries.

Important clues on the formation and evolution of planetary systems can be inferred from the stellar obliquity $\psi$. We study the distribution of obliquities using the California-Kepler Survey and the TEPCat Catalog of Rossiter-MacLaughlin (RM) measurements, from which we extract, respectively, 275 and 118 targets. We infer a "best fit" obliquity distribution in $\psi$ with a single parameter $\kappa$. Large values of $\kappa$ imply that $\psi$ is distributed narrowly around zero, while small values imply approximate isotropy. Our findings are: (1) The distribution of $\psi$ in Kepler systems is narrower than found by previous studies and consistent with $\kappa\sim15$ (mean $\langle\psi\rangle\sim19^\circ$ and spread $\psi\sim10^\circ$). (2) The value of $\kappa$ in Kepler systems does not depend, at a statistically significant level, on planet multiplicity, stellar multiplicity or stellar age; on the other hand, metal-rich hosts, small planet hosts and long-period planet hosts tend to be more oblique than the general sample (at a $\sim$2.5-$\sigma$ significance level). (3) Hot Jupiter (HJ) systems with RM measurements are consistent with $\kappa\sim2$, more broadly distributed than the general Kepler population. (4) A separation of the RM sample into cooler ($T_{\rm eff}\lesssim$6250 K) and hotter ($T_{\rm eff}\gtrsim$6250 K) HJ hosts results in two distinct distributions, $\kappa_{\rm cooler}\sim4$ and $\kappa_{\rm hotter}\sim1$ (4-$\sigma$ significance), both more oblique than the Kepler sample. We hypothesize that the total mass in planets may be behind the increasing obliquity with metallicity and planet radius, and that the period dependence could be due to primordial disk alignment rather than realignment of stellar spin.

Small bipolar magnetic features are observed to appear in the interior of individual granules in the quiet Sun, signaling the emergence of tiny magnetic loops from the solar interior. We study the origin of those features as part of the magnetoconvection process in the top layers of the convection zone. Two quiet-Sun magnetoconvection models, calculated with the radiation-magnetohydrodynamic (MHD) Bifrost code and with domain stretching from the top layers of the convection zone to the corona, are analyzed. Using 3D visualization as well as a posteriori spectral synthesis of Stokes parameters, we detect the repeated emergence of small magnetic elements in the interior of granules, as in the observations. Additionally, we identify the formation of organized horizontal magnetic sheets covering whole granules. Our approach is twofold, calculating statistical properties of the system, like joint probability density functions (JPDFs), and pursuing individual events via visualization tools. We conclude that the small magnetic loops surfacing within individual granules in the observations may originate from sites at or near the downflows in the granular and mesogranular levels, probably in the first 1 or 1.5 Mm below the surface. We also document the creation of granule-covering magnetic sheet-like structures through the sideways expansion of a small subphotospheric magnetic concentration picked up, and pulled out of the interior, by a nascent granule. The sheet-like structures we found in the models may match the recent observations of Centeno et al. (2017).

We present two-dimensional stellar and gaseous kinematics of the inner 0.7 $\times$ 1.2 kpc$^{2}$ of the Seyfert galaxy ESO 362-G18, derived from optical spectra obtained with the GMOS/IFU on the Gemini South telescope at a spatial resolution of $\approx$170 pc and spectral resolution of 36 km s$^{-1}$. ESO 362-G18 is a strongly perturbed galaxy of morphological type Sa or S0/a, with a minor merger approaching along the NE direction. Previous studies have shown that the [OIII] emission shows a fan-shaped extension of $\approx$ 10\arcsec\ to the SE. We detect the [OIII] doublet, [NII] and H${\alpha}$ emission lines throughout our field of view. The stellar kinematics is dominated by circular motions in the galaxy plane, with a kinematic position angle of $\approx$137$^{\circ}$. The gas kinematics is also dominated by rotation, with kinematic position angles ranging from 122$^{\circ}$ to 139$^{\circ}$. A double-Gaussian fit to the [OIII]$\lambda$5007 and H${\alpha}$ lines, which have the highest signal to noise ratios of the emission lines, reveal two kinematic components: (1) a component at lower radial velocities which we interpret as gas rotating in the galactic disk; and (2) a component with line of sight velocities 100-250 km s$^{-1}$ higher than the systemic velocity, interpreted as originating in the outflowing gas within the AGN ionization cone. We estimate a mass outflow rate of 7.4 $\times$ 10$^{-2}$ M$_{\odot}$ yr$^{-1}$ in the SE ionization cone (this rate doubles if we assume a biconical configuration), and a mass accretion rate on the supermassive black hole (SMBH) of 2.2 $\times$ 10$^{-2}$ M$_{\odot}$ yr$^{-1}$. The total ionized gas mass within $\sim$84 pc of the nucleus is 3.3 $\times$ 10$^{5}$ M$_{\odot}$; infall velocities of $\sim$34 km s$^{-1}$ in this gas would be required to feed both the outflow and SMBH accretion.

We carry out numerical simulations of circumbinary discs, solving the viscous hydrodynamics equations on a polar grid covering an extended disc outside the binary co-orbital region. We use carefully controlled outer boundary conditions and long-term integrations to ensure that the disc reaches a quasi-steady state, in which the time-averaged mass accretion rate onto the binary, $\langle\dot{M}\rangle$, matches the mass supply rate at the outer disc. We focus on binaries with comparable masses and a wide range of eccentricities ($e_\mathrm{B}$). For $e_\mathrm{B} \lesssim 0.05$, the mass accretion rate of the binary is modulated at about $5$ times the binary period; otherwise it is modulated at the binary period. The inner part of the circumbinary disc ($r \lesssim 6 a_\mathrm{B}$) generally becomes coherently eccentric. For low and high $e_\mathrm{B}$, the disc line of apsides precesses around the binary, but for intermediate $e_\mathrm{B}$ ($0.2 - 0.4$), it instead becomes locked with that of the binary. By considering the balance of angular momentum transport through the disc by advection, viscous stress, and gravitational torque, we determine the time-averaged net angular momentum transfer rate to the binary, $\langle\dot{J}\rangle$. The specific angular momentum, $l_0 = \langle\dot{J}\rangle/\langle\dot{M}\rangle$, depends non-monotonically on $e_\mathrm{B}$. Contrary to previous claims, we find that $l_0$ is positive for most $e_\mathrm{B}$, implying that the binary receives net angular momentum, which may cause its separation to grow with time. The minimum $l_0$ occurs at intermediate $e_\mathrm{B}$ ($0.2 - 0.4$), corresponding to the regime where the inner eccentric disc is apsidally aligned with the binary.

The presence of a planetary system can shield a planetesimal disk from the secular gravitational perturbations due to distant outer massive objects (planets or stellar companions). As the host star evolves off the main sequence to become a white dwarf, these planets can be engulfed, triggering secular instabilities and leading to the tidal disruptions of small rocky bodies. These disrupted bodies can feed the white dwarfs with rocky material and possibly explain the high-metallicity material in their atmospheres. We illustrate how this mechanism can operate when the gravitational perturbations are due to the Kozai-Lidov mechanism from a stellar binary companion. We show that this mechanism can explain the observed levels of accretion if: (1) the planetary engulfment happens fast compared to the secular timescale, which is generally the case for wide binaries ($>100$ AU) and planetary engulfment during the Asymptotic Giant Branch; (2) the planetesimal disk has a total mass of $\sim10^{-4}-10^{-2}M_\oplus$. We show that this new mechanism can provide a steady supply of material throughout the entire life of the white dwarfs for all cooling ages and can account for a large fraction (up to nearly half) of the observed polluted WDs.

We present numerical simulations of circumbinary accretion onto eccentric and circular binaries using the moving-mesh code AREPO. This is the first set of simulations to tackle the problem of binary accretion using a finite-volume scheme on a freely moving mesh, which allows for accurate measurements of accretion onto individual stars for arbitrary binary eccentricity. While accretion onto a circular binary shows bursts with period of ~5 times the binary period P_b,accretion onto an eccentric binary is predominantly modulated at the period ~1P_b. For an equal-mass circular binary, the accretion rates onto individual stars are quite similar to each other, following the same variable pattern in time. By contrast, for eccentric binaries, one of the stars can accrete at a rate 10-20 times larger than its companion. This "symmetry breaking" between the stars, however, alternates over timescales of order 200 P_b, and can be attributed to a slowly precessing, eccentric circumbinary disk. Over longer timescales, the net accretion rates onto individual stars are the same, reaching a quasi-steady state with the circumbinary disk. These results have important implications for the accretion behavior of binary T-Tauri stars and supermassive binary black holes.

Lidov-Kozai oscillations of planets in stellar binaries, combined with tidal dissipation, can lead to the formation of hot Jupiters (HJs) or tidal disruption of planets. Recent population synthesis studies have found that the fraction of systems resulting in HJs (F_HJ) depends strongly on the planet mass, host stellar type and tidal dissipation strength, while the total migration fraction F_mig = F_ HJ + F_dis (including both HJ formation and tidal disruption) exhibits much weaker dependence. We present an analytical method for calculating F_HJ and F_mig in the Lidov-Kozai migration scenario. The key ingredient of our method is to determine the critical initial planet-binary inclination angle that drives the planet to reach sufficiently large eccentricity for efficient tidal dissipation or disruption. This calculation includes the effects of the octupole potential and short-range forces on the planet. Our analytical method reproduces the resulting planet migration/disruption fractions from population synthesis, and can be easily implemented for various planet, stellar/companion types, and for different distributions of initial planetary semi-major axes, binary separations and eccentricities. We extend our calculations to planets in the super-Earth mass range and discuss the conditions for such planets to survive Lidov-Kozai migration and form close-in rocky planets.

The discovery of transiting circumbinary planets by the Kepler mission suggests that planets can form efficiently around binary stars. None of the stellar binaries currently known to host planets has a period shorter than 7 days, despite the large number of eclipsing binaries found in the Kepler target list with periods shorter than a few days. These compact binaries are believed to have evolved from wider orbits into their current configurations via the so-called Lidov-Kozai migration mechanism, in which gravitational perturbations from a distant tertiary companion induce large-amplitude eccentricity oscillations in the binary, followed by orbital decay and circularization due to tidal dissipation in the stars. Here we explore the orbital evolution of planets around binaries undergoing orbital decay by this mechanism. We show that planets may survive and become misaligned from their host binary, or may develop erratic behavior in eccentricity, resulting in their consumption by the stars or ejection from the system as the binary decays. Our results suggest that circumbinary planets around compact binaries could still exist, and we offer predictions as to what their orbital configurations should be like.

The Lidov-Kozai (LK) mechanism plays an important role in the secular evolution of many hierarchical triple systems. The standard LK mechanism consists of large-amplitude oscillations in eccentricity and inclination of a binary subject to the quadrupole potential from an outer perturber. Recent work has shown that when the octupole terms are included in the potential, the inner binary can reach more extreme eccentricities as well as undergo orientation flips. It is known that pericenter precessions due to short-range effects, such as General Relativity and tidal and rotational distortions, can limit the growth of eccentricity and even suppress standard (quadrupolar) LK oscillations, but their effect on the octupole-level LK mechanism has not been fully explored. In this paper, we systematically study how these short-range forces affect the extreme orbital behaviour found in octupole LK cycles. In general, the influence of the octupole potential is confined to a range of initial mutual inclinations Itot centered around 90deg (when the inner binary mass ratio is <<1), with this range expanding with increasing octupole strength. We find that, while the short-range forces do not change the width and location of this "window of influence", they impose a strict upper limit on the maximum achievable eccentricity. This limiting eccentricity can be calculated analytically, and its value holds even for strong octupole potential and for the general case of three comparable masses. Short-range forces also affect orbital flips, progressively reducing the range of Itot within which flips are possible as the intensity of these forces increases.

The formation and early evolution of circumstellar discs often occurs within dense, newborn stellar clusters. For the first time, we apply the moving-mesh code AREPO, to circumstellar discs in 3-D, focusing on disc-disc interactions that result from stellar fly-bys. Although a small fraction of stars are expected to undergo close approaches, the outcomes of the most violent encounters might leave an imprint on the discs and host stars that will influence both their orbits and their ability to form planets. We first construct well-behaved 3-D models of self-gravitating discs, and then create a suite of numerical experiments of parabolic encounters, exploring the effects of pericenter separation r_p, disc orientation and disc-star mass ratio (M_d/M_*) on the orbital evolution of the host stars. Close encounters (2r_p<~ disc radius) can truncate discs on very short time scales. If discs are massive, close encounters facilitate enough orbital angular momentum extraction to induce stellar capture. We find that for realistic primordial disc masses M_d<~0.1M_*, non-colliding encounters induce minor orbital changes, which is consistent with analytic calculations of encounters in the linear regime. The same disc masses produce entirely different results for grazing/colliding encounters. In the latter case, rapidly cooling discs lose orbital energy by radiating away the energy excess of the shock-heated gas, thus causing capture of the host stars into a bound orbit. In rare cases, a tight binary with a circumbinary disc forms as a result of this encounter.

In this study, we present high-resolution millimeter observations of the dust and gas disk of the T Tauri star AS 205 N and its companion, AS 205 S, obtained with the Atacama Large Millimeter Array. The gas disk around AS 205 N, for which infrared emission spectroscopy demonstrates significant deviations from Keplerian motion that has been interpreted as evidence for a disk wind (Pontoppidan et al. 2011; Bast et al. 2011), also displays significant deviations from Keplerian disk emission in the observations presented here. Detections near both AS 205 N and S are obtained in 1.3 mm continuum, 12CO 2-1, 13CO 2-1 and C18O 2-1. The 12CO emission is extended up to 2 arcsec from AS 205N, and both 12CO and 13CO display deviations from Keplerian rotation at all angular scales. Two possible explanations for these observations hold up best to close scrutiny - tidal interaction with AS 205 S or disk winds (or a combination of the two), and we discuss these possibilities in some detail.

Numerous formulations of finite volume schemes for the Euler and Navier-Stokes equations exist, but in the majority of cases they have been developed for structured and stationary meshes. In many applications, more flexible mesh geometries that can dynamically adjust to the problem at hand and move with the flow in a (quasi) Lagrangian fashion would, however, be highly desirable, as this can allow a significant reduction of advection errors and an accurate realization of curved and moving boundary conditions. Here we describe a novel formulation of viscous continuum hydrodynamics that solves the equations of motion on a Voronoi mesh created by a set of mesh-generating points. The points can move in an arbitrary manner, but the most natural motion is that given by the fluid velocity itself, such that the mesh dynamically adjusts to the flow. Owing to the mathematical properties of the Voronoi tessellation, pathological mesh-twisting effects are avoided. Our implementation considers the full Navier-Stokes equations and has been realized in the AREPO code both in 2D and 3D. We propose a new approach to compute accurate viscous fluxes for a dynamic Voronoi mesh, and use this to formulate a finite volume solver of the Navier-Stokes equations. Through a number of test problems, including circular Couette flow and flow past a cylindrical obstacle, we show that our new scheme combines good accuracy with geometric flexibility, and hence promises to be competitive with other highly refined Eulerian methods. This will in particular allow astrophysical applications of the AREPO code where physical viscosity is important, such as in the hot plasma in galaxy clusters, or for viscous accretion disk models.

We report the first detections of circularly polarized emission at submillimeter wavelengths from the compact radio source and supermassive black hole candidate Sgr A* at a level of 1.2\pm0.3% at 1.3 mm wavelength (230 GHz) and 1.6\pm0.3% at 860 microns (345 GHz) with the same handedness as observed at lower frequencies (1.4-15 GHz). The observations, taken with the Submillimeter Array in multiple epochs, also show simultaneous linear polarization (LP) at both wavelengths of about 6%. These properties differ sharply from those at wavelengths longer than 1 cm (frequencies below 30 GHz), where weak circular polarization (CP) (~ 0.5%) dominates over LP, which is not detected at similar fractional limits. We describe an extensive set of tests to ensure the accuracy of our measurements. We find no CP in any other source, including the bright quasar 1924-292, which traces the same path on the sky as Sgr A* and therefore should be subject to identical systematic errors originating in the instrument frame. Since a relativistic synchrotron plasma is expected to produce little CP, the observed CP is probably generated close to the event horizon by the Faraday conversion process. We use a simple model to show that the phase shift associated with Faraday conversion can be nearly independent of frequency, a sufficient condition to make the handedness of CP independent of frequency. Because the size of the tau=1-surface changes by more than an order of magnitude between 1.4 and 345 GHz, the magnetic field must be coherent over such scales to consistently produce left CP.