Showing 1–50 of 107 results for author: Park, M

Metals in Star-Forming Galaxies with KCWI. I. Methodology and First Results on the Abundances of Iron, Magnesium, and Oxygen

Authors: Zhuyun Zhuang, Evan N. Kirby, Charles C. Steidel, Mithi A. C. de los Reyes, Nikolaus Z. Prusinski, N. Leethochawalit, Minjung Park, Charlie Conroy, Evan H. Nuñez

Abstract: Understanding the chemical enrichment of different elements is crucial to gaining a complete picture of galaxy chemical evolution. In this study, we present a new sample of 46 low-redshift, low-mass star-forming galaxies at $M_*\sim 10^{8-10}M_{\odot}$ along with two quiescent galaxies at $M_*\sim 10^{8.8}M_{\odot}$ observed with the Keck Cosmic Web Imager (KCWI), aiming to investigate the chemica… ▽ More Understanding the chemical enrichment of different elements is crucial to gaining a complete picture of galaxy chemical evolution. In this study, we present a new sample of 46 low-redshift, low-mass star-forming galaxies at $M_*\sim 10^{8-10}M_{\odot}$ along with two quiescent galaxies at $M_*\sim 10^{8.8}M_{\odot}$ observed with the Keck Cosmic Web Imager (KCWI), aiming to investigate the chemical evolution of galaxies in the transition zone between Local Group satellites and massive field galaxies. We develop a novel method to simultaneously determine stellar abundances of iron and magnesium in star-forming galaxies. With the gas-phase oxygen abundance (O/H)$_{\rm g}$ measured using the strong line method, we are able to make the first-ever apples-to-apples comparison of $α$ elements in the stars and the ISM. We find that the [Mg/H]$_*$-[O/H]$_{\rm g}$ relation is much tighter than the [Fe/H]$_*$-[O/H]$_{\rm g}$ relation, which can be explained by the similar production processes of $α$ elements. Most galaxies in our sample exhibit higher [O/H]$_{\rm g}$ than [Fe/H]$_*$ and [Mg/H]$_*$. In addition, we construct mass-metallicity relations (MZRs) measured as three different elements (Fe$_*$, Mg$_*$, O$_{\rm g}$). Compared to the gas O-MZR, the stellar Fe- and Mg-MZRs show larger scatter driven by variations in specific star formation rates (sSFR), with star-forming galaxies exhibiting higher sSFR and lower stellar abundances at fixed mass. The excess of [O/H]$_{\rm g}$ compared to stellar abundances as well as the anti-correlation between sSFR and stellar abundance suggests that galaxy quenching of intermediate-mass galaxies at $M_*\sim 10^{8-10}M_{\odot}$ is primarily driven by starvation. △ Less

Submitted 5 July, 2024; originally announced July 2024.

Comments: 36 pages, 16 figures, 4 tables; accepted for publication in ApJ. Main results in Figure 7, 10 and 11

AGN Feedback in Quiescent Galaxies at Cosmic Noon Traced by Ionized Gas Emission

Authors: Letizia Bugiani, Sirio Belli, Minjung Park, Rebecca L. Davies, J. Trevor Mendel, Benjamin D. Johnson, Amir H. Khoram, Chloë Benton, Andrea Cimatti, Charlie Conroy, Razieh Emami, Joel Leja, Yijia Li, Gabriel Maheson, Elijah P. Mathews, Rohan P. Naidu, Erica J. Nelson, Sandro Tacchella, Bryan A. Terrazas, Rainer Weinberger

Abstract: We analyze ionized gas emission lines in deep rest-frame optical spectra of 16 quiescent galaxies at redshift $1.7<z<3.5$ observed with JWST/NIRSpec by the Blue Jay survey. Robust detection of emission lines in $75\%$ of the sample indicates the presence of ongoing ionizing sources in this passive population. The H$α$ line luminosities confirm that the population is quiescent, with star formation… ▽ More We analyze ionized gas emission lines in deep rest-frame optical spectra of 16 quiescent galaxies at redshift $1.7<z<3.5$ observed with JWST/NIRSpec by the Blue Jay survey. Robust detection of emission lines in $75\%$ of the sample indicates the presence of ongoing ionizing sources in this passive population. The H$α$ line luminosities confirm that the population is quiescent, with star formation rates that are at least ten times lower than the main sequence of star formation. The quiescent sample is clearly separate from the star-forming population in line diagnostic diagrams, and occupies a region usually populated by active galactic nuclei (AGN). Analysis of the observed line ratios, equivalent widths, and velocity dispersions leads us to conclude that in most cases the gas is ionized by AGN activity, despite the lack of X-ray detections. A subset of the sample also hosts ionized and/or neutral outflows. Our results show, for the first time using a representative sample, that low luminosity AGN are extremely common among quiescent galaxies at high redshift. These low luminosity AGN may play a key role in quenching star formation and in maintaining massive galaxies quiescent from Cosmic Noon to $z\sim0$. △ Less

Submitted 12 June, 2024; originally announced June 2024.

Comments: 20 pages, 18 figures. Submitted to ApJ

Accelerated protons produced by magnetic Penrose process in Sgr A*

Authors: Myeonghwan Oh, Myeong-Gu Park

Abstract: Typical mechanisms to extract energies from a rotating black hole are the Blandford-Znajek process and the Penrose process. The Penrose process requires a special condition that is difficult to occur in common astrophysical situations. However, the magnetic Penrose process (MPP) does not require such a special condition, and can produce ultra-high energy cosmic rays. When neutrons decay near a rot… ▽ More Typical mechanisms to extract energies from a rotating black hole are the Blandford-Znajek process and the Penrose process. The Penrose process requires a special condition that is difficult to occur in common astrophysical situations. However, the magnetic Penrose process (MPP) does not require such a special condition, and can produce ultra-high energy cosmic rays. When neutrons decay near a rotating black hole, the MPP efficiency of the produced proton is maximized. The supermassive black hole in Sagittarius A* (Sgr A*) is likely to have a radiatively inefficient accretion flow that is hot enough to produce neutrons by nuclear reactions, which can be subsequently accelerated to high-energy by the MPP. We calculate the production rate of accelerated protons from the Sgr A* to estimated the gamma-ray flux at Earth produced by these accelerated protons and the flux of the accelerated protons themselves transported from Sgr A* to Earth. We find that these very high-energy gamma rays ($E_γ\gtrsim10\,\mathrm{TeV}$) amount to a significant fraction of the flux of the gamma-ray from the HESS J1745-290 and the central molecular zone around $100\,\mathrm{TeV}$. The accelerated proton flux, when the dimensionless spin parameter $a_{*}=0.5$ and the magnetic field strength in the vicinity of the black hole $B_{0}=100\,\mathrm{G}$, is about $1.6-4.1\%$ of the cosmic ray proton flux from KASCADE experiment at about $1\,\mathrm{PeV}$. Due to the finite decay time of neutrons which need to be transported from the accretion flow to the acceleration zone, our acceleration model can operate only around black holes with mass not much greater than $\sim10^8\,M_\odot$. △ Less

Submitted 11 July, 2024; v1 submitted 7 June, 2024; originally announced June 2024.

Comments: Published in Physical Review D

Journal ref: Phys. Rev. D 110, 023015 (2024)

Discrepancies Between JWST Observations and Simulations of Quenched Massive Galaxies at $z > 3$: A Comparative Study With IllustrisTNG and ASTRID

Authors: Emma Jane Weller, Fabio Pacucci, Yueying Ni, Lars Hernquist, Minjung Park

Abstract: Recent JWST observations have uncovered an unexpectedly large population of massive quiescent galaxies at $z>3$. Using the cosmological simulations IllustrisTNG and ASTRID, we identify analogous galaxies and investigate their abundance, formation, quenching mechanisms, and post-quenching evolution for stellar masses $9.5 < \log_{10}{(M_\star/{\rm M}_\odot)} < 12$. We apply three different quenchin… ▽ More Recent JWST observations have uncovered an unexpectedly large population of massive quiescent galaxies at $z>3$. Using the cosmological simulations IllustrisTNG and ASTRID, we identify analogous galaxies and investigate their abundance, formation, quenching mechanisms, and post-quenching evolution for stellar masses $9.5 < \log_{10}{(M_\star/{\rm M}_\odot)} < 12$. We apply three different quenching definitions and find that both simulations significantly underestimate the comoving number density of quenched massive galaxies at $z \gtrsim 3$ compared to JWST observations by up to $\sim 2$ dex. This fact highlights the necessity for improved physical models of AGN feedback in galaxy formation simulations. In both simulations, the high-$z$ quenched massive galaxies often host overmassive central black holes above the standard $M_{BH}-M_\star$ relation, implying that the AGN feedback plays a crucial role in quenching galaxies in the early Universe. The typical quenching timescales for these galaxies are $\sim 200-600$ Myr. IllustrisTNG primarily employs AGN kinetic feedback, while ASTRID relies on AGN thermal feedback, which is less effective and has a longer quenching timescale. We also study the post-quenching evolution of the high-$z$ massive quiescent galaxies and find that many experience subsequent reactivation of star formation, evolving into primary progenitors of $z=0$ brightest cluster galaxies. △ Less

Submitted 4 June, 2024; originally announced June 2024.

Comments: Submitted to The Astrophysical Journal. 13 pages, 12 figures

Cue: A Fast and Flexible Photoionization Emulator for Modeling Nebular Emission Powered By Almost Any Ionizing Source

Authors: Yijia Li, Joel Leja, Benjamin D. Johnson, Sandro Tacchella, Rebecca Davies, Sirio Belli, Minjung Park, Razieh Emami

Abstract: The complex physics governing nebular emission in galaxies, particularly in the early universe, often defy simple low-dimensional models. This has proven to be a significant barrier in understanding the (often diverse) ionizing sources powering this emission. We present Cue, a highly flexible tool for interpreting nebular emission across a wide range of abundances and ionizing conditions of galaxi… ▽ More The complex physics governing nebular emission in galaxies, particularly in the early universe, often defy simple low-dimensional models. This has proven to be a significant barrier in understanding the (often diverse) ionizing sources powering this emission. We present Cue, a highly flexible tool for interpreting nebular emission across a wide range of abundances and ionizing conditions of galaxies at different redshifts. Unlike typical nebular models used to interpret extragalactic nebular emission, our model does not require a specific ionizing spectrum as a source, instead approximating the ionizing spectrum with a 4-part piece-wise power-law. We train a neural net emulator based on the CLOUDY photoionization modeling code and make self-consistent nebular continuum and line emission predictions. Along with the flexible ionizing spectra, we allow freedom in [O/H], [N/O], [C/O], gas density, and total ionizing photon budget. This flexibility allows us to either marginalize over or directly measure the incident ionizing radiation, thereby directly interrogating the source of the ionizing photons in distant galaxies via their nebular emission. Our emulator demonstrates a high accuracy, with $\sim$1% uncertainty in predicting the nebular continuum and $\sim$5% uncertainty in the emission lines. Mock tests suggest Cue is well-calibrated and produces useful constraints on the ionizing spectra when $S/N (\mathrm{H}_α) \gtrsim 10$, and furthermore capable of distinguishing between the ionizing spectra predicted by single and binary stellar models. The compute efficiency of neural networks facilitates future applications of Cue for rapid modeling of the nebular emission in large samples and Monte Carlo sampling techniques. △ Less

Submitted 7 May, 2024; originally announced May 2024.

Comments: submitted to ApJ

Impacts of bar-driven shear and shocks on star formation

Authors: Taehyun Kim, Dimitri A. Gadotti, Miguel Querejeta, Isabel Pérez, Almudena Zurita, Justus Neumann, Glenn van de Ven, Jairo Méndez-Abreu, Adriana de Lorenzo-Cáceres, Patricia Sánchez-Blázquez, Francesca Fragkoudi, Lucimara P. Martins, Luiz A. Silva-Lima, Woong-Tae Kim, Myeong-gu Park

Abstract: Bars drive gas inflow. As the gas flows inwards, shocks and shear occur along the bar dust lanes. Such shocks and shear can affect the star formation and change the gas properties. For four barred galaxies, we present Hα velocity gradient maps that highlight bar-driven shocks and shear using data from the PHANGS-MUSE and PHANGS-ALMA surveys which allow us to study bar kinematics in unprecedented d… ▽ More Bars drive gas inflow. As the gas flows inwards, shocks and shear occur along the bar dust lanes. Such shocks and shear can affect the star formation and change the gas properties. For four barred galaxies, we present Hα velocity gradient maps that highlight bar-driven shocks and shear using data from the PHANGS-MUSE and PHANGS-ALMA surveys which allow us to study bar kinematics in unprecedented detail. Velocity gradients are enhanced along the bar dust lanes, where shocks and shear are shown to occur in numerical simulations. Velocity gradient maps also efficiently pick up expanding shells around HII regions. We put pseudo slits on the regions where velocity gradients are enhanced and find that Hα and CO velocities jump up to ~170 km/s, even after removing the effects of circular motions due to the galaxy rotation. Enhanced velocity gradients either coincide with the peak of CO intensity along the bar dust lanes or are slightly offset from CO intensity peaks, depending on the objects. Using the BPT diagnostic, we identify the source of ionization on each spaxel and find that star formation is inhibited in the high velocity gradient regions of the bar, and the majority of those regions are classified as LINER or composite. This implies that star formation is inhibited where bar-driven shear and shocks are strong. Our results are consistent with the results from the numerical simulations that show star formation is inhibited in the bar where shear force is strong. △ Less

