Cold Dark Stars made of self-gravitating fermions in the degenerate limit are constructed in General Relativity and in R-square gravity, $f(R)=R+\alpha R^2$. The properties of the resulting Cold Dark Stars in both theories of gravity are studied. It is found that the same gravitational potential is generated for different election of the parameters of the model, such as the mass of the fermion, the self-interacting strength or the value of $\alpha$, thus, a possible confusion in the determination of the dark matter properties and the favored theory of gravity might arise.

The equation of state for a degenerate gas of fermions at zero temperature in the non relativistic case is a polytrope, i.e. $p=\gamma \rho^{5/3}/m_F^{8/3}$. If dark matter is modelled by such non interacting fermion, this dependence in the mass of the fermion $m_F$ explains why if dark matter is very heavy the effective pressure of dark matter is negligible. Nevertheless, if the mass of the dark matter is very small, the effective pressure can be very large, and thus, a system of self-gravitating fermions can be formed. In this work we model the dark matter halo of the Milky-Way by solving the Tolman-Oppenheimer-Volkoff equations, with the equation of state for a partially degenerate ultralight non interacting fermion. It is found that in order to fit its rotational velocity curve of the Milky Way, the mass of the fermion should be in the range $29 ~\mbox{eV} < m_F < 33~$eV. Moreover, the central density is constrained to be in the range of $46 < \rho_0<61$ GeV/cm$^3$. The fermionic dark matter halo has a very different profile as compared with the standard Navarro-Frenk-White profile, thus, the possible indirect signals for annihilating dark matter may change by orders of magnitude. We found bounds for the annihilation cross section in this case by using the Saggitarius A* spectral energy distribution. Those limits are very strong confirming the idea that the lighter the dark matter particle is, the darkest it becomes.

The Zombie Vortex Instability (ZVI) occurs in the dead zones of protoplanetary disks (PPDs) where perturbations excite baroclinic critical layers, generating "zombie" vortices and turbulence. In this work, we investigate ZVI with nonuniform vertical stratification; while ZVI is triggered in the stratified regions away from the midplane, the subsequent turbulence propagates into and fills the midplane. ZVI turbulence alters the background Keplerian shear flow, creating a steady-state zonal flow. Intermittency is observed, where the flow cycles through near-laminar phases of zonal flow punctuated by chaotic bursts of new vortices. ZVI persists in the presence of radiative damping, as long as the thermal relaxation timescale is more than a few orbital periods. We refute the premature claim by Lesur & Latter (2016) that radiative damping inhibits ZVI for disk radii r>0.3 au. Their conclusions were based on unrealistically short cooling times using opacities with virtually no grain growth. We explore different grain growth and vertical settling scenarios, and find that the gas and dust in off-midplane regions are not necessarily in local thermodynamic equilibrium (LTE) with each other. In such cases, thermal relaxation timescales can be orders of magnitude longer than the optically thin cooling times assuming LTE because of the finite time for energy to be exchanged between gas and dust grains via collisions. With minimal amounts of grain growth and dust settling, the off-midplane regions of disks are susceptible to ZVI and much of the planet-forming regions can be filled with zombie vortices and turbulence.

We present new, fully nonlinear numerical solutions to the static, spherically symmetric Einstein-Klein-Gordon system for a collection of an arbitrary odd number $N$ of complex scalar fields with an internal $U(N)$ symmetry and no self-interactions. These solutions, which we dub $\ell$-boson stars, are parametrized by an angular momentum number $\ell=(N-1)/2$, an excitation number $n$, and a continuous parameter representing the amplitude of the fields. They are regular at every point and possess a finite total mass. For $\ell = 0$ the standard spherically symmetric boson stars are recovered. We determine their generalizations for $\ell > 0$, and show that they give rise to a large class of new static configurations which might have a much larger compactness ratio than $\ell=0$ stars.

