Showing 1–50 of 119 results for author: Naoz, S

Once a Triple, Not Always a Triple: The Evolution of Hierarchical Triples that Yield Merged Inner Binaries

Authors: Cheyanne Shariat, Smadar Naoz, Kareem El-Badry, Antonio C. Rodriguez, Bradley M. S. Hansen, Isabel Angelo, Alexander P. Stephan

Abstract: More than half of all main-sequence (MS) stars have one or more companions, and many of those with initial masses <8 M$_\odot$ are born in hierarchical triples. These systems feature two stars in a close orbit (the inner binary) while a tertiary star orbits them on a wider orbit (the outer binary). In hierarchical triples, three-body dynamics combined with stellar evolution drives interactions and… ▽ More More than half of all main-sequence (MS) stars have one or more companions, and many of those with initial masses <8 M$_\odot$ are born in hierarchical triples. These systems feature two stars in a close orbit (the inner binary) while a tertiary star orbits them on a wider orbit (the outer binary). In hierarchical triples, three-body dynamics combined with stellar evolution drives interactions and, in many cases, merges the inner binary entirely to create a renovated `Post-Merger Binary' (PMB). By leveraging dynamical simulations and tracking binary interactions, we explore the outcomes of merged triples and investigate whether PMBs preserve signatures of their three-body history. Our findings indicate that in 26-54% of wide double WD binaries (s>100 au), the more massive white dwarf (WD) is a merger product, implying that these DWD binaries were previously triples. Overall, we estimate that $44\pm14\%$ of observed wide DWDs originated in triple star systems and thereby have rich dynamical histories. Additionally, our results suggest that the separations of inner and outer binaries are uncorrelated at birth, providing insights into stellar formation processes. We also examine MS+MS and MS+Red Giant mergers manifesting as Blue Straggler stars (BSSs). These PMBs have orbital configurations and ages similar to most observed BSS binaries. While the triple+merger formation channel can explain the observed chemical abundances, moderate eccentricities, and companion masses in BSS binaries, it likely only accounts for $\sim$20-25% of BSSs. Meanwhile, we predict that the majority of observed single BSSs formed as collisions in triples and harbor long-period (>10 yr) companions. Furthermore, both BSS binaries and DWDs exhibit signatures of WD birth kicks. △ Less

Submitted 8 July, 2024; originally announced July 2024.

Comments: Submitted to ApJ. Comments are welcome

More Likely Than You Think: Inclination-Driving Secular Resonances are Common in Known Exoplanet Systems

Authors: Thea Faridani, Smadar Naoz, Gongjie Li, Malena Rice, Nicholas Inzunza

Abstract: Multi-planet systems face significant challenges to detection. For example, further orbiting planets have reduced signal-to-noise ratio in radial velocity detection methods, and small mutual inclinations between planets can prevent them from all transiting. One mechanism to excite mutual inclination between planets is secular resonance, where the nodal precession frequencies of the planets align s… ▽ More Multi-planet systems face significant challenges to detection. For example, further orbiting planets have reduced signal-to-noise ratio in radial velocity detection methods, and small mutual inclinations between planets can prevent them from all transiting. One mechanism to excite mutual inclination between planets is secular resonance, where the nodal precession frequencies of the planets align such as to greatly increase the efficiency of angular momentum transport between planets. These resonances can significantly misalign planets from one another, hindering detection, and typically can only occur when there are three or more planets in the system. Naively, systems can only be in resonance for particular combinations of planet semimajor axes and masses; however, effects that alter the nodal precession frequencies of the planets, such as the decay of stellar oblateness, can significantly expand the region of parameter space where resonances occur. In this work, we explore known three-planet systems, determine whether they are in (or were in) secular resonance due to evolving stellar oblateness, and demonstrate the implications of resonance on their detectability and stability. We show that about 20% of a sample of three planet transiting systems seem to undergo these resonances early in their lives. △ Less

Submitted 13 June, 2024; originally announced June 2024.

Comments: 23 pages, 11 figures, 5 tables

Analytical models for secular descents in hierarchical triple systems

Authors: Grant C. Weldon, Smadar Naoz, Bradley M. S. Hansen

Abstract: Triple body systems are prevalent in nature, from planetary to stellar to supermassive black hole scales. In a hierarchical triple system, oscillations of the inner orbit's eccentricity and inclination can be induced on secular timescales. Over many cycles, the octupole-level terms in the secular equations of motion can drive the system to extremely high eccentricities via the Eccentric Kozai-Lido… ▽ More Triple body systems are prevalent in nature, from planetary to stellar to supermassive black hole scales. In a hierarchical triple system, oscillations of the inner orbit's eccentricity and inclination can be induced on secular timescales. Over many cycles, the octupole-level terms in the secular equations of motion can drive the system to extremely high eccentricities via the Eccentric Kozai-Lidov (EKL) mechanism. The overall decrease in the inner orbit's pericenter distance has potentially dramatic effects for realistic systems, such as tidal disruption. We present an analytical approximation in the test particle limit to describe the step-wise eccentricity evolution of the inner orbit. We also integrate the equations of motion and present approximations in both the test particle limit and in the general case to describe the overall octupole-level time evolution of the eccentricity. The analytical approximations are compared to numerical simulations to show that the models accurately describe the form and timescale of the secular descent from large distances to a close-encounter distance (e.g., the Roche limit). By circumventing the need for numerical simulations to obtain the long-term behavior, these approximations can be used to readily estimate descent timescales for populations of systems. We demonstrate by calculating rates of EKL-driven migration for Hot Jupiters in stellar binaries. △ Less

Submitted 30 May, 2024; originally announced May 2024.

Comments: 16 pages, 8 figures

The Supersonic Project: Early Star Formation with the Streaming Velocity

Authors: William Lake, Claire E. Williams, Smadar Naoz, Federico Marinacci, Blakesley Burkhart, Mark Vogelsberger, Naoki Yoshida, Gen Chiaki, Avi Chen, Yeou S. Chiou

Abstract: At high redshifts ($z\gtrsim12$), the relative velocity between baryons and dark matter (the so-called streaming velocity) significantly affects star formation in low-mass objects. Streaming substantially reduces the abundance of low-mass gas objects while simultaneously allowing for the formation of supersonically-induced gas objects (SIGOs) and their associated star clusters outside of dark matt… ▽ More At high redshifts ($z\gtrsim12$), the relative velocity between baryons and dark matter (the so-called streaming velocity) significantly affects star formation in low-mass objects. Streaming substantially reduces the abundance of low-mass gas objects while simultaneously allowing for the formation of supersonically-induced gas objects (SIGOs) and their associated star clusters outside of dark matter halos. Here, we present a study of the population-level effects of streaming on star formation within both halos and SIGOs in a set of simulations with and without streaming. Notably, we find that streaming actually enhances star formation within individual halos of all masses at redshifts between $z=12$ and $z=20$. This is demonstrated both as an increased star formation rate per object as well as an enhancement of the Kennicutt-Schmidt relation for objects with streaming. We find that our simulations are consistent with some observations at high redshift, but on a population level, they continue to under-predict star formation relative to the majority of observations. However, simulations of overdense regions (both with and without streaming) agree with observations, suggesting a strategy for extracting information about the overdensity and streaming velocity in a given survey volume in future observations. △ Less

Submitted 23 May, 2024; originally announced May 2024.

Comments: 16 pages, 8 figures

Detecting gravitational-wave bursts from black hole binaries in the Galactic Center with LISA

Authors: Alan M. Knee, Jess McIver, Smadar Naoz, Isobel M. Romero-Shaw, Bao-Minh Hoang

Abstract: Stellar-mass black hole binaries (BHBs) in galactic nuclei are gravitationally perturbed by the central supermassive black hole (SMBH) of the host galaxy, potentially inducing strong eccentricity oscillations through the eccentric Kozai-Lidov (EKL) mechanism. These highly eccentric binaries emit a train of gravitational-wave (GW) bursts detectable by the Laser Interferometer Space Antenna (LISA) -… ▽ More Stellar-mass black hole binaries (BHBs) in galactic nuclei are gravitationally perturbed by the central supermassive black hole (SMBH) of the host galaxy, potentially inducing strong eccentricity oscillations through the eccentric Kozai-Lidov (EKL) mechanism. These highly eccentric binaries emit a train of gravitational-wave (GW) bursts detectable by the Laser Interferometer Space Antenna (LISA) -- a planned space-based GW detector -- with signal-to-noise ratios (SNRs) up to ${\sim}100$ per burst. In this work, we study the GW signature of BHBs orbiting our galaxy's SMBH, Sgr A$^*$, which are consequently driven to very high eccentricities. We demonstrate that an unmodeled approach using a wavelet decomposition of the data effectively yields the time-frequency properties of each burst, provided that the GW frequency peaks between $10^{-3}\,\,\mathrm{Hz}$--$10^{-1}\,\,\mathrm{Hz}$. The wavelet parameters may be used to infer the eccentricity of the binary, measuring $\log_{10}(1-e)$ within an error of $20\%$. Our proposed search method can thus constrain the parameter space to be sampled by complementary Bayesian inference methods, which use waveform templates or orthogonal wavelets to reconstruct and subtract the signal from LISA data. △ Less

Submitted 18 April, 2024; originally announced April 2024.

Comments: 10 pages, 5 figures. To be submitted to ApJL

An Extreme Ultra-Compact X-ray Binary in a Globular Cluster: Multi-Wavelength Observations of RZ2109 Explored in a Triple System Framework

Authors: Kristen C. Dage, Arash Bahramian, Smadar Naoz, Alexey Bobrick, Wasundara Athukoralalage, McKinley C. Brumback, Daryl Haggard, Arunav Kundu, Stephen E. Zepf

Abstract: The globular cluster ultraluminous X-ray source, RZ2109, is a complex and unique system which has been detected at X-ray, ultra-violet, and optical wavelengths. Based on almost 20 years of Chandra and XMM-Newton observations, the X-ray luminosity exhibits order-of-magnitude variability, with the peak flux lasting on the order of a few hours. We perform robust time series analysis on the archival X… ▽ More The globular cluster ultraluminous X-ray source, RZ2109, is a complex and unique system which has been detected at X-ray, ultra-violet, and optical wavelengths. Based on almost 20 years of Chandra and XMM-Newton observations, the X-ray luminosity exhibits order-of-magnitude variability, with the peak flux lasting on the order of a few hours. We perform robust time series analysis on the archival X-ray observations and find that this variability is periodic on a timescale of 1.3 $\pm 0.04$ days. The source also demonstrates broad [OIII] 5007 Angstrom emission, which has been observed since 2004, suggesting a white dwarf donor and therefore an ultra-compact X-ray binary. We present new spectra from 2020 and 2022, marking eighteen years of observed [OIII] emission from this source. Meanwhile, we find that the globular cluster counterpart is unusually bright in the NUV/UVW2 band. Finally, we discuss RZ2109 in the context of the eccentric Kozai Lidov mechanism and show that the observed 1.3 day periodicity can be used to place constraints on the tertiary configuration, ranging from 20 minutes (for a 0.1 ${\rm M}_\odot$ companion) to approximately 95 minutes (for a 1 ${\rm M}_\odot$ companion), if the eccentric Kozai Lidov mechanism is at the origin of the periodic variability. △ Less

Submitted 14 March, 2024; originally announced March 2024.

Comments: Accepted to MNRAS

Journal ref: 2024MNRAS.tmp..625D

Stochastic Gravitational Wave Background from Highly-Eccentric Stellar-Mass Binaries in the Milli-hertz Band

Authors: Zeyuan Xuan, Smadar Naoz, Bence Kocsis, Erez Michaely

Abstract: Many gravitational wave (GW) sources are expected to have non-negligible eccentricity in the millihertz band. These highly eccentric compact object binaries may commonly serve as a progenitor stage of GW mergers, particularly in dynamical channels where environmental perturbations bring a binary with large initial orbital separation into a close pericenter passage, leading to efficient GW emission… ▽ More Many gravitational wave (GW) sources are expected to have non-negligible eccentricity in the millihertz band. These highly eccentric compact object binaries may commonly serve as a progenitor stage of GW mergers, particularly in dynamical channels where environmental perturbations bring a binary with large initial orbital separation into a close pericenter passage, leading to efficient GW emission and a final merger. This work examines the stochastic GW background from highly eccentric ($e\gtrsim 0.9$), stellar-mass sources in the mHz band. Our findings suggest that these binaries can contribute a substantial GW power spectrum, potentially exceeding the LISA instrumental noise at $\sim 3-7$~mHz. This stochastic background is likely to be dominated by eccentric sources within the Milky Way, thus introducing anisotropy and time dependence in LISA's detection. However, given efficient search strategies to identify GW transients from highly eccentric binaries, the unresolvable noise level can be substantially lower, approaching $\sim 2$ orders of magnitude below the LISA noise curve. Therefore, we highlight the importance of characterizing stellar-mass GW sources with extreme eccentricity, especially their transient GW signals in the millihertz band. △ Less

Submitted 12 June, 2024; v1 submitted 7 March, 2024; originally announced March 2024.

Comments: 14 pages, 4 figures, accepted by PRD

An Estimate of the Binary Star Fraction Among Young Stars at the Galactic Center: Possible Evidence of a Radial Dependence

Authors: Abhimat K. Gautam, Tuan Do, Andrea M. Ghez, Devin S. Chu, Matthew W. Hosek Jr., Shoko Sakai, Smadar Naoz, Mark R. Morris, Anna Ciurlo, Zoe Haggard, Jessica R. Lu

Abstract: We present the first estimate of the intrinsic binary fraction of young stars across the central $\approx$ 0.4 pc surrounding the supermassive black hole (SMBH) at the Milky Way Galactic center (GC). This experiment searched for photometric variability in 102 young stars, using 119 nights of 10"-wide adaptive optics imaging observations taken at Keck Observatory over 16 years in the K'- and H-band… ▽ More We present the first estimate of the intrinsic binary fraction of young stars across the central $\approx$ 0.4 pc surrounding the supermassive black hole (SMBH) at the Milky Way Galactic center (GC). This experiment searched for photometric variability in 102 young stars, using 119 nights of 10"-wide adaptive optics imaging observations taken at Keck Observatory over 16 years in the K'- and H-bands. We photometrically detected three binary stars, all of which are situated more than 1" (0.04 pc) from the SMBH and one of which, S2-36, is newly reported here with spectroscopic confirmation. To convert the observed binary fraction into an estimate of the underlying binary fraction, we determined experiment sensitivity through detailed light curve simulations, incorporating photometric effects of eclipses, irradiation, and tidal distortion in binaries. The simulations assumed a population of young binaries, with stellar ages (4 Myr) and masses matched to the most probable values measured for the GC young star population and underlying binary system parameters similar to those of local massive stars. The detections and simulations imply young, massive stars in the GC have a stellar binary fraction $\geq$ 71% (68% confidence), or $\geq$ 42% (95% confidence). This inferred GC young star binary fraction is consistent with that typically seen in young stellar populations in the solar neighborhood. Furthermore, our measured binary fraction is significantly higher than that recently reported by Chu et al. (2023) based on RV measurements of young stars <~1" of the SMBH. Constrained with these two studies, the probability that the same underlying young binary fraction extends across the entire region is <1.4%. This tension provides support for a radial dependence of the binary star fraction and, therefore, for the dynamical predictions of binary merger and evaporation events close to the SMBH. △ Less

Submitted 7 February, 2024; v1 submitted 23 January, 2024; originally announced January 2024.

Comments: 51 pages, 27 figures, 7 tables. Accepted for publication in The Astrophysical Journal. Abstract revised for final manuscript

Natal Kicks from the Galactic Center and Implications on their Environment and the Roman Space Telescope

Authors: Carlos Jurado, Smadar Naoz, Casey Y. Lam, Bao-Minh Hoang

Abstract: Most galaxies, including the Milky Way, harbor a central supermassive black hole (SMBH) weighing millions to billions of solar masses. Surrounding these SMBHs are dense regions of stars and stellar remnants, such as neutron stars and black holes. Neutron stars and possibly black holes receive large natal kicks at birth on the order of hundreds of km s$^{-1}$. The natal kicks that occur in the vici… ▽ More Most galaxies, including the Milky Way, harbor a central supermassive black hole (SMBH) weighing millions to billions of solar masses. Surrounding these SMBHs are dense regions of stars and stellar remnants, such as neutron stars and black holes. Neutron stars and possibly black holes receive large natal kicks at birth on the order of hundreds of km s$^{-1}$. The natal kicks that occur in the vicinity of an SMBH may redistribute the orbital configuration of the compact objects and alter their underlying density distribution. We model the effects of natal kicks on a Galactic Center (GC) population of massive stars and stellar binaries with different initial density distributions. Using observational constraints from stellar orbits near the GC, we place an upper limit on the steepness of the initial stellar profile and find it to be core-like. In addition, we predict that $30-70 \%$ of compact objects become unbound from the SMBH due to their kicks and will migrate throughout the galaxy. Different black hole kick prescriptions lead to distinct spatial and kinematic distributions. We suggest that the Roman Space Telescope may be able to distinguish between these distributions and thus be able to differentiate natal kick mechanisms. △ Less

Submitted 26 October, 2023; originally announced October 2023.

Comments: 18 pages, 11 Figures

The Supersonic Project: Lighting up the faint end of the JWST UV luminosity function

Authors: Claire E. Williams, William Lake, Smadar Naoz, Blakesley Burkhart, Tommaso Treu, Federico Marinacci, Yurina Nakazato, Mark Vogelsberger, Naoki Yoshida, Gen Chiaki, Yeou S. Chiou, Avi Chen

Abstract: The James Webb Space Telescope (JWST) is capable of probing extremely early eras of our Universe when the supersonic relative motions between dark matter and baryonic overdensities modulate structure formation ($z>\sim 10$). We study low-mass galaxy formation including this "stream velocity" using high resolution AREPO hydrodynamics simulations, and present theoretical predictions of the UV lumino… ▽ More The James Webb Space Telescope (JWST) is capable of probing extremely early eras of our Universe when the supersonic relative motions between dark matter and baryonic overdensities modulate structure formation ($z>\sim 10$). We study low-mass galaxy formation including this "stream velocity" using high resolution AREPO hydrodynamics simulations, and present theoretical predictions of the UV luminosity function (UVLF) and galaxy stellar mass function (GSMF) down to extremely faint and low mass galaxies ($M_{UV}>\sim-15$, $10^4M_\odot<=M_*<=10^8 M_\odot)$. We show that, although the stream velocity suppresses early star formation overall, it induces a short period of rapid star formation in some larger dwarfs, leading to an enhancement in the faint-end of the UVLF at $z=12$. We demonstrate that JWST observations are close to this enhanced regime, and propose that the UVLF may constitute an important probe of the stream velocity at high redshift for JWST and future observatories. △ Less

Submitted 15 December, 2023; v1 submitted 5 October, 2023; originally announced October 2023.

Comments: 13 pages, 7 figures. Accepted ApJL

New Dynamical Channel: Wide Binaries in the Galactic Center as a Source of Binary Interactions

Authors: Erez Michaely, Smadar Naoz

Abstract: The inner $500\rm pc$ in the galactic center is dense with stars and a dynamically hot environment. Here, we focus on wide binaries as a source of tidally or collisional interactions. These binaries were previously ignored as sources of binary interaction because they are too wide to have a close pericenter passage, or they will quickly become unbound (ionized) due to gravitational interactions wi… ▽ More The inner $500\rm pc$ in the galactic center is dense with stars and a dynamically hot environment. Here, we focus on wide binaries as a source of tidally or collisional interactions. These binaries were previously ignored as sources of binary interaction because they are too wide to have a close pericenter passage, or they will quickly become unbound (ionized) due to gravitational interactions with passing neighbors. However, we show that wide binaries tend to interact more frequently with neighboring stars due to their larger cross-section for gravitational impulse interactions. These interactions mainly torque the wide system, causing it to change its eccentricity. As a result, the eccentricity might be excited to sufficiently high values, causing the binary to interact at the pericenter. As a proof of concept, we present four channels of such interactions: binary main-sequence (MS), white-dwarf (WD) - MS, black hole - MS, and lastly, WD-WD. During Galaxy's lifetime, we predict tens of thousands of MS-MS interacting binaries that may form G2-like objects later appear younger than their environment. X-ray signatures and, perhaps, supernovae may result from thousands of WD-MS and WD-WD interacting binaries from this channel. Lastly, we predict a few hundred BH-MS interacting binaries at the inner $500$~pc. △ Less

Submitted 3 October, 2023; originally announced October 2023.

Detecting Gravitational Wave Bursts From Stellar-Mass Binaries in the Milli-hertz Band

Authors: Zeyuan Xuan, Smadar Naoz, Bence Kocsis, Erez Michaely

Abstract: The dynamical formation channels of gravitational wave (GW) sources typically involve a stage when the compact object binary source interacts with the environment, which may excite its eccentricity, yielding efficient GW emission. For the wide eccentric compact object binaries, the GW emission happens mostly near the pericenter passage, creating a unique, burst-like signature in the waveform. This… ▽ More The dynamical formation channels of gravitational wave (GW) sources typically involve a stage when the compact object binary source interacts with the environment, which may excite its eccentricity, yielding efficient GW emission. For the wide eccentric compact object binaries, the GW emission happens mostly near the pericenter passage, creating a unique, burst-like signature in the waveform. This work examines the possibility of stellar-mass bursting sources in the millihertz band for future LISA detections. Because of their long lifetime ($\sim 10^{7}\rm\, yr$) and promising detectability, the number of millihertz bursting sources can be large in the local universe. For example, based on our estimates, there will be $\sim 3 - 45$ bursting binary black holes in the Milky Way, with $\sim 10^{2} - 10^{4}$ bursts detected during the LISA mission. Moreover, we find that the number of bursting sources strongly depends on their formation history. If certain regions undergo active formation of compact object binaries in the recent few million years, there will be a significantly higher bursting source fraction. Thus, the detection of millihertz GW bursts not only serves as a clue for distinguishing different formation channels, but also helps us understand the star formation history in different regions of the Milky Way. △ Less

Submitted 21 February, 2024; v1 submitted 29 September, 2023; originally announced October 2023.

Comments: 14+7 pages, 7+2 figures, accepted by ApJ

$G$ Objects and primordial black holes

Authors: Marcos M. Flores, Alexander Kusenko, Andrea M. Ghez, Smadar Naoz

Abstract: We suggest that "$G$ objects" recently discovered in the Galactic Center may be clouds of gas bound by the gravitational field of stellar-mass black holes produced in the interactions of sublunar primordial black holes with neutron stars. If dark matter is composed of primordial black holes with masses $(10^{-16} - 10^{-10}) M_\odot$, these black holes can be captured by neutron stars in the Galac… ▽ More We suggest that "$G$ objects" recently discovered in the Galactic Center may be clouds of gas bound by the gravitational field of stellar-mass black holes produced in the interactions of sublunar primordial black holes with neutron stars. If dark matter is composed of primordial black holes with masses $(10^{-16} - 10^{-10}) M_\odot$, these black holes can be captured by neutron stars in the Galactic Center, where the dark matter density is high. After the capture, the neutron star is consumed by the black hole, resulting in a population of $(1-2) M_\odot $ black holes. These stellar-mass black holes, accompanied by gaseous atmospheres, can account for the observed properties of the $G$ objects, including their resilience to tidal disruption by the supermassive black hole in the Galactic Center while also producing emission consistent with inferred luminosities. △ Less

Submitted 16 August, 2023; originally announced August 2023.

Comments: 8 pages, 2 figures, Accepted into Phys. Rev. D

Report number: IPMU23-0026

Journal ref: Phys. Rev. D 108, L061301 (2023)

Constraining a companion of the galactic center black hole, Sgr A*

Authors: Clifford M. Will, Smadar Naoz, Aurélien Hees, Alexandria Tucker, Eric Zhang, Tuan Do, Andrea Ghez

Abstract: We use 23 years of astrometric and radial velocity data on the orbit of the star S0-2 to constrain a hypothetical intermediate-mass black hole orbiting the massive black hole Sgr A* at the Galactic center. The data place upper limits on variations of the orientation of the stellar orbit (inclination, nodal angle, and pericenter) at levels between 0.02 and 0.07 degrees per year. We use a combinatio… ▽ More We use 23 years of astrometric and radial velocity data on the orbit of the star S0-2 to constrain a hypothetical intermediate-mass black hole orbiting the massive black hole Sgr A* at the Galactic center. The data place upper limits on variations of the orientation of the stellar orbit (inclination, nodal angle, and pericenter) at levels between 0.02 and 0.07 degrees per year. We use a combination of analytic estimates and full numerical integrations of the orbit of S0-2 in the presence of a black-hole binary. For a companion IMBH whose semi-major axis $a_c$ is larger than that of S0-2 (1020 a.u.), we find that in the region between 1000 and 4000 a.u., a companion black hole with mass $m_c$ between $10^3$ and $10^5 M_\odot$ is excluded, with a boundary behaving as $a_c \sim m_c^{1/3}$. For a companion with $a_c < 1020$ a.u., we find that a black hole with mass between $10^3$ and $10^5 \, M_\odot$ is again excluded, with a boundary behaving as $a_c \sim m_c^{-1/2}$. These bounds arise from quadrupolar perturbations of the orbit of S0-2. However, significantly stronger bounds on the mass of an inner companion arise from the fact that the location of S0-2 is measured relative to the bright emission of Sgr A*. As a consequence, that separation is perturbed by the ``wobble'' of Sgr A* about the center of mass between it and the companion, leading to ``apparent'' perturbations of S0-2's orbit that also include a dipole component. The result is a set of bounds as small as $400 \, M_\odot$ at 200 a.u.; the numerical simulations suggest a bound from these effects varying as $a_c \sim m_c^{-1}$. We compare and contrast our results with those from a recent analysis by the GRAVITY collaboration. △ Less

Submitted 11 December, 2023; v1 submitted 31 July, 2023; originally announced July 2023.

Comments: 9 pages, 4 figures; minor revisions to agree with published version

Journal ref: Astrophys. J. 959, 58 (2023)

The Enhanced Population of Extreme Mass-Ratio Inspirals in the LISA Band from Supermassive Black Hole Binaries

Authors: Smadar Naoz, Zoltan Haiman

Abstract: Extreme mass ratio inspirals (EMRIs) take place when a stellar-mass black hole (BH) merges with a supermassive black hole (SMBH). The gravitational wave emission from such an event is expected to be detectable by the future Laser Interferometer Space Antenna (LISA) and other mHz detectors. It was recently suggested that the EMRI rate in SMBH binary systems is orders of magnitude higher than the EM… ▽ More Extreme mass ratio inspirals (EMRIs) take place when a stellar-mass black hole (BH) merges with a supermassive black hole (SMBH). The gravitational wave emission from such an event is expected to be detectable by the future Laser Interferometer Space Antenna (LISA) and other mHz detectors. It was recently suggested that the EMRI rate in SMBH binary systems is orders of magnitude higher than the EMRI rate around a single SMBH with the same total mass. Here we show that this high rate can produce thousands of SMBH-BH sources at a redshift of unity. We predict that LISA may detect a few hundred of these EMRIs with signal-to-noise ratio above SNR>=8 within a four-year mission lifetime. The remaining sub-threshold sources will contribute to a large confusion noise, which is approximately an order of magnitude above LISA's sensitivity level. Finally, we suggest that the individually detectable systems, as well as the background noise from the sub-threshold EMRIs, can be used to constrain the SMBH binary fraction in the low-redshift Universe. △ Less

Submitted 18 September, 2023; v1 submitted 20 July, 2023; originally announced July 2023.

Comments: 9 pages, 5 figures, accepted to ApJ letter

Dynamical Evolution of White Dwarfs in Triples in the Era of Gaia

Authors: Cheyanne Shariat, Smadar Naoz, Bradley M. S. Hansen, Isabel Angelo, Erez Michaely, Alexander P. Stephan

Abstract: The $Gaia$ mission has detected many white dwarfs (WDs) in binary and triple configurations, and while observations suggest that triple stellar systems are common in our Galaxy, not much attention was devoted to WDs in triples. For stability reasons, these triples must have hierarchical configurations, i.e., two stars are on a tight orbit (the inner binary), with the third companion on a wider orb… ▽ More The $Gaia$ mission has detected many white dwarfs (WDs) in binary and triple configurations, and while observations suggest that triple stellar systems are common in our Galaxy, not much attention was devoted to WDs in triples. For stability reasons, these triples must have hierarchical configurations, i.e., two stars are on a tight orbit (the inner binary), with the third companion on a wider orbit about the inner binary. In such a system, the two orbits torque each other via the eccentric Kozai-Lidov mechanism (EKL), which can alter the orbital configuration of the inner binary. We simulate thousands of triple stellar systems for over $10$ Gyr, tracking gravitational interactions, tides, general relativity, and stellar evolution up to their WD fate. As demonstrated here, three-body dynamics coupled with stellar evolution is a critical channel to form tight WD binaries or merge a WD binary. Amongst these triples, we explore their manifestations as cataclysmic variables, Type Ia supernovae, and gravitational-wave events. The simulated systems are then compared to a sample of WD triples selected from the $Gaia$ catalog. We find that including the effect of mass loss-induced kicks is crucial for producing a distribution of the inner binary-tertiary separations that is consistent with $Gaia$ observations. Lastly, we leverage this consistency to estimate that, at minimum, $30\%$ of solar-type stars in the local $200$ parsecs were born in triples. △ Less

Submitted 7 September, 2023; v1 submitted 22 June, 2023; originally announced June 2023.

Comments: 18 pages, 7 figures, revised version accepted for publication in ApJ

Roman CCS White Paper: Characterizing the Galactic population of isolated black holes

Authors: Casey Y. Lam, Natasha Abrams, Jeff Andrews, Etienne Bachelet, Arash Bahramian, David Bennett, Valerio Bozza, Floor Broekgaarden, Sukanya Chakrabarti, William Dawson, Kareem El-Badry, Maya Fishbach, Giacomo Fragione, Scott Gaudi, Abhimat Gautam, Ryosuke Hirai, Daniel Holz, Matthew Hosek Jr., Macy Huston, Tharindu Jayasinghe, Samson Johnson, Daisuke Kawata, Naoki Koshimoto, Jessica R. Lu, Ilya Mandel , et al. (12 additional authors not shown)

Abstract: Although there are estimated to be 100 million isolated black holes (BHs) in the Milky Way, only one has been found so far, resulting in significant uncertainty about their properties. The Galactic Bulge Time Domain Survey provides the only opportunity in the coming decades to grow this catalog by order(s) of magnitude. This can be achieved if 1) Roman's astrometric potential is fully realized in… ▽ More Although there are estimated to be 100 million isolated black holes (BHs) in the Milky Way, only one has been found so far, resulting in significant uncertainty about their properties. The Galactic Bulge Time Domain Survey provides the only opportunity in the coming decades to grow this catalog by order(s) of magnitude. This can be achieved if 1) Roman's astrometric potential is fully realized in the observation strategy and software pipelines, 2) Roman's observational gaps of the Bulge are minimized, and 3) observations with ground-based facilities are taken of the Bulge to fill in gaps during non-Bulge seasons. A large sample of isolated BHs will enable a broad range of astrophysical questions to be answered, such as massive stellar evolution, origin of gravitational wave sources, supernova physics, and the growth of supermassive BHs, maximizing Roman's scientific return. △ Less

Submitted 29 June, 2023; v1 submitted 21 June, 2023; originally announced June 2023.

Comments: 20 pages. Submitted in response to Nancy Grace Roman Space Telescope white paper call: https://roman.gsfc.nasa.gov/science/ccs_white_papers.html. v2 fixes a typo in Figure 5 axis label (days --> year)

Uncovering Hidden Massive Black Hole Companions with Tidal Disruption Events

Authors: Brenna Mockler, Denyz Melchor, Smadar Naoz, Enrico Ramirez-Ruiz

Abstract: Dynamical perturbations from supermassive black hole (SMBH) binaries can increase the rates of tidal disruption events (TDEs). However, most previous work focuses on TDEs from the heavy black hole in the SMBH binary (SMBHB) system. In this work, we focus on the lighter black holes in SMBHB systems and show that they can experience a similarly dramatic increase in their TDE rate due to perturbation… ▽ More Dynamical perturbations from supermassive black hole (SMBH) binaries can increase the rates of tidal disruption events (TDEs). However, most previous work focuses on TDEs from the heavy black hole in the SMBH binary (SMBHB) system. In this work, we focus on the lighter black holes in SMBHB systems and show that they can experience a similarly dramatic increase in their TDE rate due to perturbations from a more massive companion. While the increase in TDEs around the more massive black hole is mostly due to chaotic orbital perturbations, we find that, around the smaller black hole, the eccentric Kozai-Lidov (EKL) mechanism is dominant and capable of producing a comparably large number of TDEs. In this instance, the mass derived from the light curve and spectra of TDEs triggered by the lighter SMBH companion are expected to be significant smaller than the SMBH mass estimated from galaxy scaling relations, which are dominated by the more massive companion. This apparent inconsistency can help find SMBHB candidates that are not currently accreting as active galactic nuclei (AGN) and that are at separations too small to be resolved as two distinct sources. TDEs thus provide us with an exciting opportunity to study SMBHB, in particular when they might be detected from galaxies hosting SMBHs that are too massive to disrupt sun-like stars. △ Less

Submitted 12 June, 2023; v1 submitted 8 June, 2023; originally announced June 2023.

Comments: 14 pages, 7 figures, 1 table, submitted to ApJ, corrected ref. links

Tidal Disruption Events from the Combined Effects of Two-Body Relaxation and the Eccentric Kozai-Lidov Mechanism

Authors: Denyz Melchor, Brenna Mockler, Smadar Naoz, Sanaea Rose, Enrico Ramirez-Ruiz

Abstract: Tidal disruption events (TDEs) take place when a star ventures too close to a supermassive black hole (SMBH) and becomes ruptured. One of the leading proposed physical mechanisms often invoked in the literature involves weak two-body interactions experienced by the population of stars within the host SMBH's sphere of influence, commonly referred to as two-body relaxation. This process can alter th… ▽ More Tidal disruption events (TDEs) take place when a star ventures too close to a supermassive black hole (SMBH) and becomes ruptured. One of the leading proposed physical mechanisms often invoked in the literature involves weak two-body interactions experienced by the population of stars within the host SMBH's sphere of influence, commonly referred to as two-body relaxation. This process can alter the angular momentum of stars at large distances and place them into nearly radial orbits, thus driving them to disruption. On the other hand, gravitational perturbations from an SMBH companion via the eccentric Kozai-Lidov (EKL) mechanism have also been proposed as a promising stellar disruption channel. Here we demonstrate that the combination of EKL and two-body relaxation in SMBH binaries is imperative for building a comprehensive picture of the rates of TDEs. Here we explore how the density profile of the surrounding stellar distribution and the binary orbital parameters of an SMBH companion influence the rate of TDEs. We show that this combined channel naturally produces disruptions at a rate that is consistent with observations and also naturally forms repeated TDEs, where a bound star is partially disrupted on multiple orbits. Recent observations show stars being disrupted in short-period orbits, which is challenging to explain when these mechanisms are considered independently. However, the diffusive effect of two-body relaxation, combined with the secular nature of the eccentricity excitations from EKL, is found to drive stars on short eccentric orbits at a much higher rate. △ Less

Submitted 13 June, 2023; v1 submitted 8 June, 2023; originally announced June 2023.

Comments: Fixed affiliation

The Supersonic Project: Star Formation in Early Star Clusters without Dark Matter

Authors: William Lake, Smadar Naoz, Federico Marinacci, Blakesley Burkhart, Mark Vogelsberger, Claire E. Williams, Yeou S. Chiou, Gen Chiaki, Yurina Nakazato, Naoki Yoshida

Abstract: The formation mechanism of globular clusters (GCs) has long been debated by astronomers. It was recently proposed that Supersonically Induced Gas Objects (SIGOs), which formed in the early Universe due to the supersonic relative motion of baryons and dark matter at recombination, could be the progenitors of early globular clusters. In order to become GCs, SIGOs must form stars relatively efficient… ▽ More The formation mechanism of globular clusters (GCs) has long been debated by astronomers. It was recently proposed that Supersonically Induced Gas Objects (SIGOs), which formed in the early Universe due to the supersonic relative motion of baryons and dark matter at recombination, could be the progenitors of early globular clusters. In order to become GCs, SIGOs must form stars relatively efficiently despite forming outside of dark matter halos. We investigate the potential for star formation in SIGOs using cosmological hydrodynamic simulations, including the aforementioned relative motions of baryons and dark matter, molecular hydrogen cooling in primordial gas clouds, and including explicit star formation. We find that SIGOs do form stars and that the nascent star clusters formed through this process are accreted by dark matter halos on short timescales (a few hundreds of Myr). Thus, SIGOs may be found as intact substructures within these halos, analogous to many present-day GCs. From this result, we conclude that SIGOs are capable of forming star clusters with similar properties to globular clusters in the early Universe and we discuss their detectablity by upcoming JWST surveys. △ Less

Submitted 18 September, 2023; v1 submitted 1 June, 2023; originally announced June 2023.

Comments: 12 pages, 5 figures

Let's Sweep: The Effect of Evolving $J_2$ on the Resonant Structure of a Three-Planet System

Authors: Thea H. Faridani, Smadar Naoz, Gongjie Li, Nicholas Inzunza

Abstract: Short and ultra-short planets are a peculiar type of exoplanets with periods as short as a few days or less. Although it is challenging to detect them, already several are observed with many additional candidates. If these planets have formation pathways to their longer period counterparts, they are predicted to reside in multi-planet systems. Thus, gravitational perturbation from potential planet… ▽ More Short and ultra-short planets are a peculiar type of exoplanets with periods as short as a few days or less. Although it is challenging to detect them, already several are observed with many additional candidates. If these planets have formation pathways to their longer period counterparts, they are predicted to reside in multi-planet systems. Thus, gravitational perturbation from potential planetary neighbors may affect their orbital configuration. However, due to their close proximity to their host star, they are also subjects to general relativity precession and torques from the stellar spin quadrupole moment ($J_2$). Here we show that an evolving $J_2$ due to magnetic braking, affects the magnitude and location of secular resonances of the short period planet in a multi-planet system. Thus, driving the short period planet into and out of a secular resonance, exciting the planet's eccentricity and inclination. The high inclination can hinder transit observation, and, in some cases, the high eccentricity may result in an unstable configuration. We propose that evolving $J_2$ in a multi-planet system can be critical in understanding the detectability and stability of short-period planets. △ Less

Submitted 4 May, 2023; originally announced May 2023.

Comments: 15 pages, 7 figures

Stellar Collisions in the Galactic Center: Massive Stars, Collision Remnants, and Missing Red Giants

Authors: Sanaea C. Rose, Smadar Naoz, Re'em Sari, Itai Linial

Abstract: Like most galaxies, the Milky Way harbors a supermassive black hole (SMBH) at its center, surrounded by a nuclear star cluster. In this dense star cluster, direct collisions can occur between stars before they evolve off the main-sequence. Using a statistical approach, we characterize the outcomes of these stellar collisions within the inner parsec of the Galactic Center (GC). Close to the SMBH, w… ▽ More Like most galaxies, the Milky Way harbors a supermassive black hole (SMBH) at its center, surrounded by a nuclear star cluster. In this dense star cluster, direct collisions can occur between stars before they evolve off the main-sequence. Using a statistical approach, we characterize the outcomes of these stellar collisions within the inner parsec of the Galactic Center (GC). Close to the SMBH, where the velocity dispersion is larger than the escape speed from a Sun-like star, collisions lead to mass loss. We find that the stellar population within $0.01$ pc is halved within about a Gyr because of destructive collisions. Additionally, we predict a diffuse population of peculiar low-mass stars in the GC. These stars have been divested of their outer layers in the inner $0.01$ pc before migrating to larger distances from the SMBH. Between $0.01$ and $0.1$ pc from the SMBH, collisions can result in mergers. Our results suggest that repeated collisions between lower mass stars can produce massive ($\gtrsim 10$ M$_\odot$) stars, and there may be $\sim 100$ of them residing in this region. We provide predictions on the number of G objects, dust and gas enshrouded stellar objects, that may result from main-sequence stellar collisions. Lastly, we comment on uncertainties in our model and possible connections between stellar collisions and the missing red giants in the GC. △ Less

Submitted 20 April, 2023; originally announced April 2023.

Comments: Submitted to ApJ. Comments welcome. 24 pages, 13 figures

Evidence of a Decreased Binary Fraction for Massive Stars Within 20 Milliparsecs of the Supermassive Black Hole at the Galactic Center

Authors: Devin S. Chu, Tuan Do, Andrea Ghez, Abhimat K. Gautam, Anna Ciurlo, Kelly Kosmo O'Neil, Matthew W. Hosek Jr., Aurélien Hees, Smadar Naoz, Shoko Sakai, Jessica R. Lu, Zhuo Chen, Rory O. Bentley, Eric E. Becklin, Keith Matthews

Abstract: We present the results of the first systematic search for spectroscopic binaries within the central 2 x 3 arcsec$^2$ around the supermassive black hole at the center of the Milky Way galaxy. This survey is based primarily on over a decade of adaptive optics-fed integral-field spectroscopy (R$\sim$4000), obtained as part of the Galactic Center Orbits Initiative at Keck Observatory, and has a limiti… ▽ More We present the results of the first systematic search for spectroscopic binaries within the central 2 x 3 arcsec$^2$ around the supermassive black hole at the center of the Milky Way galaxy. This survey is based primarily on over a decade of adaptive optics-fed integral-field spectroscopy (R$\sim$4000), obtained as part of the Galactic Center Orbits Initiative at Keck Observatory, and has a limiting $K$'-band magnitude of 15.8, which is at least 4 magnitudes deeper than previous spectroscopic searches for binaries at larger radii within the central nuclear star cluster. From this primary dataset, over 600 new radial velocities are extracted and reported, increasing by a factor of 3 the number of such measurements. We find no significant periodic signals in our sample of 28 stars, of which 16 are massive, young (main-sequence B) stars and 12 are low-mass, old (M and K giant) stars. Using Monte Carlo simulations, we derive upper limits on the intrinsic binary star fraction for the young star population at 47% (at 95% confidence) located $\sim$20 mpc from the black hole. The young star binary fraction is significantly lower than that observed in the field (70%). This result is consistent with a scenario in which the central supermassive black hole drives nearby stellar binaries to merge or be disrupted and may have important implications for the production of gravitational waves and hypervelocity stars. △ Less

Submitted 29 March, 2023; originally announced March 2023.

Comments: Accepted by ApJ. 22 pages, 14 figures

Quiescent and active galactic nuclei as factories of merging compact objects in the era of gravitational-wave astronomy

Authors: Manuel Arca Sedda, Smadar Naoz, Bence Kocsis

Abstract: Galactic nuclei harbouring a central supermassive black hole (SMBH), possibly surrounded by a dense nuclear cluster (NC), represent extreme environments which house a complex interplay of many physical processes that uniquely affect stellar formation, evolution, and dynamics. The discovery of gravitational waves (GW) emitted by merging black holes (BHs) and neutron stars (NSs), funnelled a huge am… ▽ More Galactic nuclei harbouring a central supermassive black hole (SMBH), possibly surrounded by a dense nuclear cluster (NC), represent extreme environments which house a complex interplay of many physical processes that uniquely affect stellar formation, evolution, and dynamics. The discovery of gravitational waves (GW) emitted by merging black holes (BHs) and neutron stars (NSs), funnelled a huge amount of work focused on understanding how compact object binaries (COBs) can pair-up and merge together. Here, we review from a theoretical standpoint how different mechanisms concur to the formation, evolution, and merger of COBs around quiescent SMBHs and active galactic nuclei (AGNs), summarizing the main predictions for current and future (GW) detections and outlining the possible features that can clearly mark a galactic nuclei origin. △ Less

Submitted 27 February, 2023; originally announced February 2023.

Comments: 73 pages, 16 Figures, 2 Tables. Invited review to be published in the Special Issue "Binary Evolution in Galactic Nuclei", eds. Roberto Capuzzo-Dolcetta and Adrian Hamers, accepted for publication in the journal "Universe"

A Stability Timescale for Non-Hierarchical Three-Body Systems

Authors: Eric Zhang, Smadar Naoz, Clifford M. Will

Abstract: The gravitational three-body problem is a fundamental problem in physics and has significant applications to astronomy. Three-body configurations are often considered stable as long the system is hierarchical; that is, the two orbital distances are well-separated. However, instability, which is often associated with significant energy exchange between orbits, takes time to develop. Assuming two ma… ▽ More The gravitational three-body problem is a fundamental problem in physics and has significant applications to astronomy. Three-body configurations are often considered stable as long the system is hierarchical; that is, the two orbital distances are well-separated. However, instability, which is often associated with significant energy exchange between orbits, takes time to develop. Assuming two massive objects in a circular orbit and a test particle in an eccentric orbit, we develop an analytical formula estimating the time it takes for the test particle's orbital energy to change by an order of itself. We show its consistency with results from N-body simulations. For eccentric orbits in particular, the instability is primarily driven not by close encounters of the test particle with one of the other bodies, but by the fundamental susceptibility of eccentric orbits to exchange energy at their periapsis. Motivated by recent suggestions that the galactic center may host an intermediate-mass black hole (IMBH) as a companion to the massive black hole Sgr A*, we use our timescale to explore the parameter space that could harbor an IMBH for the lifetime of the S-cluster of stars surrounding Sgr A*. Furthermore, we show that the orbit of an S-star can be stable for long timescales in the presence of other orbital crossing stars, thus suggesting that the S-cluster may be stable for the lifetimes of its member stars. △ Less

Submitted 2 May, 2023; v1 submitted 19 January, 2023; originally announced January 2023.

Comments: 16 pages, 8 figures, Submitted ApJ - Comments welcome

The Swansong of the Galactic Center Source X7: An Extreme Example of Tidal Evolution near the Supermassive Black Hole

Authors: Anna Ciurlo, Randall D. Campbell, Mark R. Morris, Tuan Do, Andrea M. Ghez, Eric E. Becklin, Rory O. Bentley, Devin S. Chu, Abhimat K. Gautam, Yash A. Gursahani, Aurelien Hees, Kelly Kosmo O'Neil, Jessica R. Lu, Gregory D. Martinez, Smadar Naoz, Shoko Sakai, Rainer Schoedel

Abstract: We present two decades of new high-angular-resolution near-infrared data from the W. M. Keck Observatory that reveal extreme evolution in X7, an elongated dust and gas feature, presently located half an arcsecond from the Galactic Center supermassive black hole. With both spectro-imaging observations of Br-γ line-emission and Lp (3.8 μm) imaging data, we provide the first estimate of its orbital p… ▽ More We present two decades of new high-angular-resolution near-infrared data from the W. M. Keck Observatory that reveal extreme evolution in X7, an elongated dust and gas feature, presently located half an arcsecond from the Galactic Center supermassive black hole. With both spectro-imaging observations of Br-γ line-emission and Lp (3.8 μm) imaging data, we provide the first estimate of its orbital parameters and quantitative characterization of the evolution of its morphology and mass. We find that the leading edge of X7 appears to be on a mildly eccentric (e~0.3), relatively short-period (170 years) orbit and is headed towards periapse passage, estimated to occur in ~2036. Furthermore, our kinematic measurements rule out the earlier suggestion that X7 is associated with the stellar source S0-73 or with any other point source that has overlapped with X7 during our monitoring period. Over the course of our observations, X7 has (1) become more elongated, with a current length-to-width ratio of 9, (2) maintained a very consistent long-axis orientation (position angle of 50 deg), (3) inverted its radial velocity differential from tip to tail from -50 to +80 km/sec, and (4) sustained its total brightness (12.8 Lp magnitudes at the leading edge) and color temperature (425 K), which suggest a constant mass of ~50 MEarth. We present a simple model showing that these results are compatible with the expected effect of tidal forces exerted on it by the central black hole and we propose that X7 is the gas and dust recently ejected from a grazing collision in a binary system. △ Less

Submitted 16 January, 2023; originally announced January 2023.

Comments: accepted by ApJ, 20 pages, 16 figures

The Supersonic Project: The eccentricity and rotational support of SIGOs and DM GHOSts

Authors: Claire E. Williams, Smadar Naoz, William Lake, Yeou S. Chiou, Blakesley Burkhart, Federico Marinacci, Mark Vogelsberger, Gen Chiaki, Yurina Nakazato, Naoki Yoshida

Abstract: A supersonic relative velocity between dark matter (DM) and baryons (the stream velocity) at the time of recombination induces the formation of low mass objects with anomalous properties in the early Universe. We widen the scope of the `Supersonic Project' paper series to include objects we term Dark Matter + Gas Halos Offset by Streaming (DM GHOSts)--diffuse, DM-enriched structures formed because… ▽ More A supersonic relative velocity between dark matter (DM) and baryons (the stream velocity) at the time of recombination induces the formation of low mass objects with anomalous properties in the early Universe. We widen the scope of the `Supersonic Project' paper series to include objects we term Dark Matter + Gas Halos Offset by Streaming (DM GHOSts)--diffuse, DM-enriched structures formed because of a physical offset between the centers of mass of DM and baryonic overdensities. We present an updated numerical investigation of DM GHOSts and Supersonically Induced Gas Objects (SIGOs), including the effects of molecular cooling, in high resolution hydrodynamic simulations using the AREPO code. Supplemented by an analytical understanding of their ellipsoidal gravitational potentials, we study the population-level properties of these objects, characterizing their morphology, spin, radial mass, and velocity distributions in comparison to classical structures in non-streaming regions. The stream velocity causes deviations from sphericity in both the gas and DM components and lends greater rotational support to the gas. Low mass ($<\sim 10^{5.5}$ M$_\odot$) objects in regions of streaming demonstrate core-like rotation and mass profiles. Anomalies in the rotation and morphology of DM GHOSts could represent an early Universe analogue to observed ultra-faint dwarf galaxies with variations in DM content and unusual rotation curves. △ Less