Submitted 30 April, 2024; originally announced May 2024.

Comments: 26 pages, Accepted for publication in ApJ

Widespread rapid quenching at cosmic noon revealed by JWST deep spectroscopy

Authors: Minjung Park, Sirio Belli, Charlie Conroy, Benjamin D. Johnson, Rebecca L. Davies, Joel Leja, Sandro Tacchella, J. Trevor Mendel, Chloë Benton, Letizia Bugiani, Razieh Emami, Amirhossein Khoram, Yijia Li, Gabriel Maheson, Elijah P. Mathews, Rohan P. Naidu, Erica J. Nelson, Bryan A. Terrazas, Rainer Weinberger

Abstract: Massive quiescent galaxies in the young universe are expected to be quenched rapidly, but it is unclear whether they all experience starbursts before quenching and what physical mechanism drives rapid quenching. We study 16 massive quiescent galaxies ($\log(M_\star/M_\odot) > 10$) at $z\sim2$ selected from a representative sample of the Blue Jay survey. We reconstruct their star formation historie… ▽ More Massive quiescent galaxies in the young universe are expected to be quenched rapidly, but it is unclear whether they all experience starbursts before quenching and what physical mechanism drives rapid quenching. We study 16 massive quiescent galaxies ($\log(M_\star/M_\odot) > 10$) at $z\sim2$ selected from a representative sample of the Blue Jay survey. We reconstruct their star formation histories by fitting spectral energy distribution models to the JWST/NIRSpec $R\sim1000$ spectra. We find that massive quiescent galaxies can be split into three categories with roughly equal numbers of galaxies according to their SFHs: 1) Relatively old galaxies quenched at early epochs; 2) Galaxies that are rapidly and recently quenched after a flat or bursty formation history (depending on the assumed prior); 3) Galaxies that are rapidly and recently quenched after a major starburst. Most recently quenched galaxies show neutral gas outflows, probed by blueshifted $\rm Na\,I\,D$ absorption, and ionized gas emission, with line ratios consistent with active galactic nucleus (AGN) diagnostics. This suggests that AGN activity drives multi-phase gas outflows, leading to rapid quenching. By tracing back the SFHs of the entire sample, we predict the number density of massive quiescent galaxies at $z=4-6$: $n=3.0\pm1.4\times10^{-5}\,\rm Mpc^{-3}$. The two oldest massive quiescent galaxies in our sample appear to have extremely early formation and quenching ($z\gtrsim6$), possibly descendants of early post-starbursts at $z>3$. These galaxies still show neutral gas reservoirs and low-level star formation, consistent with weak H$α$ emission, perhaps because the ejective AGN feedback that caused rapid quenching has weakened over time. △ Less

Submitted 27 April, 2024; originally announced April 2024.

Comments: Submitted to ApJ; 21 pages, 9 figures, 3 tables (excluding appendix)

Cosmic evolution of black hole-spin and galaxy orientations: clues from the NewHorizon and Galactica simulations

Authors: Sebastien Peirani, Yasushi Suto, Ricarda S. Beckmann, Marta Volonteri, Yen-Ting Lin, Yohan Dubois, Sukyoung K. Yi, Christophe Pichon, Katarina Kraljic, Minjung Park, Julien Devriendt, San Han, Wei-Huai Chen

Abstract: (Reduced) Using the recent cosmological high-resolution zoom-in simulations, NewHorizon and Galactica, in which the evolution of black hole spin is followed on the fly, we have tracked the cosmic history of a hundred of black holes (BHs) with a mass greater than 2x10^4 Ms. For each of them, we have studied the variations of the three dimensional angle (Psi) subtended between the BH spins and the a… ▽ More (Reduced) Using the recent cosmological high-resolution zoom-in simulations, NewHorizon and Galactica, in which the evolution of black hole spin is followed on the fly, we have tracked the cosmic history of a hundred of black holes (BHs) with a mass greater than 2x10^4 Ms. For each of them, we have studied the variations of the three dimensional angle (Psi) subtended between the BH spins and the angular momentum vectors of their host galaxies. The analysis of the individual evolution of the most massive BHs suggests that they are generally passing by three different regimes. First, for a short period after their birth, low mass BHs (<3x10^4 Ms) are rapidly spun up by gas accretion and their spin tends to be aligned with their host galaxy spin. Then follows a second phase in which the accretion of gas onto low mass BHs (<10^5 Ms) is quite chaotic and inefficient, reflecting the complex and disturbed morphologies of forming proto-galaxies at high redshifts. The variations of Psi are rather erratic during this phase and are mainly driven by the rapid changes of the direction of the galaxy angular momentum. Then, in a third and long phase, BHs are generally well settled in the center of galaxies around which the gas accretion becomes much more coherent (>10^5 Ms). In this case, the BH spins tend to be well aligned with the angular momentum of their host galaxy and this configuration is generally stable even though BH merger episodes can temporally induce misalignment. We have also derived the distributions of cos(Psi) at different redshifts and found that BHs and galaxy spins are generally aligned. Finally, based on a Monte Carlo method, we also predict statistics for the 2-d projected spin-orbit angles lambda. In particular, the distribution of lambda traces well the alignment tendency in the 3-d analysis. Such predictions provide an interesting background for future observational analyses. △ Less

Submitted 25 March, 2024; v1 submitted 8 January, 2024; originally announced January 2024.

Comments: Accepted for publication in A&A (12/03/2024)

Journal ref: A&A 686, A233 (2024)

Observational Signatures of AGN Feedback in the Morphology and the Ionization States of Milky Way-like Galaxies

Authors: Nadia Qutob, Razieh Emami, Kung-Yi Su, Randall Smith, Lars Hernquist, Dian P. Triani, Cameron Hummels, Drummond Fielding, Philip F. Hopkins, Rachel S. Somerville, David R. Ballantyne, Mark Vogelsberger, Grant Tremblay, James F. Steiner, Douglas Finkbeiner, Ramesh Narayan, Minjung Park, Josh Grindlay, Priyamvada Natarajan, Christopher C. Hayward, Dušan Kereš, Sam B. Ponnada, Sirio Belli, Rebecca Davies, Gabriel Maheson , et al. (2 additional authors not shown)

Abstract: We make an in-depth analysis of different AGN jet models' signatures, inducing quiescence in galaxies with a halo mass of $10^{12} M_\odot$. Three jet models, including cosmic ray-dominant, hot thermal, and precessing kinetic jets, are studied at two energy flux levels each, compared to a jet-free, stellar feedback-only simulation. We examine the distribution of Mg II, O VI, and O VIII ions, along… ▽ More We make an in-depth analysis of different AGN jet models' signatures, inducing quiescence in galaxies with a halo mass of $10^{12} M_\odot$. Three jet models, including cosmic ray-dominant, hot thermal, and precessing kinetic jets, are studied at two energy flux levels each, compared to a jet-free, stellar feedback-only simulation. We examine the distribution of Mg II, O VI, and O VIII ions, alongside gas temperature and density profiles. Low-energy ions, like Mg II, concentrate in the ISM, while higher energy ions, e.g., O VIII, prevail at the AGN jet cocoon's edge. High-energy flux jets display an isotropic ion distribution with lower overall density. High-energy thermal or cosmic ray jets pressurize at smaller radii, significantly suppressing core density. The cosmic ray jet provides extra pressure support, extending cool and warm gas distribution. A break in the ion-to-mass ratio slope in O VI and O VIII is demonstrated in the ISM-to-CGM transition (between 10-30 kpc), growing smoothly towards the CGM at greater distances. △ Less

Submitted 22 December, 2023; originally announced December 2023.

Comments: 17 pages, 13 figures

JWST Reveals Widespread AGN-Driven Neutral Gas Outflows in Massive z ~ 2 Galaxies

Authors: Rebecca L. Davies, Sirio Belli, Minjung Park, J. Trevor Mendel, Benjamin D. Johnson, Charlie Conroy, Chloë Benton, Letizia Bugiani, Razieh Emami, Joel Leja, Yijia Li, Gabriel Maheson, Elijah P. Mathews, Rohan P. Naidu, Erica J. Nelson, Sandro Tacchella, Bryan A. Terrazas, Rainer Weinberger

Abstract: We use deep JWST/NIRSpec R~1000 slit spectra of 113 galaxies at 1.7 < z < 3.5, selected from the mass-complete Blue Jay survey, to investigate the prevalence and typical properties of neutral gas outflows at cosmic noon. We detect excess Na I D absorption (beyond the stellar contribution) in 46% of massive galaxies ($\log$ M$_*$/M$_\odot >$ 10), with similar incidence rates in star-forming and que… ▽ More We use deep JWST/NIRSpec R~1000 slit spectra of 113 galaxies at 1.7 < z < 3.5, selected from the mass-complete Blue Jay survey, to investigate the prevalence and typical properties of neutral gas outflows at cosmic noon. We detect excess Na I D absorption (beyond the stellar contribution) in 46% of massive galaxies ($\log$ M$_*$/M$_\odot >$ 10), with similar incidence rates in star-forming and quenching systems. Half of the absorption profiles are blueshifted by at least 100 km/s, providing unambiguous evidence for neutral gas outflows. Galaxies with strong Na I D absorption are distinguished by enhanced emission line ratios consistent with AGN ionization. We conservatively measure mass outflow rates of 3 - 100 $M_\odot$ yr$^{-1}$; comparable to or exceeding ionized gas outflow rates measured for galaxies at similar stellar mass and redshift. The outflows from the quenching systems (log(sSFR)[yr$^{-1}$] $\lesssim$ -10) have mass loading factors of 4 - 360, and the energy and momentum outflow rates exceed the expected injection rates from supernova explosions, suggesting that these galaxies could possibly be caught in a rapid blowout phase powered by the AGN. Our findings suggest that AGN-driven ejection of cold gas may be a dominant mechanism for fast quenching of star formation at z~2. △ Less

Submitted 30 January, 2024; v1 submitted 27 October, 2023; originally announced October 2023.

Comments: 17 pages, 8 figures, accepted for publication in MNRAS

Detection of Accretion Shelves Out to the Virial Radius of a Low-Mass Galaxy with JWST

Authors: Charlie Conroy, Benjamin D. Johnson, Pieter van Dokkum, Alis Deason, Sandro Tacchella, Sirio Belli, William P. Bowman, Rohan P. Naidu, Minjung Park, Roberto Abraham, Razieh Emami

Abstract: We report the serendipitous discovery of an extended stellar halo surrounding the low-mass galaxy Ark 227 ($M_\ast=5\times10^9 M_\odot$; d=35 Mpc) in deep JWST NIRCam imaging from the Blue Jay Survey. The F200W-F444W color provides robust star-galaxy separation, enabling the identification of stars at very low density. By combining resolved stars at large galactocentric distances with diffuse emis… ▽ More We report the serendipitous discovery of an extended stellar halo surrounding the low-mass galaxy Ark 227 ($M_\ast=5\times10^9 M_\odot$; d=35 Mpc) in deep JWST NIRCam imaging from the Blue Jay Survey. The F200W-F444W color provides robust star-galaxy separation, enabling the identification of stars at very low density. By combining resolved stars at large galactocentric distances with diffuse emission from NIRCam and Dragonfly imaging at smaller distances, we trace the surface brightness and color profiles of this galaxy over the entire extent of its predicted dark matter halo, from 0.1-100 kpc. Controlled N-body simulations have predicted that minor mergers create "accretion shelves" in the surface brightness profile at large radius. We observe such a feature in Ark 227 at 10-20 kpc, which, according to models, could be caused by a merger with total mass ratio 1:10. The metallicity declines over this radial range, further supporting the minor merger scenario. There is tentative evidence of a second shelf at $μ_V\approx 35$ mag arcsec$^{-2}$ extending from 50-100 kpc, along with a corresponding drop in metallicity. The stellar mass in this outermost envelope is $\approx10^7M_\odot$. These results suggest that Ark 227 experienced multiple mergers with a spectrum of lower-mass galaxies -- a scenario that is broadly consistent with the hierarchical growth of structure in a cold dark matter-dominated universe. Finally, we identify an ultra-faint dwarf associated with Ark 227 with $M_\ast\approx10^5 M_\odot$ and $μ_{V,e}=28.1$ mag arcsec$^{-2}$, demonstrating that JWST is capable of detecting very low-mass dwarfs to distances of at least ~30 Mpc. △ Less

