Eccentric gaseous discs are unstable to a parametric instability involving the resonant interaction between inertial-gravity waves and the eccentric mode in the disc. We present 3D global hydrodynamical simulations of inviscid circumbinary discs that form an inner cavity and become eccentric through interaction with the central binary. The parametric instability grows and generates turbulence that transports angular momentum with stress parameter $\alpha \sim 5 \times 10^{-3}$ at distances $\lesssim 7 \;a_{bin} $, where $a_{bin}$ is the binary semi-major axis. Vertical turbulent diffusion occurs at a rate corresponding to $\alpha_{diff}\sim 1-2\times 10^{-3}$. We examine the impact of turbulent diffusion on the vertical settling of pebbles, and on the rate of pebble accretion by embedded planets. In steady state, dust particles with Stokes numbers ${\it St} \lesssim 0.1$ form a layer of finite thickness $H_d \gtrsim 0.1 H$, where $H$ is the gas scale height. Pebble accretion efficiency is then reduced by a factor $r_{acc}/H_d$, where $r_{acc}$ is the accretion radius, compared to the rate in a laminar disc. For accreting core masses with $m_p \lesssim 0.1\; M_\oplus$, pebble accretion for particles with ${\it St} \gtrsim 0.5$ is also reduced because of velocity kicks induced by the turbulence. These effects combine to make the time needed by a Ceres-mass object to grow to the pebble isolation mass, when significant gas accretion can occur, longer than typical disc lifetimes. Hence, the origins of circumbinary planets orbiting close to their central binary systems, as discovered by the Kepler mission, are difficult to explain using an in situ model that invokes a combination of the streaming instability and pebble accretion.

We present simulations of low-mass planet--disc interactions in inviscid three-dimensional discs. We show that a wind-driven laminar accretion flow through the surface layers of the disc does not significantly modify the migration torque experienced by embedded planets. More importantly, we find that 3D effects lead to a dramatic change in the behaviour of the dynamical corotation torque compared to earlier 2D theory and simulations. Although it was previously shown that the dynamical corotation torque could act to slow and essentially stall the inward migration of a low-mass planet, our results in 3D show that the dynamical corotation torque has the complete opposite effect and speeds up inward migration. Our numerical experiments implicate buoyancy resonances as the cause. These have two effects: (i) they exert a direct torque on the planet, whose magnitude relative to the Lindblad torque is measured in our simulations to be small; (ii) they torque the gas librating on horseshoe orbits in the corotation region and drive evolution of its vortensity, leading to the negative dynamical corotation torque. This indicates that at low turbulent viscosity, the detailed vertical thermal structure of the protoplanetary disc plays an important role in determining the migration behaviour of embedded planets. If this result holds up under a more refined treatment of disc thermal evolution, then it has important implications for understanding the formation and early evolution of planetary systems.

Convergent migration involving multiple planets embedded in a viscous protoplanetary disc is expected to produce a chain of planets in mean motion resonances, but the multiplanet systems observed by the Kepler spacecraft are generally not in resonance. We demonstrate that under equivalent conditions, where in a viscous disc convergent migration will form a long-term stable system of planets in a chain of mean motion resonances, migration in an inviscid disc often produces a system which is highly dynamically unstable. In particular, if planets are massive enough to significantly perturb the disc surface density and drive vortex formation, the smooth capture of planets into mean motion resonances is disrupted. As planets pile up in close orbits, not protected by resonances, close encounters increase the probability of planet-planet collisions, even while the gas disc is still present. While inviscid discs often produce unstable non-resonant systems, stable, closely packed, non-resonant systems can also be formed. Thus, when examining the expectation for planet migration to produce planetary systems in mean motion resonances, the effective turbulent viscosity of the protoplanetary disc is a key parameter.

