Hyper-accretion occurs when the gas inflow rate onto a black hole (BH) is so high that the radiative feedback cannot reverse the accretion flow. This extreme process is a promising mechanism for the rapid growth of seed BHs in the early universe, which can explain high-redshift quasars powered by billion solar mass BHs. In theoretical models, spherical symmetry is commonly adopted for hyper-accretion flow, however its sustainability on timescales corresponding to the BH growth has not been addressed yet. Here we show that stochastic interactions between the ionizing radiation from the BH and non-uniform accretion flow can lead to a formation of rotating gas disk around the BH. Once the disk is formed, the supply of gas to the BH preferentially occurs via biconical dominated accretion flow (BDAF) perpendicular to the disk, avoiding the centrifugal barrier of the disk. BDAFs from opposite directions collide in the vicinity of the BH supplying high-density, low angular momentum gas to the BH, whereas most of the gas with non-negligible angular momentum is deflected to the rotationally supported outflowing decretion disk. The disk becomes reinforced progressively as more mass from the biconical flow transfers to the disk and some of the outflowing gas from the disk is redirected to the biconical accretion funnels through a meridional structure. This axisymmetric hydrodynamic structure of a BDAF and decretion disk continues to provide uninterrupted flow of high-density gas to the BH.

Recent observations provide evidence that some cool-core clusters (CCCs) host quasars in their brightest cluster galaxies (BCGs). Motivated by these findings we use 3D radiation-hydrodynamic simulations with the code Enzo to explore the joint role of the kinetic and radiative feedback from supermassive black holes (SMBHs) in BCGs. We implement kinetic feedback as sub-relativistic plasma outflows and model radiative feedback using the ray-tracing radiative transfer or thermal energy injection. In our simulations the central SMBH transitions between the radiatively efficient and radiatively inefficient states on timescales of a few Gyr, as a function of its accretion rate. The timescale for this transition depends primarily on the fraction of power allocated to each feedback mode, and to a lesser degree on the overall feedback luminosity of the active galactic nucleus (AGN). Specifically, we find that (a) kinetic feedback must be present at both low and high accretion rates in order to prevent the cooling catastrophe, and (b) its contribution likely accounts for > 10% of the total AGN feedback power, since below this threshold simulated BCGs tend to host radio-loud quasars most of the time, in apparent contrast with observations. We also find a positive correlation between the AGN feedback power and the mass of the cold gas filaments in the cluster core, indicating that observations of H$\alpha$ filaments can be used as a measure of AGN feedback.

This paper describes the liquid scintillator veto system for the COSINE-100 dark matter experiment and its performance. The COSINE-100 detector consists of 8 NaI(Tl) crystals immersed in 2,200~L of linear alkylbenzene based liquid scintillator. The liquid scintillator tags between 65 and 75\% of the internal $^{40}$K background in the 2--6~keV energy region. We also describe the background model for the liquid scintillator, which is primarily used to assess its energy calibration and threshold.

Dynamical friction is thought to be a principal mechanism responsible for orbital evolution of massive black holes (MBHs) in the aftermath of galactic mergers and an important channel for formation of gravitationally bound MBH binaries. We use 2D radiative hydrodynamic simulations to investigate the efficiency of dynamical friction in the presence of radiative feedback from an MBH moving through a uniform density gas. We find that ionizing radiation that emerges from the innermost parts of the MBH's accretion flow strongly affects the dynamical friction wake and renders dynamical friction inefficient for a range of physical scenarios. MBHs in this regime tend to experience positive net acceleration, meaning that they speed up, contrary to the expectations for gaseous dynamical friction in absence of radiative feedback. The magnitude of this acceleration is however negligibly small and should not significantly alter the velocity of MBHs over relevant physical timescales. Our results suggest that suppression of dynamical friction is more severe at the lower mass end of the MBH spectrum which, compounded with inefficiency of the gas drag for lower mass objects in general, implies that $< 10^7$ solar mass MBHs have fewer means to reach the centers of merged galaxies. These findings provide formulation for a sub-resolution model of dynamical friction in presence of MBH radiative feedback that can be easily implemented in large scale simulations.

The alpha effect is known to be an indispensable energy source of the poloidal magnetic field in the sun or planet. However, the alpha effect is quenched as the magnetic field grows due to the conservation of magnetic helicity. With these conventional understanding, what indeed generates and sustains the observed poloidal magnetic field remains a mystery. To solve this contradiction between theory and real nature, we derived a semi-analytic representation of alpha and beta using large scale magnetic helicity and energy. Applying the simulation data to alpha and beta, we found that the negative beta effect is a promising substitution of the quenched alpha effect. However, since the negative beta effect contradicts the conventional dynamo theory, we derived the new beta expression referring to the field structure model. This analytic result with the field relation between velocity U and magnetic field B shows that the beta effect in the helical system is not a fixed one. Rather, it plays a variable and dynamic role in dynamo depending on the interaction between the poloidal velocity field and relative strength of the large scale magnetic field.