Submitted 13 February, 2023; v1 submitted 3 November, 2022; originally announced November 2022.

Comments: 27 pages, 20 figures. Accepted ApJ

Detecting Accelerating Eccentric Binaries in the LISA Band

Authors: Zeyuan Xuan, Smadar Naoz, Xian Chen

Abstract: Many gravitational wave (GW) sources in the LISA band are expected to have non-negligible eccentricity. Furthermore, many of them can undergo acceleration because they reside in the presence of a tertiary. Here we develop analytical and numerical methods to quantify how the compact binary's eccentricity enhances the detection of its peculiar acceleration. We show that the general relativistic prec… ▽ More Many gravitational wave (GW) sources in the LISA band are expected to have non-negligible eccentricity. Furthermore, many of them can undergo acceleration because they reside in the presence of a tertiary. Here we develop analytical and numerical methods to quantify how the compact binary's eccentricity enhances the detection of its peculiar acceleration. We show that the general relativistic precession pattern can disentangle the binary's acceleration-induced frequency shift from the chirp-mass-induced frequency shift in GW template fitting, thus relaxing the signal-to-noise ratio requirement for distinguishing the acceleration by a factor of $10\sim100$. Moreover, by adopting the GW templates of the accelerating eccentric compact binaries, we can enhance the acceleration measurement accuracy by a factor of $\sim100$, compared to the zero-eccentricity case, and detect the source's acceleration even if it does not change during the observational time. For example, a stellar-mass binary black hole (BBH) with moderate eccentricity in the LISA band yields an error of the acceleration measurement $\sim10^{-7}m\cdot s^{-2}$ for $\rm{SNR}=20$ and observational time of $4$ yrs. In this example, we can measure the BBHs' peculiar acceleration even when it is $\sim1\rm pc$ away from a $4\times 10^{6}\rm M_{\odot}$ SMBH. Our results highlight the importance of eccentricity to the LISA-band sources and show the necessity of developing GW templates for accelerating eccentric compact binaries. △ Less