The future detection and measurement of the diffuse neutrino supernova background will shed light on the rate of supernovae events in the Universe, the star formation rate and the neutrino spectrum from each supernova. Little has been said about what those measurements will tell us about the expansion history of the universe. The purpose of this article is to show that the detection of the diffuse supernova neutrino background will be a complementary tool for the study and possible discrimination of cosmological models. In particular, we study three different cosmological models: the $\Lambda$ Cold Dark Matter model, the Logotropic universe and a bulk viscous matter-dominated universe. By fitting the free parameters of each model with the supernova Ia probe, we found that the predicted number of events computed with the best fit parameters for the $\Lambda$-Cold dark matter model and with the Logotropic model are the same, while a bulk viscous matter-dominated cosmological model predicts $\sim 3$ times more events. We show that the current limit set by Super-Kamiokande on the diffuse supernova neutrino background flux gives complementary constraints on the free parameters of a bulk viscous matter-dominated universe. Furthermore, this limit implies, within a $\Lambda$ Cold Dark Matter model, that the universe should be expanding with $H_0 > 21.5 ~\rm{Km/sec/Mpc}$ independently of the content of dark matter $\Omega_m$.

It has long been known that no static, spherically symmetric, asymptotically flat Klein-Gordon scalar field configuration surrounding a nonrotating black hole can exist in general relativity. In a series of previous papers we proved that, at the effective level, this no-hair theorem can be circumvented by relaxing the staticity assumption: for appropriate model parameters there are quasi-bound scalar field configurations living on a fixed Schwarzschild background which, although not being strictly static, have a larger lifetime than the age of the universe. This situation arises when the mass of the scalar field distribution is much smaller than the black hole mass, and following the analogies with the hair in the literature we dubbed these long-lived field configurations wigs. Here we extend our previous work to include the gravitational backreaction produced by the scalar wigs. We derive new approximate solutions of the spherically symmetric Einstein-Klein-Gordon system which represent self-gravitating scalar wigs surrounding black holes. These configurations interpolate between boson star configurations and Schwarzschild black holes dressed with the long-lived scalar test field distributions discussed in previous papers. Nonlinear numerical evolutions of initial data sets extracted from our approximate solutions support the validity of our approach. Arbitrarily large lifetimes are still possible, although for the parameter space that we analyze in this paper they seem to decay faster than the quasi-bound states. Finally, we speculate about the possibility that these configurations could describe the innermost regions of dark matter halos.

In Zombie Vortex Instability (ZVI), perturbations excite critical layers in stratified, rotating shear flow (as in protoplanetary disks), causing them to generate vortex layers, which roll-up into anticyclonic zombie vortices and cyclonic vortex sheets. The process is self-sustaining as zombie vortices perturb new critical layers, spawning a next generation of zombie vortices. Here, we focus on two issues: the minimum threshold of perturbations that trigger self-sustaining vortex generation, and the properties of the late-time zombie turbulence on large and small scales. The critical parameter that determines whether ZVI is triggered is the magnitude of the vorticity on the small scales (and not velocity), the minimum Rossby number needed for instability is $Ro_{crit}\sim0.2$ for $\beta\equiv N/\Omega = 2$, where $N$ is the Brunt-V\"ais\"al\"a frequency. While the threshold is set by vorticity, it is useful to infer a criterion on the Mach number, for Kolmogorov noise, the critical Mach number scales with Reynolds number: $Ma_{crit}\sim Ro_{crit}Re^{-1/2}$. In protoplanetary disks, this is $Ma_{crit}\sim10^{-6}$. On large scales, zombie turbulence is characterized by anticyclones and cyclonic sheets with typical Rossby number $\sim$0.3. The spacing of the cyclonic sheets and anticyclones appears to have a "memory" of the spacing of the critical layers. On the small scales, zombie turbulence has no memory of the initial conditions and has a Kolmogorov-like energy spectrum. While our earlier work was in the limit of uniform stratification, we have demonstrated that ZVI works for non-uniform Brunt-V\"ais\"al\"a frequency profiles that may be found in protoplanetary disks.