There is a long standing debate about whether or not the annual modulation signal reported by the DAMA/LIBRA collaboration is induced by Weakly Interacting Massive Particles~(WIMP) in the galaxy's dark matter halo scattering from nuclides in their NaI(Tl) crystal target/detector. This is because regions of WIMP-mass vs. WIMP-nucleon cross-section parameter space that can accommodate the DAMA/LIBRA-phase1 modulation signal in the context of the standard WIMP dark matter galactic halo and isospin-conserving~(canonical), spin-independent~(SI) WIMP-nucleon interactions have been excluded by many of other dark matter search experiments including COSINE-100, which uses the same NaI(Tl) target/detector material. Moreover, the recently released DAMA/LIBRA-phase2 results are inconsistent with an interpretation as WIMP-nuclide scattering via the canonical SI interaction and prefer, instead, isospin-violating or spin-dependent interactions. Dark matter interpretations of the DAMA/LIBRA signal are sensitive to the NaI(Tl) scintillation efficiency for nuclear recoils, which is characterized by so-called quenching factors~(QF), and the QF values used in previous studies differ significantly from recently reported measurements, which may have led to incorrect interpretations of the DAMA/LIBRA signal. In this article, the compatibility of the DAMA/LIBRA and COSINE-100 results, in light of the new QF measurements is examined for different possible types of WIMP-nucleon interactions. The resulting allowed parameter space regions associated with the DAMA/LIBRA signal are explicitly compared with 90\% confidence level upper limits from the initial 59.5~day COSINE-100 exposure. With the newly measured QF values, the allowed 3$\sigma$ regions from the DAMA/LIBRA data are still generally excluded by the COSINE-100 data.

COSINE-100 is a direct detection dark matter search experiment that uses a 106 kg array of eight NaI(Tl) crystals that are kept underground at the Yangyang Underground Laboratory to avoid cosmogenic activation of radioisotopes by cosmic rays. Even though the cosmogenic activity is declining with time, there are still significant background rates from the remnant nuclides. In this paper, we report measurements of cosmogenic isotope contaminations with less than one year half-lives that are based on extrapolations of the time dependent activities of their characteristic energy peaks to activity rates at the time the crystals were deployed underground. For longer-lived $^{109}$Cd ($T_{1/2}=1.6$ y) and $^{22}$Na ($T_{1/2}=2.6$ y), we investigate time correlations of characteristic $\gamma$/X-ray peaks. The inferred sea-level production rates are compared with caluclations based on the ACTIVIA and MENDL-2 model calculations and experimental data. For $^{3}$H, which has a long, 12.3 year half-life, we evaluated the activity levels from the exposure times and determined a cosmogenic activation rate that is consistent with other measurements.

We present new constraints on the dark matter-induced annual modulation signal using 1.7 years, of COSINE-100 data with a total exposure of 97.7 kg$\cdot$years. The COSINE-100 experiment, consisting of 106 kg of NaI(Tl) target material, is designed to carry out a model-independent test of DAMA/LIBRA's claim of WIMP discovery by searching for the same annual modulation signal using the same NaI(Tl) target. The crystal data show a 2.7 cpd/kg/keV background rate on average in the 2--6 keV energy region of interest. Using a $\chi$-squared minimization method we observe best fit values for modulation amplitude and phase of 0.0092$\pm$0.0067 cpd/kg/keV and 127.2$\pm$45 d, respectively.

We present results from a search for solar axions with the COSINE-100 experiment. We find no evidence of solar axion events from a data-set of 6,303.9 kg$\cdot$days exposure and set a 90\,\% confidence level upper limit on the axion-electron coupling, $g_{ae}$, at 1.70~$\times$~$10^{-11}$ for an axion mass less than 1\,keV/c$^2$. This limit excludes QCD axions heavier than 0.59\,eV/c$^2$ in the DFSZ model and 168.1\,eV/c$^2$ in the KSVZ model.

Assuming a standard Maxwellian for the WIMP velocity distribution, we obtain the bounds from null WIMP search results of 59.5 days of COSINE-100 data on the DAMA/LIBRA-phase2 modulation effect within the context of the non-relativistic effective theory of WIMP-nucleus scattering. Here, we systematically assume that one of the effective operators allowed by Galilean invariance dominates in the effective Hamiltonian of a spin-1/2 dark matter (DM) particle. We find that, although DAMA/LIBRA and COSINE-100 use the same sodium-iodide target, the comparison of the two results still depends on the particle-physics model. This is mainly due to two reasons: i) the WIMP signal spectral shape; ii) the expected modulation fractions, when the upper bound on the time-averaged rate in COSINE-100 is converted into a constraint on the annual modulation component in DAMA/LIBRA. We find that the latter effect is the dominant one. For several effective operators the expected modulation fractions are larger than in the standard spin-independent or spin-dependent interaction cases. As a consequence, compatibility between the modulation effect observed in DAMA/LIBRA and the null result from COSINE-100 is still possible for several non-relativistic operators. At low WIMP masses such relatively high values of the modulation fractions arise because COSINE-100 is mainly sensitive to WIMP-sodium scattering events, due to the higher threshold compared to DAMA/LIBRA. A next COSINE analysis is expected to have a full sensitivity for the 5$\sigma$ region of DAMA/LIBRA.