Submitted 11 January, 2023; v1 submitted 6 October, 2022; originally announced October 2022.

Comments: 15 pages (+references), 9 figures. Accepted by PRD

Formation History of HD106906 and the Vertical Warping of Debris Disks by an External Inclined Companion

Authors: Nathaniel W. H. Moore, Gongjie Li, Lee Hassenzahl, Erika R. Nesvold, Smadar Naoz, Fred C. Adams

Abstract: HD106906 is a planetary system that hosts a wide-orbit companion, as well as an eccentric and flat debris disk, which hold important constraints on its formation and subsequent evolution. The recent observations of the companion constrain its orbit to be eccentric and inclined relative to the plane of the debris disk. Here, we show that, in the presence of the inclined companion, the debris disk q… ▽ More HD106906 is a planetary system that hosts a wide-orbit companion, as well as an eccentric and flat debris disk, which hold important constraints on its formation and subsequent evolution. The recent observations of the companion constrain its orbit to be eccentric and inclined relative to the plane of the debris disk. Here, we show that, in the presence of the inclined companion, the debris disk quickly ($\lesssim5$ Myr) becomes warped and puffy. This suggests that the current configuration of the system is relatively recent. We explore the possibility that a recent close encounter with a free floating planet could produce a companion with orbital parameters that agree with observations of HD106906b. We find that this scenario is able to recreate the structure of the debris disk while producing a companion in agreement with observation. △ Less