There is considerable interest in hydrodynamic instabilities in dead zones of protoplanetary disks as a mechanism for driving angular momentum transport and as a source of particle-trapping vortices to mix chondrules and incubate planetesimal formation. We present simulations with a pseudo-spectral anelastic code and with the compressible code Athena, showing that stably stratified flows in a shearing, rotating box are violently unstable and produce space-filling, sustained turbulence dominated by large vortices with Rossby numbers of order 0.2-0.3. This Zombie Vortex Instability (ZVI) is observed in both codes and is triggered by Kolmogorov turbulence with Mach numbers less than 0.01. It is a common view that if a given constant density flow is stable, then stable vertical stratification should make the flow even more stable. Yet, we show that sufficient vertical stratification can be unstable to ZVI. ZVI is robust and requires no special tuning of boundary conditions, or initial radial entropy or vortensity gradients (though we have studied ZVI only in the limit of infinite cooling time). The resolution of this paradox is that stable stratification allows for a new avenue to instability: baroclinic critical layers. ZVI has not been seen in previous studies of flows in rotating, shearing boxes because those calculations frequently lacked vertical density stratification and/or sufficient numerical resolution. Although we do not expect appreciable angular momentum transport from ZVI in the small domains in this study, we hypothesize that ZVI in larger domains with compressible equations may lead to angular transport via spiral density waves.

In this work we model galactic halos describing the dark matter as a non zero pressure fluid and derive, not impose, a dark matter equation of state by using observational data of the rotation curves of galaxies. In order to reach hydrostatic equilibrium, as expected for the halo, it is mandatory that dark fluid's pressure should not be zero. The equation of state is obtained by solving the matter-geometry system of equations assuming different dark matter density or rotational velocity profiles. The resulting equations of state are, in general, different to a barotropic equation of state. The free parameters of the equation of state are fixed by fitting the observed rotational velocities of a set of galaxies.

If dark matter is mainly composed of axions, the density distribution can be nonuniformly distributed, being clumpy instead. By solving the Einstein-Klein-Gordon system of a scalar field with the potential energy density of an axionlike particle, we obtain the maximum mass of the self-gravitating system made of axions, called axion stars. The collision of axion stars with neutron stars may release the energy of axions due to the conversion of axions into photons in the presence of the neutron star's magnetic field. We estimate the energy release and show that it should be much less than previous estimates.Future data from femtolensing should strongly constrain this scenario.

Motivated by the stringent flux limits for UHE neutrinos coming from gamma ray burst or active galactic nuclei, we explore the possibility that the active neutrinos generated in such astrophysical objects could oscillate to sterile right handed states due to a neutrino magnetic moment mu_nu. We find that a value as small as mu_nu ~1E-15 mu_B could produce such a transition thanks to the intense magnetic fields that are expected in these objects.

We study the evolution of a massive scalar field surrounding a Schwarzschild black hole and find configurations that can survive for arbitrarily long times, provided the black hole or the scalar field mass is small enough. In particular, both ultra-light scalar field dark matter around supermassive black holes and axion-like scalar fields around primordial black holes can survive for cosmological times. Moreover, these results are quite generic, in the sense that fairly arbitrary initial data evolves, at late times, as a combination of those long-lived configurations.

We study the possible suppression of the extragalactic neutrino flux due to a nonstandard interaction during its propagation. In particular, we study neutrino interaction with an ultra-light scalar field dark matter. It is shown that the extragalactic neutrino flux may be suppressed by such an interaction, leading to a new mechanism to reduce the ultra-high energy neutrino flux. We study both the cases of non-self-conjugate as well as self-conjugate dark matter. In the first case, the suppression is independent of the neutrino and dark matter masses. We conclude that care must be taken when explaining limits on the neutrino flux through source acceleration mechanisms only, since there could be other mechanisms for the reduction of the neutrino flux.

Constraints to the mass of a scalar field and the strength of its self-interacting coupling constant are obtained. This was done using observations of stellar dynamics at the center of our galaxy and by assuming that the dark compact object responsible of such dynamics is a boson star and not a supermassive black hole. We show that if such scalar field represents a spin-zero particle with cross section high enough to be considered collisional dark matter, there is a region of parameters compatible with both conditions: that the scalar field play the role of collisional dark matter and that it can form objects with the mass and compactness compatible with stellar kinematics.