Observations of galaxies and primordial radiation suggest that the Universe is made mostly of non-luminous dark matter. Several types of new fundamental particles have been proposed as candidates for dark matter such as weakly interacting massive particles (WIMPs) but no definitive signal has been seen despite concerted efforts by many collaborations. One exception is the much-debated claim by the DAMA collaboration of a statistically significant annual modulation in the event rate of their experiment with a period and phase consistent with that expected from WIMP dark matter. Several groups have been working to develop experiments with the aim of reproducing DAMA's results using the same target medium. Here we report results from the initial operation of the COSINE-100 experiment. COSINE-100 uses sodium iodide as the target medium-the same medium as DAMA-and is designed to carry out a model-independent test of DAMA's claim. Initial data based on the first 59.5 days indicate that there is no excess of events over the expected background, confirming that DAMA's annual modulation signal is in severe tension with results from other experiments under the assumption of dark matter having spin independent interactions and the Standard Halo Model. COSINE-100 is now taking data to study the presence of dark matter-induced annual modulation in the event rate of the sodium iodide detectors.

Intensive Monitoring Survey of Nearby Galaxies (IMSNG) is a high cadence observation program monitoring nearby galaxies with high probabilities of hosting supernovae (SNe). IMSNG aims to constrain the SN explosion mechanism by inferring sizes of SN progenitor systems through the detection of the shock-heated emission that lasts less than a few days after the SN explosion. To catch the signal, IMSNG utilizes a network of 0.5-m to 1-m class telescopes around the world and monitors the images of 60 nearby galaxies at distances D < 50 Mpc to a cadence as short as a few hours. The target galaxies are bright in near-ultraviolet (NUV) with M_NUV < -18.4 AB mag and have high probabilities of hosting SNe (0.06 SN/yr per galaxy). With this strategy, we expect to detect the early light curves of 3.4 SNe per year to a depth of R ~ 19.5 mag, enabling us to detect the shock-heated emission from a progenitor star with a radius as small as 0.1 R_sun. The accumulated data will be also useful for studying faint features around the target galaxies and other science projects. So far, 18 SNe have occurred in our target fields (16 in IMSNG galaxies) over 5 years, confirming our SN rate estimate of 0.06 SN/yr per galaxy.

A search for inelastic boosted dark matter (iBDM) using the COSINE-100 detector with 59.5 days of data is presented. This relativistic dark matter is theorized to interact with the target material through inelastic scattering with electrons, creating a heavier state that subsequently produces standard model particles, such as an electron-positron pair. In this study, we search for this electron-positron pair in coincidence with the initially scattered electron as a signature for an iBDM interaction. No excess over the predicted background event rate is observed. Therefore, we present limits on iBDM interactions under various hypotheses, one of which allows us to explore an area of the experimental search for iBDM using a terrestrial detector.

We explain the (non)helical dynamo process using a field-structure model based on magnetic induction equation in an intuitive way. We show how nonhelical kinetic energy converts into magnetic energy and cascades toward smaller eddies in a mechanically forced plasma system. Also, we show how helical magnetic energy is inversely cascaded ($\alpha$ effect) toward large scale magnetic eddies in a mechanically or magnetically forced system. We, then, compare the simulation results with the model qualitatively for the verification of the model. In addition to these intuitive and numerical approaches, we show how to get $\alpha$ and $\beta$ coefficient semi-analytically from the temporally evolving large scale magnetic energy and magnetic helicity.

COSINE-100 is a dark matter direct detection experiment designed to test the annual modulation signal observed by the DAMA/LIBRA experiment. COSINE-100 consists of 8 NaI(Tl) crystals with a total mass of 106 kg, a 2200 L liquid scintillator veto, and 37 muon detector panels. We present details of the data acquisition system of COSINE-100, including waveform storage using flash analog-to-digital converters for crystal events and integrated charge storage using charge-sensitive analog-to-digital converters for liquid scintillator and plastic scintillator muon veto events. We also discuss several trigger conditions developed in order to distinguish signal events from photomultiplier noise events. The total trigger rate observed for the crystal/liquid scintillator (plastic scintillator) detector is 15 Hz (24 Hz).

The COSINE-100 dark matter search experiment is an array of NaI(Tl) crystal detectors located in the Yangyang Underground Laboratory (Y2L). To understand measured backgrounds in the NaI(Tl) crystals we have performed Monte Carlo simulations using the Geant4 toolkit and developed background models for each crystal that consider contributions from both internal and external sources, including cosmogenic nuclides. The background models are based on comparisons of measurement data with Monte Carlo simulations that are guided by a campaign of material assays and are used to evaluate backgrounds and identify their sources. The average background level for the six crystals (70 kg total mass) that are studied is 3.5 counts/day/keV/kg in the (2-6) keV energy interval. The dominant contributors in this energy region are found to be $^{210}$Pb and $^3$H.

COSINE is a dark matter search experiment based on an array of low background NaI(Tl) crystals located at the Yangyang underground laboratory. The assembly of COSINE-100 was completed in the summer of 2016 and the detector is currently collecting physics quality data aimed at reproducing the DAMA/LIBRA experiment that reported an annual modulation signal. Stable operation has been achieved and will continue for at least two years. Here, we describe the design of COSINE-100, including the shielding arrangement, the configuration of the NaI(Tl) crystal detection elements, the veto systems, and the associated operational systems, and we show the current performance of the experiment.