Submitted 28 November, 2022; v1 submitted 23 September, 2022; originally announced September 2022.

Comments: 13 pages, 7 figures. Accepted for publication in ApJ

The Supersonic Project: The Early Evolutionary Path of SIGOs

Authors: William Lake, Smadar Naoz, Blakesley Burkhart, Federico Marinacci, Mark Vogelsberger, Gen Chiaki, Yeou S. Chiou, Naoki Yoshida, Yurina Nakazato, Claire E. Williams

Abstract: Supersonically Induced Gas Objects (SIGOs) are a class of early Universe objects that have gained attention as a potential formation route for globular clusters. SIGOs have only recently begun to be studied in the context of molecular hydrogen cooling, which is key to characterizing their structure and evolution. Studying the population-level properties of SIGOs with molecular cooling is important… ▽ More Supersonically Induced Gas Objects (SIGOs) are a class of early Universe objects that have gained attention as a potential formation route for globular clusters. SIGOs have only recently begun to be studied in the context of molecular hydrogen cooling, which is key to characterizing their structure and evolution. Studying the population-level properties of SIGOs with molecular cooling is important for understanding their potential for collapse and star formation, and central for addressing whether SIGOs can survive to the present epoch. Here, we investigate the evolution of SIGOs before they form stars, using a combination of numerical and analytical analysis. For example, we study various timescales important to the evolution of SIGOs at a population level in the presence of molecular cooling. Revising the previous formulation for the critical density of collapse for SIGOs allows us to show that their prolateness tends to act as an inhibiting factor to collapse. We find that simulated SIGOs are limited by artificial two-body relaxation effects that tend to disperse them, an effect of their limited resolution. We expect that SIGOs in nature will be longer-lived compared to our simulations. Further, the fall-back timescale on which SIGOs fall into nearby dark matter halos, potentially producing a globular-cluster-like system, is frequently longer than their cooling timescale and the collapse timescale on which they shrink through gravity. Therefore, some SIGOs have time to cool and collapse outside of halos despite initially failing to exceed the critical density, even without considering metal line cooling. From this analysis we conclude that SIGOs should form stars outside of halos in non-negligible stream velocity patches in the Universe. △ Less

