A third object in the vicinity of a binary system causes variations in the eccentricity and the inclination of the binary through the Kozai-Lidov effect. We examine if such variations leave a detectable imprint on the gravitational waves of a binary consisting of intermediate mass black holes and stellar mass objects. As a proof of concept, we present an example where LISA may detect the Kozai-Lidov modulated gravitational wave signals of such sources from at least a distance of 1Mpc if the perturbation is caused by a supermassive black hole tertiary. Although the quick pericenter precession induced by general relativity significantly reduces the appropriate parameter space for this effect by quenching the Kozai-Lidov oscillations, we still find reasonable parameters where the Kozai-Lidov effect may be detected with high signal-to-noise ratios.

The central 0.1 parsecs of the Milky Way host a supermassive black hole identified with the position of the radio and infrared source Sagittarius A*, a cluster of young, massive stars (the S stars) and various gaseous features. Recently, two unusual objects have been found to be closely orbiting Sagittarius A*: the so-called G sources, G1 and G2. These objects are unresolved (having a size of the order of 100 astronomical units, except at periapse, where the tidal interaction with the black hole stretches them along the orbit) and they show both thermal dust emission and line emission from ionized gas. G1 and G2 have generated attention because they appear to be tidally interacting with the supermassive Galactic black hole, possibly enhancing its accretion activity. No broad consensus has yet been reached concerning their nature: the G objects show the characteristics of gas and dust clouds but display the dynamical properties of stellar-mass objects. Here we report observations of four additional G objects, all lying within 0.04 parsecs of the black hole and forming a class that is probably unique to this environment. The widely varying orbits derived for the six G objects demonstrate that they were commonly but separately formed.

The hierarchical nature of galaxy formation suggests that a supermassive black hole binary could exist in our galactic center. We propose a new approach to constraining the possible orbital configuration of such a binary companion to the galactic center black hole Sgr A* through the measurement of stellar orbits. Focusing on the star S0-2, we show that requiring its orbital stability in the presence of a companion to Sgr A* yields stringent constraints on the possible configurations of such a companion. Furthermore, we show that precise measurements of {\it time variations} in the orbital parameters of S0-2 could yield stronger constraints. Using existing data on S0-2 we derive upper limits on the binary black hole separation as a function of the companion mass. For the case of a circular orbit, we can rule out a 10^5 M_sun companion with a semimajor axis greater than 170 astronomical units or 0.8 mpc. This is already more stringent than bounds obtained from studies of the proper motion of Sgr A*. Including other stars orbiting the galactic center should yield stronger constraints that could help uncover the presence of a companion to Sgr A*. We show that a companion can also affect the accretion process, resulting in a variability which may be consistent with the measured infrared flaring timescales and amplitudes. Finally, if such a companion exists, it will emit gravitational wave radiation, potentially detectable with LISA.

A natural consequence of the galaxy formation paradigm is the existence of supermassive black hole (SMBH) binaries. Gravitational perturbations from a far-away SMBH companion can induce high orbital eccentricities on dark matter particles orbiting the primary SMBH via the eccentric Kozai-Lidov mechanism. This process yields an influx of dark matter particles into the primary SMBH ergosphere, where test particles linger for long timescales. This influx results in high self-gravitating densities, forming a dark matter clump that is extremely close to the SMBH. In such a situation, the gravitational wave emission between the dark matter clump and the SMBH is potentially detectable by LISA. If dark matter self-annihilates, the high densities of the clump will result in a unique co-detection of gravitational wave emission and high energy electromagnetic signatures.

The binary star Par 1802 in the Orion Nebula presents an interesting puzzle in the field of stellar dynamics and evolution. Binary systems such as Par 1802 are thought to form from the same natal material and thus the stellar members are expected to have very similar physical attributes. However, Par 1802's stars have significantly different temperatures despite their identical (within $3\%$) masses of about $0.39$ solar mass. The leading proof-of-concept idea is that a third companion gravitationally induced the two stars to orbit closer than their Roche-limit, which facilitated heating through tidal effects. Here we expand on this idea and study the three-body dynamical evolution of such a system, including tidal and pre-main-sequence evolution. We also include tidal heating and mass transfer at the onset of Roche-limit Crossing. We show, as a proof-of-concept, that mass transfer combined with tidal heating can naturally explain the observed temperature discrepancy. We also predict the orbital configuration of the possible tertiary companion. Finally we suggest that the dynamical evolution of such a system has pervasive consequences. We expect an abundance of systems to undergo mass transfer during their pre-main-sequence time, which can cause temperature differences.

The Central Molecular Zone, which spans up to about 500 pc around the center of our Galaxy, consists of a turbulent and complex gaseous environment predominantly in the form of relatively dense molecular gas. Here we explore the effect of turbulence on the gas dynamics of this region. We use Smoothed Particle Hydrodynamics simulations, which include self-gravity and a recipe for turbulence driving, to follow the long-term evolution of the gas in the rotating bar potential of the Galaxy. We find that before the system reaches a steady state, a twisted ring is formed, similar to the observed one. As expected, this feature later dissipates leaving a smooth gaseous disk. We also find that the average densities resulting from our simulations are, in the absence of major perturbations, consistent with the observed ones. Our results suggest that the current gas in the Central Molecular Zone is not in a steady state.

Exoplanets have been observed around stars at all stages of stellar evolution, in many cases orbiting in configurations that will eventually lead to the planets being engulfed or consumed by their host stars, such as Hot Jupiters or ultra-short period planets. Furthermore, objects such as polluted white dwarfs provide strong evidence that the consumption of planets by stars is a common phenomenon. This consumption causes several significant changes in the stellar properties, such as changes to the stellar spin, luminosity, chemical composition, or mass loss processes. Here, we explore this wide variety of effects for a comprehensive range of stellar and planetary masses and stages of stellar evolution, from the main sequence over red giants to the white dwarfs. We determine that planet consumption can cause transient luminosity features that last on the order of centuries to millennia, and that the post-consumption stellar spins can often reach break-up speeds. Furthermore, stellar moss loss can be caused by this spin-up, as well as through surface grazing interactions, leading to to the formation of unusual planetary nebula shapes or collimated stellar gas ejections. Our results highlight several observable stellar features by which the presence or previous existence of a planet around a given star can be deduced. This will provide future observational campaigns with the tools to better constrain exoplanet demographics, as well as planetary formation and evolution histories.

General Relativity predicts that a star passing close to a supermassive black hole should exhibit a relativistic redshift. We test this using observations of the Galactic center star S0-2. We combine existing spectroscopic and astrometric measurements from 1995-2017, which cover S0-2's 16-year orbit, with measurements in 2018 March to September which cover three events during its closest approach to the black hole. We detect the combination of special relativistic- and gravitational-redshift, quantified using a redshift parameter, $\Upsilon$. Our result, $\Upsilon=0.88 \pm 0.17$, is consistent with General Relativity ($\Upsilon=1$) and excludes a Newtonian model ($\Upsilon=0$ ) with a statistical significance of 5 $\sigma$.

The current disparity in computational knowledge is a critical hindrance to the diversity and success of the field. Recommendations are outlined for policies and funding models to enable the growth and retention of a new generation of computational researchers that reflect the demographics of the undergraduate population in Astronomy and Physics.

A recent study using $Hubble$ $Space$ $Telescope$ observations found periodic, high-speed, collimated ejections (or "bullets") from the star V Hya. The authors of that study proposed a model associating these bullets with the periastron passage of an unseen, substellar companion in an eccentric orbit and with an orbital period of $\sim8$ yrs. Here we propose that V Hya is part of a triple system, with a substellar companion having an orbital period of $\sim8$ yrs, and a tertiary object on a much wider orbit. In this model, the more distant object causes high-eccentricity excitations on the substellar companion's orbit via the Eccentric Kozai-Lidov mechanism. These eccentricities can reach such high values that they lead to Roche-lobe crossing, producing the observed bullet ejections via a strongly enhanced accretion episode. For example, we find that a ballistic bullet ejection mechanism can be produced by a brown-dwarf-mass companion, while magnetically driven outflows are consistent with a Jovian-mass companion. Finally, we suggest that the distant companion may reside at few a hundred AU on an eccentric orbit.

Supersonically induced gas objects (SIGOs) with little to no dark matter component are predicted to exist in patches of the Universe with non-negligible relative velocity between baryons and the dark matter at the time of recombination. Using {\sc arepo} hydrodynamic simulations we find that the gas densities inside these objects are high enough to allow stars to form. An estimate of the luminosity of the first star clusters formed within these SIGOs suggests that they may be observed at high redshift using future HST and JWST observations. Furthermore, our simulations indicate that SIGOs lie in a distinct place in the luminosity-radius parameter space, which can be used observationally to distinguish SIGOs from dark-matter hosting gas systems. Finally, as a proof-of-concept, we model star formation before reionization and evolve these systems to current times. We find that SIGOs occupy a similar part of the magnitude-radius parameter space as globular clusters. These results suggest that SIGOs may be linked with present-day metal-poor local globular clusters. Since the relative velocity between the baryons and dark matter is coherent over a few Mpc scales, we predict that if this is the dominant mechanism for the formation of globular clusters, their abundance should vary significantly over these scales.

