Comprehensive Broadband X-Ray and Multiwavelength Study of Active Galactic Nuclei in 57 Local Luminous and Ultraluminous Infrared Galaxies Observed with NuSTAR and/or Swift/BAT

We perform a systematic X-ray spectroscopic analysis of 57 local luminous and ultraluminous infrared galaxy systems (containing 84 individual galaxies) observed with the Nuclear Spectroscopic Telescope Array and/or Swift/BAT. Combining soft X-ray data obtained with Chandra, XMM-Newton, Suzaku, and/or Swift/XRT, we identify 40 hard (>10 keV) X-ray-detected active galactic nuclei (AGNs) and constrain their torus parameters with the X-ray clumpy torus model XCLUMPY. Among the AGNs at z < 0.03, for which sample biases are minimized, the fraction of Compton-thick (N _H ≥ 10²⁴ cm^-2) AGNs reaches ${64}_{-15}^{+14}$ % (6/9 sources) in late mergers, while it is ${24}_{-10}^{+12}$ % (3/14 sources) in early mergers, consistent with the tendency reported by Ricci et al. We find that the bolometric AGN luminosities derived from the infrared data increase but the X-ray to bolometric luminosity ratios decrease with merger stage. The X-ray-weak AGNs in late mergers ubiquitously show massive outflows at subparsec to kiloparsec scales. Among them, the most luminous AGNs (L _bol,AGN ~ 10⁴⁶ erg s^-1) have relatively small column densities of ≲10²³ cm^-2 and almost super-Eddington ratios (λ _Edd ~ 1.0). Their torus covering factors (C _T ⁽²²⁾ ~ 0.6) are larger than those of Swift/BAT-selected AGNs with similarly high Eddington ratios. These results suggest a scenario where, in the final stage of mergers, multiphase strong outflows are produced due to chaotic quasi-spherical inflows, and the AGN becomes extremely X-ray weak and deeply buried due to obscuration by inflowing and/or outflowing material.