Submitted 9 January, 2023; v1 submitted 11 August, 2022; originally announced August 2022.

Comments: 18 pages, 11 figures

Tidal disruption events from eccentric orbits and lessons learned from the noteworthy ASASSN-14ko

Authors: Chang Liu, Brenna Mockler, Enrico Ramirez-Ruiz, Ricardo Yarza, Jamie A. P. Law-Smith, Smadar Naoz, Denyz Melchor, Sanaea Rose

Abstract: Stars grazing supermassive black holes (SMBHs) on bound orbits may survive tidal disruption, causing periodic flares. Inspired by the recent discovery of the periodic nuclear transient ASASSN-14ko, a promising candidate for a repeating tidal disruption event (TDE), we study the tidal deformation of stars approaching SMBHs on eccentric orbits. With both analytical and hydrodynamics methods, we show… ▽ More Stars grazing supermassive black holes (SMBHs) on bound orbits may survive tidal disruption, causing periodic flares. Inspired by the recent discovery of the periodic nuclear transient ASASSN-14ko, a promising candidate for a repeating tidal disruption event (TDE), we study the tidal deformation of stars approaching SMBHs on eccentric orbits. With both analytical and hydrodynamics methods, we show the overall tidal deformation of a star is similar to that in a parabolic orbit provided that the eccentricity is above a critical value. This allows one to make use of existing simulation libraries from parabolic encounters to calculate the mass fallback rate in eccentric TDEs. We find the flare structures of eccentric TDEs show a complicated dependence on both the SMBH mass and the orbital period. For stars orbiting SMBHs with relatively short periods, we predict significantly shorter-lived duration flares than those in parabolic TDEs, which can be used to predict repeating events if the mass of the SMBH can be independently measured. Using an adiabatic mass loss model, we study the flare evolution over multiple passages, and show the evolved stars can survive many more passages than main sequence stars. We apply this theoretical framework to the repeating TDE candidate ASASSN-14ko and suggest that its recurrent flares originate from a moderately massive ($M\gtrsim 1 \mathrm{M_\odot}$), extended ($\approx$ a few $\mathrm{R_\odot}$), evolved star on a grazing, bound orbit around the SMBH. Future hydrodynamics simulations of multiple tidal interactions will enable realistic models on the individual flare structure and the evolution over multiple flares. △ Less