The Galactic Center (GC) is dominated by the gravity of a super-massive black hole (SMBH), Sagittarius A$^*$, and is suspected to contain a sizable population of binary stars. Such binaries form hierarchical triples with the SMBH, undergoing Eccentric Kozai-Lidov (EKL) evolution, which can lead to high eccentricity excitations for the binary companions' mutual orbit. This effect can lead to stellar collisions or Roche-lobe crossings, as well as orbital shrinking due to tidal dissipation. In this work we investigate the dynamical and stellar evolution of such binary systems, especially with regards to the binaries' post-main-sequence evolution. We find that the majority of binaries (~75%) is eventually separated into single stars, while the remaining binaries (~25%) undergo phases of common-envelope evolution and/or stellar mergers. These objects can produce a number of different exotic outcomes, including rejuvenated stars, G2-like infrared-excess objects, stripped giant stars, Type Ia supernovae (SNe), cataclysmic variables (CVs), symbiotic binaries (SBs), or compact object binaries. We estimate that, within a sphere of 250 Mpc radius, about 7.5 to 15 Type Ia SNe per year should occur in galactic nuclei due to this mechanism, potentially detectable by ZTF and ASAS-SN. Likewise we estimate that, within a sphere of 1 Gpc$^3$ volume, about 10 to 20 compact object binaries form per year that could become gravitational wave sources. Based on results of EKL-driven compact object binary mergers in galactic nuclei by Hoang at al. (2018), this compact object binary formation rate translates to about 15 to 30 events per year detectable by Advanced LIGO.

We analyze the behavior of the argument of pericenter $\omega_2$ of an outer particle in the elliptical restricted three-body problem, focusing on the inverse Lidov-Kozai resonance. First, we calculate the contribution of the terms of quadrupole, octupole, and hexadecapolar order of the secular approximation of the potential to the outer particle's $\omega_2$ precession rate $(d\omega_2/d\tau)$. Then, we derive analytical criteria that determine the vanishing of the $\omega_2$ quadrupole precession rate $(d\omega_2/d\tau)_{\text{quad}}$ for different values of the inner perturber's eccentricity $e_1$. Finally, we use such analytical considerations and describe the behavior of $\omega_2$ of outer particles extracted from N-body simulations. Our analytical study indicates that the values of the inclination $i_2$ and the ascending node longitude $\Omega_2$ associated with the outer particle that vanish $(d\omega_2/d\tau)_{\text{quad}}$ strongly depend on the eccentricity $e_1$ of the inner perturber. In fact, if $e_1 <$ 0.25 ($>$ 0.40825), $(d\omega_2/d\tau)_{\text{quad}}$ is only vanished for particles whose $\Omega_2$ circulates (librates). For $e_1$ between 0.25 and 0.40825, $(d\omega_2/d\tau)_{\text{quad}}$ can be vanished for any particle for a suitable selection of pairs ($\Omega_2$, $i_2$). Our analysis of the N-body simulations shows that the inverse Lidov-Kozai resonance is possible for small, moderate and high values of $e_1$. Moreover, such a resonance produces distinctive features in the evolution of a particle in the ($\Omega_2$, $i_2$) plane. In fact, if $\omega_2$ librates and $\Omega_2$ circulates, the extremes of $i_2$ at $\Omega_2 =$ 90$^{\circ}$ and 270$^{\circ}$ do not reach the same value, while if $\omega_2$ and $\Omega_2$ librate, the evolutionary trajectory of the particle in the ($\Omega_2$, $i_2$) plane evidences an asymmetry respect to $i_2 =$ 90$^{\circ}$.

Stellar-mass black hole binaries (BHBs) near supermassive black holes (SMBH) in galactic nuclei undergo eccentricity oscillations due to gravitational perturbations from the SMBH. Previous works have shown that this channel can contribute to the overall BHB merger rate detected by the Laser Interferometer Gravitational-Wave Observatory (LIGO) and Virgo Interferometer. Significantly, the SMBH gravitational perturbations on the binary's orbit may produce eccentric BHBs which are expected to be visible using the upcoming Laser Interferometer Space Antenna (LISA) for a large fraction of their lifetime before they merge in the LIGO/Virgo band. For a proof-of-concept, we show that the eccentricity oscillations of these binaries can be detected with LISA for BHBs in the local universe up to a few Mpcs, with observation periods shorter than the mission lifetime, thereby disentangling this merger channel from others. The approach presented here is straightforward to apply to a wide variety of compact object binaries with a tertiary companion.

At least $70\%$ of massive OBA-type stars reside in binary or higher-order systems. The dynamical evolution of these systems can lend insight into the origins of extreme phenomena such as X-ray binaries and gravitational wave sources. In one such dynamical process, the Eccentric Kozai-Lidov (EKL) Mechanism, a third companion star alters the secular evolution of a binary system. For dynamical stability, these triple systems must have a hierarchical configuration. We explore the effects of a distant third companion's gravitational perturbations on a massive binary's orbital configuration before significant stellar evolution has taken place ($\leq 10$ Myr). We include tidal dissipation and general relativistic precession. With large ($38,000$ total) Monte-Carlo realizations of massive hierarchical triples, we characterize imprints of the birth conditions on the final orbital distributions. Specifically, we find that the final eccentricity distribution over the range $0.1-0.7$ is an excellent indicator of its birth distribution. Furthermore, we find that the period distributions have a similar mapping for wide orbits. Finally, we demonstrate that the observed period distribution for approximately $10$ Myr-old massive stars is consistent with EKL evolution.

As the closest example of a galactic nucleus, the Galactic center (GC) presents an exquisite laboratory for learning about supermassive black holes (SMBH) and their environment. We describe several exciting new research directions that, over the next 10 years, hold the potential to answer some of the biggest scientific questions raised in recent decades: Is General Relativity (GR) the correct description for supermassive black holes? What is the nature of star formation in extreme environments? How do stars and compact objects dynamically interact with the supermassive black hole? What physical processes drive gas accretion in low-luminosity black holes? We describe how the high sensitivity, angular resolution, and astrometric precision offered by the next generation of large ground-based telescopes with adaptive optics will help us answer these questions. First, it will be possible to obtain precision measurements of stellar orbits in the Galaxy's central potential, providing both tests of GR in the unexplored regime near a SMBH and measurements of the extended dark matter distribution that is predicted to exist at the GC. Second, we will probe stellar populations at the GC to significantly lower masses than are possible today, down to brown dwarfs. Their structure and dynamics will provide an unprecedented view of the stellar cusp around the SMBH and will distinguish between models of star formation in this extreme environment. This increase in depth will also allow us to measure the currently unknown population of compact remnants at the GC by observing their effects on luminous sources. Third, uncertainties on the mass of and distance to the SMBH can be improved by a factor of $\sim$10. Finally, we can also study the near-infrared accretion onto the black hole at unprecedented sensitivity and time resolution, which can reveal the underlying physics of black hole accretion.

Most massive stars, if not all, are in binary configuration or higher multiples. These massive stars undergo supernova explosions and end their lives as either black holes or neutron stars. Recent observations have suggested that neutron stars and perhaps even black holes receive large velocity kicks at birth. Such natal kicks and the sudden mass loss can significantly alter the orbital configuration of the system. Here we derive general analytical expressions that describe the effects of natal kicks in binaries on hierarchical triple systems. We explore several proof-of-concept applications such as black hole and neutron stars binaries and X-ray binaries with either stellar or Supermassive Black Hole (SMBH) companions on a wide orbit. Kicks can disrupt the hierarchical configuration, although it is harder to escape the potential well of an SMBH. Some binary systems do escape the SMBH system resulting in hyper-velocity binary system. Furthermore, kicks can result in increasing or decreasing the orbital separations. Decreasing the orbital separation may have significant consequences in these astrophysical systems. For example, shrinking the separation post-supernova kick can lead to the shrinking of an inner compact binary that then may merge via Gravitational Wave (GW) emission. This process yields a supernova that is shortly followed by a possible GW-LIGO event. Interestingly, we find that the natal kick can result in shrinking the outer orbit, and the binary may cross the tertiary Roche limit, breaking up the inner binary. Thus, in the case of SMBH companion, this process can lead to either a tidal disruption event or a GW-LISA detection event with a supernova precursor.