Modern astronomy has finally been able to observe protoplanetary disks in reasonable resolution and detail, unveiling the processes happening during planet formation. These observed processes are understood under the framework of disk-planet interaction, a process studied analytically and modeled numerically for over 40 years. Long a theoreticians' game, the wealth of observational data has been allowing for increasingly stringent tests of the theoretical models. Modeling efforts are crucial to support the interpretation of direct imaging analyses, not just for potential detections but also to put meaningful upper limits on mass accretion rates and other physical quantities in current and future large-scale surveys. This white paper addresses the questions of what efforts on the computational side are required in the next decade to advance our theoretical understanding, explain the observational data, and guide new observations. We identified the nature of accretion, ab initio planet formation, early evolution, and circumplanetary disks as major fields of interest in computational planet formation. We recommend that modelers relax the approximations of alpha-viscosity and isothermal equations of state, on the grounds that these models use flawed assumptions, even if they give good visual qualitative agreement with observations. We similarly recommend that population synthesis move away from 1D hydrodynamics. The computational resources to reach these goals should be developed during the next decade, through improvements in algorithms and the hardware for hybrid CPU/GPU clusters. Coupled with high angular resolution and great line sensitivity in ground based interferometers, ELTs and JWST, these advances in computational efforts should allow for large strides in the field in the next decade.

We present the highest resolution study to date of super-Earths migrating in inviscid and low-viscosity discs, motivated by the connection to laminar, wind-driven models of protoplanetary discs. Our models unveil the critical role of vortices in determining the migration behaviour for partial gap-opening planets. Vortices form in pressure maxima at gap edges, and prevent the disc-feedback stopping of migration for intermediate planets in low-viscosity and inviscid discs, contrary to the concept of the `inertial limit' or `disc feedback' halting predicted from analytical models. Vortices may also form in the corotation region, and can dramatically modify migration behaviour through direct gravitational interaction with the planet. These features become apparent at high resolution, and for all but the highest viscosities there exist significant difficulties in obtaining numerically converged results. The migration of partial gap-opening planets, however, clearly becomes chaotic for sufficiently low viscosities. At moderate viscosity, a smooth disc-feedback regime is found in which migration can slow substantially, and the migration time-scale observed corresponds to migration being driven by diffusive relaxation of the gap edges. At high viscosity classical Type I migration is recovered. For Jupiter-analogue planets in inviscid discs, a wide, deep gap is formed. Transient Type II migration occurs over radial length-scales corresponding to the gap width, beyond which migration can stall. Finally, we examine the particle trapping driven by structures left in inviscid discs by a migrating planet, and find that particle traps in the form of multiple rings and vortices can persist long after the planet has passed. In this case, the observation of particle traps by submillimetre interferometers such as ALMA cannot be used to infer the current presence of an adjacent planet.

Recent developments in non-ideal magnetohydrodynamic simulations of protoplanetary disks suggest that instead of being traditional turbulent (viscous) accretion disks, they have a largely laminar flow with accretion driven by large-scale wind torques. These disks are possibly threaded by Hall-effect generated large-scale horizontal magnetic fields. We have examined the dynamics of the corotation region of a low mass planet embedded in such a disk and the evolution of the associated migration torque. These disks lack strong turbulence and associated turbulent diffusion, and the presence of a magnetic field and radial gas flow presents a situation outside the applicability of previous corotation torque theory. We summarize the analytical analysis of the corotation torque, give details on the numerical methods used, and in particular the relative merits of different numerical schemes for the inviscid problem.

We examine the migration of low mass planets in laminar protoplanetary discs, threaded by large scale magnetic fields in the dead zone that drive radial gas flows. As shown in Paper I, a dynamical corotation torque arises due to the flow-induced asymmetric distortion of the corotation region and the evolving vortensity contrast between the librating horseshoe material and background disc flow. Using simulations of laminar torqued discs containing migrating planets, we demonstrate the existence of the four distinct migration regimes predicted in Paper I. In two regimes, the migration is approximately locked to the inward or outward radial gas flow, and in the other regimes the planet undergoes outward runaway migration that eventually settles to fast steady migration. In addition, we demonstrate torque and migration reversals induced by midplane magnetic stresses, with a bifurcation dependent on the disc surface density. We develop a model for fast migration, and show why the outward runaway saturates to a steady speed, and examine phenomenologically its termination due to changing local disc conditions. We also develop an analytical model for the corotation torque at late times that includes viscosity, for application to discs that sustain modest turbulence. Finally, we use the simulation results to develop torque prescriptions for inclusion in population synthesis models of planet formation.