Submitted 19 October, 2023; originally announced October 2023.

Comments: 10 pages, 9 figures, submitted to ApJ

A Search for exoplanets around northern circumpolar stars VIII. filter out a planet cycle from the multi-period radial velocity variations in M giant HD 3638

Authors: Byeong-Cheol Lee, Gwanghui Jeong, Jae-Rim Koo, Beomdu Lim, Myeong-Gu Park, Tae-Yang Bang, Yeon-Ho Choi, Hyeong-Ill Oh, Inwoo Han

Abstract: This paper is written as a follow-up observations to reinterpret the radial velocity (RV) of HD 36384, where the existence of planetary systems is known to be ambiguous. In giants, it is, in general, difficult to distinguish the signals of planetary companions from those of stellar activities. Thus, known exoplanetary giant hosts are relatively rare. We, for many years, have obtained RV data in ev… ▽ More This paper is written as a follow-up observations to reinterpret the radial velocity (RV) of HD 36384, where the existence of planetary systems is known to be ambiguous. In giants, it is, in general, difficult to distinguish the signals of planetary companions from those of stellar activities. Thus, known exoplanetary giant hosts are relatively rare. We, for many years, have obtained RV data in evolved stars using the high-resolution, fiber-fed Bohyunsan Observatory Echelle Spectrograph (BOES) at the Bohyunsan Optical Astronomy Observatory (BOAO). Here, we report the results of RV variations in the M giant HD 36384. We have found two significant periods of 586d and 490d. Considering the orbital stability, it is impossible to have two planets at so close orbits. To determine the nature of the RV variability variations, we analyze the HIPPARCOS photometric data, some indicators of stellar activities, and line profiles. A significant period of 580d was revealed in the HIPPARCOS photometry. Hα EW variations also show a meaningful period of 582d. Thus, the period of 586d may be closely related to the rotational modulations and/or stellar pulsations. On the other hand, the other significant period of 490d is interpreted as the result of the orbiting companion. Our orbital fit suggests that the companion was a planetary mass of 6.6 MJ and is located at 1.3 AU from the host. △ Less

Submitted 11 August, 2023; originally announced August 2023.

Comments: 5 pages, 2 figures, Accepted for publication in JKAS

Star Formation Shut Down by Multiphase Gas Outflow in a Galaxy at a Redshift of 2.45

Authors: Sirio Belli, Minjung Park, Rebecca L. Davies, J. Trevor Mendel, Benjamin D. Johnson, Charlie Conroy, Chloë Benton, Letizia Bugiani, Razieh Emami, Joel Leja, Yijia Li, Gabriel Maheson, Elijah P. Mathews, Rohan P. Naidu, Erica J. Nelson, Sandro Tacchella, Bryan A. Terrazas, Rainer Weinberger

Abstract: Large-scale outflows driven by supermassive black holes are thought to play a fundamental role in suppressing star formation in massive galaxies. However, direct observational evidence for this hypothesis is still lacking, particularly in the young universe where star formation quenching is remarkably rapid, thus requiring effective removal of gas as opposed to slow gas heating. While outflows of… ▽ More Large-scale outflows driven by supermassive black holes are thought to play a fundamental role in suppressing star formation in massive galaxies. However, direct observational evidence for this hypothesis is still lacking, particularly in the young universe where star formation quenching is remarkably rapid, thus requiring effective removal of gas as opposed to slow gas heating. While outflows of ionized gas are commonly detected in massive distant galaxies, the amount of ejected mass is too small to be able to suppress star formation. Gas ejection is expected to be more efficient in the neutral and molecular phases, but at high redshift these have only been observed in starbursts and quasars. Here we report JWST spectroscopy of a massive galaxy experiencing rapid quenching at redshift z=2.445. We detect a weak outflow of ionized gas and a powerful outflow of neutral gas, with a mass outflow rate that is sufficient to quench the star formation. Neither X-ray or radio activity are detected; however, the presence of a supermassive black hole is suggested by the properties of the ionized gas emission lines. We thus conclude that supermassive black holes are able to rapidly suppress star formation in massive galaxies by efficiently ejecting neutral gas. △ Less

Submitted 10 May, 2024; v1 submitted 10 August, 2023; originally announced August 2023.

Comments: 21 pages, 9 figures. Published in Nature

On the Significance of the Thick Disks of Disk Galaxies

Authors: Sukyoung K. Yi, J. K. Jang, Julien Devriendt, Yohan Dubois, San Han, Taysun Kimm, Katarina Kraljic, Minjung Park, Sebastien Peirani, Christophe Pichon, Jinsu Rhee

Abstract: Thick disks are a prevalent feature observed in numerous disk galaxies including our own Milky Way. Their significance has been reported to vary widely, ranging from a few to 100% of the disk mass, depending on the galaxy and the measurement method. We use the NewHorizon simulation which has high spatial and stellar mass resolutions to investigate the issue of thick disk mass fraction. We also use… ▽ More Thick disks are a prevalent feature observed in numerous disk galaxies including our own Milky Way. Their significance has been reported to vary widely, ranging from a few to 100% of the disk mass, depending on the galaxy and the measurement method. We use the NewHorizon simulation which has high spatial and stellar mass resolutions to investigate the issue of thick disk mass fraction. We also use the NewHorizon2 simulation that was run on the same initial conditions but additionally traced nine chemical elements. Based on a sample of 27 massive disk galaxies with M* > 10^10 M_{\odot} in NewHorizon, the contribution of the thick disk was found to be 34 \pm 15% in r-band luminosity or 48 \pm 13% in mass to the overall galactic disk, which seems in agreement with observational data. The vertical profiles of 0, 22, and 5 galaxies are best fitted by 1, 2, or 3 sech2 components, respectively. The NewHorizon2 data show that the selection of thick disk stars based on a single [α/Fe] cut is severely contaminated by stars of different kinematic properties while missing a bulk of kinematically thick disk stars. Vertical luminosity profile fits recover the key properties of thick disks reasonably well. The majority of stars are born near the galactic mid-plane with high circularity and get heated with time via fluctuation in the force field. Depending on the star formation and merger histories, galaxies may naturally develop thick disks with significantly different properties. △ Less

Submitted 7 August, 2023; originally announced August 2023.

Comments: 28 pages, 24 figures. Submitted to ApJS

Long-period radial velocity variations of nine M red giants: The detection of sub-stellar companions around HD 6860 and HD 112300

Authors: Byeong-Cheol Lee, Hee-Jin Do, Myeong-Gu Park, Beomdu Lim, Yeon-Ho Choi, Jae-Rim Koo, Tae-Yang Bang, Hyeong-Ill Oh, Inwoo Han, Heon-Young Chang

Abstract: Context. Certain periodic variations of radial velocities (RV) of wobbling giants originate from exoplanets. Indeed, a number of exoplanets have been discovered around giant stars. Aims. The purpose of our study is to find low-amplitude and long-period RV variations around bright M (super) giants in the RGB (or AGB) stage, which are long-period variables (LPVs) or high-proper-motion (HPM) stars.… ▽ More Context. Certain periodic variations of radial velocities (RV) of wobbling giants originate from exoplanets. Indeed, a number of exoplanets have been discovered around giant stars. Aims. The purpose of our study is to find low-amplitude and long-period RV variations around bright M (super) giants in the RGB (or AGB) stage, which are long-period variables (LPVs) or high-proper-motion (HPM) stars. Methods. High-resolution, fiber-fed Bohyunsan Observatory Echelle Spectrograph (BOES) at the Bohyunsan Optical Astronomy Observatory (BOAO) was used to record numerous spectra of nine giants. The observation period for the targets spans 16 years, from 2005 to 2022. Results. We found from the precise RV observations of nine M giants two sub-stellar companions, one with a 28.26$^{+2.05}_{-2.17}$ $M_{J}$ orbiting period of 663.87$^{+4.61}_{-4.31}$ days at a distance of 2.03$^{+0.01}_{-0.01}$ AU (HD 6860) and the other, with a 15.83$^{+2.33}_{-2.74}$ $M_{J}$ orbiting period of 466.63 $^{+1.47}_{-1.28}$ days at a distance of 1.33 $^{+0.08}_{-0.11}$ AU (HD 112300). Our estimate of the stellar parameters for HD 6860 makes it currently the largest star with a sub-stellar companion. We also found RV variations mimicking a planetary companion in HD 18884 and confirmed LPVs in two stars, HD 39801 and HD 42995. The RV variations of some stars seem to be associated with stellar activities rather than reflex orbital motion due to their companions. Such variations are also detected even for HD 6860 and HD 112300, hosting sub-stellar companions. △ Less

Submitted 29 July, 2023; originally announced July 2023.

Comments: 19 pages, 15 figures, Accepted for publication in A&A

The First Quiescent Galaxies in TNG300

Authors: Abigail I. Hartley, Erica J. Nelson, Katherine A. Suess, Alex M. Garcia, Minjung Park, Lars Hernquist, Rachel Bezanson, Rebecca Nevin, Annalisa Pillepich, Aimee L. Schechter, Bryan A. Terrazas, Paul Torrey, Sarah Wellons, Katherine E. Whitaker, Christina C. Williams

Abstract: We identify the first quiescent galaxies in TNG300, the largest volume of the IllustrisTNG cosmological simulation suite, and explore their quenching processes and time evolution to z=0. We find that the first quiescent galaxies with stellar masses M_* > 3 x 10^{10} M_sun and specific star formation rates sSFR < 10^{-11} yr^{-1} emerge at z~4.2 in TNG300. Suppression of star formation in these gal… ▽ More We identify the first quiescent galaxies in TNG300, the largest volume of the IllustrisTNG cosmological simulation suite, and explore their quenching processes and time evolution to z=0. We find that the first quiescent galaxies with stellar masses M_* > 3 x 10^{10} M_sun and specific star formation rates sSFR < 10^{-11} yr^{-1} emerge at z~4.2 in TNG300. Suppression of star formation in these galaxies begins with a thermal mode of AGN feedback at z~6, and a kinetic feedback mode acts in each galaxy by z~4.7 to complete the quenching process, which occurs on a time-scale of ~0.35 Gyr. Surprisingly, we find that the majority of these galaxies are not the main progenitors of their z=0 descendants; instead, four of the five galaxies fall into more massive galaxies in subsequent mergers at a range of redshifts 2.5 < z < 0.2. By z=0, these descendants are the centres of galaxy clusters with average stellar masses of 8 x 10^{11} M_sun. We make predictions for the first quenched galaxies to be located by the James Webb Space Telescope (JWST). △ Less

Submitted 18 April, 2023; originally announced April 2023.

Comments: 6 pages, 4 figures

A search for exoplanets around north circumpolar stars. VII. Detection of planetary companion orbiting the largest host star HD 18438

Authors: Byeong-Cheol Lee, Jae-Rim Koo, Gwanghui Jeong, Myeong-Gu Park, Inwoo Han, Yeon-Ho Choi

Abstract: We have been conducting a exoplanet search survey using Bohyunsan Observatory Echelle Spectrograph (BOES) for the last 18 years. We present the detection of exoplanet candidate in orbit around HD 18438 from high-precision radial velocity (RV) mesurements. The target was already reported in 2018 (Bang et al. 2018). They conclude that the RV variations with a period of 719 days are likely to be caus… ▽ More We have been conducting a exoplanet search survey using Bohyunsan Observatory Echelle Spectrograph (BOES) for the last 18 years. We present the detection of exoplanet candidate in orbit around HD 18438 from high-precision radial velocity (RV) mesurements. The target was already reported in 2018 (Bang et al. 2018). They conclude that the RV variations with a period of 719 days are likely to be caused by the pulsations because the Lomb-Scargle periodogram of HIPPARCOS photometric and Ha EW variations for HD 18438 show peaks with periods close to that of RV variations and there were no correlations between bisectors and RV measurements. However, the data were not sufficient to reach a firm conclusion. We obtained more RV data for four years. The longer time baseline yields a more accurate determination with a revised period of 803 +/- 5 days and the planetary origin of RV variations with a minimum planetary companion mass of 21 +/- 1 MJup. Our current estimate of the stellar parameters for HD 18438 makes it currently the largest star with a planetary companion. △ Less

Submitted 15 March, 2023; originally announced March 2023.

Comments: 6 pages, 4 figures, Accept to the Journal of the Korean Astronomical Society

Primordial cosmic complexity and effects of reheating

Authors: Pankaj Saha, Myeonghun Park