The COSINE-100 dark matter search experiment has started taking physics data with the goal of performing an independent measurement of the annual modulation signal observed by DAMA/LIBRA. A muon detector was constructed by using plastic scintillator panels in the outermost layer of the shield surrounding the COSINE-100 detector. It is used to detect cosmic ray muons in order to understand the impact of the muon annual modulation on dark matter analysis. Assembly and initial performance test of each module have been performed at a ground laboratory. The installation of the detector in Yangyang Underground Laboratory (Y2L) was completed in the summer of 2016. Using three months of data, the muon underground flux was measured to be 328 $\pm$ 1(stat.)$\pm$ 10(syst.) muons/m$^2$/day. In this report, the assembly of the muon detector and the results from the analysis are presented.

Accretion of gas and interaction of matter and radiation are at the heart of many questions pertaining to black hole (BH) growth and coevolution of massive BHs and their host galaxies. To answer them it is critical to quantify how the ionizing radiation that emanates from the innermost regions of the BH accretion flow couples to the surrounding medium and how it regulates the BH fueling. In this work we use high resolution 3-dimensional (3D) radiation-hydrodynamic simulations with the code Enzo, equipped with adaptive ray tracing module Moray, to investigate radiation-regulated BH accretion of cold gas. Our simulations reproduce findings from an earlier generation of 1D/2D simulations: the accretion powered UV and X-ray radiation forms a highly ionized bubble, which leads to suppression of BH accretion rate characterized by quasi-periodic outbursts. A new feature revealed by the 3D simulations is the highly turbulent nature of the gas flow in vicinity of the ionization front. During quiescent periods between accretion outbursts, the ionized bubble shrinks in size and the gas density that precedes the ionization front increases. Consequently, the 3D simulations show oscillations in the accretion rate of only ~2-3 orders of magnitude, significantly smaller than 1D/2D models. We calculate the energy budget of the gas flow and find that turbulence is the main contributor to the kinetic energy of the gas but corresponds to less than 10% of its thermal energy and thus does not contribute significantly to the pressure support of the gas.

We have developed ultra-low-background NaI(Tl) crystals to reproduce the DAMA results with the ultimate goal of achieving purity levels that are comparable to or better than those of the DAMA/LIBRA crystals. Even though the achieved background level does not approach that of DAMA/LIBRA, it is crucial to have a quantitative understanding of the backgrounds. We have studied background simulations toward a deeper understanding of the backgrounds and developed background models for a 9.16-kg NaI(Tl) crystal used in the test arrangement. In this paper we describe the contributions of background sources quantitatively by performing Geant4 Monte Carlo simulations that are fitted to the measured data to quantify the unknown fractions of the background compositions. In the fitted results, the overall simulated background spectrum well describes the measured data with a 9.16-kg NaI(Tl) crystal and shows that the background sources are dominated by surface $^{210}$Pb and internal $^{40}$K in the 2 to 6-keV energy interval, which produce 2.4 counts/day/keV/kg (dru) and 0.5 dru, respectively.

Evidence for dust around supermassive black holes (SMBHs) in the early Universe is strongly suggested by recent observations. However, the accretion mechanism of SMBHs in dusty gas is not well understood yet. We investigate the growth of intermediate-mass black-holes (IMBHs) of $\sim 10^{5}~M_{\odot}$ in dusty clouds by using one-dimensional radiative-hydrodynamics simulations. We find that the accretion of dusty gas onto IMBHs proceeds gently with small fluctuations of the accretion rate, whereas that of pristine gas causes more violent periodic bursts. At dust-to-gas mass ratios similar to the solar neighborhood, the time averaged luminosity becomes smaller than that for primordial gas by one order of magnitude and the time-averaged Eddington ratio ranges from $\sim 10^{-4}$ to $\sim 10^{-2}$ in clouds with initial gas densities of $n_{\rm H} = 10 - 1000~\rm cm^{-3}$. Our calculations show that the effect of dust opacity alone is secondary compared to the radiation pressure on dust in regulating the BH growth. We also derive spectral energy distributions at IR bands by calculating dust thermal emission and show that the flux ratio between $\lambda \lesssim 20~\rm \mu m$ and $\gtrsim 100~\rm \mu m$ is closely related to the Eddington ratio. Thermal emission from hot dust near the BH dominates only during the high accretion phase, producing higher flux density at $\lesssim 20~\rm \mu m$. Therefore, we suggest that the combinations of MIR observations by JWST and FIR observation by ALMA or Spitzer can be used to estimate the Eddington ratio of massive BHs.

We report the performance of an active veto system using a liquid scintillator with NaI(Tl) crystals for use in a dark matter search experiment. When a NaI(Tl) crystal is immersed in the prototype detector, the detector tags 48% of the internal K-40 background in the 0-10 keV energy region. We also determined the tagging efficiency for events at 6-20 keV as 26.5 +/- 1.7% of the total events, which corresponds to 0.76 +/- 0.04 events/keV/kg/day. According to a simulation, approximately 60% of the background events from U, Th, and K radioisotopes in photomultiplier tubes are tagged at energies of 0-10 keV. Full shielding with a 40-cm-thick liquid scintillator can increase the tagging efficiency for both the internal K-40 and external background to approximately 80%.