We show that the inclusion of an axion-like effective potential in the construction of a self-gravitating system made of scalar fields leads to a decrease on its compactness when the value of the self-interaction coupling constant is increased. By including the current values for the axion mass m and decay constant f_a, we have computed the mass and the radius for self-gravitating systems made of axion particles. It is found that such objects will have asteroid-size masses and radius of few meters, then, the self-gravitating system made of axions could play the role of scalar mini-machos that are mimicking a cold dark matter model for the galactic halo.

The existence and detection of scalar fields could provide solutions to long-standing puzzles about the nature of dark matter, the dark compact objects at the centre of most galaxies, and other phenomena. Yet, self-interacting scalar fields are very poorly constrained by astronomical observations, leading to great uncertainties in estimates of the mass $m_\phi$ and the self-interacting coupling constant $\lambda$ of these fields. To counter this, we have systematically employed available astronomical observations to develop new constraints, considerably restricting this parameter space. In particular, by exploiting precise observations of stellar dynamics at the centre of our Galaxy and assuming that these dynamics can be explained by a single boson star, we determine an upper limit for the boson star compactness and impose significant limits on the values of the properties of possible scalar fields. Requiring the scalar field particle to follow a collisional dark matter model further narrows these constraints. Most importantly, we find that if a scalar dark matter particle does exist, then it cannot account for both the dark-matter halos and the existence of dark compact objects in galactic nuclei

Gravitational instability (GI) of a dust-rich layer at the midplane of a gaseous circumstellar disk is one proposed mechanism to form planetesimals, the building blocks of rocky planets and gas giant cores. Self-gravity competes against the Kelvin-Helmholtz instability (KHI): gradients in dust content drive a vertical shear which risks overturning the dusty subdisk and forestalling GI. To understand the conditions under which the disk can resist the KHI, we perform 3D simulations of stratified subdisks in the limit that dust particles are small and aerodynamically well coupled to gas. This limit screens out the streaming instability and isolates the KHI. Each subdisk is assumed to have a vertical density profile given by a spatially constant Richardson number Ri. We vary Ri and the midplane dust-to-gas ratio mu and find that the critical Richardson number dividing KH-unstable from KH-stable flows is not unique; rather Ri_crit grows nearly linearly with mu for mu=0.3-10. Only for disks of bulk solar metallicity is Ri_crit ~ 0.2, close to the classical value. Our results suggest that a dusty sublayer can gravitationally fragment and presumably spawn planetesimals if embedded within a solar metallicity gas disk ~4x more massive than the minimum-mass solar nebula; or a minimum-mass disk having ~3x solar metallicity; or some intermediate combination of these two possibilities. Gravitational instability seems possible without resorting to the streaming instability or to turbulent concentration of particles.

Dust at the midplane of a circumstellar disk can become gravitationally unstable and fragment into planetesimals if the local dust-to-gas density ratio mu is sufficiently high. We simulate how dust settles in passive disks and ask how high mu can become. We settle the dust using a 1D code and test for dynamical stability using a 3D shearing box code. This scheme allows us to explore the behavior of small particles having short but non-zero stopping times in gas: 0 < t_stop << the orbital period. The streaming instability is thereby filtered out. Dust settles until shearing instabilities in the edges of the dust layer threaten to overturn the entire layer. In this state of marginal stability, mu=2.9 for a disk whose bulk (height-integrated) metallicity is solar. For a disk whose bulk metallicity is 4x solar, mu reaches 26.4. These maximum values of mu, which depend on the background radial pressure gradient, are so large that gravitational instability of small particles is viable in disks whose bulk metallicities are just a few (<4) times solar. Earlier studies assumed that dust settles until the Richardson number Ri is spatially constant. Our simulations are free of this assumption but provide support for it within the dust layer's edges, with the proviso that Ri increases with bulk metallicity in the same way that we found in Paper I. Only modest enhancements in bulk metallicity are needed to spawn planetesimals directly from small particles.