Abstract: We study the effects of the reheating phase on the evolution of complexities for the primordial curvature perturbation using the squeezed formalism. We examine the evolution of the out-of-time correlator, the quantum discord, and circuit complexity, starting from the inflationary epoch to the radiation-dominated epoch with different reheating scenarios. We find that for a mode that reenters the ho… ▽ More We study the effects of the reheating phase on the evolution of complexities for the primordial curvature perturbation using the squeezed formalism. We examine the evolution of the out-of-time correlator, the quantum discord, and circuit complexity, starting from the inflationary epoch to the radiation-dominated epoch with different reheating scenarios. We find that for a mode that reenters the horizon after reheating, the effect of a finite reheating epoch on the characteristic \textit{freeze-in} amplitude of these primordial complexities can only be distinguished up to three different classes depending on whether the equation of state parameter: $(i)$ $w_\mathrm{re}=1/3$ $(ii)$ $w_\mathrm{re}<1/3$, or, (iii) $w_\mathrm{re}>1/3$. For reheating with different EOS within these classes, the final amplitude will be the same -- hence, the detailed signature of reheating with a class on the complexity measures will be lost. Taking the central value of the scalar spectral index ($n_s=0.9649$) from Planck and the equation of state during reheating $w_\mathrm{re}=0.25$ as benchmark values, we found that the behavior of the complexities for all modes smaller than $1.27\times10^{16}\mathrm{Mpc^{-1}}$ can be classified as above. However, for the small-scale modes reentering the horizon during reheating, the signature of EOS on the evolution of these two complexities will be embedded in each of the cases separately. △ Less

Submitted 21 October, 2023; v1 submitted 28 December, 2022; originally announced December 2022.

Comments: 17 pages, 8 figures; replaced to match published version

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

A New Test of Dynamical Dark Energy Models and Cosmic Tensions in Hořava Gravity

Authors: Eleonora Di Valentino, Nils A. Nilsson, Mu-In Park

Abstract: Horava gravity has been proposed as a renormalizable, higher-derivative, Lorentz-violating quantum gravity model without ghost problems. A Horava gravity based dark energy (HDE) model for dynamical dark energy has been also proposed earlier by identifying all the extra (gravitational) contributions from the Lorentz-violating terms as an effective energy-momentum tensor in Einstein equation. We con… ▽ More Horava gravity has been proposed as a renormalizable, higher-derivative, Lorentz-violating quantum gravity model without ghost problems. A Horava gravity based dark energy (HDE) model for dynamical dark energy has been also proposed earlier by identifying all the extra (gravitational) contributions from the Lorentz-violating terms as an effective energy-momentum tensor in Einstein equation. We consider a complete CMB, BAO, and SNe Ia data test of the HDE model by considering general perturbations over the background perfect HDE fluid. Except from BAO, we obtain the preference of non-flat universes for all other data-set combinations. We obtain a positive result on the cosmic tensions between the Hubble constant H0 and the cosmic shear S8, because we have a shift of H0 towards a higher value, though not enough for resolving the H0 tension, but the value of S8 is unaltered. This is in contrast to a rather decreasing H0 but increasing S8 in a non-flat LCDM. For all other parameters, like Omega_m and Omega_Lambda, we obtain quite comparable results with those of LCDM for all data sets, especially with BAO, so that our results are close to a cosmic concordance between the datasets, contrary to the standard non-flat LCDM. We also obtain some undesirable features, like an almost null result on Omegak, which gives back the flat LCDM, if we do not predetermine the sign of Omegak, but we propose several promising ways for improvements by generalizing our analysis. △ Less

Submitted 1 January, 2023; v1 submitted 15 December, 2022; originally announced December 2022.

Comments: 22 pages, 10 figures, Added comments on (1) projectable vs. non-projectable models (footnotes No.1, 7) (2) RG flows (Discussion No.4), Added references, Corrected typos (Eqs. (10),(11)), Accepted in MNRAS

Rapid Quenching of Galaxies at Cosmic Noon

Authors: Minjung Park, Sirio Belli, Charlie Conroy, Sandro Tacchella, Joel Leja, Sam E. Cutler, Benjamin D. Johnson, Erica J. Nelson, Razieh Emami

Abstract: The existence of massive quiescent galaxies at high redshift seems to require rapid quenching, but it is unclear whether all quiescent galaxies have gone through this phase and what physical mechanisms are involved. To study rapid quenching, we use rest-frame colors to select 12 young quiescent galaxies at $z \sim 1.5$. From spectral energy distribution fitting, we find that they all experienced i… ▽ More The existence of massive quiescent galaxies at high redshift seems to require rapid quenching, but it is unclear whether all quiescent galaxies have gone through this phase and what physical mechanisms are involved. To study rapid quenching, we use rest-frame colors to select 12 young quiescent galaxies at $z \sim 1.5$. From spectral energy distribution fitting, we find that they all experienced intense starbursts prior to rapid quenching. We confirm this with deep Magellan/FIRE spectroscopic observations for a subset of seven galaxies. Broad emission lines are detected for two galaxies and are most likely caused by AGN activity. The other five galaxies do not show any emission features, suggesting that gas has already been removed or depleted. Most of the rapidly quenched galaxies are more compact than normal quiescent galaxies, providing evidence for a central starburst in the recent past. We estimate an average transition time of $300\,\rm Myr$ for the rapid quenching phase. Approximately $4\%$ of quiescent galaxies at $z=1.5$ have gone through rapid quenching; this fraction increases to $23\%$ at $z=2.2$. We identify analogs in the TNG100 simulation and find that rapid quenching for these galaxies is driven by AGN, and for half of the cases, gas-rich major mergers seem to trigger the starburst. We conclude that these massive quiescent galaxies are not just rapidly quenched but also rapidly formed through a major starburst. We speculate that mergers drive gas inflow towards the central regions and grow supermassive black holes, leading to rapid quenching by AGN feedback. △ Less

Submitted 7 October, 2022; originally announced October 2022.

Comments: Submitted to ApJ. Comments are welcome

The Atacama Cosmology Telescope: The Persistence of Neutrino Self-Interaction in Cosmological Measurements

Authors: Christina D. Kreisch, Minsu Park, Erminia Calabrese, Francis-Yan Cyr-Racine, Rui An, J. Richard Bond, Olivier Dore, Jo Dunkley, Patricio Gallardo, Vera Gluscevic, J. Colin Hill, Adam D. Hincks, Mathew S. Madhavacheril, Jeff McMahon, Kavilan Moodley, Thomas W. Morris, Federico Nati, Lyman A. Page, Bruce Partridge, Maria Salatino, Cristobal Sifon, David N. Spergel, Cristian Vargas, Edward J. Wollack

Abstract: We use data from the Atacama Cosmology Telescope (ACT) DR4 to search for the presence of neutrino self-interaction in the cosmic microwave background. Consistent with prior works, the posterior distributions we find are bimodal, with one mode consistent with $Λ$CDM and one where neutrinos strongly self-interact. By combining ACT data with large-scale information from WMAP, we find that a delayed o… ▽ More We use data from the Atacama Cosmology Telescope (ACT) DR4 to search for the presence of neutrino self-interaction in the cosmic microwave background. Consistent with prior works, the posterior distributions we find are bimodal, with one mode consistent with $Λ$CDM and one where neutrinos strongly self-interact. By combining ACT data with large-scale information from WMAP, we find that a delayed onset of neutrino free streaming caused by significantly strong neutrino self-interaction is compatible with these data at the $2-3σ$ level. As seen in the past, the preference shifts to $Λ$CDM with the inclusion of Planck data. We determine that the preference for strong neutrino self-interaction is largely driven by angular scales corresponding to $700 \lesssim \ell \lesssim 1000$ in the ACT E-mode polarization data. This region is expected to be key to discriminate between neutrino self-interacting modes and will soon be probed with more sensitive data. △ Less

Submitted 7 August, 2022; v1 submitted 7 July, 2022; originally announced July 2022.

Comments: added additional references, 9+10 pages, 4+8 figures

Beyond dark energy Fisher forecasts: how DESI will constrain LCDM and quintessence models

Authors: Samuel Goldstein, Minsu Park, Marco Raveri, Bhuvnesh Jain, Lado Samushia

Abstract: We baseline with current cosmological observations to forecast the power of the Dark Energy Spectroscopic Instrument (DESI) in two ways: 1. the gain in constraining power of parameter combinations in the standard $Λ$CDM model, and 2. the reconstruction of quintessence models of dark energy. For the former task we use a recently developed formalism to extract the leading parameter combinations cons… ▽ More We baseline with current cosmological observations to forecast the power of the Dark Energy Spectroscopic Instrument (DESI) in two ways: 1. the gain in constraining power of parameter combinations in the standard $Λ$CDM model, and 2. the reconstruction of quintessence models of dark energy. For the former task we use a recently developed formalism to extract the leading parameter combinations constrained by different combinations of cosmological survey data. For the latter, we perform a non-parametric reconstruction of quintessence using the Effective Field Theory of Dark Energy. Using mock DESI observations of the Hubble parameter, angular diameter distance, and growth rate, we find that DESI will provide significant improvements over current datasets on $Λ$CDM and quintessence constraints. Including DESI mocks in our $Λ$CDM analysis improves constraints on $Ω_m$, $H_0$, and $σ_8$ by a factor of two, where the improvement results almost entirely from the angular diameter distance and growth of structure measurements. Our quintessence reconstruction suggests that DESI will considerably improve constraints on a range of quintessence properties, such as the reconstructed potential, scalar field excursion, and the dark energy equation of state. The angular diameter distance measurements are particularly constraining in the presence of a non-$Λ$CDM signal in which the potential cannot be accounted for by shifts in $H_0$ and $Ω_m$. △ Less

Submitted 4 July, 2022; originally announced July 2022.

No Evolution in the Half-mass Radius of Milky Way-type Galaxies over the Last 10 Gyr

Authors: Maryam Hasheminia, Moein Mosleh, Sandro Tacchella, S. Zahra Hosseini-ShahiSavandi, Minjung Park, Rohan P. Naidu

Abstract: The Milky Way (MW) galaxy is in focus, thanks to new observational data. Here we shed new light on the MW's past by studying the structural evolution of MW progenitors, which we identify from extragalactic surveys. Specifically, we constrain the stellar-mass growth history (SMGH) of the MW with two methods: ($i$) direct measurement of the MW's star formation history, and ($ii$) assuming the MW is… ▽ More The Milky Way (MW) galaxy is in focus, thanks to new observational data. Here we shed new light on the MW's past by studying the structural evolution of MW progenitors, which we identify from extragalactic surveys. Specifically, we constrain the stellar-mass growth history (SMGH) of the MW with two methods: ($i$) direct measurement of the MW's star formation history, and ($ii$) assuming the MW is a typical star-forming galaxy that remains on the star-forming main sequence. We select MW progenitors based on these two SMGHs at $z=0.2-2.0$ from the CANDELS/3D-HST data. We estimate the structural parameters (including half-mass radius $r_{50}$ and Sérsic index) from the stellar-mass profiles. Our key finding is that the progenitors of the MW galaxy grow self-similarly on spatially resolved scales with roughly a constant half-mass radius ($\sim2-3$ kpc) over the past 10 Gyr, while their stellar masses increase by about 1 dex, implying little-to-no inside-out growth. We discover that the radius containing $20\%$ of the stellar mass ($r_{20}$) decreases by $60\%$ between redshifts of $z=2.0$ and $z=0.7$, while the central stellar-mass density ($Σ_1$) increases by a factor of 1.3 dex over the same time and the Sérsic index changes as $n\propto(1+z)^{-1.41\pm0.19}$. This is consistent with an early ($z>1$) formation of a thick disk, followed by the formation of a bar that led to an increase in the mass in the core. The formation and evolution of the thin disk had only little impact on the overall half-mass size. We also show that the constant-size evolution of the MW progenitors challenges semiempirical approaches and numerical simulations. △ Less

Submitted 27 June, 2022; v1 submitted 9 June, 2022; originally announced June 2022.

Comments: Published in ApJL; 8 pages, 5 figures; doi:10.3847/2041-8213/ac76c8

Journal ref: The Astrophysical Journal Letters, 932:L23 (6pp), 2022 June 20

On the formation of massive quiescent galaxies with diverse morphologies in the TNG50 simulation

Authors: Minjung Park, Sandro Tacchella, Erica J. Nelson, Lars Hernquist, Rainer Weinberger, Benedikt Diemer, Dylan Nelson, Annalisa Pillepich, Federico Marinacci, Mark Vogelsberger