The Korea Invisible Mass Search (KIMS) collaboration has developed low-background NaI(Tl) crystals that are suitable for the direct detection of WIMP dark matter. With experience built on the KIMS-CsI programs, the KIMS-NaI experiment will consist of a 200~kg NaI(Tl) crystal array surrounded by layers of shielding structures and will be operated at the Yangyang underground laboratory. The goal is to provide an unambiguous test of the DAMA/LIBRA's annual modulation signature. Measurements of six prototype crystals show progress in the reduction of internal contaminations of radioisotopes. Based on our understanding of these measurements, we expect to achieve a background level in the final detector configuration that is less than 1~count/day/keV/kg for recoil energies around 2~keV. The annual modulation sensitivity for the KIMS-NaI experiment shows that an unambiguous 7$\sigma$ test of the DAMA/LIBRA signature would be possible with a 600~kg$\cdot$year exposure with this system.

The inverse cascade of magnetic energy occurs when helicity or rotational instability exists in the magnetohydrodynamic (MHD) system. This well known phenomenon has been considered as a basis for the large scale magnetic field in universe. However nonhelical magnetic energy in a decaying MHD system also migrates toward the large scale, which holds vital clues to the origin of large scale magnetic field in a quiescent astrophysical system. Zeldovich's rope dynamo model is considered as a basic and symbolistic model of magnetic field amplification. However, the rope model assuming specific external forces like buoyancy or Coriolis force is not appropriate for a decaying turbulent system without any external force. So we suggest a new dynamo model based on magnetic induction equation excluding a forcing source. This model shows the expansion and growth of magnetic field (flux) is basically the redistribution of energy in the system. The transfer of magnetic energy is in fact a successive induction of magnetic field resulted from the interaction between the fluid motion and seed magnetic field. We also discuss about an analytic theorem based on the scaling invariant MHD equation.

We examine radiation-regulated accretion onto intermediate-mass and massive black holes (BHs) embedded in a bulge component. Using spherically symmetric one-dimensional radiation-hydrodynamics simulations, we track the growth of BHs accreting from a cold, neutral gas reservoir with temperature T=10^4 K. We find that the accretion rate of BHs embedded in bulges is proportional to r_{B,eff}/r_B, where r_{B,eff} is the increased effective Bondi radius that includes the gravitational potential of the bulge, and r_B is the Bondi radius of the BH. The radiative feedback from the BH suppresses the cold accretion rate to ~1 percent of the Bondi rate when a bulge is not considered. However, we find that the BH fueling rate increases rapidly when the bulge mass M_bulge is greater than the critical value of 10^6 M_sun and is proportional to M_bulge. Since the critical bulge mass is independent of the central BH mass M_{BH}, the growth rate of BHs with masses of 10^2, 10^4, and 10^6 M_sun exhibits distinct dependencies on the bulge-to-BH mass ratio. Our results imply that light seed BHs (<= 10^2 M_sun) which might be the remnants of the Pop III stars, cannot grow through accretion coevally with the early assembly of the bulge of the host galaxies until the bulge reaches the critical mass. However, massive BH seeds (>= 10^5 M_sun) that may form via direct collapse, are more likely to be embedded in a supercritical bulge and thus can grow efficiently coupling to the host galaxies and driving the early evolution of the M_{BH}-$\sigma$ relationship.

We report on the response of a high light-output NaI(Tl) crystal to nuclear recoils induced by neutrons from an Am-Be source and compare the results with the response to electron recoils produced by Compton scattered 662 keV $\gamma$-rays from a $^{137}$Cs source. The measured pulse-shape discrimination (PSD) power of the NaI(Tl) crystal is found to be significantly improved because of the high light output of the NaI(Tl) detector. We quantify the PSD power with a quality factor and estimate the sensitivity to the interaction rate for weakly interacting massive particles (WIMPs) with nucleons, and the result is compared with the annual modulation amplitude observed by the DAMA/LIBRA experiment. The sensitivity to spin-independent WIMP-nucleon interactions based on 100 kg$\cdot$year of data from NaI detectors is estimated with simulated experiments, using the standard halo model.

We examine the luminosity function (LF) of [OII] emission-line galaxies in the high-resolution cosmological simulation MassiveBlack-II (MBII). From the spectral energy distribution of each galaxy, we select a sub-sample of star-forming galaxies at $0.06 \le z \le 3.0$ using the [OII] emission line luminosity L([OII]). We confirm that the specific star formation rate matches that in the GAMA survey. We show that the [OII] LF at z=1.0 from the MBII shows a good agreement with the LFs from several surveys below L([OII])=$10^{43.0}$ erg/s while the low redshifts ($z \le 0.3$) show an excess in the prediction of bright [OII] galaxies, but still displaying a good match with observations below L([OII])=$10^{41.6}$ erg/s. Based on the validity in reproducing the properties of [OII] galaxies at low redshift ($z \le 1$), we forecast the evolution of the [OII] LF at high redshift ($z \le 3$), which can be tested by upcoming surveys such as the HETDEX and DESI. The slopes of the LFs at bright and faint ends range from -3 to -2 showing minima at z=2. The slope of the bright end evolves approximately as 1/(z+1) at z=2 while the faint end evolves as ~3/(z+1) at $0.6 \le z \le 2$. In addition, a similar analysis is applied for the evolution of [OIII] LFs, which is to be explored in the forthcoming survey WFIRST-AFTA. Finally, we show that the auto-correlation function of [OII] and [OIII] emitting galaxies shows a rapid evolution from z=2 to 1.