We introduce the "Supersonic Project," aimed at investigating the effects of the supersonic relative velocity between dark matter (DM) and baryons at high redshift using a combination of analytical calculations and cosmological simulations. In this paper, we study the effect of this stream velocity on the angular momentum of the first structures in the early Universe using simulations. We focus on DM haloes and their gas component as well as the recently predicted supersonically-induced gas objects (SIGOs) that arise as a result of the stream velocity phase shift. We find that the spin parameter of the gas component in these first haloes is increased with the stream velocity. Moreover, we find that when the stream velocity is taken into account, the angular momentum vectors of the DM component and the gas component are typically misaligned and this misalignment angle has a nearly isotropic distribution. The spin parameter value of the gas component is higher than in the no stream velocity case, which even in the absence of cooling, may result in more prolate objects. We also generalize the spin parameter to the SIGOs and find that they typically have a larger spin parameter with respect to their dark matter counterparts and that there is no correlation of the spin parameter and the prolateness of such structures. We speculate that SIGOs may be observed as very low luminosity objects in the early Universe and may serve as potential progenitors of Little Blue Dot-like systems.

Hot Jupiters (HJs), gas giant planets orbiting their host stars with periods on the order of days, commonly occur in the Galaxy, including relatively massive ($1.6-2.4$ $M_\odot$, i.e., A-type main sequence stars) and evolved stars. The majority of A-type main sequence stars have stellar binary companions, which can strongly affect the dynamical evolution of planets around either star. In this work, we investigate the effects of gravitational perturbations by a far away stellar companion on the orbital evolution of gas giant planets orbiting A-type stars, the so-called Eccentric Kozai-Lidov (EKL) mechanism, including the effects of general relativity, post-main sequence stellar evolution, and tides. We find that only $0.15~\%$ of A-type stars will host HJs during their main sequence lifetime. However, we also find a new class of planets, Temporary Hot Jupiters (THJs), that form during the post-main sequence lifetime of about $3.7~\%$ of former A-type main sequence stars. These THJs orbit on periods of tens to a hundred days and only exist for a few $100{,}000$ years before they are engulfed, but they reach similar temperatures as `classical' HJs due to the increased stellar luminosities. THJs' spin-orbit angles will mostly be misaligned. THJ effects on the host stars' evolution could also be observable for longer than a few $100{,}000$ years. Overall, we find that approximately $70~\%$ of all gas giant planets orbiting A-type stars will eventually be destroyed or engulfed by their star, about $25~\%$ during the main sequence lifetime, about $45~\%$ during post-main sequence evolution.

One of the predictions of high eccentricity planetary migration is that many planets will end up plunging into their host stars. We investigate the consequence of planetary mergers on their stellar hosts' spin-period. Energy and angular momentum conservation yield that a planet consumption by a star will spin-up the star. We find that our calculations align with the observed bifurcation in the stellar spin-period in young clusters. For example, after a Sun-like star has eaten a Jupiter-mass planet it will spin up by ~60% (i.e., spin-period is reduced by ~60%), causing an apparent gap in the stellar spin period, between stars that consumed a planet and those that did not. The spun-up star will later spin down due to magnetic braking, consistent with the disappearance of this bifurcation in clusters (>300Myr). The agreement between the calculations presented here, and the observed spin-period color diagram of stars in young clusters provides circumstantial evidence that planetary accretion onto their host stars is a generic feature of planetary-system evolution.

Recent studies have analyzed the dynamical evolution of outer small body populations under the effects of an eccentric inner massive perturber. Such outer reservoirs are composed of particles on prograde and retrograde orbits, as well as particles whose orbit flips from prograde to retrograde and back again showing a coupling between the inclination i and the ascending node longitude \Omega (Type-F particles). We analyze the role of the General Relativity (GR) on the dynamics of outer particles under the influence of an inner eccentric Jupiter-mass planet produced by a planetary scattering event. In particular, we study how the GR affects the dynamical evolution of the outer Type-F particles, which experience an eccentric Lidov-Kozai mechanism. We carry out N-body simulations with and without GR effects. When the GR is included, the extreme values of \Omega are obtained for retrograde inclinations, while the minimum and maximum inclinations allowed for Type-F particles increase in comparison with that derived without GR effects. According to this, if the GR is included in the simulations, the range of prograde (retrograde) inclinations of the libration region is reduced (increased) respect to that obtained in absence of GR. We find two new class of particles when the GR effects are included in the simulations. On the one hand, particles whose orbital plane flips from prograde to retrograde and back again without experiencing a coupling between i and \Omega. On the other hand, retrograde particles that show a strong coupling between i and \Omega. We infer that the GR may significantly modify the dynamical properties of the outer reservoirs that evolve under the effects of an eccentric inner perturber.

Nuclear star clusters around massive black holes are expected to be abundant in stellar mass black holes and black hole binaries. These binaries form a hierarchical triple system with the massive black hole at the center. Gravitational perturbations from the massive black hole can cause high eccentricity excitation. During this process, the eccentricity may approach unity, and the pericenter distance may become sufficiently small that gravitational wave emission drives the binary to merge. In this paper, we consider a simple proof of concept and explore the effect of the eccentric Kozai-Lidov mechanism for unequal mass binaries. We perform a set of Monte Carlo simulations on BH-BH binaries in galactic nuclei with quadrupole and octupole-level secular perturbations, general relativistic precession, and gravitational wave emission. For a nominal number of steady-state BH-BH binaries, our model gives a total merger rate $\sim 1 - 3$$Gpc^{-3} yr^{-1}$, depending on the assumed density profile. Thus, our model potentially competes with other dynamical mechanisms, such as the dynamical formations and mergers of BH binaries in globular clusters or dense nuclear clusters without a massive black hole. We provide predictions for the distributions of these LIGO sources in galactic nuclei.

Recent observational campaigns have shown that multi-planet systems seem to be abundant in our Galaxy. Moreover, it seems that these systems might have distant companions, either planets, brown-dwarfs or other stellar objects. These companions might be inclined with respect to the inner planets, and could potentially excite the eccentricities of the inner planets through the Eccentric Kozai-Lidov mechanism. These eccentricity excitations could perhaps cause the inner orbits to cross, disrupting the inner system. We study the stability of two-planet systems in the presence of a distant, inclined, giant planet. Specifically, we derive a stability criterion, which depends on the companion's separation and eccentricity. We show that our analytic criterion agrees with the results obtained from numerically integrating an ensemble of systems. Finally, as a potential proof-of-concept, we provide a set of predictions for the parameter space that allows the existence of planetary companions for the Kepler-419, Kepler-56, Kepler-448, Kepler-88, Kepler-109, and Kepler-36 systems.

The population of young stars near the Supermassive black hole (SMBH) in the Galactic Center (GC) has presented an unexpected challenge to theories of star formation. Kinematics measurements of these stars have revealed a stellar disk structure (with an apparent 20% disk membership) that has provided important clues to the origin of these mysterious young stars. However many of the apparent disk properties are difficult to explain, including the low disk membership fraction and the high eccentricities, given the youth of this population. Thus far, all efforts to derive the properties of this disk have made the simplifying assumption that stars at the GC are single stars. Nevertheless, stellar binaries are prevalent in our Galaxy, and recent investigations suggested that they may also be abundant in the Galactic Center. Here we show that binaries in the disk can largely alter the apparent orbital properties of the disk. The motion of binary members around each other adds a velocity component, which can be comparable to the magnitude of the velocity around the SMBH in the GC. Thus neglecting the contribution of binaries can significantly vary the inferred stars' orbital properties. While the disk orientation is unaffected the apparent disk's 2D width is increased to about 11.2deg, similar to the observed width. For a population of stars orbiting the SMBH with zero eccentricity, unaccounted for binaries will create a wide apparent eccentricity distribution with an average of 0.23.This is consistent with the observed average eccentricity of the stars' in the disk. We suggest that this high eccentricity value, which poses a theoretical challenge, may be an artifact of binary stars. Finally our results suggest that the actual disk membership might be significantly higher than the one inferred by observations that ignore the contribution of binaries, alleviating another theoretical challenge.

Compact stellar binaries are expected to survive in the dense environment of the Galactic Center. The stable binaries may undergo Kozai-Lidov oscillations due to perturbations from the central supermassive black hole (Sgr A*), yet the General Relativistic precession can suppress the Kozai-Lidov oscillations and keep the stellar binaries from merging. However, it is challenging to resolve the binary sources and distinguish them from single stars. The close separation of the stable binaries allow higher eclipse probabilities. Here, we consider the massive star SO2 as an example and calculate the probability of detecting eclipses, assuming it is a binary. We find that the eclipse probability is ~ 30-50%, reaching higher values when the stellar binary is more eccentric or highly inclined relative to its orbit around Sgr A*.

The star S0-2, which orbits the supermassive black hole (SMBH) in our Galaxy with a period of 16 years, provides the strongest constraint on both the mass of the SMBH and the distance to the Galactic center. S0-2 will soon provide the first measurement of relativistic effects near a SMBH. We report the first limits on the binarity of S0-2 from radial velocity monitoring, which has implications for both understanding its origin and robustness as a probe of the central gravitational field. With 87 radial velocity measurements, which include 12 new observations presented, we have the data set to look for radial velocity variations from S0-2's orbital model. Using a Lomb-Scargle analysis and orbit fitting for potential binaries, we detect no radial velocity variation beyond S0-2's orbital motion and do not find any significant periodic signal. The lack of a binary companion does not currently distinguish between different formation scenarios for S0-2. The upper limit on the mass of a companion star ($M_{\text{comp}}$) still allowed by our results has a median upper limit of $M_{\text{comp}}$ $\sin i \leq$ 1.6 M$_{\odot}$ for periods between 1 and 150 days, the longest period to avoid tidal break up of the binary. We also investigate the impact of the remaining allowed binary system on the measurement of the relativistic redshift at S0-2's closest approach in 2018. While binary star systems are important to consider for this experiment, we find plausible binaries for S0-2 will not alter a 5$\sigma$ detection of the relativistic redshift.

We investigate the nature and evolution of large-scale structure within the SSA22 protocluster region at $z=3.09$ using cosmological simulations. A redshift histogram constructed from current spectroscopic observations of the SSA22 protocluster reveals two separate peaks at $z = 3.065$ (blue) and $z = 3.095$ (red). Based on these data, we report updated overdensity and mass calculations for the SSA22 protocluster. We find $\delta_{b,gal}=4.8 \pm 1.8$, $\delta_{r,gal}=9.5 \pm 2.0$ for the blue and red peaks, respectively, and $\delta_{t,gal}=7.6\pm 1.4$ for the entire region. These overdensities correspond to masses of $M_b = (0.76 \pm 0.17) \times 10^{15} h^{-1} M_{\odot}$, $M_r = (2.15 \pm 0.32) \times 10^{15} h^{-1} M_{\odot}$, and $M_t=(3.19 \pm 0.40) \times 10^{15} h^{-1} M_{\odot}$ for the red, blue, and total peaks, respectively. We use the Small MultiDark Planck (SMDPL) simulation to identify comparably massive $z\sim 3$ protoclusters, and uncover the underlying structure and ultimate fate of the SSA22 protocluster. For this analysis, we construct mock redshift histograms for each simulated $z\sim 3$ protocluster, quantitatively comparing them with the observed SSA22 data. We find that the observed double-peaked structure in the SSA22 redshift histogram corresponds not to a single coalescing cluster, but rather the proximity of a $\sim 10^{15}h^{-1} M_{\odot}$ protocluster and at least one $>10^{14} h^{-1} M_{\odot}$ cluster progenitor. Such associations in the SMDPL simulation are easily understood within the framework of hierarchical clustering of dark matter halos. We finally find that the opportunity to observe such a phenomenon is incredibly rare, with an occurrence rate of $7.4h^3 \mbox{ Gpc}^{-3}$.

White dwarfs have atmospheres that are expected to consist nearly entirely of hydrogen and helium, since heavier elements will sink out of sight on short timescales. However, observations have revealed atmospheric pollution by heavier elements in about a quarter to a half of all white dwarfs. While most of the pollution can be accounted for with asteroidal or dwarf planetary material, recent observations indicate that larger planetary bodies, as well as icy and volatile material from Kuiper belt analog objects, are also viable sources of pollution. The commonly accepted pollution mechanisms, namely scattering interactions between planetary bodies orbiting the white dwarfs, can hardly account for pollution by objects with large masses or long-period orbits. Here we report on a mechanism that naturally leads to the emergence of massive body and icy and volatile material pollution. This mechanism occurs in wide binary stellar systems, where the mass loss of the planets' host stars during post main sequence stellar evolution can trigger the Eccentric Kozai-Lidov mechanism. This mechanism leads to large eccentricity excitations, which can bring massive and long-period objects close enough to the white dwarfs to be accreted. We find that this mechanism readily explains and is consistent with observations.

We analyze the process of planetary scattering around M0-type stars. To do this, we carry out N-body simulations with three Jupiter-mass planets close to their instability limit together with an outer planetesimal disk. This paper focuses on the analysis of systems in which a single Jupiter-mass planet survives after the dynamical instability event. The small body reservoirs show different dynamical behaviors. In fact, our simulations produce particles on prograde and retrograde orbits, as well as particles whose orbital plane flips from prograde to retrograde and back again along their evolution. Such particles are called "Type-F particles". We find strong correlations between the inclination $i$ and the ascending node longitude $\Omega$ of such particles. First, $\Omega$ librates around 90$^{\circ}$ or/and 270$^{\circ}$. This property is very important since it represents a necessary and sufficient condition for the flipping of an orbit. Moreover, the libration periods of $i$ and $\Omega$ are equal and they are out to phase by a quarter period. We also remark that the larger the libration amplitude of $i$, the larger the libration amplitude of $\Omega$. Finally, we analyze the initial conditions of Type-F particles of all our simulations immediately after the dynamical instability event, when a single Jupiter-mass planet survives in the system. We carry out this study with the goal to determine the parameter space that lead to the flipping of an orbit. Our results suggest that the orbit of a test particle can flip for any value of its initial eccentricity, although we found only two Type-F particles with initial inclinations $i <$ 17$^{\circ}$. Moreover, our study indicates that the minimum value of the inclination of the Type-F particles in a given system decreases with an increase in the eccentricity of the giant planet.

We present new Adaptive Optics (AO) imaging and spectroscopic measurements of Galactic Center source G1 from W. M. Keck Observatory. Our goal is to understand its nature and relationship to G2, which is the first example of a spatially-resolved object interacting with the supermassive black hole (SMBH). Both objects have been monitored with AO for the past decade (2003 - 2014) and are comparatively close to the black hole ($a_{\rm{min}} \sim$200-300 AU) on very eccentric orbits ($e_{\rm{G1}}\sim$0.99; $e_{\rm{G2}}\sim$0.96). While G2 has been tracked before and during periapse passage ($T_{0} \sim$ 2014.2), G1 has been followed since soon after emerging from periapse ($T_{0} \sim$ 2001.3). Our observations of G1 double the previously reported observational time baseline, which improves its orbital parameter determinations. G1's orbital trajectory appears to be in the same plane as that of G2, but with a significantly different argument of periapse ($\Delta\omega$ = 21$\pm$4 degrees). This suggests that G1 is an independent object and not part of a gas stream containing G2 as has been proposed. Furthermore, we show for the first time that: (1) G1 is extended in the epochs closest to periapse along the direction of orbital motion and (2) G1 becomes significantly smaller over time, (450 AU in 2004 to less than 170 AU in 2009). Based on these observations, G1 appears to be the second example of an object tidally interacting with a SMBH. G1's existence 14 years after periapse, along with its compactness in epochs further from the time of periapse, suggest that this source is stellar in nature.

Many giant exoplanets are found near their Roche limit and in mildly eccentric orbits. In this study we examine the fate of such planets through Roche-lobe overflow as a function of the physical properties of the binary components, including the eccentricity and the asynchronicity of the rotating planet. We use a direct three-body integrator to compute the trajectories of the lost mass in the ballistic limit and investigate the possible outcomes. We find three different outcomes for the mass transferred through the Lagrangian point $L_{1}$: (i) self-accretion by the planet, (ii) direct impact on the stellar surface, (iii) disk formation around the star. We explore the parameter space of the three different regimes and find that at low eccentricities, $e\lesssim 0.2$, mass overflow leads to disk formation for most systems, while for higher eccentricities or retrograde orbits self-accretion is the only possible outcome. We conclude that the assumption often made in previous work that when a planet overflows its Roche lobe it is quickly disrupted and accreted by the star is not always valid.

Motion of short-period stars orbiting the supermassive black hole in our Galactic Center has been monitored for more than 20 years. These observations are currently offering a new way to test the gravitational theory in an unexplored regime: in a strong gravitational field, around a supermassive black hole. In this proceeding, we present three results: (i) a constraint on a hypothetical fifth force obtained by using 19 years of observations of the two best measured short-period stars S0-2 and S0-38 ; (ii) an upper limit on the secular advance of the argument of the periastron for the star S0-2 ; (iii) a sensitivity analysis showing that the relativistic redshift of S0-2 will be measured after its closest approach to the black hole in 2018.

In this Letter, we demonstrate that short-period stars orbiting around the supermassive black hole in our Galactic Center can successfully be used to probe the gravitational theory in a strong regime. We use 19 years of observations of the two best measured short-period stars orbiting our Galactic Center to constrain a hypothetical fifth force that arises in various scenarios motivated by the development of a unification theory or in some models of dark matter and dark energy. No deviation from General Relativity is reported and the fifth force strength is restricted to an upper 95% confidence limit of $\left|\alpha\right| < 0.016$ at a length scale of $\lambda=$ 150 astronomical units. We also derive a 95% confidence upper limit on a linear drift of the argument of periastron of the short-period star S0-2 of $\left|\dot \omega_\textrm{S0-2} \right|< 1.6 \times 10^{-3}$ rad/yr, which can be used to constrain various gravitational and astrophysical theories. This analysis provides the first fully self-consistent test of the gravitational theory using orbital dynamic in a strong gravitational regime, that of a supermassive black hole. A sensitivity analysis for future measurements is also presented.

The secular approximation of the hierarchical three body systems has been proven to be very useful in addressing many astrophysical systems, from planets, stars to black holes. In such a system two objects are on a tight orbit, and the tertiary is on a much wider orbit. Here we study the dynamics of a system by taking the tertiary mass to zero and solve the hierarchical three body system up to the octupole level of approximation. We find a rich dynamics that the outer orbit undergoes due to gravitational perturbations from the inner binary. The nominal result of the precession of the nodes is mostly limited for the lowest order of approximation, however, when the octupole-level of approximation is introduced the system becomes chaotic, as expected, and the tertiary oscillates below and above 90deg, similarly to the non-test particle flip behavior (e.g., Naoz 2016). We provide the Hamiltonian of the system and investigate the dynamics of the system from the quadrupole to the octupole level of approximations. We also analyze the chaotic and quasi-periodic orbital evolution by studying the surfaces of sections. Furthermore, including general relativity, we show case the long term evolution of individual debris disk particles under the influence of a far away interior eccentric planet. We show that this dynamics can naturally result in retrograde objects and a puffy disk after a long timescale evolution (few Gyr) for initially aligned configuration.

Models of debris disk morphology are often focused on the effects of a planet orbiting interior to or within the disk. Nonetheless, an exterior planetary-mass perturber can also excite eccentricities in a debris disk, via Laplace-Lagrange secular perturbations in the coplanar case or Kozai-Lidov perturbations for mutually inclined companions and disks. HD 106906 is an ideal example of such a system, as it harbors a confirmed exterior 11 M_Jup companion at a projected separation of 650 au outside a resolved, asymmetric disk. We use collisional and dynamical simulations to investigate the interactions between the disk and the companion, and to use the disk's observed morphology to place constraints on the companion's orbit. We conclude that the disk's observed morphology is consistent with perturbations from the observed exterior companion. Generalizing this result, we suggest that exterior perturbers, as well as interior planets, should be considered when investigating the cause of observed asymmetries in a debris disk.

Most, if not all, stars in the field are born in binary configurations or higher multiplicity systems. In dense stellar environment such as the Galactic Center (GC), many stars are expected to be in binary configurations as well. These binaries form hierarchical triple body systems, with the massive black hole (MBH) as the third, distant object. The stellar binaries are expected to undergo large amplitude eccentricity and inclination oscillations via the so-called "eccentric Kozai-Lidov" (EKL) mechanism. These eccentricity excitations, combined with post main sequence stellar evolution, can drive the inner stellar binaries to merge. We study the mergers of stellar binaries in the inner 0.1 pc of the GC caused by gravitational perturbations due to the MBH. We run a large set of Monte Carlo simulations that include the secular evolution of the orbits, general relativistic precession, tides, and post-main-sequence stellar evolution. We find that about 13 % of the initial binary population will have merged after a few million years and about 29 % after a few billion years. These expected merged systems represent a new class of objects at the GC and we speculate that they are connected to G2-like objects and the young stellar population.

The first indication of the presence of a circumstellar debris disk is usually the detection of excess infrared emission from the population of small dust grains orbiting the star. This dust is short-lived, requiring continual replenishment, and indicating that the disk must be excited by an unseen perturber. Previous theoretical studies have demonstrated that an eccentric planet orbiting interior to the disk will stir the larger bodies in the belt and produce dust via interparticle collisions. However, motivated by recent observations, we explore another possible mechanism for heating a debris disk: a stellar-mass perturber orbiting exterior to and inclined to the disk and exciting the disk particles' eccentricities and inclinations via the Kozai-Lidov mechanism. We explore the consequences of an exterior perturber on the evolution of a debris disk using secular analysis and collisional N-body simulations. We demonstrate that a Kozai-Lidov excited disk can generate a dust disk via collisions and we compare the results of the Kozai-Lidov excited disk with a simulated disk perturbed by an interior eccentric planet. Finally, we propose two observational tests of a dust disk that can distinguish whether the dust was produced by an exterior brown dwarf or stellar companion or an interior eccentric planet.

The hierarchical triple body approximation has useful applications to a variety of systems from planetary and stellar scales to supermassive black holes. In this approximation, the energy of each orbit is separately conserved and therefore the two semi-major axes are constants. On timescales much larger than the orbital periods, the orbits exchange angular momentum which leads to eccentricity and orientation (i.e., inclination) oscillations. The orbits' eccentricity can reach extreme values leading to a nearly radial motion, which can further evolve into short orbit periods and merging binaries. Furthermore, the orbits' mutual inclination may change dramatically from pure prograde to pure retrograde leading to misalignment and a wide range of inclinations. This dynamical behavior is coined as the "eccentric Kozai-Lidov" mechanism. The behavior of such a system is exciting, rich and chaotic in nature. Furthermore, these dynamics are accessible from a large part of the triple body parameter space and can be applied to diverse range of astrophysical settings and used to gain insights to many puzzles.

Using adiabatic high-resolution numerical simulations we quantify the effect of the streaming motion of baryons with respect to dark matter at the time of recombination on structure formation and evolution. Formally a second order effect, the baryonic stream velocity has proven to have significant impact on dark matter halo abundance, as well as on the gas content and morphology of small galaxy clusters. In this work, we study the impact of stream velocity on the formation and gas content of haloes with masses up to $10^9 M_{\odot}$, an order of magnitude larger than previous studies. We find that the non-zero stream velocity has a sizable impact on the number density of haloes with masses $\lesssim$ few $\times 10^7 M_{\odot}$ up to $z=10$, the final redshift of our simulations. Furthermore, the gas stream velocity induces a suppression of the gas fraction in haloes, which at z=10 is $\sim 10 \%$ for objects with $M\sim10^7M_{\odot}$, as well as a flattening of the gas density profiles in the inner regions of haloes. We further identify and study the formation, in the context of a non-zero stream velocity, of moderately long lived gas dominated structures at intermediate redshifts $10 < z < 20$, which Naoz and Narayan have recently proposed as potential progenitors of globular clusters.

The formation of Black Hole (BH) Low-Mass X-ray Binaries (LMXB) poses a theoretical challenge, as low-mass companions are not expected to survive the common-envelope scenario with the BH progenitor. Here we propose a formation mechanism that skips the common-envelope scenario and relies on triple-body dynamics. We study the evolution of hierarchical triples, following the secular dynamical evolution up to the octupole-level of approximation, including general relativity, tidal effects and post-main-sequence evolution, such as mass loss, changes to stellar radii and supernovae. During the dynamical evolution of the triple system, the "eccentric Kozai-Lidov" mechanism can cause large eccentricity excitations in the LMXB progenitor, resulting in three main BH-LMXB formation channels. Here we define BH-LMXB candidates as systems where the inner BH companion star crosses its Roche limit. In the "eccentric" channel (~ 81% of the LMXBs in our simulations), the donor star crosses its Roche limit during an extreme eccentricity excitation, while still on a wide orbit. Second, we find a "giant" LMXB channel (~ 11%), where a system undergoes only moderate eccentricity excitations, but the donor star fills its Roche lobe after evolving toward the giant branch. Third, we identify a "classical" channel (~8%), where tidal forces and magnetic braking shrink and circularize the orbit to short periods, triggering mass transfer. Finally, for the giant channel, we predict an eccentric ($\sim 0.3-0.6$), preferably inclined (~40, ~140 degreed) tertiary, typically on a wide enough orbit (~10^4AU), to potentially become unbound later in the triple evolution.

An enhanced rate of stellar tidal disruption events (TDEs) may be an important characteristic of supermassive black hole (SMBH) binaries at close separations. Here we study the evolution of the distribution of stars around a SMBH binary due to the eccentric Kozai-Lidov (EKL) mechanism, including octupole effects and apsidal precession caused by the stellar mass distribution and general relativity. We identify a region around one of the SMBHs in the binary where the EKL mechanism drives stars to high eccentricities, which ultimately causes the stars to either scatter off the second SMBH or get disrupted. For SMBH masses 10^7 Msun and 10^8 Msun, the TDE rate can reach 10^{-2} yr and deplete a region of the stellar cusp around the secondary SMBH in ~0.5 Myr. As a result, the final geometry of the stellar distribution between 0.01 and 0.1 pc around the secondary SMBH is a torus. These effects may be even more prominent in nuclear stellar clusters hosting a supermassive and an intermediate mass black hole.

We report new observations of the Galactic Center source G2 from the W. M. Keck Observatory. G2 is a dusty red object associated with gas that shows tidal interactions as it nears closest approach with the Galaxy's central black hole. Our observations, conducted as G2 passed through periapse, were designed to test the proposal that G2 is a 3 earth mass gas cloud. Such a cloud should be tidally disrupted during periapse passage. The data were obtained using the Keck II laser guide star adaptive optics system (LGSAO) and the facility near-infrared camera (NIRC2) through the K' [2.1 $\mu$m] and L' [3.8 $\mu$m] broadband filters. Several results emerge from these observations: 1) G2 has survived its closest approach to the black hole as a compact, unresolved source at L'; 2) G2's L' brightness measurements are consistent with those over the last decade; 3) G2's motion continues to be consistent with a Keplerian model. These results rule out G2 as a pure gas cloud and imply that G2 has a central star. This star has a luminosity of $\sim$30 $L_{\odot} $ and is surrounded by a large ($\sim$2.6 AU) optically thick dust shell. The differences between the L' and Br-$\gamma$ observations can be understood with a model in which L' and Br-$\gamma$ emission arises primarily from internal and external heating, respectively. We suggest that G2 is a binary star merger product and will ultimately appear similar to the B-stars that are tightly clustered around the black hole (the so-called S-star cluster).

We explore the effects of long term secular perturbations on the distribution of dark matter particles around Supermassive Black Hole (BH) binaries. We show that in the hierarchical (in separation) three-body problem, one of the BHs and a dark matter particle form an inner binary. Gravitational perturbations from the BH companion, on a much wider orbit, can cause the dark matter particle to reach extremely high eccentricities and even get accreted onto the BH, by what is known as the Eccentric Kozai-Lidov (EKL) mechanism. We show that this may produce a torus-like configuration for the dark matter distribution around the less massive member of the BH binary. We first consider an Intermediate BH (IMBH) in the vicinity of our Galactic Center, which may be a relic of a past minor merger. We show that if the IMBH is close enough (i.e., near the stellar disk) the EKL mechanism is very efficient in exciting the eccentricity of dark matter particles in near-polar configurations to extremely high values where they are accreted by the IMBH. We show that this mechanism is even more effective if the central BH grows in mass, where we have assumed adiabatic growth. Since near-polar configurations are disrupted, a torus-like shape is formed. We also show that this behavior is also likely to be relevant for Supermassive BH binaries. We suggest that if the BHs are spinning, the accreted dark matter particles may linger in the ergosphere and thereby may generate self-annihilations and produce an indirect signature of potential interest.

Many close stellar binaries are accompanied by a far-away star. The "eccentric Kozai-Lidov" (EKL) mechanism can cause dramatic inclination and eccentricity fluctuations, resulting in tidal tightening of inner binaries of triple stars. We run a large set of Monte-Carlo simulations including the secular evolution of the orbits, general relativistic precession and tides, and we determine the semimajor axis, eccentricity, inclination and spin-orbit angle distributions of the final configurations. We find that the efficiency of forming tight binaries (<~16 d) when taking the EKL mechanism into account is ~ 21%, and about 4% of all simulated systems ended up in a merger event. These merger events can lead to the formation of blue-stragglers. Furthermore, we find that the spin-orbit angle distribution of the inner binaries carries a signature of the initial setup of the system, thus observations can be used to disentangle close binaries' birth configuration. The resulting inner and outer final orbits' period distributions, and their estimated fraction, suggests secular dynamics may be a significant channel for the formation of close binaries in triples and even blue stragglers.

Kepler-56 is a multi-planet system containing two coplanar inner planets that are in orbits misaligned with respect to the spin axis of the host star, and an outer planet. Various mechanisms have been proposed to explain the broad distribution of spin-orbit angles among exoplanets, and these theories fall under two broad categories. The first is based on dynamical interactions in a multi-body system, while the other assumes that disk migration is the driving mechanism in planetary configuration and that the star (or disk) is titled with respect to the planetary plane. Here we show that the large observed obliquity of Kepler-56 system is consistent with a dynamical origin. In addition, we use observations by Huber et al. (2013) to derive the obliquity's probability distribution function, thus improving the constrained lower limit. The outer planet may be the cause of the inner planets' large obliquities, and we give the probability distribution function of its inclination, which depends on the initial orbital configuration of the planetary system. We show that even in the presence of precise measurement of the true obliquity, one cannot distinguish the initial configurations. Finally we consider the fate of the system as the star continues to evolve beyond the main sequence, and we find that the obliquity of the system will not undergo major variations as the star climbs the red giant branch. We follow the evolution of the system and find that the innermost planet will be engulfed in ~129 Myr. Furthermore we put an upper limit of ~155 Myr for the engulfment of the second planet. This corresponds to ~ 3% of the current age of the star.

The formation of purely baryonic globular clusters with no gravitationally bound dark matter is still a theoretical challenge. We show that these objects might form naturally whenever there is a relative stream velocity between baryons and dark matter. The stream velocity causes a phase shift between linear modes of baryonic and dark matter perturbations, which translates to a spatial offset between the two components when they collapse. For a 2sigma (3sigma) density fluctuation, baryonic clumps with masses in the range 1e5 - 2.5e6 Msun (1e5 - 4e6 Msun) collapse outside the virial radii of their counterpart dark matter halos. These objects could survive as long-lived dark matter-free objects and might conceivably become globular clusters. In addition, their dark matter counterparts, which were deprived of gas, might become dark satellite galaxies.

The Kozai-Lidov mechanism can be applied to a vast variety of astrophysical systems involving hierarchical three-body systems. Here, we study the Kozai-Lidov mechanism systematically in the test particle limit at the octupole level of approximation. We investigate the chaotic and quasiperiodic orbital evolution by studying surfaces of section and the Lyapunov exponents. We find that the resonances introduced by the octupole level of approximation cause orbits to flip from prograde to retrograde and back as well as cause significant eccentricity excitation, and the chaotic behaviors occur when the mutual inclination between the inner and the outer binary is high. We characterize the parameter space that allows large amplitude oscillations in eccentricity and inclination.

The dynamical evolution of a hierarchical three body system is well characterized by the eccentric Kozai-Lidov mechanism, where the inner orbit can undergo large eccentricity and inclination oscillations. It was shown before that starting with a circular inner orbit, large mutual inclination (40-140 degree) can produce long timescale modulations that drives the eccentricity to extremely large value and can flip the orbit. Here, we demonstrate that starting with an almost coplanar configuration, for eccentric inner and outer orbits, the eccentricity of the inner orbit can still be excited to high values, and the orbit can flip by ~180 degree, rolling over its major axis. The ~180 degree flip criterion and the flip timescale are described by simple analytic expressions that depend on the initial orbital parameters. With tidal dissipation, this mechanism can produce counter-orbiting exo-planetary systems. In addition, we also show that this mechanism has the potential to enhance the tidal disruption or collision rates for different systems.

We study the secular, hierarchical three-body problem to first-order in a post-Newtonian expansion of General Relativity. We expand the first-order post-Newtonian Hamiltonian to leading-order in the ratio of the semi-major axis of the two orbits. In addition to the well-known terms that correspond to the GR precession of the inner and outer orbits, we find a new secular post-Newtonian interaction term that can affect the long-term evolution of the triple. We explore the parameter space for highly inclined and eccentric systems, where the Kozai-Lidov mechanism can produce large-amplitude oscillations in the eccentricities. The standard lore, i.e., that General Relativity effects suppress eccentricity, is only consistent with the parts of phase space where the General Relativity timescales are several orders of magnitude shorter than the secular Newtonian one. In other parts of phase space, however, post-Newtonian corrections combined with the three body ones, can excite eccentricities. In particular, for systems where the General Relativity timescale is comparable to the secular Newtonian timescales, the three-body interactions give rise to a resonant-like eccentricity excitation. Furthermore, for triples with a comparable-mass inner binary, where the eccentric Kozai-Lidov mechanism is suppressed, post-Newtonian corrections can further increase the eccentricity and lead to orbital flips even when the timescale of the former is much longer than the timescale of the secular Kozai-Lidov quadrupole perturbations.

Observations of exoplanets over the last two decades have revealed a new class of Jupiter-size planets with orbital periods of a few days, the so-called "hot Jupiters". Recent measurements using the Rossiter-McLaughlin effect have shown that many (~ 50%) of these planets are misaligned; furthermore, some (~ 15%) are even retrograde with respect to the stellar spin axis. Motivated by these observations, we explore the possibility of forming retrograde orbits in hierarchical triple configurations consisting of a star-planet inner pair with another giant planet, or brown dwarf, in a much wider orbit. Recently Naoz et al. (2011) showed that in such a system, the inner planet's orbit can flip back and forth from prograde to retrograde, and can also reach extremely high eccentricities. Here we map a significant part of the parameter space of dynamical outcomes for these systems. We derive strong constraints on the orbital configurations for the outer perturber that could lead to the formation of hot Jupiters with misaligned or retrograde orbits. We focus only on the secular evolution, neglecting other dynamical effects such as mean-motion resonances, as well as all dissipative forces. For example, with an inner Jupiter-like planet initially on a nearly circular orbit at 5 AU, we show that a misaligned hot Jupiter is likely to be formed in the presence of a more massive planetary companion (> 2 MJ) within 140 AU of the inner system, with mutual inclination 50 degrees and eccentricity above 0.25. This is in striking contrast to the test-particle approximation, where an almost perpendicular configuration can still cause large eccentricity excitations, but flips of an inner Jupiter-like planet are much less likely to occur. The constraints we derive can be used to guide future observations, and, in particular, searches for more distant companions in systems containing a hot Jupiter.

X-ray photons, because of their long mean-free paths, can easily escape the galactic environments where they are produced, and interact at long distances with the inter-galactic medium, potentially having a significant contribution to the heating and reionization of the early Universe. The two most important sources of X-ray photons in the Universe are active galactic nuclei (AGN) and X-ray binaries (XRBs). In this Letter we use results from detailed, large scale population synthesis simulations to study the energy feedback of XRBs, from the first galaxies (z~ 20) until today. We estimate that X-ray emission from XRBs dominates over AGN at z>6-8. The shape of the spectral energy distribution of the emission from XRBs shows little change with redshift, in contrast to its normalization which evolves by ~4 orders of magnitude, primarily due to the evolution of the cosmic star-formation rate. However, the metallicity and the mean stellar age of a given XRB population affect significantly its X-ray output. Specifically, the X-ray luminosity from high-mass XRBs per unit of star-formation rate varies an order of magnitude going from solar metallicity to less than 10% solar, and the X-ray luminosity from low-mass XRBs per unit of stellar mass peaks at an age of ~300 Myr and then decreases gradually at later times, showing little variation for mean stellar ages > 3 Gyr. Finally, we provide analytical and tabulated prescriptions for the energy output of XRBs, that can be directly incorporated in cosmological simulations.

Magnetic fields appear to be present in all galaxies and galaxy clusters. Recent measurements indicate that a weak magnetic field may be present even in the smooth low density intergalactic medium. One explanation for these observations is that a seed magnetic field was generated by some unknown mechanism early in the life of the Universe, and was later amplified by various dynamos in nonlinear objects like galaxies and clusters. We show that a primordial magnetic field is expected to be generated in the early Universe on purely linear scales through vorticity induced by scale-dependent temperature fluctuations or equivalently, a spatially varying speed of sound of the gas. Residual free electrons left over after recombination tap into this vorticity to generate magnetic field via the Biermann battery process. Although the battery operates even in the absence of any relative velocity between dark matter and gas at the time of recombination, the presence of such a relative velocity modifies the predicted spatial power spectrum of the magnetic field. At redshifts of order a few tens, we estimate a root mean square field strength of order 1e-25 -1e-24 G on comoving scales~10 kpc. This field, which is generated purely from linear perturbations, is expected to be amplified significantly after reionization, and to be further boosted by dynamo processes during nonlinear structure formation.

The secular approximation for the evolution of hierarchical triple configurations has proven to be very useful in many astrophysical contexts, from planetary to triple-star systems. In this approximation the orbits may change shape and orientation, on time scales longer than the orbital time scales, but the semi major axes are constant. For example, for highly inclined triple systems, the Kozai-Lidov mechanism can produce large-amplitude oscillations of the eccentricities and inclinations. Here we revisit the secular dynamics of hierarchical triple systems. We derive the secular evolution equations to octupole order in Hamiltonian perturbation theory. Our derivation corrects an error in some previous treatments of the problem that implicitly assumed a conservation of the z-component of the angular momentum of the inner orbit (i.e., parallel to the total angular momentum of the system). Already to quadrupole order, our results show new behaviors including the possibility for a system to oscillate from prograde to retrograde orbits. At the octupole order, for an eccentric outer orbit, the inner orbit can reach extremely high eccentricities and undergo chaotic flips in its orientation. We discuss applications to a variety of astrophysical systems, from stellar triples to merging compact binaries and planetary systems. Our results agree with those of previous studies done to quadrupole order only in the limit in which one of the inner two bodies is a massless test particle and the outer orbit is circular;our results agree with previous studies at octupole order for the eccentricity evolution, but not for the inclination evolution.

Understanding the gas content of high redshift halos is crucial for studying the formation of the first generation of galaxies and reionization. Recently, Tseliakhovich & Hirata showed that the relative "stream" velocity between the dark matter and baryons at the time of recombination - formally a second order effect, but an unusually large one - can influence the later structure formation history of the Universe. We quantify the effect of the stream velocity on the so-called "characteristic mass" - the minimum mass of a dark matter halo capable of retaining most of its baryons throughout its formation epoch - using three different high-resolution sets of cosmological simulations (with separate transfer functions for baryons and dark matter) that vary in box size, particle number, and the value of the relative velocity between the dark matter and baryons. In order to understand this effect theoretically, we generalize the linear theory filtering mass to properly account for the difference between the dark matter and baryonic density fluctuation evolution induced by the stream velocity. We show that the new filtering mass provides an accurate estimate for the characteristic mass, while other theoretical ansatzes for the characteristic mass are substantially less precise.

We study the production of Hot Jupiters (HJs) in stellar binaries. We show that the "eccentric Kozai-Lidov" (EKL) mechanism can play a key role in the dynamical evolution of a star-planet-star triple system. We run a large set of Monte Carlo simulations including the secular evolution of the orbits, general relativistic precession, and tides, and we determine the semi-major axis, eccentricity, inclination and spin-orbit angle distributions of the HJs that are produced. We explore the effect of different tidal friction parameters on the results. We find that the efficiency of forming HJs when taking the EKL mechanism into account is higher then previously estimated. Accounting for the frequency of stellar binaries, we find that this production mechanism can account for about 30% of the observed HJ population. Current observations of spin-orbit angles are consistent with this mechanism producing \sim 30% of all HJs, and up to 100% of the misaligned systems. Based on the properties of binaries without a HJ in our simulations, we predict the existence of many Jupiter-like planets with moderately eccentric and inclined orbits and semi-major axes of several AU.

It has been recently shown that the relative velocity between the dark matter and the baryons (vbc) at the time of recombination can affect the structure formation in the early universe (Tseliakhovich & Hirata 2010). We statistically quantify this effect using large cosmological simulations. We use three different high resolution sets of simulations (with separate transfer functions for baryons and dark matter) that vary in box size, particle number, and the value of the relative velocity between dark matter and baryons. We show that the total number density of halos is suppressed by ~ 20% at z = 25 for vbc = 1{\sigma}(vbc), where {\sigma}(vbc) is the variance of the relative velocity, while for vbc = 3.4{\sigma}(vbc) the relative suppression at the same redshift reaches 50%, remaining at or above the 30% level all the way to z = 11. We also find high abundance of "empty halos", i.e., halos that have gas fraction below half of the cosmic mean baryonic fraction fb. Specifically we find that for vbc = 1{\sigma}(vbc) all halos below 10^5M\odot are empty at z \geq 19. The high abundance of empty halos results in significant delay in the formation of gas rich mini-halos and the first galaxies.

We study the dynamical evolution of a test particle that orbits a star in the presence of an exterior massive planet, considering octupole-order secular interactions. In the standard Kozai mechanism (SKM), the planet's orbit is circular, and so the particle conserves vertical angular momentum. As a result, the particle's orbit oscillates periodically, exchanging eccentricity for inclination. However, when the planet's orbit is eccentric, the particle's vertical angular momentum varies and its Kozai oscillations are modulated on longer timescales---we call this the eccentric Kozai mechanism (EKM). The EKM can lead to behavior that is dramatically different from the SKM. In particular, the particle's orbit can flip from prograde to retrograde and back again, and it can reach arbitrarily high eccentricities given enough time. We map out the conditions under which this dramatic behavior (flipping and extreme eccentricities) occurs, and show that when the planet's eccentricity is sufficiently high, it occurs quite generically. For example, when the planet's eccentricity exceeds a few percent of the ratio of semi-major axes (outer to inner), around half of randomly oriented test particle orbits will flip and reach extreme eccentricities. The SKM has often been invoked for bringing pairs of astronomical bodies (star-star, planet-star, compact-object pairs) close together. Including the effect of the EKM will enhance the rate at which such matchmaking occurs.

About 25 per cent of `hot Jupiters' (extrasolar Jovian-mass planets with close-in orbits) are actually orbiting counter to the spin direction of the star. Perturbations from a distant binary star companion can produce high inclinations, but cannot explain orbits that are retrograde with respect to the total angular momentum of the system. Such orbits in a stellar context can be produced through secular (that is, long term) perturbations in hierarchical triple-star systems. Here we report a similar analysis of planetary bodies, including both octupole-order effects and tidal friction, and find that we can produce hot Jupiters in orbits that are retrograde with respect to the total angular momentum. With distant stellar mass perturbers, such an outcome is not possible. With planetary perturbers, the inner orbit's angular momentum component parallel to the total angular momentum need not be constant. In fact, as we show here, it can even change sign, leading to a retrograde orbit. A brief excursion to very high eccentricity during the chaotic evolution of the inner orbit allows planet-star tidal interactions to rapidly circularize that orbit, decoupling the planets and forming a retrograde hot Jupiter.