Orbital advection is a significant bottleneck in disk simulations, and a particularly tricky one when used in connection with magnetohydrodynamics. We have developed an orbital advection algorithm suitable for the induction equation with magnetic potential. The electromotive force is split into advection and shear term, and we find that we do not need an advective gauge as solving the orbital advection implicitly precludes the shear term from canceling the advection term. We prove and demonstrate the third order in time accuracy of the scheme. The algorithm is also suited for non-magnetic problems. Benchmarked results of (hydrodynamical) planet-disk interaction and the of the magnetorotational instability are reproduced. We include detailed descriptions of the construction and selection of stabilizing dissipations (or high frequency filters) needed to generate practical results. The scheme is self-consistent, accurate, and elegant in its simplicity, making it particularly efficiently in straightforward finite-difference methods. As a result of the work, the algorithm is incorporated in the public version of the Pencil Code, where it can be used by the community.

Motivated by models suggesting that the inner planet forming regions of protoplanetary discs are predominantly lacking in viscosity-inducing turbulence, and are possibly threaded by Hall-effect generated large-scale horizontal magnetic fields, we examine the dynamics of the corotation region of a low-mass planet in such an environment. The corotation torque in an inviscid, isothermal, dead zone ought to saturate, with the libration region becoming both symmetrical and of a uniform vortensity, leading to fast inward migration driven by the Lindblad torques alone. However, in such a low viscosity situation, the material on librating streamlines essentially preserves its vortensity. If there is relative radial motion between the disc gas and the planet, the librating streamlines will no longer be symmetrical. Hence, if the gas is torqued by a large scale magnetic field so that it undergoes a net inflow or outflow past the planet, driving evolution of the vortensity and inducing asymmetry of the corotation region, the corotation torque can grow, leading to a positive torque. In this paper we treat this effect by applying a symmetry argument to the previously studied case of a migrating planet in an inviscid disc. Our results show that the corotation torque due to a laminar Hall-induced magnetic field in a dead zone behaves quite differently from that studied previously for a viscous disc. Furthermore, the magnetic field induced corotation torque and the dynamical corotation torque in a low viscosity disc can be regarded as one unified effect.

In protoplanetary disks, the differential gravity-driven settling of dust grains with respect to gas and with respect to grains of varying sizes determines the observability of grains, and sets the conditions for grain growth and eventually planet formation. In this work we explore the effect of photophoresis on the settling of large dust grains in the inner regions of actively accreting protoplanetary disks. Photophoretic forces on dust grains result from the collision of gas molecules with differentially heated grains. We undertake one dimensional dust settling calculations to determine the equilibrium vertical distribution of dust grains in each column of the disk. In the process we introduce a new treatment of the photophoresis force which is consistent at all optical depths with the representation of the radiative intensity field in a two-stream radiative transfer approximation. The levitation of large dust grains creates a photophoretic dust trap several scale heights above the mid-plane in the inner regions of the disk where the dissipation of accretion energy is significant. We find that differential settling of dust grains is radically altered in these regions of the disk, with large dust grains trapped in a layer below the stellar irradiation surface in where the dust to gas mass ratio can be enhanced by a factor of a hundred for the relevant particles. The photophoretic trapping effect has a strong dependence on particle size and porosity.

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.

Disks with a barotropic equilibrium structure, for which the pressure is only a function of the density, rotate on cylinders in the presence of a gravitational potential, so that the angular frequency of such a disk is independent of height. Such disks with barotropic equilibria can be approximately modeled using the shearing box framework, representing a small disk volume with height-independent angular frequency. If the disk is in baroclinic equilibrium, the angular frequency does generally depend on height, and it is thus necessary to go beyond the standard shearing box approach. In this paper, we show that given a global disk model, it is possible to develop approximate models that are local in horizontal planes without an expansion in height with shearing-periodic boundary conditions. We refer to the resulting framework as the vertically global shearing box (VGSB). These models can be non-axisymmetric for globally barotropic equilibria but should be axisymmetric for globally baroclinic equilibria. We provide explicit equations for this VGSB which can be implemented in standard magnetohydrodynamic codes by generalizing the shearing-periodic boundary conditions to allow for a height-dependent angular frequency and shear rate. We also discuss the limitations that result from the radial approximations that are needed in order to impose height-dependent shearing periodic boundary conditions. We illustrate the potential of this framework by studying a vertical shear instability and examining the modes associated with the magnetorotational instability.