Abstract: Observations have shown that the star-formation activity and the morphology of galaxies are closely related, but the underlying physical connection is not well understood. Using the TNG50 simulation, we explore the quenching and the morphological evolution of the 102 massive quiescent galaxies in the mass range of $10.5<\log(M_{\rm stellar}/M_{\odot})<11.5$ selected at $z=0$. The morphology of gal… ▽ More Observations have shown that the star-formation activity and the morphology of galaxies are closely related, but the underlying physical connection is not well understood. Using the TNG50 simulation, we explore the quenching and the morphological evolution of the 102 massive quiescent galaxies in the mass range of $10.5<\log(M_{\rm stellar}/M_{\odot})<11.5$ selected at $z=0$. The morphology of galaxies is quantified based on their kinematics, and we measure the quenching timescale of individual galaxies directly from star formation history. We show that galaxies tend to be quenched more rapidly if they: (i) are satellites in massive halos, (ii) have lower star-forming gas fractions, or (iii) inject a larger amount of black hole kinetic feedback energy. By following the global evolutionary pathways, we conclude that quiescent discs are mainly disc galaxies that are recently and slowly quenched. Approximately half of the quiescent ellipticals at $z=0$ are rapidly quenched at higher redshifts while still disc-like. While being quiescent, they gradually become more elliptical mostly by disc heating, yet these ellipticals still retain some degree of rotation. The other half of quiescent ellipticals with the most random motion-dominated kinematics build up large spheroidal components before quenching primarily by mergers, or in some cases, misaligned gas accretion. However, the mergers that contribute to morphological transformation do not immediately quench galaxies in many cases. In summary, we find that quenching and morphological transformation are decoupled. We conclude that the TNG black hole feedback -- in combination with the stochastic merger history of galaxies -- leads to a large diversity of quenching timescales and a rich morphological landscape. △ Less

Submitted 14 December, 2021; originally announced December 2021.

Comments: Submitted to MNRAS. 9 Figures 17 Pages. Comments are very welcome!

What does a cosmological experiment really measure? Covariant posterior decomposition with normalizing flows

Authors: Tara Dacunha, Marco Raveri, Minsu Park, Cyrille Doux, Bhuvnesh Jain

Abstract: We present methods to rigorously extract parameter combinations that are constrained by data from posterior distributions. The standard approach uses linear methods that apply to Gaussian distributions. We show the limitations of the linear methods for current surveys, and develop non-linear methods that can be used with non-Gaussian distributions, and are independent of the parameter basis. These… ▽ More We present methods to rigorously extract parameter combinations that are constrained by data from posterior distributions. The standard approach uses linear methods that apply to Gaussian distributions. We show the limitations of the linear methods for current surveys, and develop non-linear methods that can be used with non-Gaussian distributions, and are independent of the parameter basis. These are made possible by the use of machine-learning models, normalizing flows, to learn posterior distributions from their samples. These models allow us to obtain the local covariance of the posterior at all positions in parameter space and use its inverse, the Fisher matrix, as a local metric over parameter space. The posterior distribution can then be non-linearly decomposed into the leading constrained parameter combinations via parallel transport in the metric space. We test our methods on two non-Gaussian, benchmark examples, and then apply them to the parameter posteriors of the Dark Energy Survey and Planck CMB lensing. We illustrate how our method automatically learns the survey-specific, best constrained effective amplitude parameter $S_8$ for cosmic shear alone, cosmic shear and galaxy clustering, and CMB lensing. We also identify constrained parameter combinations in the full parameter space, and as an application we estimate the Hubble constant, $H_0$, from large-structure data alone. △ Less

Submitted 10 December, 2021; originally announced December 2021.

Comments: 26 pages, 21 figures

Search for Exoplanets around Northern Circumpolar Stars V. Three likely planetary companions to the giant stars HD 19615, HD 150010, and HD 174205

Authors: G. Jeong, B. C. Lee, M. G. Park, T. Y. Bang, I. Han

Abstract: Aims. We report the detection of long-period radial velocity (RV) variations in three giant stars, HD 19615, HD 150010, and HD 174205, using precise RV measurements. Methods. These detections are part of the Search for Exoplanets around Northern Circumpolar Stars (SENS) survey being conducted at the Bohyunsan Optical Astronomy Observatory (BOAO). The nature of the RV variations was investigated by… ▽ More Aims. We report the detection of long-period radial velocity (RV) variations in three giant stars, HD 19615, HD 150010, and HD 174205, using precise RV measurements. Methods. These detections are part of the Search for Exoplanets around Northern Circumpolar Stars (SENS) survey being conducted at the Bohyunsan Optical Astronomy Observatory (BOAO). The nature of the RV variations was investigated by analyzing the photometric and line shape variations. We found no variability with the RV period in these quantities and conclude that the RV variations are most likely caused by planetary companions. Results. Orbital solutions for the three stars yield orbital periods of 402 d, 562 d, and 582 d and minimum masses of 8.5 MJ , 2.4 MJ , and 4.2 MJ , respectively. These masses and periods are typical for planets around intermediate-mass stars, although some unclear interpretations and recent studies may being calling some planet convictions into question. Nevertheless, the SENS program is contributing to our knowledge of giant planets around intermediate-mass stars. △ Less

Submitted 16 November, 2021; originally announced November 2021.

Comments: 6 pages, 4 figures

Journal ref: A&A 662, A12 (2022)

Properties of Fast and Slow Bars Classified by Epicyclic Frequency Curves from Photometry of Barred Galaxies

Authors: Yun Hee Lee, Myeong-Gu Park, Ho Seong Hwang, Hong Bae Ann, Haeun Chung, Taehyun Kim

Abstract: We test the idea that bar pattern speeds decrease with time owing to angular momentum exchange with a dark matter halo. If this process actually occurs, then the locations of the corotation resonance and other resonances should generally increase with time. We therefore derive the angular velocity $Ω$ and epicyclic frequency $κ$ as functions of galactocentric radius for 85 barred galaxies using ph… ▽ More We test the idea that bar pattern speeds decrease with time owing to angular momentum exchange with a dark matter halo. If this process actually occurs, then the locations of the corotation resonance and other resonances should generally increase with time. We therefore derive the angular velocity $Ω$ and epicyclic frequency $κ$ as functions of galactocentric radius for 85 barred galaxies using photometric data. Mass maps are constructed by assuming a dynamical mass-to-light ratio and then solving the Poisson equation for the gravitatonal potential. The location of Lindblad resonances and the corotation resonance radius are then derived using the standard precession frequency curves in conjunction with bar pattern speeds recently estimated from the Tremaine-Weinberg method as applied to Integral Field Spectroscopy (IFS) data. Correlations between physical properties of bars and their host galaxies indicate that bar {\it length} and the corotation radius depend on the disk circular velocity while bar {\it strength} and pattern speed do not. As the bar pattern speed decreases, bar strength, length, and corotation radius incease, but when bars are subclassified into fast, medium, and slow domains, no significant change in bar length is found. Only a hint of an increase of bar strength from fast to slow bars is found. These results suggest that bar length in galaxies undergoes little evolution, being instead determined mainly by the size of their host galaxy. △ Less

Submitted 16 November, 2021; originally announced November 2021.

Comments: 24 pages, 11 figures

Cosmic Evolution of Barred Galaxies up to z ~ 0.84

Authors: Taehyun Kim, E. Athanassoula, Kartik Sheth, Albert Bosma, Myeong-Gu Park, Yun Hee Lee, Hong Bae Ann

Abstract: We explore the cosmic evolution of the bar length, strength, and light deficit around the bar for 379 barred galaxies at 0.2 < z $\leq$ 0.835 using F814W images from the COSMOS survey. Our sample covers galaxies with stellar mass 10.0 $\leq$ log(M*/Msun) $\leq$ 11.4 and various Hubble types. The bar length is strongly related to the galaxy mass, the disk scale length (h), R50, and R90, where the l… ▽ More We explore the cosmic evolution of the bar length, strength, and light deficit around the bar for 379 barred galaxies at 0.2 < z $\leq$ 0.835 using F814W images from the COSMOS survey. Our sample covers galaxies with stellar mass 10.0 $\leq$ log(M*/Msun) $\leq$ 11.4 and various Hubble types. The bar length is strongly related to the galaxy mass, the disk scale length (h), R50, and R90, where the last two are the radii containing 50 and 90% of total stellar mass, respectively. Bar length remains almost constant, suggesting little or no evolution in bar length over the last 7 Gyrs. The normalized bar lengths (Rbar/h, Rbar/R50, and Rbar/R90) do not show any clear cosmic evolution. Also, the bar strength (A2 and Qb) and the light deficit around the bar reveal little or no cosmic evolution. The constancy of the normalized bar lengths over cosmic time implies that the evolution of bars and of disks is strongly linked over all times. We discuss our results in the framework of predictions from numerical simulations. We conclude there is no strong disagreement between our results and up-to-date simulations. △ Less

Submitted 8 September, 2021; originally announced September 2021.

Comments: 16 pages, 8 figures, Accepted for publication in ApJ

Tests of Standard Cosmology in Horava Gravity and Bayesian Evidence for a Closed Universe, and the Hubble Tension

Authors: Nils A. Nilsson, Mu-In Park

Abstract: We consider some background tests of standard cosmology in the context of Horava gravity with different scaling dimensions for space and time, which has been proposed as a renormalizable, higher-derivative, Lorentz-violating quantum gravity model without ghost problems. We obtain the "very strong" and "strong" Bayesian evidence for our two cosmology models A and B, respectively, depending on the c… ▽ More We consider some background tests of standard cosmology in the context of Horava gravity with different scaling dimensions for space and time, which has been proposed as a renormalizable, higher-derivative, Lorentz-violating quantum gravity model without ghost problems. We obtain the "very strong" and "strong" Bayesian evidence for our two cosmology models A and B, respectively, depending on the choice of parametrization based on Horava gravity, against the standard, spatially-flat, LCDM cosmology model based on general relativity. An MCMC analysis with observational data, including BAO, shows (a) preference of a closed universe with the curvature density parameter Omega_k=-0.005+- 0.0007, -0.004+0.003-0.001 and (b) reduction of the Hubble tension with the Hubble constant H_0=71.4+1.2-0.9, 69.5+1.6-0.9 km s^{-1} Mpc^{-1} for the models A, B. We comment on some possible further improvements for the "cosmic-tension problem" by considering more complete early universe physics, based on the Lorentz-violating standard model with anisotropic space-time scaling, consistently with Horava gravity, as well as the observational data which are properly adopted for the closed universe. △ Less

Submitted 16 September, 2022; v1 submitted 18 August, 2021; originally announced August 2021.

Comments: Changed title, Corrected typos, Added Ref; 23 pages, 6 figures, 2 tables

The NewHorizon Simulation -- To Bar Or Not To Bar

Authors: J. Reddish, K. Kraljic, M. S. Petersen, K. Tep, Y. Dubois, C. Pichon, S. Peirani, F. Bournaud, H. Choi, J. Devriendt, R. Jackson, G. Martin, M. J. Park, M. Volonteri, S. K. Yi

Abstract: We use the NewHorizon simulation to study the redshift evolution of bar properties and fractions within galaxies in the stellar masses range $M_{\star} = 10^{7.25} - 10^{11.4} \ \rm{M}_{\odot}$ over the redshift range $z = 0.25 - 1.3$. We select disc galaxies using stellar kinematics as a proxy for galaxy morphology. We employ two different automated bar detection methods, coupled with visual insp… ▽ More We use the NewHorizon simulation to study the redshift evolution of bar properties and fractions within galaxies in the stellar masses range $M_{\star} = 10^{7.25} - 10^{11.4} \ \rm{M}_{\odot}$ over the redshift range $z = 0.25 - 1.3$. We select disc galaxies using stellar kinematics as a proxy for galaxy morphology. We employ two different automated bar detection methods, coupled with visual inspection, resulting in observable bar fractions of $f_{\rm bar} = 0.070_{-0.012}^{+0.018}$ at $z\sim$ 1.3, decreasing to $f_{\rm bar} = 0.011_{-0.003}^{+0.014}$ at $z\sim$ 0.25. Only one galaxy is visually confirmed as strongly barred in our sample. This bar is hosted by the most massive disk and only survives from $z=1.3$ down to $z=0.7$. Such a low bar fraction, in particular amongst Milky Way-like progenitors, highlights a missing bars problem, shared by literally all cosmological simulations with spatial resolution $<$100 pc to date. The analysis of linear growth rates, rotation curves and derived summary statistics of the stellar, gas and dark matter components suggest that galaxies with stellar masses below $10^{9.5}-10^{10} \ \rm{M}_{\odot}$ in NewHorizon appear to be too dominated by dark matter to form a bar, while more massive galaxies typically have formed large bulges that prevent bar persistence at low redshift. This investigation confirms that the evolution of the bar fraction puts stringent constraints on the assembly history of baryons and dark matter onto galaxies. △ Less

Submitted 24 February, 2022; v1 submitted 4 June, 2021; originally announced June 2021.

Comments: 25 pages, 13 figures, 5 tables, accepted by Monthly Notices of Royal Astronomical Society (MNRAS) on 16/02/22

Reconstructing Quintessence

Authors: Minsu Park, Marco Raveri, Bhuvnesh Jain