Submitted 24 January, 2023; v1 submitted 27 June, 2022; originally announced June 2022.

Comments: 19 pages, 9 figures. Updated to accepted version (ApJ)

Journal ref: ApJ, 944, 184 (2023)

Ultra wide black-hole - neutron star binaries as a possible source for gravitational waves and short gamma ray bursts

Authors: Erez Michaely, Smadar Naoz

Abstract: The third observing run of the LIGO/Virgo/KARGA collaboration reported a few neutron star - black hole (NSBH) merger events. While NSBH mergers have yet to receive extensive theoretical attention, they may have a promising electromagnetic signature in the form of short gamma - ray bursts. Here we show that NSBH dynamical mergers can naturally form from ultra - wide binaries in the field. Flyby gra… ▽ More The third observing run of the LIGO/Virgo/KARGA collaboration reported a few neutron star - black hole (NSBH) merger events. While NSBH mergers have yet to receive extensive theoretical attention, they may have a promising electromagnetic signature in the form of short gamma - ray bursts. Here we show that NSBH dynamical mergers can naturally form from ultra - wide binaries in the field. Flyby gravitational interactions with other neighbors in the galaxy in these ultra - wide systems may result in high eccentricity that drives the binary into a merger. We show that this process can result in a merger rate at the order of $\sim 10$~Gpc$^{-3}$~yr$^{-1}$ ($\sim 5$~Gpc$^{-3}$~yr$^{-1}$) for elliptical (spiral) galaxies. This channel predicts higher merger rate with higher velocity dispersion of the host - galaxy, delay time distribution which shallower than uniform but steeper that $1/t$, higher merger rate for lower BH to NS mass ratio. △ Less

Submitted 30 May, 2022; originally announced May 2022.

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

Authors: Bao-Minh Hoang, Smadar Naoz, Melodie Sloneker

Abstract: 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… ▽ More 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. △ Less

Submitted 9 June, 2022; v1 submitted 7 April, 2022; originally announced April 2022.

Comments: Accepted to ApJ

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

Authors: Isabel Angelo, Smadar Naoz, Erik Petigura, Mason MacDougall, Alexander Stephan, Howard Isaacson, Andrew W. Howard

Abstract: 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 o… ▽ More 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. △ Less

Submitted 28 April, 2022; v1 submitted 31 March, 2022; originally announced April 2022.