Intra Cluster Media (ICMs) located at galaxy clusters is in the state of hot, tenuous, magnetized, and highly ionized X-ray emitting plasmas. This overall collisionless, viscous, and conductive magnetohydrodynamic (MHD) turbulence in ICM is simulated using hyper and physical magnetic diffusivity. The results show that fluctuating random plasma motion amplifies the magnetic field, which cascades toward the diffusivity scale passing through the viscous scale. The kinetic eddies in the subviscous scale are driven and constrained by the magnetic tension which finally gets balanced with the highly damping effect of the kinetic eddies. However, the saturated kinetic energy spectrum is deeper than that of the incompressible or compressible hydrodynamics fluid. To explain this unusual field profile we set up two simultaneous differential equations for the kinetic and magnetic energy spectrum using an Eddy Damped Quasi Normal Markovianized (EDQNM) approximation. The analytic solution tells us that the magnetic energy in addition to the viscous damping effect constrains the plasma motion leading to the power spectra: kinetic energy spectrum $E_V^k\sim k^{-3}$ and corresponding representative magnetic energy spectrum $E_M^k\sim k^{-1/2}$. Also the comparison of simulation results with different resolutions and magnetic diffusivities implies the role of small scale magnetic energy in dynamo.

We investigate the role of Compton heating in radiation-regulated accretion on to black holes from a neutral dense medium using 1D radiation-hydrodynamic simulations. We focus on the relative effects of Compton-heating and photo-heating as a function of the spectral slope {\alpha}, assuming a power-law spectrum in the energy range of 13.6 eV--100 keV. While Compton heating is dominant only close to the black hole, it can reduce the accretion rate to 0.1 % ($l \propto \dot{m}^2$ model)--0.01 % ($l \propto \dot{m}$ model) of the Bondi accretion rate when the BH radiation is hard ({\alpha} ~ 1), where $l$ and $\dot{m}$ are the luminosity and accretion rate normalised by Eddington rates, respectively. The oscillatory behaviour otherwise typically seen in simulations with {\alpha} > 1, become suppressed when {\alpha} ~ 1 only for the $l \propto \dot{m}$ model. The relative importance of the Compton heating over photo-heating decreases and the oscillatory behaviour becomes stronger as the spectrum softens. When the spectrum is soft ({\alpha} > 1.5), photo-heating prevails regardless of models making the effect of Compton heating negligible. On the scale of the ionization front, where the gas supply into the Str\"omgren sphere from large scale is regulated, photo-heating dominates. Our simulations show consistent results with the advection-dominated accretion flow ($l \propto \dot{m}^2$) where the accretion is inefficient and the spectrum is hard ({\alpha} ~ 1).

The influence of initially given small scale magnetic energy($E_M(0)$) and helicity($H_M(0)$) on the magnetohydrodynamics(MHD) dynamo was investigated. Equations for $E_M$(t), $H_M$(t), and electromotive force($\langle {\bf v}\times {\bf b}\rangle$, $EMF$) were derived and solved. The solutions indicate small scale magnetic field(${\bf b}_i$) caused by $E_M$(0) modifies $EMF$ and generates additional terms of which effect depends on magnetic diffusivity $\eta$, position of initial conditions($IC$s) $k_f$, and time ($\sim e^{-\eta k_f^2 t}$). ${\bf b}_i$ increases the inverse cascade of energy resulting in the enhanced growth of large scale magnetic field($\overline{{\bf B}}$). Simulation data show that $E_M$(0) in small scale boosts the growth rate, which also proportionally depends on $H_M(0)$. If $E_M$(0) is the same, positive $H_M(0)$ is more effective for MHD dynamo than negative $H_M(0)$ is. It was discussed why large scale magnetic helicity should have the opposite sign of the injected kinetic helicity.

We examine the role of ionizing radiation emitted from black holes (BHs) in suppressing the growth of the Rayleigh-Taylor instability (RTI) across the ionization front (I-front) that forms when the gas fueling the BH is neutral. We use radiation-hydrodynamic simulations to show that the RTI is suppressed for non-accelerating fronts on all scales resolved in our simulations. A necessary condition for the stability of the I-front is that the radius of the Str\"omgren sphere is larger than the Bondi radius. When this condition is violated the I-front collapses producing an accretion luminosity burst. Transient growth of the RTI occurs only during the accretion burst when the effective acceleration in the frame of reference of the I-front increases significantly due to the rapid expansion of the Str\"omgren sphere.

To investigate the effect of energy and helicity on the growth of magnetic field, helical kinetic forcing was applied to the magnetohydrodynamic(MHD) system that had a specific distribution of energy and helicity as initial conditions. Simulation results show the saturation of a system is not influenced by the initial conditions, but the growth rate of large scale magnetic field is proportionally dependent on the initial large scale magnetic energy and helicity. It is already known that the helical component of small scale magnetic field(i.e., current helicity $<{\bf j}\cdot {\bf b}>$) quenches the growth of large scale magnetic field. However, $<{\bf j}\cdot {\bf b}>$ can also boost the growth of large scale magnetic field by changing its sign and magnitude. In addition, simulation shows the nonhelical magnetic field can suppress the velocity field through Lorentz force. Comparison of the profiles of evolving magnetic and kinetic energy indicates that kinetic energy migrates backward when the external energy flows into the three dimensional MHD system, which means the velocity field may play a preceding role in the very early MHD dynamo stage.

