Source author record

Smadar Naoz

Smadar Naoz appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

40works
7topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

40 published item(s)

preprint2026arXiv

Possible Dynamical Pathways to the Misalignment of the VHS 1256-1257 System

Circumbinary planets (CBPs) provide a unique window into planet formation and dynamical evolution in complex gravitational environments. Their orbits are shaped not only by the protoplanetary disk but also by the perturbations from two stellar hosts, making them sensitive probes of both early- and late-stage dynamical processes. In this work, we investigate the unusual architecture of the VHS J125601.92-125723.9 system, where a retrograde, nearly polar tertiary orbits an extremely low-mass substellar binary in a hierarchical triple configuration. We find that triple body dynamics can naturally reproduce the observed high eccentricity of the inner binary and the tertiary's near-polar obliquity. However, this configuration alone cannot account for the observed mutual inclination, which is both near-polar and retrograde. This tension suggests two possible formation pathways: either the planet formed in an aligned, protoplanetary disk-like configuration and was later tilted by an additional, undetected fourth companion (below current Gaia limits), or the system formed close to its current state. Stellar flybys, in contrast, are unlikely due to their long timescales. Our results highlight both the explanatory power and the limitations of triple dynamics, and the potential role of hidden companions in shaping extreme planetary architectures.

preprint2022arXiv

$\mathrm{H_2}$ cooling and gravitational collapse of supersonically induced gas objects

We study the formation and gravitational collapse of supersonically induced gas objects (SIGOs) in the early universe. We run cosmological hydrodynamics simulations of SIGOs, including relative streaming motions between baryons and dark matter. Our simulations also follow nonequilibrium chemistry and molecular hydrogen cooling in primordial gas clouds. A number of SIGOs are formed in the run with fast-streaming motions of 2 times the rms of the cosmological velocity fluctuations. We identify a particular gas cloud that condensates by H$_2$ cooling without being hosted by a dark matter halo. The SIGO remains outside the virial radius of its closest halo, and it becomes Jeans unstable when the central gas-particle density reaches $\sim 100~{\rm cm}^{-3}$ with a temperature of $\sim$ 200 K. The corresponding Jeans mass is $\sim 10^5 M_{\odot}$, and thus the formation of primordial stars or a star cluster is expected in the SIGO.

preprint2022arXiv

Binary Natal Kicks in the Galactic Center: X-ray Binaries, Hypervelocity Stars, and Gravitational Waves

Theoretical and observational studies suggest that stellar binaries exist in large numbers in galactic nuclei like our own Galactic Center. Neutron stars (NSs), and debatedly, black holes (BHs) and white dwarfs (WDs), receive natal kicks at birth. In this work we study the effect of two successive natal kicks on a population of stellar binaries orbiting the massive black hole (MBH) in our Galactic Center. These natal kicks can significantly alter the binary orbit in a variety of ways, and also the orbit of the binary around the MBH. We found a variety of dynamical outcomes resulting from these kicks, including a steeper cusp of single NSs relative to the initial binary distribution. Furthermore, hypervelocity star and binary candidates, including hypervelocity X-ray binaries, are a common outcome of natal kicks. In addition, we show that the population of X-ray binaries in the Galactic Center can be used as a diagnostic for the BH natal kick distribution. Finally, we estimate the rate of gravitational wave (GW) events triggered by natal kicks, including binary mergers and EMRIs.

preprint2022arXiv

Hiding Planets Near and Far: The Parameter Space of Hidden Companions for Known Planetary Systems

Recent ground and space-based observations show that stars with multiple planets are common in the galaxy. Most of these observational methods are biased toward detecting large planets near to their host stars. Because of these observational biases, these systems can hide small, close-in planets or far-orbiting (big or small) companions. These planets can still exert dynamical influence on known planets and have such influence exerted upon them in turn. In certain configurations, this influence can destabilize the system; in others, the star's gravitational influence can instead further stabilize the system. For example, in systems with planets close to the host star, effects arising from general relativity can help to stabilize the configuration. We derive criteria for hidden planets orbiting both beyond and within known planets that quantify how strongly general relativistic effects can stabilize systems that would otherwise be unstable. As a proof-of-concept, we investigate the several planets in a system based on Kepler 56, and show that the outermost planet will not disrupt the system even at high eccentricities, and show that an Earth-radius planet could be stable within this system if it orbits below 0.08 au. Furthermore, we provide specific predictions to known observed systems by constraining the parameter space of possible hidden planets.

preprint2022arXiv

Kepler-1656b's Extreme Eccentricity: Signature of a Gentle Giant