Comments: Submitted to AJ;22 pages, 14 figures

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

Authors: Smadar Naoz, Sanaea C. Rose, Erez Michaely, Denyz Melchor, Enrico Ramirez-Ruiz, Brenna Mockler, Jeremy D. Schnittman

Abstract: 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 proce… ▽ More 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. △ Less

Submitted 21 March, 2022; v1 submitted 24 February, 2022; originally announced February 2022.

Comments: 11 pages, 6 figures, accepted to ApJ-Letters

The Formation of Intermediate Mass Black Holes in Galactic Nuclei

Authors: Sanaea C. Rose, Smadar Naoz, Re'em Sari, Itai Linial

Abstract: 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… ▽ More 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). △ Less

Submitted 6 July, 2022; v1 submitted 31 December, 2021; originally announced January 2022.

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

Authors: Yurina Nakazato, Gen Chiaki, Naoki Yoshida, Smadar Naoz, William Lake, Yeou S. Chiou

Abstract: 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… ▽ More 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. △ Less

Submitted 2 March, 2022; v1 submitted 19 November, 2021; originally announced November 2021.

Comments: 8 pages, 4 figures, Accepted for publication in ApJL

Raining Rocks: An analytical formulation for collision timescales in planetary systems

Authors: Santiago Torres, Smadar Naoz, Gongjie Li, Sanaea C. Rose

Abstract: The dynamical interaction of minor bodies (such as comets or asteroids) with planets plays an essential role in the planetary system's architecture and evolution. As a result of these interactions, structures like the Kuiper belt and the Oort cloud can be created. In particular, the collision of minor bodies with planets can drastically change the planet's internal and orbital evolution. We presen… ▽ More The dynamical interaction of minor bodies (such as comets or asteroids) with planets plays an essential role in the planetary system's architecture and evolution. As a result of these interactions, structures like the Kuiper belt and the Oort cloud can be created. In particular, the collision of minor bodies with planets can drastically change the planet's internal and orbital evolution. We present an analytic formulation to determine the collision timescale for a minor body to impact a planet for arbitrary geometry. By comparing with a suite of detailed N-body simulations and an analytical method for collision timescales in the solar system, we confirmed the accuracy of our analytic formulation. As a proof of concept, we focused on the collision rate of minor bodies randomly distributed around a Jupiter-like planet, emulating a Kuiper belt-like disk. We show that our analytical method yields in good agreement with the numerical simulations. The formalism presented here thus provides a succinct and accurate alternative to numerical calculations. △ Less

Submitted 13 June, 2022; v1 submitted 5 October, 2021; originally announced October 2021.

Comments: 10 pages, 7 figures

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

Authors: Thea Faridani, Smadar Naoz, Lingfeng Wei, Will M. Farr

Abstract: 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 p… ▽ More 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. △ Less

Submitted 8 June, 2022; v1 submitted 15 July, 2021; originally announced July 2021.

Comments: 14 pages, 6 figures, submitted to AAS Journals

Relativistic Dynamical Stability Criterion of Multi-Planet Systems with a Distant Companion

Authors: Lingfeng Wei, Smadar Naoz, Thea Faridani, Will M. Farr

Abstract: Multi-planetary systems are prevalent in our Galaxy. The long-term stability of such systems may be disrupted if a distant inclined companion excites the eccentricity and inclination of the inner planets via the eccentric Kozai-Lidov mechanism. However, the star-planet and the planet-planet interactions can help stabilize the system. In this work, we extend the previous stability criterion that on… ▽ More Multi-planetary systems are prevalent in our Galaxy. The long-term stability of such systems may be disrupted if a distant inclined companion excites the eccentricity and inclination of the inner planets via the eccentric Kozai-Lidov mechanism. However, the star-planet and the planet-planet interactions can help stabilize the system. In this work, we extend the previous stability criterion that only considered the companion-planet and planet-planet interactions by also accounting for short-range forces or effects, specifically, relativistic precession induced by the host star. A general analytical stability criterion is developed for planetary systems with $N$ inner planets and a relatively distant inclined perturber by comparing precession rates of relevant dynamical effects. Furthermore, we demonstrate as examples that in systems with $2$ and $3$ inner planets, the analytical criterion is consistent with numerical simulations using a combination of Gauss's averaging method and direct N-body integration. Finally, the criterion is applied to observed systems, constraining the orbital parameter space of a possible undiscovered companion. This new stability criterion extends the parameter space in which an inclined companion of multi-planet systems can inhabit. △ Less

Submitted 7 December, 2023; v1 submitted 4 June, 2021; originally announced June 2021.

Comments: 16 pages, 10 figures, 3 tables, accepted by the Astrophysical Journal

Journal ref: ApJ 923 (2021) 118

The Supersonic Project: SIGOs, a Proposed Progenitor to Globular Clusters, and their Connections to Gravitational Wave Anisotropies

Authors: William Lake, Smadar Naoz, Yeou S. Chiou, Blakesley Burkhart, Federico Marinacci, Mark Vogelsberger, Kyle Kremer

Abstract: Supersonically Induced Gas Objects (SIGOs), are structures with little to no dark matter component predicted to exist in regions of the Universe with large relative velocities between baryons and dark matter at the time of recombination. They have been suggested to be the progenitors of present-day globular clusters. Using simulations, SIGOs have been studied on small scales (around 2 Mpc), where… ▽ More Supersonically Induced Gas Objects (SIGOs), are structures with little to no dark matter component predicted to exist in regions of the Universe with large relative velocities between baryons and dark matter at the time of recombination. They have been suggested to be the progenitors of present-day globular clusters. Using simulations, SIGOs have been studied on small scales (around 2 Mpc), where these relative velocities are coherent. However, it is challenging to study SIGOs using simulations on large scales due to the varying relative velocities at scales larger than a few Mpc. Here, we study SIGO abundances semi-analytically: using perturbation theory, we predict the number density of SIGOs analytically, and compare these results to small-box numerical simulations. We use the agreement between the numerical and analytic calculations to extrapolate the large-scale variation of SIGO abundances over different stream velocities. As a result, we predict similar large-scale variations of objects with high gas densities before reionization that could possibly be observed by JWST. If indeed SIGOs are progenitors of globular clusters, then we expect a similar variation of globular cluster abundances over large scales. Significantly, we find that the expected number density of SIGOs is consistent with observed globular cluster number densities. As a proof-of-concept, and because globular clusters were proposed to be natural formation sites for gravitational wave sources from binary black hole (BBH) mergers, we show that SIGOs should imprint an anisotropy on the gravitational wave signal on the sky, consistent with SIGOs' distribution. △ Less

Submitted 24 August, 2021; v1 submitted 22 April, 2021; originally announced April 2021.

Comments: 17 pages, 8 figures, Accepted to ApJ

Effects of turbulence in the Circumnuclear Disk

Authors: Cuc K. Dinh, Jesus M. Salas, Mark R. Morris, Smadar Naoz

Abstract: A Circumnuclear Disk (CND) of molecular gas occupies the central few parsecs of the Galactic Center. It is likely subject to turbulent disruptions from violent events in its surrounding environment, but the effect of such perturbations has not yet been investigated in detail. Here we perform 3D, N-body/smoothed particle hydrodynamic (SPH) simulations with an adapted general turbulence driving meth… ▽ More A Circumnuclear Disk (CND) of molecular gas occupies the central few parsecs of the Galactic Center. It is likely subject to turbulent disruptions from violent events in its surrounding environment, but the effect of such perturbations has not yet been investigated in detail. Here we perform 3D, N-body/smoothed particle hydrodynamic (SPH) simulations with an adapted general turbulence driving method to investigate the CND's structural evolution, in particular its reaction to varied scales of injected turbulence. We find that, because of shear flow in the disk, transient arcs of gas (streams) naturally arise when turbulence is driven on large scales (up to ~4 pc), as might occur when a supernova blast wave encounters the CND. Because energetic events arise naturally and often in the central parsecs of our Galaxy, this result suggests that the transient structures that characterize the CND do not imply that the CND itself is a transient structure. We also note that features similar to the density concentrations, or "clumps", detailed in literature emerge when we account for the observed orientation of the disk and for the spatial resolution of observations. As such, clumps could be an artifact of observational limitations. △ Less

Submitted 22 September, 2021; v1 submitted 18 February, 2021; originally announced February 2021.

Comments: 19 pages, 9 figures. Accepted for publication by ApJ

The Stationary Points of the Hierarchical Three Body Problem

Authors: Bradley M. S. Hansen, Smadar Naoz

Abstract: We study the stationary points of the hierarchical three body problem in the planetary limit (m_2, m_3 << m_1) at both the quadrupole and octupole orders. We demonstrate that the extension to octupole order preserves the principal stationary points of the quadrupole solution in the limit of small outer eccentricity e_2 but that new families of stable fixed points occur in both prograde and retrogr… ▽ More We study the stationary points of the hierarchical three body problem in the planetary limit (m_2, m_3 << m_1) at both the quadrupole and octupole orders. We demonstrate that the extension to octupole order preserves the principal stationary points of the quadrupole solution in the limit of small outer eccentricity e_2 but that new families of stable fixed points occur in both prograde and retrograde cases. The most important new equilibria are those that branch off from the quadrupolar solutions and extend to large e_2. The apsidal alignment of these families is a function of mass and inner planet eccentricity, and is determined by the relative directions of precession of omega_1 and omega_2 at the quadrupole level. These new equilibria are also the most resilient to the destabilizing effects of relativistic precession. We find additional equilibria that enable libration of the inner planet argument of pericentre in the limit of radial orbits and recover the non-linear analogue of the Laplace-Lagrange solutions in the coplanar limit. Finally, we show that the chaotic diffusion and orbital flips identified with the Eccentric Kozai Lidov mechanism and its variants can be understood in terms of the stationary points discussed here. △ Less

Submitted 13 November, 2020; originally announced November 2020.

Comments: 27 pages, 36 figures

Journal ref: MNRAS, 499, 1682 (2020)

Gravitational-Wave Signatures from Compact Object Binaries in the Galactic Center

Authors: Huiyi Wang, Alexander P. Stephan, Smadar Naoz, Bao-Minh Hoang, Katelyn Breivik