Protoplanetary disks are believed to accrete onto their central T Tauri star because of magnetic stresses. Recently published shearing box simulations indicate that Ohmic resistivity, ambipolar diffusion and the Hall effect all play important roles in disk evolution. In the presence of a vertical magnetic field, the disk remains laminar between 1-5au, and a magnetocentrifugal disk wind forms that provides an important mechanism for removing angular momentum. Questions remain, however, about the establishment of a true physical wind solution in the shearing box simulations because of the symmetries inherent in the local approximation. We present global MHD simulations of protoplanetary disks that include Ohmic resistivity and ambipolar diffusion, where the time-dependent gas-phase electron and ion fractions are computed under FUV and X-ray ionization with a simplified recombination chemistry. Our results show that the disk remains laminar, and that a physical wind solution arises naturally in global disk models. The wind is sufficiently efficient to explain the observed accretion rates. Furthermore, the ionization fraction at intermediate disk heights is large enough for magneto-rotational channel modes to grow and subsequently develop into belts of horizontal field. Depending on the ionization fraction, these can remain quasi-global, or break-up into discrete islands of coherent field polarity. The disk models we present here show a dramatic departure from our earlier models including Ohmic resistivity only. It will be important to examine how the Hall effect modifies the evolution, and to explore the influence this has on the observational appearance of such systems, and on planet formation and migration.

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.

Rossby wave instability (RWI) at dead zone boundaries may play an important role in planet formation. Viscous hydrodynamics results suggest RWI is excited only when the viscosity changes over a radial distance less than two density scale heights. However in the disks around Solar-mass T Tauri stars, it is not viscosity but magnetic forces that provide the accretion stress beyond about 10 AU, where surface densities are low enough so stellar X-rays and interstellar cosmic rays can penetrate. Here we aim to explore the conditions for RWI in the smooth transition with increasing distance, from resistive and magnetically-dead to conducting and magnetically-active. We perform 3D unstratified MHD simulations with the Pencil Code, using static resistivity profiles. As a result, we find that in MHD, contrary to viscous models, the RWI is triggered even with a gradual change in resistivity extending from 10 to 40 AU (i.e., spanning 15 scale heights for aspect ratio 0.1). This is because magneto-rotational turbulence sets in abruptly when the resistivity reaches a threshold level. At higher resistivities the longest unstable wavelength is quenched, resulting in a sharp decline of the Maxwell stress towards the star. The sharp gradient in the magnetic forces leads to a localized density bump, that is in turn Rossby wave unstable. We conclude that even weak gradients in the resistivity can lead to sharp transitions in the Maxwell stress. The upshot is that the RWI is more easily activated in the outer disk than previously thought. Rossby vortices at the outer dead zone boundary thus could underlie the dust asymmetries seen in the outer reaches of transition disks.

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.

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.

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.

Recently, there has been a significant level of discussion of the correct treatment of Kelvin-Helmholtz instability in the astrophysical community. This discussion relies largely on how the KHI test is posed and analyzed. We pose a stringent test of the initial growth of the instability. The goal is to provide a rigorous methodology for verifying a code on two dimensional Kelvin-Helmholtz instability. We ran the problem in the Pencil Code, Athena, Enzo, NDSPHMHD, and Phurbas. A strict comparison, judgment, or ranking, between codes is beyond the scope of this work, though this work provides the mathematical framework needed for such a study. Nonetheless, how the test is posed circumvents the issues raised by tests starting from a sharp contact discontinuity yet it still shows the poor performance of Smoothed Particle Hydrodynamics. We then comment on the connection between this behavior to the underlying lack of zeroth-order consistency in Smoothed Particle Hydrodynamics interpolation. We comment on the tendency of some methods, particularly those with very low numerical diffusion, to produce secondary Kelvin-Helmholtz billows on similar tests. Though the lack of a fixed, physical diffusive scale in the Euler equations lies at the root of the issue, we suggest that in some methods an extra diffusion operator should be used to damp the growth of instabilities arising from grid noise. This statement applies particularly to moving-mesh tessellation codes, but also to fixed-grid Godunov schemes.