Abstract: We present an Effective Field Theory based reconstruction of quintessence models of dark energy directly from cosmological data. We show that current cosmological data possess enough constraining power to test several quintessence model properties for redshifts $z\in [0,1.5]$ with no assumptions about the behavior of the scalar field potential. We use measurements of the cosmic microwave backgroun… ▽ More We present an Effective Field Theory based reconstruction of quintessence models of dark energy directly from cosmological data. We show that current cosmological data possess enough constraining power to test several quintessence model properties for redshifts $z\in [0,1.5]$ with no assumptions about the behavior of the scalar field potential. We use measurements of the cosmic microwave background, supernovae distances, and the clustering and lensing of galaxies to constrain the evolution of the dark energy equation of state, Swampland Conjectures, the shape of the scalar field reconstructed potential, and the structure of its phase space. The standard cosmological model still remains favored by data and, within quintessence models, deviations from its expansion history are bounded to be below the 10% level at 95% confidence at any redshift below $z=1.5$. △ Less

Submitted 12 January, 2021; originally announced January 2021.

Comments: 10 pages, 8 figures

Journal ref: Phys. Rev. D 103, 103530 (2021)

Exploring the origin of thick disks using the NewHorizon and Galactica simulations

Authors: Minjung J. Park, Sukyoung K. Yi, Sebastien Peirani, Christophe Pichon, Yohan Dubois, Hoseung Choi, Julien Devriendt, Sugata Kaviraj, Taysun Kimm, Katarina Kraljic, Marta Volonteri

Abstract: Ever since a thick disk was proposed to explain the vertical distribution of the Milky Way disk stars, its origin has been a recurrent question. We aim to answer this question by inspecting 19 disk galaxies with stellar mass greater than $10^{10}\,\rm M_\odot$ in recent cosmological high-resolution zoom-in simulations: Galactica and NewHorizon. The thin and thick disks are reasonably reproduced by… ▽ More Ever since a thick disk was proposed to explain the vertical distribution of the Milky Way disk stars, its origin has been a recurrent question. We aim to answer this question by inspecting 19 disk galaxies with stellar mass greater than $10^{10}\,\rm M_\odot$ in recent cosmological high-resolution zoom-in simulations: Galactica and NewHorizon. The thin and thick disks are reasonably reproduced by the simulations with scale heights and luminosity ratios as observed. We then spatially classify the thin and thick disks and find that the thick disk stars are older, metal-poorer, kinematically-hotter, and higher in accreted star fraction, while both disks are dominated by the stars formed in situ. Half of the in-situ stars in the thick disks are formed before the galaxies develop their disks, and the rest are formed in spatially and kinematically thinner disks and then thickened with time by heating. However, the 19 galaxies have various properties and evolutionary routes, highlighting the need for statistically-large samples to draw general conclusions. We conclude from our simulations that the thin and thick disk components are not entirely distinct in terms of formation processes, but rather markers of the evolution of galactic disks. Moreover, as the combined result of the thickening of the existing disk stars and the continued formation of young thin-disk stars, the vertical distribution of stars does not change much after the disks settle, pointing to the modulation of both orbital diffusion and star formation by the same confounding factor: the proximity of galaxies to marginal stability. △ Less

Submitted 25 February, 2021; v1 submitted 25 September, 2020; originally announced September 2020.

Comments: Accepted for publication in ApJS

Introducing the NewHorizon simulation: Galaxy properties with resolved internal dynamics across cosmic time

Authors: Yohan Dubois, Ricarda Beckmann, Frédéric Bournaud, Hoseung Choi, Julien Devriendt, Ryan Jackson, Sugata Kaviraj, Taysun Kimm, Katarina Kraljic, Clotilde Laigle, Garreth Martin, Min-Jung Park, Sébastien Peirani, Christophe Pichon, Marta Volonteri, Sukyoung K. Yi

Abstract: Hydrodynamical cosmological simulations are increasing their level of realism by considering more physical processes and having greater resolution or larger statistics. However, usually either the statistical power of such simulations or the resolution reached within galaxies are sacrificed. Here, we introduce the NewHorizon project in which we simulate at high resolution a zoom-in region of… ▽ More Hydrodynamical cosmological simulations are increasing their level of realism by considering more physical processes and having greater resolution or larger statistics. However, usually either the statistical power of such simulations or the resolution reached within galaxies are sacrificed. Here, we introduce the NewHorizon project in which we simulate at high resolution a zoom-in region of $\sim(16\,\rm Mpc)^3$ that is larger than a standard zoom-in region around a single halo and is embedded in a larger box. A resolution of up to 34 pc is reached within galaxies; this allows the simulation to capture the multi-phase nature of the interstellar medium and the clumpy nature of the star formation process in galaxies. In this introductory paper, we present several key fundamental properties of galaxies and their black holes, including the galaxy mass function, cosmic star formation rate, galactic metallicities, the Kennicutt-Schmidt relation, the stellar-to-halo mass relation, galaxy sizes, stellar kinematics and morphology, gas content within galaxies and its kinematics, and the black hole mass and spin properties over time. The various scaling relations are broadly reproduced by NewHorizon with some differences with the standard observables. Owing to its exquisite spatial resolution, NewHorizon captures the inefficient process of star formation in galaxies, which evolve over time from being more turbulent, gas rich, and star bursting at high redshift. These high-redshift galaxies are also more compact, and they are more elliptical and clumpier until the level of internal gas turbulence decays enough to allow for the formation of discs. The NewHorizon simulation gives access to a broad range of galaxy formation and evolution physics at low-to-intermediate stellar masses, which is a regime that will become accessible in the near future through surveys such as the LSST. △ Less

Submitted 28 June, 2021; v1 submitted 22 September, 2020; originally announced September 2020.

Comments: 29 pages, 28 figures, A&A accepted

Rotating Black Holes in Three-Dimensional Horava Gravity Revisited

Authors: Mu-In Park

Abstract: I revisit rotating black hole solutions in three-dimensional Horava gravity with z = 2 as a simpler set-up of the renormalizable quantum gravity `a la Lifshitz and DeWitt. The solutions have a curvature singularity at the origin for a non-vanishing rotation parameter J, unlike the black holes in three-dimensional Einstein gravity. For anti-de Sitter space, there are black hole event horizons as us… ▽ More I revisit rotating black hole solutions in three-dimensional Horava gravity with z = 2 as a simpler set-up of the renormalizable quantum gravity `a la Lifshitz and DeWitt. The solutions have a curvature singularity at the origin for a non-vanishing rotation parameter J, unlike the black holes in three-dimensional Einstein gravity. For anti-de Sitter space, there are black hole event horizons as usual and the singularity is not naked, in agreement with the cosmic censorship. On the other hand, for flat or de Sitter space, the earlier solution has also a cosmic-censorship problem because there are no conventional black hole horizons as in Einstein gravity, other than the usual cosmological horizon for the latter case, so that the singularity could be naked in Horava gravity. However, with the help of recent corrections, I show that the solutions have a peculiar black hole horizon at the origin so that the singularity is not naked even without the conventional black hole horizons in flat or de Sitter case, due to the Lorentz-violating higher-derivative terms. On the other hand, I note also that a new ``cosmological" horizon exists even for the flat case, contrary to the usual wisdom, due to combined effects of the higher derivatives and the angular-momentum barrier. I study an unified treatment of their unusual black hole thermodynamics for the flat and de Sitter spaces, as well as the anti-de Sitter space, which might be due to lack of the absolute horizons in the Lorentz-violating gravity. △ Less

Submitted 20 September, 2022; v1 submitted 14 August, 2020; originally announced August 2020.

Comments: Changed title, Added comments, 16 pages, 7 figures

Bar Classification based on the Potential Map

Authors: Yun Hee Lee, Myeong-Gu Park, Hong Bae Ann, Taehyun Kim, Woo-Young Seo

Abstract: We introduce a new approach to classify barred galaxies that utilizes the transverse-to-radial force ratio map (ratio map, hereafter) in a different manner from previous studies. When we display the ratio map in polar coordinates, barred galaxies appear as four aligned, horizontal thick slabs. This characteristic feature enables us to successfully classify barred and nonbarred galaxies with an acc… ▽ More We introduce a new approach to classify barred galaxies that utilizes the transverse-to-radial force ratio map (ratio map, hereafter) in a different manner from previous studies. When we display the ratio map in polar coordinates, barred galaxies appear as four aligned, horizontal thick slabs. This characteristic feature enables us to successfully classify barred and nonbarred galaxies with an accuracy of 87%. It yields the bar fraction of 53%, including both SBs and SABs, when applied to 884 nearby (z < 0.01) spiral galaxies from the Sloan Digital Sky Survey/DR7. It also provides the bar strength and length measurements, in particular, separated from the spiral arms. They show good correlations with the measures estimated from ellipse fitting and Fourier analysis. However, we find different tendencies of the bar strength measurements in terms of the Hubble sequence: as the Hubble sequence increases (towards late-type), the bar strength and bar ellipticity increase, whereas the dipole Fourier amplitude decreases. We show that the bulge affects the estimation of the bar strength differently, depending on the classification methods. The bulge causes the bar length to be overestimated in all three methods. Meanwhile, we find that barred galaxies show two types of radial profiles of the angle-averaged force ratio: one has a maximum peak (type M) and the other a plateau (type P). Comparison with numerical simulations suggests that type M bars are more mature than type P bars in terms of evolutionary stage. △ Less

Submitted 11 July, 2020; v1 submitted 8 July, 2020; originally announced July 2020.

Comments: 29 pages, 21 figures

Hybrid star HD 81817 accompanied by brown dwarf or substellar companion

Authors: Tae-Yang Bang, Byeong-Cheol Lee, V. Perdelwitz, Gwanghui Jeong, Inwoo Han, Hyeong-Il Oh, Myeong-Gu Park

Abstract: HD 81817 is known as a hybrid star. Hybrid stars have both cool stellar wind properties and Ultraviolet (UV) or even X-ray emission features of highly ionized atoms in their spectra. A white dwarf companion has been suggested as the source of UV or X-ray features. HD 81817 has been observed since 2004 as a part of a radial velocity (RV) survey program to search for exoplanets around K giant stars… ▽ More HD 81817 is known as a hybrid star. Hybrid stars have both cool stellar wind properties and Ultraviolet (UV) or even X-ray emission features of highly ionized atoms in their spectra. A white dwarf companion has been suggested as the source of UV or X-ray features. HD 81817 has been observed since 2004 as a part of a radial velocity (RV) survey program to search for exoplanets around K giant stars using the Bohyunsan Observatory Echelle Spectrograph at the 1.8 m telescope of Bohyunsan Optical Astronomy Observatory in Korea. We obtained 85 RV measurements between 2004 and 2019 for HD 81817 and found two periodic RV variations. The amplitudes of RV variations are around 200 m s^-1, which are significantly lower than that expected from a closely orbiting white dwarf companion. Photometric data and relevant spectral lines were also analyzed to help determine the origin of the periodic RV variations. We conclude that 627.4-day RV variations are caused by intrinsic stellar activities such as long-term pulsations or rotational modulations of surface activities based on Hα equivalent width (EW) variations of a similar period. On the other hand, 1047.1-day periodic RV variations are likely to be caused by a brown dwarf or substellar companion, which is corroborated by a recent GAIA proper motion anomaly for HD 81817. The Keplerian fit yields a minimum mass of 27.1 M_Jup, a semimajor axis of 3.3 AU, and an eccentricity of 0.17 for the stellar mass of 4.3 M_sun for HD 81817. The inferred mass puts HD 81817 b in the brown dwarf desert. △ Less

Submitted 19 May, 2020; originally announced May 2020.

Journal ref: A&A 638, A148 (2020)

Stellar Interferometry for Gravitational Waves

Authors: I. H. Park, K. -Y. Choi, J. Hwang, S. Jung, D. H. Kim, M. H. Kim, C. -H. Lee, K. H. Lee, S. H. Oh, M. -G. Park, S. C. Park, A. Pozanenko, C. D. Rho, N. Vedenkin, E. Won

Abstract: We propose a new method to detect gravitational waves, based on spatial coherence interferometry with stellar light, as opposed to the conventional temporal coherence interferometry with laser sources. The proposed method detects gravitational waves by using two coherent beams of light from a single distant star measured at separate space-based detectors with a long baseline. This method can be ap… ▽ More We propose a new method to detect gravitational waves, based on spatial coherence interferometry with stellar light, as opposed to the conventional temporal coherence interferometry with laser sources. The proposed method detects gravitational waves by using two coherent beams of light from a single distant star measured at separate space-based detectors with a long baseline. This method can be applied to either the amplitude or intensity interferometry. This experiment allows for the search of gravitational waves in the lower frequency range of $10^{-6}$ to $10^{-4}$ Hz. In this work, we present the detection sensitivity of the proposed stellar interferometer by taking the detector response and shot and acceleration noises into account. Furthermore, the proposed experimental setup is capable of searching for primordial black holes and studying the size of the target neutron star, which are also discussed in the paper. △ Less

Submitted 6 November, 2021; v1 submitted 14 June, 2019; originally announced June 2019.