Motivated by the increasing interest in models which consider scalar fields as viable dark matter candidates, we have constructed a generalization of relativistic Boson Stars (BS) composed of two coexisting states of the scalar field, the ground state and the first excited state. We have studied the dynamical evolution of these Multi-state Boson Stars (MSBS) under radial perturbations, using numerical techniques. We show that stable MSBS can be constructed, when the number of particles in the first excited state, N2, is smaller than the number of particles in the ground state, N1. On the other hand, when N2 > N1, the configurations are initially unstable. However, they evolve and settle down into stable configurations. In the stabilization process, the initially ground state is excited and ends in a first excited state, whereas the initially first excited state ends in a ground state. During this process, both states emit scalar field radiation, decreasing their number of particles. This behavior shows that even though BS in the first excited state are intrinsically unstable under finite perturbations, the configuration resulting from the combination of this state with the ground state produces stable objects. Finally we show in a qualitative way, that stable MSBS could be realistic models of dark matter galactic halos, as they produce rotation curves that are flatter at large radii than the rotation curves produced by BS with only one state.

The current understanding of the turbulent dissipation in stellar convective zones is based on the assumption that the turbulence follows Kolmogorov scaling. This assumption is valid for some cases in which the time frequency of the external shear is high (e.g., solar p modes). However, for many cases of astrophysical interest (e.g., binary orbits, stellar pulsations, etc.), the timescales of interest lie outside the regime of applicability of Kolmogorov scaling. We present direct calculations of the dissipation efficiency of the turbulent convective flow in this regime, using simulations of anelastic convection with external forcing. We show that the effects of the turbulent flow are well represented by an effective viscosity coefficient, we provide the values of the effective viscosity as a function of the perturbation frequency and compare our results to the perturbative method for finding the effective viscosity of Penev et al. that can be applied to actual simulations of the surface convective zones of stars.

We have adapted the anelastic spectral code of Barranco & Marcus (2006) to simulate a turbulent convective layer with the intention of studying the effectiveness of turbulent eddies in dissipating external shear (e.g. tides). We derive the anelastic equations, show the time integration scheme we use to evolve these equations and present the tests we ran to confirm that our code does what we expect. Further we apply a perturbative approach to find an approximate scaling of the effective eddy viscosity with frequency, and find that it is in general agreement with an estimate obtained by applying the same procedure to a realistic simulation of the upper layers of the solar convective zone.

We present the results of a deep (15 ~< r ~< 23), 20 night survey for transiting planets in the intermediate age open cluster M37 (NGC 2099) using the Megacam wide-field mosaic CCD camera on the 6.5m MMT. We do not detect any transiting planets among the ~1450 observed cluster members. We do, however, identify a ~ 1 R_J candidate planet transiting a ~ 0.8 Msun Galactic field star with a period of 0.77 days. The source is faint (V = 19.85 mag) and has an expected velocity semi-amplitude of K ~ 220 m/s (M/M_J). We conduct Monte Carlo transit injection and recovery simulations to calculate the 95% confidence upper limit on the fraction of cluster members and field stars with planets as a function of planetary radius and orbital period. Assuming a uniform logarithmic distribution in orbital period, we find that < 1.1%, < 2.7% and < 8.3% of cluster members have 1.0 R_J planets within Extremely Hot Jupiter (EHJ, 0.4 < T < 1.0 day), Very Hot Jupiter (VHJ, 1.0 < T < 3.0 days) and Hot Jupiter (HJ, 3.0 < T < 5.0 days) period ranges respectively. For 0.5 R_J planets the limits are < 3.2%, and < 21% for EHJ and VHJ period ranges, while for 0.35 R_J planets we can only place an upper limit of < 25% on the EHJ period range. For a sample of 7814 Galactic field stars, consisting primarily of FGKM dwarfs, we place 95% upper limits of < 0.3%, < 0.8% and < 2.7% on the fraction of stars with 1.0 R_J EHJ, VHJ and HJ assuming the candidate planet is not genuine. If the candidate is genuine, the frequency of ~ 1.0 R_J planets in the EHJ period range is 0.002% < f_EHJ < 0.5% with 95% confidence. We place limits of < 1.4%, < 8.8% and < 47% for 0.5 R_J planets, and a limit of < 16% on 0.3 R_J planets in the EHJ period range. This is the first transit survey to place limits on the fraction of stars with planets as small as Neptune.

