Using cosmological particle hydrodynamical simulations and uniform ultraviolet backgrounds, we compare Lyman-$\alpha$ forest flux spectra predicted by the conventional cold dark matter (CDM) model, the free-particle wave dark matter (FP$\psi$DM) model and extreme-axion wave dark matter (EA$\psi$DM) models of different initial axion field angles against the BOSS Lyman-$\alpha$ forest absorption spectra with a fixed boson mass $m_b\sim 10^{-22}$eV. We recover results reported previously (Ir\v{s}i\v{c} et al. 2017b; Armengaud et al. 2017) that the CDM model agrees better with the BOSS data than the FP$\psi$DM model by a large margin, and we find the difference of total $\chi^2$'s is $120$ for $420$ data bins. These previous results demand a larger boson mass by a factor $>10$ to be consistent with the date and are in tension with the favoured value determined from local satellite galaxies. We however find that such tension is removed as some EA$\psi$DM models predict Lyman-$\alpha$ flux spectra agreeing better with the BOSS data than the CDM model, and the difference of total $\chi^2$'s can be as large as $24$ for the same bin number. This finding arises with no surprise since EA$\psi$DM models have unique spectral shapes with spectral bumps in excess of the CDM power near the small-scale cutoff typical of $\psi$DM linear matter power spectra as well as more extended cutoffs than FP$\psi$DM (Zhang & Chiueh 2017a,b).

GAMER, a parallel Graphic-processing-unit-accelerated Adaptive-MEsh-Refinement hydrodynamic code, has been extended to support magnetohydrodynamics (MHD) with both the corner-transport-upwind (CTU) and MUSCL-Hancock schemes and the constraint transport (CT) technique. The divergent preserving operator for adaptive mesh refinement (AMR) has been applied to reinforce the divergence-free constraint on the magnetic field. GAMER-MHD has fully exploited the concurrent executions between the GPU MHD solver and other CPU computation pertinent to AMR. We perform various standard tests to demonstrate that GAMER-MHD is both second-order accurate and robust, producing results as accurate as those given by high-resolution uniform-grid runs. We also explore a new 3D MHD test, where the magnetic field assumes the Arnold-Beltrami-Childress (ABC) configuration, temporarily becomes turbulent with current sheets and finally settles to a lowest-energy equilibrium state. This 3D problem is adopted for the performance test of GAMER-MHD. The single-GPU performance reaches $1.2\times 10^8$ and $5.5\times 10^7$ cell-updates/sec for the single- and double-precision calculations, respectively, on Tesla P100. We also demonstrate a parallel efficiency of $\sim 70\%$ for both weak and strong scaling using $1,024$ XK nodes on the Blue Waters supercomputers.

Linear perturbations of the wave dark matter, or $\psi$ dark matter ($\psi$DM), of particle mass $\sim 10^{-22}$eV in the radiation-dominant era are analyzed, and the matter power spectrum at the photon-matter equality is obtained. We identify four phases of evolution for $\psi$DM perturbations, where the dynamics can be vastly different from the counterparts of cold dark matter (CDM). While in late stages after mass oscillation long-wave $\psi$DM perturbations are almost identical to CDM perturbations, some subtle differences remain, let alone intermediate-to-short waves that bear no resemblance with those of CDM throughout the whole evolutionary history. The dissimilarity is due to quantum mechanical effects which lead to severe mode suppression. We also discuss the axion model with a cosine field potential. The power spectrum of axion models are generally almost identical to those of $\psi$DM, but in the extreme case when the initial axion angle is near the field potential top, this axion model predict a power excess over a range of wave number and a higher spectral cutoff than $\psi$DM as if $\psi$DM had a higher particle mass.

Sub-horizon perturbations under the extreme initial condition of the axion model are investigated, where initial axion angles start near the potential maximum. This work focuses on a few new features found in the extreme axion model but absent in the free-particle model. A particularly novel new feature is the spectral excess relative to the CDM model in some wave number range, where the excess may be so large that landscapes of high-redshift universe beyond $z=10$ can be significantly altered. For axions of particle mass $10^{-22}$ eV, this range of wave number corresponds to first galaxies of few times $10^9-10^{10} M_\odot$. We demonstrate that sub-horizon perturbations are accurately described by Mathieu's equation and subject to parametric instability, which explains this novel feature. Actually the axion model is not a special one; perturbations in a wide range of scalar field models can share the similar characteristic.

Inspired by Barnard 68, a Bok globule, that undergoes stable oscillations, we perform multi-phase hydrodynamic simulations to analyze the stability of Bok globules. We show that a high-density soft molecular core, with an adiabatic index $\gamma$ = 0.7 embedded in a warm isothermal diffuse gas, must have a small density gradient to retain the stability. Despite being stable, the molecular core can still collapse spontaneously as it will relax to develop a sufficiently large density gradient after tens of oscillations, or a few $10^7$ years. However, during its relaxation, the core may abruptly collapse triggered by the impingement of small-amplitude, long-wavelength ($\sim$ 6 $-$ 36 pc) sound waves in the warm gas. This triggered collapse mechanism is similar to a sonoluminescence phenomenon, where underwater ultrasounds can drive air bubble coalescence. The collapse configuration is found to be different from both inside-out and outside-in models of low-mass star formation; nonetheless the mass flux is close to the prediction of the inside-out model. The condition and the efficiency for this core collapse mechanism are identified. Generally speaking, a broad-band resonance condition must be met, where the core oscillation frequency and the wave frequency should match each other within a factor of several. A consequence of our findings predicts the possibility of propagating low-mass star formation, for which collapse of cores, within a mass range short of one order of magnitude, takes place sequentially tracing the wave front across a region of few tens of pc over $10^7$ years.

We present the implementation and performance of a class of directionally unsplit Riemann-solver-based hydrodynamic schemes on Graphic Processing Units (GPU). These schemes, including the MUSCL-Hancock method, a variant of the MUSCL-Hancock method, and the corner-transport-upwind method, are embedded into the adaptive-mesh-refinement (AMR) code GAMER. Furthermore, a hybrid MPI/OpenMP model is investigated, which enables the full exploitation of the computing power in a heterogeneous CPU/GPU cluster and significantly improves the overall performance. Performance benchmarks are conducted on the Dirac GPU cluster at NERSC/LBNL using up to 32 Tesla C2050 GPUs. A single GPU achieves speed-ups of 101(25) and 84(22) for uniform-mesh and AMR simulations, respectively, as compared with the performance using one(four) CPU core(s), and the excellent performance persists in multi-GPU tests. In addition, we make a direct comparison between GAMER and the widely-adopted CPU code Athena (Stone et al. 2008) in adiabatic hydrodynamic tests and demonstrate that, with the same accuracy, GAMER is able to achieve two orders of magnitude performance speed-up.