Comments: 12 pages, 7 figures

Report number: SKKU-Astro-2019-001, KU-Cosmology-2019-001

Journal ref: Journal of Cosmology and Astroparticle Physics, 11, 008 (2021)

Structural study of analogues of Titan s haze by trapped ion mobility coupled with a Fourier transform ion cyclotron mass spectrometer

Authors: Christopher P. Rueger, Julien Maillard, Johann Le Maitre, Mark Ridgeway, Christopher J. Thompson, Isabelle Schmitz-Afonso, Thomas Gautier, Nathalie Carrasco, Melvin A. Park, Pierre Giusti, Carlos Afonso

Abstract: The aerosols present in the atmosphere of the Saturn s moon Titan are of particular planetary science interest and several spacecraft missions already allowed to gather spectroscopic data. Titan haze s analogs, so-called tholins, were produced on earth to push forward the comprehension of their formation and properties. In this study, this highly complex mixture was analyzed here for the first tim… ▽ More The aerosols present in the atmosphere of the Saturn s moon Titan are of particular planetary science interest and several spacecraft missions already allowed to gather spectroscopic data. Titan haze s analogs, so-called tholins, were produced on earth to push forward the comprehension of their formation and properties. In this study, this highly complex mixture was analyzed here for the first time by trapped ion mobility spectrometry coupled to ultra-high resolution mass spectrometry (FTICR MS). Electrospray ionization revealed the characteristic CHNx-class components, with CHN5-6 and DBE 6-7 most abundant. Deploying specialized visualization, enabled by accurate mass measurements and elemental composition assignments, the adapted Kendrick mass defect analysis highlights the C2H3N homolog series, whereas the nitrogen-modified van Krevelen diagram exhibits a clear trend towards H/C 1.5 and N/C 0.5. More interestingly, the representation of m/z versus collision cross section (CCS) allowed hypothesizing a ramified N-PAH structural motif. State-of-the-art IMS is currently not able to resolve the isomeric continuum of ultra-complex mixtures; thus peak parameter other than the CCS-value are explored. As such, analyzing the mobility peak width versus m/z shows a linear increase in isomeric diversity between m/z 170 and 350 and a near plateau in diversity at higher m/z for the N-PAH-like structure. Due to the high complexity of the sample, these structural insights are only to be revealed by TIMS-FTICR MS. △ Less

Submitted 7 June, 2019; originally announced June 2019.

Comments: Journal of The American Society for Mass Spectrometry, 2019

New Horizon: On the origin of the stellar disk and spheroid of field galaxies at $z=0.7$

Authors: Min-Jung Park, Sukyoung K. Yi, Yohan Dubois, Christophe Pichon, Taysun Kimm, Julien Devriendt, Hoseung Choi, Marta Volonteri, Sugata Kaviraj, Sebastien Peirani

Abstract: The origin of the disk and spheroid of galaxies has been a key open question in understanding their morphology. Using the high-resolution cosmological simulation, New Horizon, we explore kinematically decomposed disk and spheroidal components of 144 field galaxies with masses greater than $\rm 10^9\,M_{\odot}$ at $z=0.7$. The origins of stellar particles are classified according to their birthplac… ▽ More The origin of the disk and spheroid of galaxies has been a key open question in understanding their morphology. Using the high-resolution cosmological simulation, New Horizon, we explore kinematically decomposed disk and spheroidal components of 144 field galaxies with masses greater than $\rm 10^9\,M_{\odot}$ at $z=0.7$. The origins of stellar particles are classified according to their birthplace (in situ or ex situ) and their orbits at birth. Before disk settling, stars form mainly through chaotic mergers between proto-galaxies and become part of the spheroidal component. When disk settling starts, we find that more massive galaxies begin to form disk stars from earlier epochs; massive galaxies commence to develop their disks at $z\sim1-2$, while low-mass galaxies do after $z\sim1$. The formation of disks is affected by accretion as well, as mergers can trigger gas turbulence or induce misaligned gas infall that prevents galaxies from forming co-rotating disk stars. The importance of accreted stars is greater in more massive galaxies, especially in developing massive spheroids. A significant fraction of the spheroids comes from the disk stars that are perturbed, which becomes more important at lower redshifts. Some ($\sim12.5\%$) of our massive galaxies develop counter-rotating disks from the gas infall misaligned with the existing disk plane, which can last for more than a Gyr until they become the dominant component, and flip the angular momentum of the galaxy in the opposite direction. The final disk-to-total ratio of a galaxy needs to be understood in relation to its stellar mass and accretion history. We quantify the significance of the stars with different origins and provide them as guiding values. △ Less

Submitted 13 August, 2019; v1 submitted 6 May, 2019; originally announced May 2019.

Comments: Accepted for publication in ApJ. 23 pages with 13 figures

$Λ$CDM or self-interacting neutrinos? - how CMB data can tell the two models apart

Authors: Minsu Park, Christina D. Kreisch, Jo Dunkley, Boryana Hadzhiyska, Francis-Yan Cyr-Racine

Abstract: Of the many proposed extensions to the $Λ$CDM paradigm, a model in which neutrinos self-interact until close to the epoch of matter-radiation equality has been shown to provide a good fit to current cosmic microwave background (CMB) data, while at the same time alleviating tensions with late-time measurements of the expansion rate and matter fluctuation amplitude. Interestingly, CMB fits to this m… ▽ More Of the many proposed extensions to the $Λ$CDM paradigm, a model in which neutrinos self-interact until close to the epoch of matter-radiation equality has been shown to provide a good fit to current cosmic microwave background (CMB) data, while at the same time alleviating tensions with late-time measurements of the expansion rate and matter fluctuation amplitude. Interestingly, CMB fits to this model either pick out a specific large value of the neutrino interaction strength, or are consistent with the extremely weak neutrino interaction found in $Λ$CDM, resulting in a bimodal posterior distribution for the neutrino self-interaction cross section. In this paper, we explore why current cosmological data select this particular large neutrino self-interaction strength, and by consequence, disfavor intermediate values of the self-interaction cross section. We show how it is the $\ell \gtrsim 1000$ CMB temperature anisotropies, most recently measured by the Planck satellite, that produce this bimodality. We also establish that smaller scale temperature data, and improved polarization data measuring the temperature-polarization cross-correlation, will best constrain the neutrino self-interaction strength. We forecast that the upcoming Simons Observatory should be capable of distinguishing between the models. △ Less

Submitted 22 April, 2019; v1 submitted 4 April, 2019; originally announced April 2019.

Comments: 7 pages, 7 figures, comments welcome, references added, version submitted to PRD

Journal ref: Phys. Rev. D 100, 063524 (2019)

Bar Fraction in Early- and Late-type Spirals

Authors: Yun Hee Lee, Hong Bae Ann, Myeong-Gu Park

Abstract: Bar fractions depend on the properties of the host galaxies, which are closely related to the formation and evolution of bars. However, observational studies do not provide consistent results. We investigate the bar fraction and its dependence on the properties of the host galaxies by using three bar classification methods: visual inspection, ellipse fitting method, and Fourier analysis. Our volum… ▽ More Bar fractions depend on the properties of the host galaxies, which are closely related to the formation and evolution of bars. However, observational studies do not provide consistent results. We investigate the bar fraction and its dependence on the properties of the host galaxies by using three bar classification methods: visual inspection, ellipse fitting method, and Fourier analysis. Our volume-limited sample consists of 1,698 spiral galaxies brighter than $M_{\rm r}=-15.2$ with $z < 0.01$ from SDSS/DR7 visually classified by Ann et al. (2015). We first compare the consistency of classification among the three methods. Automatic classifications detect visually determined strongly barred galaxies with the concordance of $74\% \sim 85\%$. However, they have some difficulty in identifying bars, in particular, in bulge-dominated galaxies, which affects the distribution of bar fraction as a function of the Hubble type. We obtain, for the same sample, different bar fractions of 63%, 48%, and 36% by visual inspection, ellipse fitting, and Fourier analysis, respectively. The difference is mainly due to how many weak bars are included. Moreover, we find the different dependence of bar fraction on the Hubble type for strong versus weak bars: SBs are preponderant in early-type spirals whereas SABs in late-type spirals. This causes a contradictory dependence on host galaxy properties when different classification methods are used. We propose that strong bars and weak bars experience different processes for their formation, growth, and dissipation by interacting with different inner galactic structures of early-type and late-type spirals. △ Less

Submitted 16 January, 2019; originally announced January 2019.

Comments: 23 pages, 12 figures

A Search for Exoplanets around Northern Circumpolar Stars. IV. Six Planet Candidates to the K Giants, HD 44385, HD 97619, HD 106574, HD 118904, HD 164428, and HD 202432

Authors: Gwanghui Jeong, Inwoo Han, Myeong-Gu Park, Artie P. Hatzes, Tae-Yang Bang, Shenghong Gu, Jinming Bai, Byeong-Cheol Lee

Abstract: We report the discovery of long-period radial velocity (RV) variations in six intermediate-mass K-giant stars using precise RV measurements. These discoveries are part of the Search for Exoplanets around Northern Circumpolar Stars (SENS) survey being conducted at the Bohyunsan Optical Astronomy Observatory. The nature of the RV variations was investigated by looking for photometric and line shape… ▽ More We report the discovery of long-period radial velocity (RV) variations in six intermediate-mass K-giant stars using precise RV measurements. These discoveries are part of the Search for Exoplanets around Northern Circumpolar Stars (SENS) survey being conducted at the Bohyunsan Optical Astronomy Observatory. The nature of the RV variations was investigated by looking for photometric and line shape variations. We can find no variability with the RV period in these quantities and conclude that RV variations are most likely due to unseen sub-stellar companions. Orbital solutions for the six stars yield orbital periods in the range 418-1065 days and minimum masses in the range 1.9-8.5 MJ. These properties are typical on planets around intermediate-mass stars. Our SENS survey so far has about an 8% confirmed planet occurrence rate, and it will provide better statistics on planets around giant stars when the survey is completed. △ Less

Submitted 3 August, 2018; originally announced August 2018.

Comments: 14 figures

Search for exoplanets around northern circumpolar stars III. long-period radial velocity variations in hd 18438 and hd 158996

Authors: Tae-Yang Bang, Byeong-Cheol Lee, Gwang-Hui Jeong, Inwoo Han, Myeong-Gu Park

Abstract: Detecting exoplanets around giant stars sheds light on the later-stage evolution of planetary systems. We observed the M giant HD 18438 and the K giant HD 158996 as part of a Search for Exoplanets around Northern circumpolar Stars (SENS) and obtained 38 and 24 spectra from 2010 to 2017 using the high-resolution Bohyunsan Observatory Echelle Spectrograph (BOES) at the 1.8m telescope of Bohyunsan Op… ▽ More Detecting exoplanets around giant stars sheds light on the later-stage evolution of planetary systems. We observed the M giant HD 18438 and the K giant HD 158996 as part of a Search for Exoplanets around Northern circumpolar Stars (SENS) and obtained 38 and 24 spectra from 2010 to 2017 using the high-resolution Bohyunsan Observatory Echelle Spectrograph (BOES) at the 1.8m telescope of Bohyunsan Optical Astronomy Observatory in Korea. We obtained precise RV measurements from the spectra and found long-period radial velocity (RV) variations with period 719.0 days for HD 18438 and 820.2 days for HD 158996. We checked the chromospheric activities using Ca ii H and H_alpha lines, HIPPARCOS photometry and line bisectors to identify the origin of the observed RV variations. In the case of HD 18438, we conclude that the observed RV variations with period 719.0 days are likely to be caused by the pulsations because the periods of HIPPARCOS photometric and H alpha EW variations for HD 18438 are similar to that of RV variations in Lomb-Scargle periodogram, and there are no correlations between bisectors and RV measurements. In the case of HD 158996, on the other hand, we did not find any similarity in the respective periodograms nor any correlation between RV variations and line bisector variations. In addition, the probability that the real rotational period can be as longer than the RV period for HD 158996 is only about 4.3%. Thus we conclude that observed RV variations with a period of 820.2 days of HD 158996 are caused by a planetary companion, which has the minimum mass of 14.0 MJup, the semi-major axis of 2.1 AU, and eccentricity of 0.13 assuming the stellar mass of 1.8 M_sun. HD 158996 is so far one of the brightest and largest stars to harbor an exoplanet candidate. △ Less

Submitted 19 March, 2018; originally announced March 2018.

Comments: 9 pages, 11 figures

Magnetic field and radial velocities of the star Chi Draconis A

Authors: Byeong-Cheol Lee, D. Gadelshin, Inwoo Han, Dong-Il Kang, Kang-Min Kim, G. Valyavin, G. Galazutdinov, Gwanghui Jeong, N. Beskrovnaya, T. Burlakova, A. Grauzhanina, N. R. Ikhsanov, A. F. Kholtygin, A. Valeev, V. Bychkov, Myeong-Gu Park