We present an algorithm for simulating the equations of ideal magnetohydrodynamics and other systems of differential equations on an unstructured set of points represented by sample particles. The particles move with the fluid, so the time step is not limited by the Eulerian Courant-Friedrichs-Lewy condition. Full spatial adaptivity is required to ensure the particles fill the computational volume, and gives the algorithm substantial flexibility and power. A target resolution is specified for each point in space, with particles being added and deleted as needed to meet this target. We have parallelized the code by adapting the framework provided by GADGET-2. A set of standard test problems, including 1e-6 amplitude linear MHD waves, magnetized shock tubes, and Kelvin-Helmholtz instabilities is presented. Finally we demonstrate good agreement with analytic predictions of linear growth rates for magnetorotational instability in a cylindrical geometry. This paper documents the Phurbas algorithm as implemented in Phurbas version 1.1.

We present an algorithm for simulating the equations of ideal magnetohydrodynamics and other systems of differential equations on an unstructured set of points represented by sample particles. Local, third-order, least-squares, polynomial interpolations (Moving Least Squares interpolations) are calculated from the field values of neighboring particles to obtain field values and spatial derivatives at the particle position. Field values and particle positions are advanced in time with a second order predictor-corrector scheme. The particles move with the fluid, so the time step is not limited by the Eulerian Courant-Friedrichs-Lewy condition. Full spatial adaptivity is implemented to ensure the particles fill the computational volume, which gives the algorithm substantial flexibility and power. A target resolution is specified for each point in space, with particles being added and deleted as needed to meet this target. Particle addition and deletion is based on a local void and clump detection algorithm. Dynamic artificial viscosity fields provide stability to the integration. The resulting algorithm provides a robust solution for modeling flows that require Lagrangian or adaptive discretizations to resolve. This paper derives and documents the Phurbas algorithm as implemented in Phurbas version 1.1. A following paper presents the implementation and test problem results.

In astrophysical magnetohydrodynamics (MHD) and electrodynamics simulations, numerically enforcing the divB=0 constraint on the magnetic field has been difficult. We observe that for point-based discretization, as used in finite-difference type and pseudo-spectral methods, the divB=0 constraint can be satisfied entirely by a choice of interpolation used to define the derivatives of B. As an example we demonstrate a new class of finite-difference type derivative operators on a regular grid which has the divB=0 property. This principle clarifies the nature of divB != 0 errors. The principles and techniques demonstrated in this paper are particularly useful for the magnetic field, but can be applied to any vector field. This paper serves as a brief introduction to the method and demonstrates an implementation showing convergence.

We demonstrate a significant difference in the angular momentum transport properties of galactic disks between regions in which the interstellar medium is single phase or two phase. Our study is motivated by observations of HI in extended galactic disks which indicate velocity dispersions of nonthermal origin, suggesting that turbulence in the gas may be contributing significantly to the observed dispersion. To address this, we have implemented a shearing-box framework within the FLASH code. The new code was used to perform local simulations of galactic disks that incorporate differential rotation, self-gravity, vertical stratification, hydrodynamics and cooling. These simulations explore plausible mechanisms for driving turbulent motions via the thermal and self-gravitational instabilities coupling to differential rotation. Where a two-phase medium develops, gravitational angular momentum transporting stresses are much greater, creating a possible mechanism for transferring energy from galactic rotation to turbulence. In simulations where the disk conditions do not trigger the formation of a two-phase medium, it is found that perturbations to the flow damp without leading to a sustained mechanism for driving turbulence. The differing angular momentum transport properties of the single- and two-phase regimes of the disk suggest that a significant, dynamically motivated division can be drawn between the two, even when this division occurs far outside the star formation cutoff in a galactic disk.

We have conducted a deep multi-color imaging survey of 0.2 degrees^2 centered on the Hubble Deep Field North (HDF-N). We shall refer to this region as the Hawaii-HDF-N. Deep data were collected in U, B, V, R, I, and z' bands over the central 0.2 degrees^2 and in HK' over a smaller region covering the Chandra Deep Field North (CDF-N). The data were reduced to have accurate relative photometry and astrometry across the entire field to facilitate photometric redshifts and spectroscopic followup. We have compiled a catalog of 48,858 objects in the central 0.2 degrees^2 detected at 5 sigma significance in a 3" aperture in either R or z' band. Number counts and color-magnitude diagrams are presented and shown to be consistent with previous observations. Using color selection we have measured the density of objects at 3<z<7. Our multi-color data indicates that samples selected at z>5.5 using the Lyman break technique suffer from more contamination by low redshift objects than suggested by previous studies.