Highly eccentric orbits are one of the major surprises of exoplanets relative to the Solar System and indicate rich and tumultuous dynamical histories. One system of particular interest is Kepler-1656, which hosts a sub-Jovian planet with an eccentricity of 0.8. Sufficiently eccentric orbits will shrink in semi-major axis due to tidal dissipation of orbital energy during periastron passage. Here our goal was to assess whether Kepler-1656b is currently undergoing such high-eccentricity migration, and to further understand the system's origins and architecture. We confirm a second planet in the system with $M_{\rm c}= 0.40 \pm 0.09M_{\rm jup}$ and P$_{\rm c}= 1919\pm 27\,$days. We simulated the dynamical evolution of planet b in the presence of planet c and find a variety of possible outcomes for the system, such as tidal migration and engulfment. The system is consistent with an in situ dynamical origin of planet b followed by subsequent Eccentric Kozai Lidov (EKL) perturbations that excite Kepler-1656b's eccentricity gently, i.e. without initiating tidal migration. Thus, despite its high eccentricity, we find no evidence that planet b is or has migrated through the high-eccentricity channel. Finally, we predict the outer orbit to be mutually inclined in a nearly perpendicular configuration with respect to the inner planet orbit based on the outcomes of our simulations, and make observable predictions for the inner planet's spin-orbit angle. Our methodology can be applied to other eccentric or tidally locked planets to constrain their origins, orbital configurations and properties of a potential companion.

preprint2022arXiv

The Combined Effects of Two-Body Relaxation Processes and the Eccentric Kozai-Lidov Mechanism on the EMRI Rate

Gravitational wave (GW) emissions from extreme-mass-ratio inspirals (EMRIs) are promising sources for low-frequency GW-detectors. They result from a compact object, such as a stellar-mass black-hole (BH), captured by a supermassive black hole (SMBH). Several physical processes have been proposed to form EMRIs. In particular, weak two-body interactions over a long time scale (i.e., relaxation processes) have been proposed as a likely mechanism to drive the BH orbit to high eccentricity. Consequently, it is captured by the SMBH and becomes an EMRI. Here we demonstrate that EMRIs are naturally formed in SMBH binaries. Gravitational perturbations from an SMBH companion, known as the eccentric Kozai-Lidov (EKL) mechanism, combined with relaxation processes, yield a significantly more enhanced rate than any of these processes operating alone. Since EKL is sensitive to the orbital configuration, two-body relaxation can alter the orbital parameters, rendering the system in a more EKL-favorable regime. As SMBH binaries are expected to be prevalent in the Universe, this process predicts a substantially high EMRI rate.

preprint2022arXiv

The Formation of Intermediate Mass Black Holes in Galactic Nuclei

Most stellar evolution models predict that black holes (BHs) should not exist above approximately $50-70$ M$_\odot$, the lower limit of the pair-instability mass gap. However, recent LIGO/Virgo detections indicate the existence of BHs with masses at and above this threshold. We suggest that massive BHs, including intermediate mass black holes (IMBHs), can form in galactic nuclei through collisions between stellar-mass black holes and the surrounding main-sequence stars. Considering dynamical processes such as collisions, mass segregation, and relaxation, we find that this channel can be quite efficient, forming IMBHs as massive as $10^4$ M$_\odot$. This upper limit assumes that (1) the BHs accrete a substantial fraction of the stellar mass captured during each collision and (2) that the rate at which new stars are introduced into the region near the SMBH is high enough to offset depletion by stellar disruptions and star-star collisions. We discuss deviations from these key assumptions in the text. Our results suggest that BHs in the pair-instability mass gap and IMBHs may be ubiquitous in galactic centers. This formation channel has implications for observations. Collisions between stars and BHs can produce electromagnetic signatures, for example, from x-ray binaries and tidal disruption events. Additionally, formed through this channel, both black holes in the mass gap and IMBHs can merge with the supermassive black hole at the center of a galactic nucleus through gravitational waves. These gravitational wave events are extreme and intermediate mass ratio inspirals (EMRIs and IMRIs, respectively).

preprint2020arXiv

A population of dust-enshrouded objects orbiting the Galactic black hole

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.

preprint2020arXiv

Demographics of triple systems in dense star clusters

Depending on the stellar type, more than $\sim 50$\% and $\sim 15\%$ of stars in the field have at least one and two stellar companions, respectively. Hierarchical systems can be assembled dynamically in dense star clusters, as a result of few-body encounters among stars and/or compact remnants in the cluster core. In this paper, we present the demographics of stellar and compact-object triples formed via binary--binary encounters in the \texttt{CMC Cluster Catalog}, a suite of cluster simulations with present-day properties representative of the globular clusters (GCs) observed in the Milky Way. We show how the initial properties of the host cluster set the typical orbital parameters and formation times of the formed triples. We find that a cluster typically assembles hundreds of triples with at least one black hole (BH) in the inner binary, while only clusters with sufficiently small virial radii are efficient in producing triples with no BHs, as a result of the BH-burning process. We show that a typical GC is expected to host tens of triples with at least one luminous component at present day. We discuss how the Lidov-Kozai mechanism can drive the inner binary of the formed triples to high eccentricities, whenever it takes place before the triple is dynamically reprocessed by encountering another cluster member. Some of these systems can reach sufficiently large eccentricities to form a variety of transients and sources, such as blue stragglers, X-ray binaries, Type Ia Supernovae, Thorne-Zytkow objects, and LIGO/Virgo sources.