We run very large cosmological N-body hydrodynamical simulations in order to study statistically the baryon fractions in early dark matter halos. We critically examine how differences in the initial conditions affect the gas fraction in the redshift range z = 11--21. We test three different linear power spectra for the initial conditions: (1) A complete heating model, which is our fiducial model; this model follows the evolution of overdensities correctly, according to Naoz & Barkana (2005), in particular including the spatial variation of the speed of sound of the gas due to Compton heating from the CMB. (2) An equal-{\delta} model, which assumes that the initial baryon fluctuations are equal to those of the dark matter, while conserving sigma8 of the total matter. (3) A mean cs model, which assumes a uniform speed of sound of the gas. The latter two models are often used in the literature. We calculate the baryon fractions for a large sample of halos in our simulations. Our fiducial model implies that before reionization and significant stellar heating took place, the minimum mass needed for a minihalo to keep most of its baryons throughout its formation was ~ 3 * 10^4 Msun. However, the alternative models yield a wrong (higher by about 50%) minimum mass, since the system retains a memory of the initial conditions. We also demonstrate this using the "filtering mass" from linear theory, which accurately describes the evolution of the baryon fraction throughout the simulated redshift range.

Many binary minor planets (BMPs; both binary asteroids and binary Trans-Neptunians objects; TNOs) are known to exist in the Solar system. The currently observed orbital and physical properties of BMPs hold essential information and clues about their origin, their evolution and the conditions underwhich they evolved. Here we study the orbital properties of BMPs with currently known orbital solutions (the BMPs mutual orbits and not the BMP orbits around the sun). We find that BMPs are typically highly inclined relative to their orbit around the sun, with a distribution consistent with an isotropic distribution. BMPs not affected by tidal forces are found to have high eccentricities with non-thermal eccentricity distribution peaking at intermediate eccentricities (typically 0.4-0.6). The high inclinations and eccentricities of the BMPs suggest that BMPs evolved in a dense collisional environment, in which gravitational encounters in addition to tidal and secular Kozai affects played an important role in their orbital evolution.