Understanding large scale magnetic field growth in turbulent plasmas in the magnetohydrodynamic limit is a goal of magnetic dynamo theory. In particular, assessing how well large scale helical field growth and saturation in simulations matches that predicted by existing theories is important for progress. Using numerical simulations of isotropically forced turbulence without large scale shear with the implications, we focus on several aspects of this comparison that have not been previously tested: (1) Leading mean field dynamo theories which break the field into large and small scales predict that large scale helical field growth rates are determined by the difference between kinetic helicity and current helicity with no dependence on the non-helical energy in small scale magnetic fields. Our simulations show that the growth rate of the large scale field from fully helical forcing is indeed unaffected by the presence or absence of small scale magnetic fields amplified in a precursor non-helical dynamo. However, because the precursor non helical dynamo in our simulations produced fields that were strongly sub-equipartition with respect to the kinetic energy, we cannot yet rule out the potential influence of stronger non- helical small scale fields. (2) We have identified two features in our simulations which cannot be explained by the most minimalist versions of two-scale mean field theory: (i) fully helical small scale forcing produces significant non-helical large scale magnetic energy and (ii) the saturation of the large scale field growth is time-delayed with respect to what minimalist theory predicts. We comment on desirable generalizations to the theory in this context and future desired work.

In this third paper of a series, we study the growth and luminosity of black holes (BHs) in motion with respect to their surrounding medium. We run a large set of two-dimensional (2D) axis-symmetric simulations to explore a large parameter space of initial conditions and formulate an analytical model for the accretion. Contrary to the case without radiation feedback, the accretion rate increases with increasing BH velocity $\vbh$ reaching a maximum value at $\vbh=2c_{\rm s,in}\sim 50$ km/s, where $c_{\rm s,in}$ is the sound speed inside the "cometary-shaped" \hii region around the BH, before decreasing as $\vbh^{-3}$ when the I-front becomes {\it R}-type (rarefied) and the accretion rate approaches the classical Bondi-Hoyle-Lyttleton solution. The increase of the accretion rate with $\vbh$ is produced by the formation of a {\it D}-type (dense) ionization front (I-front) preceded by a standing bow-shock that reduces the downstream gas velocity to transonic values. There is a range of densities and velocities where the dense shell is unstable producing periodic accretion rate peaks which can significantly increase the detectability of intermediate mass BHs. We find that the mean accretion rate for a moving BH is larger than that of a stationary BH of the same mass if the medium temperature is $T_\infty<10^4$ K. This result could be important for the growth of seed BHs in the multi-phase medium of the first galaxies and for building and early X-ray background that may affect the formation of the first galaxies and the reionization process.

Progress toward understanding principles of nonlinear growth and saturation of large scale magnetic fields has emerged from comparison of theoretical models that incorporate the evolution of magnetic helicity with numerical simulations for problems that are more idealized than expected in astrophysical circumstances, but still fully non-linear. We carry out a new comparison of this sort for the magnetic field growth from forced isotropic helical turbulence in a periodic box. Previous comparisons be- tween analytic theory and simulations of this problem have shown that a two-scale model compares well with the simulations in agreeing that the driver of large scale field growth is the difference between kinetic and current helicities associated with the small scale field, and that the backreaction that slows the growth of the large scale field as the small scale current helicity grows. However, the two-scale model artificially re- stricts the small scale current helicity to reside at the same scale as the driving kinetic helicity. In addition, previous comparisons have focused on the late time saturation regime, and less on the early-time growth regime. Here we address these issues by comparing a three scale model to new simulations for both early and late time growth regimes. We find that the minimalist extension to three scales provides a better model for the field evolution than the two-scale model because the simulations show that the small scale current helicity does not reside at the same scale as that of the driving kinetic helicity. The simulations also show that the peak of the small scale current helicity migrates toward lower wave numbers as the growth evolves from the fast to saturated growth regimes.

This paper, the second of a series on radiation-regulated accretion onto black holes (BHs) from galactic scales, focuses on the effects of radiation pressure and angular momentum of the accreting gas. We simulate accretion onto intermediate-mass black holes, but we derive general scaling relationships that are solutions of the Bondi problem with radiation feedback valid for any mass of the BH $M_{\rm bh}$. Thermal pressure of the ionized sphere around the BH regulates the accretion rate producing periodic and short-lived luminosity bursts. We find that for ambient gas densities exceeding $n^{\rm cr}_{\rm H,\infty} \propto M_{\rm bh}^{-1}$, the period of the oscillations decreases rapidly and the duty cycle increases from 6%, in agreement with observations of the fraction of active galactic nuclei at $z\sim 3$, to 50%. The mean accretion rate becomes Eddington limited for $n_{\rm H,\infty}>n_{\rm H,\infty}^{\rm Edd} \simeq n_{\rm H,\infty}^{\rm cr} T_{\infty,4}^{-1}$ where $T_{\infty,4}$ is the gas temperature in units of $10^4$ K. In the sub-Eddington regime, the mean accretion rate onto BHs is about $1% T_{\infty,4}^{2.5}$ of the Bondi rate, thus is proportional to the thermal pressure of the ambient medium. The period of the oscillations coincides with depletion timescale of the gas inside the ionized bubble surrounding the BH. Gas depletion is dominated by a pressure gradient pushing the gas outward if $n_{\rm H,\infty}<n^{\rm cr}_{\rm H,\infty}$ and by accretion onto the BH otherwise. Generally, for $n_{\rm H,\infty}<n^{\rm cr}_{\rm H,\infty}$ angular momentum does not affect significantly the accretion rate and period of the oscillations.