In the course of conducting a deep (14.5 ~< r ~< 23), 20 night survey for transiting planets in the rich ~550 Myr old open cluster M37 we have measured the rotation periods of 575 stars which lie near the cluster main sequence, with masses 0.2 Msun ~< M ~< 1.3 Msun. This is the largest sample of rotation periods for a cluster older than 500 Myr. Using this rich sample we investigate a number of relations between rotation period, color and the amplitude of photometric variability. Stars with M >~ 0.8 Msun show a tight correlation between period and mass with heavier stars rotating more rapidly. There is a group of 4 stars with P > 15 days that fall well above this relation, which, if real, would present a significant challenge to theories of stellar angular momentum evolution. Below 0.8 Msun the stars continue to follow the period-mass correlation but with a broad tail of rapid rotators that expands to shorter periods with decreasing mass. We combine these results with observations of other open clusters to test the standard theory of lower-main sequence stellar angular momentum evolution. We find that the model reproduces the observations for solar mass stars, but discrepancies are apparent for stars with 0.6 ~< M ~< 1.0 Msun. We also find that for late-K through early-M dwarf stars in this cluster rapid rotators tend to be bluer than slow rotators in B-V but redder than slow rotators in V-I_{C}. This result supports the hypothesis that the significant discrepancy between the observed and predicted temperatures and radii of low-mass main sequence stars is due to stellar activity.

In this work we estimate the radius and the mass of a self-gravitating system made of axions. The quantum axion field satisfies the Klein-Gordon equation in a curved space-time and the metric components of this space-time are solutions to the Einstein equations with a source term given by the vacuum expectation value of the energy-momentum operator constructed from the axion field. As a first step towards an axion star we consider the up to the sixth term in the axion potential expansion. We found that axion stars would have masses of the order of asteroids and radius of the order of few centimeters.

We have conducted a deep ($15 \la r \la 23$), 20 night survey for transiting planets in the intermediate age ($\sim 550 {\rm Myr}$) open cluster M37 (NGC 2099) using the Megacam wide-field mosaic CCD camera on the 6.5m Multiple Mirror Telescope (MMT). In this paper we present a catalog and light curves for 1445 variable stars; 1430 (99%) of these are new discoveries. We have discovered 20 new eclipsing binaries and 31 new short period ($P < 1 {\rm day}$) pulsating stars. The bulk of the variables are most likely rapidly rotating young low-mass stars, including a substantial number ($\ga 500$) that are members of the cluster. We identify and analyze five particularly interesting individual variables including a previously identified variable which we suggest is probably a hybrid $\gamma$ Doradus/$\delta$ Scuti pulsator, two possible quiescent cataclysmic variables, a detached eclipsing binary (DEB) with at least one $\gamma$ Doradus pulsating component (only the second such variable found in an eclipsing binary), and a low mass ($M_{P} \sim M_{S} \sim 0.6 M_{\odot}$) DEB that is a possible cluster member. A preliminary determination of the physical parameters for the DEB+$\gamma$ Doradus system yields $M_{P} = 1.58 \pm 0.04 M_{\odot}$, $M_{S} = 1.58 \pm 0.04 M_{\odot}$, $R_{P} = 1.39 \pm 0.07 R_{\odot}$ and $R_{S} = 1.38 \pm 0.07 R_{\odot}$.

We have conducted a deep ($15 \la r \la 23$), 20 night survey for transiting planets in the intermediate age open cluster M37 (NGC 2099) using the Megacam wide-field mosaic CCD camera on the 6.5m Multiple Mirror Telescope (MMT). In this paper we describe the observations and data reduction procedures for the survey and analyze the stellar content and dynamical state of the cluster. By combining high resolution spectroscopy with existing $BVI_{C}K_{S}$ and new $gri$ color magnitude diagrams we determine the fundamental cluster parameters: $t = 485 \pm 28$ Myr without overshooting ($t = 550 \pm 30 {\rm Myr}$ with overshooting), $E(B-V) = 0.227 \pm 0.038$, $(m-M)_{V} = 11.57 \pm 0.13$ and $[M/H] = +0.045 \pm 0.044$ which are in good agreement with, though more precise than, previous measurements. We determine the mass function down to $0.3 M_{\odot}$ and use this to estimate the total cluster mass of $3640 \pm 170 M_{\odot}$.