Abstract: Almost every galaxy has a supermassive black hole (SMBH) residing at its center, the Milky Way included. Recent studies suggest that these unique places are expected to host a high abundance of stellar and compact object binaries. These binaries form hierarchical triple systems with the SMBH and undergo the eccentric Kozai-Lidov (EKL) mechanism. Here we estimate the detectability of potential Grav… ▽ More Almost every galaxy has a supermassive black hole (SMBH) residing at its center, the Milky Way included. Recent studies suggest that these unique places are expected to host a high abundance of stellar and compact object binaries. These binaries form hierarchical triple systems with the SMBH and undergo the eccentric Kozai-Lidov (EKL) mechanism. Here we estimate the detectability of potential Gravitational-Wave emissions from these compact objects within the frequency band of the Laser Interferometer Space Antenna (LISA) and Laser Interferometer Gravitational-Wave Observatory (LIGO) and Virgo detectors. We generate a post EKL population of stars at the onset of Roche limit crossing and follow their evolution to compact object binaries. As a proof-of-concept, we adopt two metallicities, solar metallicity ($Z = 0.02$) and $15\%$ of it ($Z = 0.003$). We demonstrate that over the observation timescale of LISA, black hole binaries (BH-BH) and white dwarf binaries provide the most prominent GW sources via the EKL assisted merger channel. Systems involving neutron stars are less observable but possibly abundant through different merger channels. Our population synthesis of BH-BH with $Z = 0.02$ ($Z = 0.003$) translate to $\sim$ $4$ ($24$) events per year with LIGO within a 1 ${\rm Gpc}^3$ sphere. We also estimated the number of binaries visible in the LISA band within the inner parsec of our galactic center (and possibly other galaxies) to be about 14 - 150 WD-WD, 0 - 2 NS-BH, 0.2 - 4 NS-NS, and 0.3 - 20 BH-BH. △ Less

Submitted 6 May, 2021; v1 submitted 29 October, 2020; originally announced October 2020.

Comments: 19 pages, 8 figures

Giant Planets, Tiny Stars: Producing Short-Period Planets around White Dwarfs with the Eccentric Kozai-Lidov Mechanism

Authors: Alexander P. Stephan, Smadar Naoz, B. Scott Gaudi

Abstract: The recent discoveries of WD J091405.30+191412.25 (WD J0914 hereafter), a white dwarf likely accreting material from an ice giant planet, and WD 1856+534 b (WD 1856 b hereafter), a Jupiter-sized planet transiting a white dwarf, are the first direct evidence of giant planets orbiting white dwarfs. However, for both systems the observations indicate that the planets' current orbital distances would… ▽ More The recent discoveries of WD J091405.30+191412.25 (WD J0914 hereafter), a white dwarf likely accreting material from an ice giant planet, and WD 1856+534 b (WD 1856 b hereafter), a Jupiter-sized planet transiting a white dwarf, are the first direct evidence of giant planets orbiting white dwarfs. However, for both systems the observations indicate that the planets' current orbital distances would have put them inside the stellar envelope during the red giant phase, implying that the planets must have migrated to their current orbits after their host stars became white dwarfs. Furthermore, WD J0914 is a very hot white dwarf with a short cooling time that indicates a fast migration mechanism. Here, we demonstrate that the Eccentric Kozai-Lidov (EKL) Mechanism, combined with stellar evolution and tidal effects, can naturally produce the observed orbital configurations, assuming that the white dwarfs have distant stellar companions. Indeed, WD 1856 is part of a stellar triple system, being a distant companion to a stellar binary. We provide constraints for the orbital and physical characteristics for the potential stellar companion of WD J0914 and determine the initial orbital parameters of the WD 1856 system. △ Less

Submitted 5 November, 2021; v1 submitted 20 October, 2020; originally announced October 2020.

Comments: Updated version, includes more details regarding observational constraints, accepted for publication in ApJ

The effects of turbulence on galactic nuclear gas rings

Authors: Jesus M. Salas, Smadar Naoz, Mark R. Morris

Abstract: The gas dynamics in the inner few kiloparsecs of barred galaxies often results in configurations that give rise to nuclear gas rings. However, the generic dynamical description of the formation of galactic nuclear rings does not take into account the effects of thermal pressure and turbulence. Here we perform 3D hydrodynamic simulations of gas in a galactic barred potential out to a radius of $2$~… ▽ More The gas dynamics in the inner few kiloparsecs of barred galaxies often results in configurations that give rise to nuclear gas rings. However, the generic dynamical description of the formation of galactic nuclear rings does not take into account the effects of thermal pressure and turbulence. Here we perform 3D hydrodynamic simulations of gas in a galactic barred potential out to a radius of $2$~kpc. We include self-gravity and a large-scale turbulence-driving module. We study how the formation of gaseous nuclear rings is affected by changing the bar pattern speed and the strength of the turbulence, and compare the results to simulations with a relatively high sound speed. We use two values for the bar pattern speed that have been discussed in the literature for our Milky Way Galaxy (40 and 63 km s$^{-1}$ kpc$^{-1}$). Our results show that turbulence produces broader and smaller nuclear rings, and enhances gas migration towards the inner few pc of the galaxy, compared to simulations without turbulence. △ Less

Submitted 26 April, 2021; v1 submitted 8 October, 2020; originally announced October 2020.

Comments: 13 pages, 3 figures, 2 tables

Black Hole Mergers from Hierarchical Triples in Dense Star Clusters

Authors: Miguel A. S. Martinez, Giacomo Fragione, Kyle Kremer, Sourav Chatterjee, Carl L. Rodriguez, Johan Samsing, Claire S. Ye, Newlin C. Weatherford, Michael Zevin, Smadar Naoz, Frederic A. Rasio

Abstract: Hierarchical triples are expected to be produced by the frequent binary-mediated interactions in the cores of globular clusters. In some of these triples, the tertiary companion can drive the inner binary to merger following large eccentricity oscillations, as a result of the eccentric Kozai-Lidov mechanism. In this paper, we study the dynamics and merger rates of black hole (BH) hierarchical trip… ▽ More Hierarchical triples are expected to be produced by the frequent binary-mediated interactions in the cores of globular clusters. In some of these triples, the tertiary companion can drive the inner binary to merger following large eccentricity oscillations, as a result of the eccentric Kozai-Lidov mechanism. In this paper, we study the dynamics and merger rates of black hole (BH) hierarchical triples, formed via binary--binary encounters in the CMC Cluster Catalog, a suite of cluster simulations with present-day properties representative of the Milky Way's globular clusters. We compare the properties of the mergers from triples to the other merger channels in dense star clusters, and show that triple systems do not produce significant differences in terms of mass and effective spin distribution. However, they represent an important pathway for forming eccentric mergers, which could be detected by LIGO--Virgo/KAGRA (LVK), and future missions such as LISA and DECIGO. We derive a conservative lower limit for the merger rate from this channel of $0.35$ Gpc$^{-3}$yr$^{-1}$ in the local Universe and up to $\sim9\%$ of these events may have a detectable eccentricity at LVK design sensitivity. Additionally, we find that triple systems could play an important role in retaining second-generation BHs, which can later merge again in the core of the host cluster. △ Less

Submitted 21 September, 2020; v1 submitted 17 September, 2020; originally announced September 2020.

Comments: 21 Pages, 11 Figures, 2 Tables, Accepted for publication in ApJ

On Socially Distant Neighbors: Using Binaries to Constrain the Density of Objects in the Galactic Center

Authors: Sanaea C. Rose, Smadar Naoz, Abhimat K. Gautam, Andrea M. Ghez, Tuan Do, Devin Chu, Eric Becklin

Abstract: Stars often reside in binary configurations. The nuclear star cluster surrounding the supermassive black hole (SMBH) in the Galactic Center (GC) is expected to include a binary population. In this dense environment, a binary frequently encounters and interacts with neighboring stars. These interactions vary from small perturbations to violent collisions. In the former case, weak gravitational inte… ▽ More Stars often reside in binary configurations. The nuclear star cluster surrounding the supermassive black hole (SMBH) in the Galactic Center (GC) is expected to include a binary population. In this dense environment, a binary frequently encounters and interacts with neighboring stars. These interactions vary from small perturbations to violent collisions. In the former case, weak gravitational interactions unbind a soft binary over the evaporation timescale, which depends on the binary properties as well as the density of surrounding objects and velocity dispersion. Similarly, collisions can also unbind a binary, and the collision rate depends on the density. Thus, the detection of a binary with known properties can constrain the density profile in the GC with implications for the number of compact objects, which are otherwise challenging to detect. We estimate the density necessary to unbind a binary within its lifetime for an orbit of arbitrary eccentricity about the SMBH. We find that the eccentricity has a minimal impact on the density constraint. In this proof of concept, we demonstrate that this procedure can probe the density in the GC using hypothetical young and old binaries as examples. Similarly, a known density profile provides constraints on the binary orbital separation. Our results highlight the need to consider multiple dynamical processes in tandem. In certain cases, often closer to the SMBH, the collision timescale rather than the evaporation timescale gives the more stringent density constraint, while other binaries farther from the SMBH provide unreliable density constraints because they migrate inwards due to mass segregation. △ Less

Submitted 2 November, 2020; v1 submitted 14 August, 2020; originally announced August 2020.

Comments: Accepted to ApJ, 20 pages, 11 figures

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

Authors: Yeou S. Chiou, Smadar Naoz, Blakesley Burkhart, Federico Marinacci, Mark Vogelsberger

Abstract: 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… ▽ More 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. △ Less

Submitted 8 November, 2020; v1 submitted 6 August, 2020; originally announced August 2020.

Comments: 12 pages, 9 figures, accepted by ApJ, expanded motivation and updated figures

Demographics of triple systems in dense star clusters

Authors: Giacomo Fragione, Miguel A. S. Martinez, Kyle Kremer, Sourav Chatterjee, Carl L. Rodriguez, Claire S. Ye, Newlin C. Weatherford, Smadar Naoz, Frederic A. Rasio

Abstract: 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 fo… ▽ More 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. △ Less

Submitted 22 July, 2020; originally announced July 2020.

Comments: 28 pages, 12 figures, 2 tables, accepted by ApJ