We study the gas content of halos in the early universe using high resolution hydrodynamical simulations. We extract from the simulations and also predict based on linear theory the halo mass for which the enclosed baryon fraction equals half of the mean cosmic fraction. We find a rough agreement between the simulations and the predictions, which suggests that during the high-redshift era before stellar heating, the minimum mass needed for a minihalo to keep most of its baryons throughout its formation was $\sim 3 \times 10^4$ M$_\odot$. We also carry out a detailed resolution analysis and show that in order to determine a halo's gas fraction even to 20% accuracy the halo must be resolved into at least 500 dark matter particles.

In recent years many binary minor planets (BMPs) have been discovered in the Solar system. Many models have been suggested for their formation, but these encounter difficulties explaining their observed characteristics. Here we show that secular perturbations by the Sun (Kozai mechanism) fundamentally change the evolution and the initial distribution of BMPs predicted by such models and lead to unique observational signatures. The Kozai mechanism can lead to a large periodic oscillations in the eccentricity and inclination of highly inclined BMP orbits, where we predict such effects to be observable with current accuracy within a few years (e.g. for the binary asteroid Huenna). In addition, the combined effects of the Kozai mechanism and tidal friction (KCTF) drives BMPs into short period circular orbits. We predict a specific inclination dependent distribution of the separation and eccentricity of BMPs, due to these effects, including a zone of avoidance at the highest inclinations. Specifically the Kozai evolution could explain the recently observed peculiar orbit of the Kuiper belt binary 2001 QW322 . Additionally, the KCTF process could lead to BMPs coalescence and serve as an important route for the formation of irregular shaped single minor planets with large axial tilts.