Models of large scale (magnetohydrodynamic) dynamos (LSD) which couple large scale field growth to total magnetic helicity evolution best predict the saturation of LSDs seen in simulations. For the simplest so called "{\alpha}2" LSDs in periodic boxes, the electromotive force driving LSD growth depends on the difference between the time-integrated kinetic and current helicity associated with fluctuations. When the system is helically kinetically forced (KF), the growth of the large scale helical field is accompanied by growth of small scale magnetic (and current) helicity which ultimately quench the LSD. Here, using both simulations and theory, we study the complementary magnetically forced(MF) case in which the system is forced with an electric field that supplies magnetic helicity. For this MF case, the kinetic helicity becomes the back-reactor that saturates the LSD. Simulations of both MF and KF cases can be approximately modeled with the same equations of magnetic helicity evolution, but with complementary initial conditions. A key difference between KF and MF cases is that the helical large scale field in the MF case grows with the same sign of injected magnetic helicity, whereas the large and small scale magnetic helicities grow with opposite sign for the KF case. The MF case can arise even when the thermal pressure is approximately smaller than the magnetic pressure, and requires only that helical small scale magnetic fluctuations dominate helical velocity fluctuations in LSD driving. We suggest that LSDs in accretion discs and Babcock models of the solar dynamo are actually MF LSDs.

We study the effect of radiative feedback on accretion onto intermediate mass black holes (IMBHs) using the hydrodynamical code ZEUS-MP with a radiative transfer algorithm. In this paper, the first of a series, we assume accretion from a uniformly dense gas with zero angular momentum and extremely low metallicity. Our 1D and 2D simulations explore how X-ray and UV radiation emitted near the black hole regulates the gas supply from large scales. Both 1D and 2D simulations show similar accretion rate and period between peaks in accretion, meaning that the hydro-instabilities that develop in 2D simulations do not affect the mean flow properties. We present a suite of simulations exploring accretion across a large parameter space, including different radiative efficiencies and radiation spectra, black hole masses, density and temperature, $T_\infty$, of the neighboring gas. In agreement with previous studies we find regular oscillatory behavior of the accretion rate, with duty cycle $\sim 6%$, mean accretion rate 3% $(T_{\infty}/10^4 {\rm K})^{2.5}$ of the Bondi rate and peak accretion $\sim 10$ times the mean for $T_{\infty}$ ranging between 3000K and 15000K. We derive parametric formulas for the period between bursts, the mean accretion rate and the peak luminosity of the bursts and thus provide a formulation of how feedback regulated accretion operates. The temperature profile of the hot ionized gas is crucial in determining the accretion rate, while the period of the bursts is proportional to the mean size of the Str\"{o}mgren sphere and we find qualitatively different modes of accretion in the high vs. low density regimes. We also find that softer spectrum of radiation produces higher mean accretion rate.

Faraday rotation (FR) is widely used to infer the orientation and strength of magnetic fields in astrophysical plasmas. Although the absence of electron-positron pairs is a plausible assumption in many astrophysical environments, the magnetospheres of pulsars and black holes and their associated jets may involve a significant pair plasma fraction. This motivates being mindful of the effect of positrons on FR. Here we derive and interpret exact expressions of FR for a neutral plasma of arbitrary composition. We focus on electron-ion-positron plasmas in which charge neutrality is maintained by an arbitrary combination of ions and positrons. Because a pure electron-positron plasma has zero FR, the greater the fraction of positrons the higher the field strength required to account for the same FR. We first obtain general formulae and then specifically consider parameters relevant to active galctic nuclei (AGN) jets to illustrate the significant differences in field strengths that FR measurements from radio frequency measurements. Complementarily, using galaxy cluster core plasmas as examples, we discuss how plasma composition can be constrained if independent measurements of the field strength and number density are available and combined with FR.

The helicity amplitudes of the electroexcitation of the Roper resonance are extracted for 1.7 < Q2 < 4.5 GeV2 from recent high precision JLab-CLAS cross section and longitudinally polarized beam asymmetry data for pi+ electroproduction on protons at W=1.15-1.69 GeV. The analysis is made using two approaches, dispersion relations and a unitary isobar model, which give consistent results. It is found that the transverse helicity amplitude A_{1/2} for the gamma* p -> N(1440)P11 transition, which is large and negative at Q2=0, becomes large and positive at Q2 ~ 2 GeV2, and then drops slowly with Q2. The longitudinal helicity amplitude S_{1/2}, which was previously found from CLAS ep -> eppi0,enpi+ data to be large and positive at Q2=0.4,0.65 GeV2, drops with Q2. Available model predictions for gamma* p -> N(1440)P11 allow us to conclude that these results provide strong evidence in favor of N(1440)P11 as a first radial excitation of the 3q ground state. The results of the present paper also confirm the conclusion of our previous analysis for Q2 < 1 GeV2 that the presentation of N(1440)P11 as a 3qG hybrid state is ruled out.