As dust settles in a protoplanetary disk, a vertical shear develops because the dust-rich gas in the midplane orbits at a rate closer to true Keplerian than the slower-moving dust-depleted gas above and below. A classical analysis (neglecting the Coriolis force and differential rotation) predicts that Kelvin-Helmholtz instability occurs when the Richardson number of the stratified shear flow is below roughly one-quarter. However, earlier numerical studies showed that the Coriolis force makes layers more unstable, whereas horizontal shear may stabilize the layers. Simulations with a 3D spectral code were used to investigate these opposing influences on the instability in order to resolve whether such layers can ever reach the dense enough conditions for the onset of gravitational instability. I confirm that the Coriolis force, in the absence of radial shear, does indeed make dust layers more unstable, however the instability sets in at high spatial wavenumber for thicker layers. When radial shear is introduced, the onset of instability depends on the amplitude of perturbations: small amplitude perturbations are sheared to high wavenumber where further growth is damped; whereas larger amplitude perturbations grow to magnitudes that disrupt the dust layer. However, this critical amplitude decreases sharply for thinner, more unstable layers. In 3D simulations of unstable layers, turbulence mixes the dust and gas, creating thicker, more stable layers. I find that layers with minimum Richardson numbers in the approximate range 0.2 -- 0.4 are stable in simulations with horizontal shear.

We are conducting a transit survey of the open cluster M37 using the Megacam instrument on the 6.5 m Multiple-Mirror Telescope. We have obtained ~4500 images of this cluster over 18.5 nights and have achieved the precision necessary to detect planets smaller than Saturn. In this presentation we provide an overview of the project, describe the ongoing data reduction/analysis and present some of our preliminary results.

We have developed a three-dimensional (3D) spectral hydrodynamic code to study vortex dynamics in rotating, shearing, stratified systems (e.g. the atmosphere of gas giant planets, protoplanetary disks around newly forming protostars). The time-independent background state is stably stratified in the vertical direction and has a unidirectional linear shear flow aligned with one horizontal axis. Superposed on this background state is an unsteady, subsonic flow that is evolved with the Euler equations subject to the anelastic approximation to filter acoustic phenomena. A Fourier-Fourier basis in a set of quasi-Lagrangian coordinates that advect with the background shear is used for spectral expansions in the two horizontal directions. For the vertical direction, two different sets of basis functions have been implemented: (1) Chebyshev polynomials on a truncated, finite domain, and (2) rational Chebyshev functions on an infinite domain. Use of this latter set is equivalent to transforming the infinite domain to a finite one with a cotangent mapping, and using cosine and sine expansions in the mapped coordinate. The nonlinear advection terms are time integrated explicitly, whereas the Coriolis force, buoyancy terms, and pressure/enthalpy gradient are integrated semi-implicitly. We show that internal gravity waves can be damped by adding new terms to the Euler equations. The code exhibits excellent parallel performance with the Message Passing Interface (MPI). As a demonstration of the code, we simulate the merger of two 3D vortices in the midplane of a protoplanetary disk.

We present the results of high-resolution, three-dimensional (3D) hydrodynamic simulations of the dynamics and formation of coherent, long-lived vortices in stably-stratified protoplanetary disks. Tall, columnar vortices that extend vertically through many scale heights in the disk are unstable to small perturbations; such vortices cannot maintain vertical alignment over more than a couple scale heights and are ripped apart by the Keplerian shear. Short, finite-height vortices that extend only one scale height above and below the midplane are also unstable, but for a different reason: we have isolated an antisymmetric (with respect to the midplane) eigenmode that grows with an e-folding time of only a few orbital periods; the nonlinear evolution of this instability leads to the destruction of the vortex. Serendipitously, we observe the formation of 3D vortices that are centered not in the midplane, but at one to three scale heights above and below. Breaking internal gravity waves create vorticity; anticyclonic regions of vorticity roll-up and coalesce into new vortices, whereas cyclonic regions shear into thin azimuthal bands. Unlike the midplane-centered vortices that were placed ad hoc in the disk and turned out to be linearly unstable, the off-midplane vortices form naturally out of perturbations in the disk, and are stable and robust for many hundreds of orbits.