New observations over the next few years of the emission of distant objects will help unfold the chapter in cosmic history around the era of the first galaxies. These observations will use the neutral hydrogen emission or absorption at a wavelength of 21-cm as a detector of the hydrogen abundance. We predict the signature on the 21-cm signal of the early generations of galaxies. We calculate the 21-cm power spectrum including two physical effects that were neglected in previous calculations. The first is the redistribution of the UV photons from the first galaxies due to their scattering off of the neutral hydrogen, which results in an enhancement of the 21-cm signal. The second is the presence of an ionized hydrogen bubble near each source, which produces a cutoff at observable scales. We show that the resulting clear signature in the 21-cm power spectrum can be used to detect and study the population of galaxies that formed just 200 million years after the Big Bang.

Large telescopes have allowed astronomers to observe galaxies that formed as early as 850 million years after the Big Bang. We predict when the first star that astronomers can observe formed in the universe, accounting for the first time for the size of the universe and for three essential ingredients: the light travel time from distant galaxies, Poisson and density fluctuations on all scales, and the effect of very early cosmic history on galaxy formation. We find that the first observable star is most likely to have formed 30 million years after the Big Bang (at redshift 65), much earlier than previously expected. Also, the first galaxy as massive as our own Milky Way likely formed when the universe was only 400 Myr old (at redshift 11). We also show that significant modifications are required in current methods of numerically simulating the formation of galaxies at redshift 20 and above.