Abstract: We present high-resolution spectropolarimetric observations of the spectroscopic binary Chi Dra. Spectral lines in the spectrum of the main component Chi Dra A show variable Zeeman displacement, which confrms earlier suggestions about the presence of a weak magnetic field on the surface of this star. Within about 2 years of time base of our observations, the longitudinal component BL of the magnet… ▽ More We present high-resolution spectropolarimetric observations of the spectroscopic binary Chi Dra. Spectral lines in the spectrum of the main component Chi Dra A show variable Zeeman displacement, which confrms earlier suggestions about the presence of a weak magnetic field on the surface of this star. Within about 2 years of time base of our observations, the longitudinal component BL of the magnetic field exhibits variation from -11.5 +/- 2.5 G to +11.1 +/- 2.1 G with a period of about 23 days. Considering the rotational velocity of Chi DraA in the literature and that newly measured in this work, this variability may be explained by the stellar rotation under the assumption that the magnetic field is globally stable. Our new measurements of the radial velocities (RV) in high-resolution I-spectra of Chi Dra A refined the orbital parameters and reveal persistent deviations of RVs from the orbital curve. We suspect that these deviations may be due to the in uence of local magnetically generated spots, pulsations, or a Jupiter-size planet orbiting the system. △ Less

Submitted 25 September, 2017; originally announced September 2017.

Comments: 5pages, 3figures, accepted for publication in Monthly Notices of the Royal Astronomical Society Letters

On Gauge Invariant Cosmological Perturbations in UV-modified Horava Gravity: A Brief Introduction

Authors: Mu-In Park

Abstract: We revisit gauge invariant cosmological perturbations in UV-modified, z = 3 Horava gravity with one scalar matter field, which has been proposed as a renormalizable gravity theory without the ghost problem in four dimensions. We confirm that there is no extra graviton modes and general relativity is recovered in IR, which achieves the consistency of the model. From the UV-modification terms which… ▽ More We revisit gauge invariant cosmological perturbations in UV-modified, z = 3 Horava gravity with one scalar matter field, which has been proposed as a renormalizable gravity theory without the ghost problem in four dimensions. We confirm that there is no extra graviton modes and general relativity is recovered in IR, which achieves the consistency of the model. From the UV-modification terms which break the detailed balance condition in UV, we obtain scale-invariant power spectrums for non-inflationary backgrounds, like the power-law expansions, without knowing the details of early expansion history of Universe. This could provide a new framework for the Big Bang cosmology. △ Less

Submitted 6 September, 2017; originally announced September 2017.

Comments: 5 pages, Prepared for Proceeding of International Joint Conference of ICGAC-XIII and IK15 (July 3-7, 2017, Seoul, Republic of Korea)

Search for Exoplanets around Northern Circumpolar Stars- II. The Detection of Radial Velocity Variations in M Giant Stars HD 36384, HD 52030, and HD 208742

Authors: Byeong-Cheol Lee, Gwanghui Jeong, Myeong-Gu Park, Inwoo Han, David E. Mkrtichian, Artie P. Hatzes, Shenghong Gu, Jinming Bai, Sang-Min Lee, Hyeong-Il Oh, Kang-Min Kim

Abstract: We present the detection of long-period RV variations in HD 36384, HD 52030, and HD 208742 by using the high-resolution, fiber-fed Bohyunsan Observatory Echelle Spectrograph (BOES) for the precise radial velocity (RV) survey of about 200 northern circumpolar stars. Analyses of RV data, chromospheric activity indicators, and bisector variations spanning about five years suggest that the RV variatio… ▽ More We present the detection of long-period RV variations in HD 36384, HD 52030, and HD 208742 by using the high-resolution, fiber-fed Bohyunsan Observatory Echelle Spectrograph (BOES) for the precise radial velocity (RV) survey of about 200 northern circumpolar stars. Analyses of RV data, chromospheric activity indicators, and bisector variations spanning about five years suggest that the RV variations are compatible with planet or brown dwarf companions in Keplerian motion. However, HD 36384 shows photometric variations with a period very close to that of RV variations as well as amplitude variations in the weighted wavelet Z-transform (WWZ) analysis, which argues that the RV variations in HD~36384 are from the stellar pulsations. Assuming that the companion hypothesis is correct, HD~52030 hosts a companion with minimum mass 13.3 M_Jup$ orbiting in 484 days at a distance of 1.2 AU. HD~208742 hosts a companion of 14.0 M_Jup at 1.5 AU with a period of 602 days. All stars are located at the asymptotic giant branch (AGB) stage on the H-R diagram after undergone the helium flash and left the giant clump.With stellar radii of 53.0 R_Sun and 57.2 R_Sun for HD 52030 and HD 208742, respectively, these stars may be the largest yet, in terms of stellar radius, found to host sub-stellar companions. However, given possible RV amplitude variations and the fact that these are highly evolved stars the planet hypothesis is not yet certain. △ Less

Submitted 28 May, 2017; originally announced May 2017.

Comments: accepted for publication in ApJ; 30 pages, 13 figures, 5 tables

On Gauge Invariant Cosmological Perturbations in UV-modified Horava Gravity

Authors: Sunyoung Shin, Mu-In Park

Abstract: We consider gauge invariant cosmological perturbations in UV-modified, z=3 Horava gravity with one scalar matter field, which has been proposed as a renormalizable gravity theory without the ghost problem in four dimensions. In order to exhibit its dynamical degrees of freedom, we consider the Hamiltonian reduction method and find that, by solving "all" the constraint equations, the degrees of fre… ▽ More We consider gauge invariant cosmological perturbations in UV-modified, z=3 Horava gravity with one scalar matter field, which has been proposed as a renormalizable gravity theory without the ghost problem in four dimensions. In order to exhibit its dynamical degrees of freedom, we consider the Hamiltonian reduction method and find that, by solving "all" the constraint equations, the degrees of freedom are the same as those of Einstein gravity: One scalar and two tensor (graviton) modes when a scalar matter field presents. However, we confirm that there is no extra graviton modes and general relativity is recovered in IR, which achieves the consistency of the model. From the UV-modification terms which break the detailed balance condition in UV, we obtain scale-invariant power spectrums for "non"-inflationary backgrounds, like the power-law expansions, without knowing the details of early expansion history of Universe. This could provide a new framework for the Big Bang cosmology. Moreover, we find that "tensor and scalar fluctuations travel differently in UV, generally". We present also some clarifying remarks about confusing points in the literatures. △ Less

Submitted 19 December, 2017; v1 submitted 13 January, 2017; originally announced January 2017.

Comments: Typos corrected, References added, Minor revisions for clarification, Accepted in JCAP

Long-period variations in the radial velocity of spectroscopic binary m giant mu ursae majoris

Authors: Byeong-Cheol Lee, Inwoo Han, Myeong-Gu Park, David E. Mkrtichian, Artie P. Hatzes, Gwanghui Jeong, Kang-Min Kim

Abstract: We report that the spectroscopic binary ? Ursae Majoris (?mu UMa) has secondary RV variations of 471.2 days in addition to those of 230.0 days already known. Keplerian orbit analysis yields stellar mass companions of 1.6 M_Sun for the 230-d period and 0.14 M_Sun for the 471-d period. However, the HIPPARCOS photometries show a period similar to the stellar rotational period, which is one-quarter of… ▽ More We report that the spectroscopic binary ? Ursae Majoris (?mu UMa) has secondary RV variations of 471.2 days in addition to those of 230.0 days already known. Keplerian orbit analysis yields stellar mass companions of 1.6 M_Sun for the 230-d period and 0.14 M_Sun for the 471-d period. However, the HIPPARCOS photometries show a period similar to the stellar rotational period, which is one-quarter of the RV period. Variations in the bisector velocity curvature show a period of 463.6 days. We also find ~?473-day variations in the equivalent width (EW) measurements of the H_alpha and H_beta lines, whose origin is probably stellar activity. We note that the nature of 471-day variations is similar to one observed in "Sequence D" of symptotic Giant Branch (AGB) pulsating stars. We therefore conclude that the RV and the EW variations in the spectroscopic binary M giant mu UMa A originate from the complex pulsations and the chromospheric activity. △ Less

Submitted 22 February, 2016; originally announced February 2016.

Comments: 20 pages, 10 figures, 3 tables, Accepted for publication in The Astronomical Journal

Search for exoplanet around northern circumpolar stars - Four planets around HD 11755, HD 12648, HD 24064, and 8 Ursae Minoris

Authors: B. -C. Lee, M. -G. Park, S. -M. Lee, G. Jeong, H. -I. Oh, I. Han, J. W. Lee, C. -U. Lee, S. -L. Kim, K. -M. Kim

Abstract: Aims. This program originated as the north pole region extension of the established exoplanet survey using 1.8 m telescope at Bohyunsan Optical Astronomy Observatory (BOAO). The aim of our paper is to find exoplanets in northern circumpolar stars with a precise radial velocity (RV) survey. Methods. We have selected about 200 northern circumpolar stars with the following criteria: Dec. > 70 degree,… ▽ More Aims. This program originated as the north pole region extension of the established exoplanet survey using 1.8 m telescope at Bohyunsan Optical Astronomy Observatory (BOAO). The aim of our paper is to find exoplanets in northern circumpolar stars with a precise radial velocity (RV) survey. Methods. We have selected about 200 northern circumpolar stars with the following criteria: Dec. > 70 degree, 0.6 < B-V < 1.6, HIPPARCOS_scat < 0.05 magnitude, and 5.0 < mv < 7.0. The high-resolution, fiber-fed Bohyunsan Observatory Echelle Spectrograph (BOES) was used for the RV survey. Chromospheric activities, the HIPPARCOS photometry, and line bisectors were analyzed to exclude other causes for the RV variations. Results. In 2010, we started to monitor the candidates and have completed initial screening for all stars for the last five years. We present the detection of four new exoplanets. Stars HD 11755, HD 12648, HD 24064, and 8 UMi all show evidence for giant planets in Keplerian motion. The companion to HD 11755 has a minimum mass of 6.5 M_Jup in a 433-day orbit with an eccentricity of 0.19. HD 12648 is orbited by a companion of minimum mass of 2.9 M_Jup having a period of 133 days and an eccentricity of 0.04. Weak surface activity was suspected in HD 24064. However, no evidence was found to be associated with the RV variations. Its companion has a minimum mass of 9.4 M_Jup, a period of 535 days, and an eccentricity of 0.35. Finally, 8 UMi has a minimum mass of 1.5 M_Jup, a period of 93 days with an eccentricity of 0.06. △ Less

Submitted 30 September, 2015; originally announced September 2015.

Comments: 9 pages, 14 figures, 6 tables, accepted for publisation in Astronomy & Astrophysics

Journal ref: A&A 584, A79 (2015)

Searches for millisecond pulsar candidates among the unidentified Fermi objects

Authors: C. Y. Hui, S. M. Park, C. P. Hu, L. C. C. Lin, K. L. Li, A. K. H. Kong, P. H. T. Tam, J. Takata, K. S. Cheng, Ruolan Jin, T. -C. Yen, Chunglee Kim

Abstract: Here we report the results of searching millisecond pulsar (MSP) candidates from the Fermi LAT second source catalog (2FGL). Seven unassociated $γ-$ray sources in this catalog are identified as promising MSP candidates based on their $γ$-ray properties. Through the X-ray analysis, we have detected possible X-ray counterparts, localized to an arcsecond accuracy. We have systematically estimated the… ▽ More Here we report the results of searching millisecond pulsar (MSP) candidates from the Fermi LAT second source catalog (2FGL). Seven unassociated $γ-$ray sources in this catalog are identified as promising MSP candidates based on their $γ$-ray properties. Through the X-ray analysis, we have detected possible X-ray counterparts, localized to an arcsecond accuracy. We have systematically estimated their X-ray fluxes and compared with the corresponding $γ$-ray fluxes. The X-ray to $γ$-ray flux ratios for 2FGL J1653.6-0159 and 2FGL J1946.4-5402 are comparable with the typical value for pulsars. For 2FGL J1625.2-0020, 2FGL J1653.6-0159 and 2FGL J1946.4-5402, their candidate X-ray counterparts are bright enough for performing a detailed spectral and temporal analysis to discriminate their thermal/non thermal nature and search for the periodic signal. We have also searched for possible optical/IR counterparts at the X-ray positions. For the optical/IR source coincident with the brightest X-ray object that associated with 2FGL J1120.0-2204, its spectral energy distribution is comparable with a late-type star. Evidence for the variability has also been found by examining its optical light curve. All the aforementioned 2FGL sources resemble a pulsar in one or more aspects, which make them as the promising targets for follow-up investigations. △ Less

Submitted 9 July, 2015; originally announced July 2015.

Comments: 37 pages, 12 figures, 6 tables, accepted for publication in ApJ