preprint2020arXiv

The Supersonic Project: To cool or not to cool Supersonically Induced Gas Objects (SIGOs)?

Supersonically Induced Gas Objects (SIGOs) primarily form in the early Universe, outside of dark matter halos due to the presence of a relative stream velocity between baryons and dark matter. These structures may be the progenitors of globular clusters. Since SIGOs are made out of pristine gas, we investigate the effect of atomic cooling on their properties. We run a suite of simulations by using the moving-mesh code {\sc arepo}, with and without baryon-dark matter relative velocity and with and without the effects of atomic cooling. We show that SIGO's density, temperature, and prolateness are determined by gravitational interactions rather than cooling. The cold gas fraction in SIGOs is much higher than that of dark matter halos. Specifically, we show that the SIGO's characteristic low temperature and extreme high gas density forges a nurturing ground for the earliest star formation sites.

preprint2019arXiv

Eating Planets for Lunch and Dinner: Signatures of Planet Consumption by Evolving Stars

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.

preprint2016arXiv

Circumstellar Debris Disks: Diagnosing the Unseen Perturber

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.

preprint2016arXiv

Merging Binaries in the Galactic Center: The eccentric Kozai-Lidov mechanism with stellar evolution

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.

preprint2016arXiv

The Eccentric Kozai-Lidov Effect and Its Applications

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.

preprint2015arXiv

Formation of Black Hole Low-Mass X-ray Binaries in Hierarchical Triple Systems

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.

preprint2015arXiv

Gas rich and gas poor structures through the stream velocity effect

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.

preprint2015arXiv

Implications of the Eccentric Kozai-Lidov Mechanism for Stars Surrounding Supermassive Black Hole Binaries

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.

preprint2014arXiv

Chaos in the Test Particle Eccentric Kozai-Lidov Mechanism

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.

preprint2014arXiv

Detection of Galactic Center source G2 at 3.8 $μ$m during periapse passage

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 $μ$m] and L' [3.8 $μ$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-$γ$ observations can be understood with a model in which L' and Br-$γ$ 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).

preprint2014arXiv

Eccentricity growth and orbit flip in coplanar hierarchical three body systems

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.

preprint2014arXiv

Formation of Dark Matter Torii Around Supermassive Black Holes Via The Eccentric Kozai-Lidov Mechanism

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.

preprint2014arXiv

Globular Clusters and Dark Satellite Galaxies through the Stream Velocity

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.

preprint2014arXiv

Mergers and Obliquities in Stellar Triples

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.

preprint2014arXiv

The Dynamics of the Multi-planet System Orbiting Kepler-56

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.

preprint2013arXiv

Energy Feedback from X-ray Binaries in the Early Universe

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.

preprint2013arXiv

Extreme orbital evolution from hierarchical secular coupling of two giant planets

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.

preprint2013arXiv

Generation of Primordial Magnetic Fields on Linear Over-density Scales

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.

preprint2013arXiv

Resonant Post-Newtonian Eccentricity Excitation in Hierarchical Three-body 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.

preprint2013arXiv

Secular Dynamics in Hierarchical Three-Body Systems

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.

preprint2012arXiv

On the Formation of Hot Jupiters in Stellar Binaries

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.

preprint2012arXiv

Simulations of Early Baryonic Structure Formation with Stream Velocity: I. Halo Abundance

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σ(vbc), where σ(vbc) is the variance of the relative velocity, while for vbc = 3.4σ(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σ(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.

preprint2012arXiv

Simulations of Early Baryonic Structure Formation with Stream Velocity: II. The Gas Fraction

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.

preprint2011arXiv

Hot Jupiters from Secular Planet--Planet Interactions

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.

preprint2011arXiv

The Eccentric Kozai Mechanism for a Test Particle

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.

preprint2011arXiv

The nonlinear evolution of baryonic overdensities in the early universe: Initial conditions of numerical simulations

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-δ 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.

preprint2010arXiv

Orbital properties of binary minor planets

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.

preprint2009arXiv

Gas in Simulations of High Redshift Galaxies and Minihalos

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.

preprint2009arXiv

Kozai cycles, tidal friction and the dynamical evolution of binary minor planets

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.

preprint2005arXiv

Growth of Linear Perturbations before the Era of the First Galaxies

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.

preprint2005arXiv

Open Cluster Birth Analysis and Multiple Spiral Arm Sets in the Milky Way

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.