We predict the redshift of the first observable (i.e., in our past light cone) Gamma Ray Burst (GRB) and calculate the GRB-rate redshift distribution of the Population III stars at very early times (z=20-60). Using the last 2 years of data from Swift we place an upper limit on the efficiency (\eta_{GRB}) of GRB production per solar mass from the first generation of stars. We find that the first observable GRB is most likely to have formed at redshift 60. The observed rate of extremely high redshift GRBs (XRGs) is a subset of a group of 15 long GRBs per year, with no associated redshift and no optical afterglow counterparts, detected by Swift. Taking this maximal rate we get that \eta_{GRB}<1.1~10^{-4} GRBs per solar mass in stars. A more realistic evaluation, e.g., taking a subgroup of 5% of the total sample of Swift gives an upper limit of \eta_{GRB}<3.2~10^{-5} GRBs per solar mass.

We investigate the suppression of the baryon density fluctuations compared to the dark matter in the linear regime. Previous calculations predict that the suppression occurs up to a characteristic mass scale of ~ 1,000,000 solar masses, which suggests that pressure has a central role in determining the properties of the first luminous objects at early times. We show that the expected characteristic mass scale is in fact substantially lower (by a factor of ~ 3-10, depending on redshift), and thus the effect of baryonic pressure on the formation of galaxies up to reionization is only moderate. This result is due to the influence on perturbation growth of the high pressure that prevailed in the period from cosmic recombination to z ~ 200, when the gas began to cool adiabatically and the pressure then dropped. At z ~ 10 the suppression of the baryon fluctuations is still sensitive to the history of pressure in this high-redshift era. We calculate the fraction of the cosmic gas that is in minihalos and find that it is substantially higher than would be expected with the previously-estimated characteristic mass. Expanding our investigation to the non-linear regime, we calculate in detail the spherical collapse of high-redshift objects in a Lambda-CDM universe. We include the gravitational contributions of the baryons and radiation and the memory of their kinematic coupling before recombination. We use our results to predict a more accurate halo mass function as a function of redshift.

We calculate the evolution of linear density and temperature perturbations in a universe with dark matter, baryons, and radiation, from cosmic recombination until the epoch of the first galaxies. In addition to gravity, the perturbations are effected by electron scattering with the radiation, by radiation pressure, and by gas pressure. We include the effect of spatial fluctuations in the baryonic sound speed and show that they induce a >10% change in the baryonic density power spectrum on small scales, and a larger change on all scales in the power spectrum of gas temperature fluctuations. A precise calculation of the growth of linear perturbations is essential since they provide the initial conditions for the formation of galaxies and they can also be probed directly via cosmological 21cm fluctuations. We also show that in general the thermal history of the cosmic gas can be measured from 21cm fluctuations using a small-scale anisotropic cutoff due to the thermal width of the 21cm line.

The Milky Way spiral arm dynamics is studied using the birth place of open clusters. We separately study the nearby spiral arms, and find evidence for multiple spiral sets. In particular, the Sagittarius-Carina arm appears to be a superposition of two sets. The first has a pattern speed of Omega_P,Carina,1=16.5^{+1.2}_{-1.4sys}+/-1.1_{stat} (km/s)/kpc, while the second has Omega_P,Carina,2=29.8^{+0.6}_{-1.3sys}+/-1.3_{stat} (km/s)/kpc. The Perseus arm located at larger galactic radii exhibits only one clear set with Omega_P,Perseus=20.0^{+1.7}_{-1.2sys}+/-1.6_{stat} (km/s)/kpc. The Perseus and "slower" Carina arms are most probably part of an m=4 set. The Orion "armlet" appears to have a pattern speed of Omega_P,Orion=28.9 ^{+1.3}_{-1.2sys}+/-0.8_{stat} (km/s)/kpc. In other words, both the "faster" Carina arm and the Orion arm are nearly corotating with the solar system.