Search

Galactic outflows produced by stellar feedback are known to be multiphase in nature. Both observations and simulations indicate that the material within several kpc of galactic disk mid-planes consists of warm clouds embedded within a hot wind. A theoretical understanding of the outflow phenomenon, including both winds and fountain flows, requires study of the interactions among thermal phases. We develop a method to quantify these interactions via measurements of mass, momentum, and energy flux exchanges using temporally and spatially averaged quantities and conservation laws. We apply this method to a star-forming ISM MHD simulation based on the TIGRESS framework, for Solar neighbourhood conditions. To evaluate the extent of interactions among the phases, we first examine the validity of the ``ballistic model,'' which predicts trajectories of the warm phase ($5050\,\rm{K}5\times10^5\,\rm{K}$) to the warm phase. The large energy flux from the hot outflow that is transferred to the warm and intermediate phases is quickly radiated away. A simple interaction model implies an effective warm cloud size in the fountain flow of a few 100~pc, showing that warm-hot flux exchange mainly involves a few large clouds rather than many small ones., Comment: 22 pages, Submitted to ApJ

Supermassive black holes (SMBHs) are thought to provide energy that prevents catastrophic cooling in the centers of massive galaxies and galaxy clusters. However, it remains unclear how this "feedback" process operates. We use high-resolution optical data to study the kinematics of multi-phase filamentary structures by measuring the velocity structure function (VSF) of the filaments over a wide range of scales in the centers of three nearby galaxy clusters: Perseus, Abell 2597 and Virgo. We find that the motions of the filaments are turbulent in all three clusters studied. There is a clear correlation between features of the VSFs and the sizes of bubbles inflated by SMBH driven jets. Our study demonstrates that SMBHs are the main driver of turbulent gas motions in the centers of galaxy clusters and suggests that this turbulence is an important channel for coupling feedback to the environment. Our measured amplitude of turbulence is in good agreement with Hitomi Doppler line broadening measurement and X-ray surface brightness fluctuation analysis, suggesting that the motion of the cold filaments is well-coupled to that of the hot gas. The smallest scales we probe are comparable to the mean free path in the intracluster medium (ICM). Our direct detection of turbulence on these scales provides the clearest evidence to date that isotropic viscosity is suppressed in the weakly-collisional, magnetized intracluster plasma., Comment: 10 pages, 5 figures. Submitted to ApJL

Feedback is indispensable in galaxy formation. However, lacking resolutions, cosmological simulations often use ad hoc feedback parameters. Conversely, small-box simulations, while better resolving the feedback, cannot capture gas evolution beyond the simulation domain. We aim to bridge the gap by implementing small-box results of supernovae-driven outflows into dark matter halo-scale simulations and studying their impact on large scales. Galactic outflows are multiphase, but small-box simulations show that the hot phase (T$\approx$ 10$^{6-7}$ K) carries the majority of energy and metals. We implement hot outflows in idealized simulations of the Milky Way halo, and examine how they impact the circumgalactic medium (CGM). In this paper, we discuss the case when the star formation surface density is low and therefore the emerging hot outflows are gravitationally bound by the halo. We find that outflows form a large-scale, metal-enriched atmosphere with fountain motions. As hot gas accumulates, the inner atmosphere becomes "saturated". Cool gas condenses, with a rate balancing the injection of the hot outflows. This balance leads to a universal density profile of the hot atmosphere, independent of mass outflow rate. The atmosphere has a radially-decreasing temperature, naturally producing the observed X-ray luminosity and column densities of O VI, O VII, O VIII. The self-regulated atmosphere has a baryon and a metal mass of $(0.5-1.2)\times 10^{10}M_\odot$ and $(0.6-1.4)\times 10^8 M_\odot$, respectively, small compared to the ``missing" baryons and metals from the halo. We conjecture that the missing materials reside at even larger radii, ejected by more powerful outflows in the past., Comment: 26 pages, 18 figures, accepted by ApJ

With the development of Internet of Things (IoT) and communication technology, the number of next-generation IoT devices has increased explosively, and the delay requirement for content requests is becoming progressively higher. Fortunately, the edge-caching scheme can satisfy users' demands for low latency of content. However, the existing caching schemes are not smart enough. To address these challenges, we propose a human-like hybrid caching architecture based on the software defined edge cloud, which simultaneously considers the content popularity and the fine-grained user characteristics. Then, an optimization problem with a caching hit ratio as an optimization objective is formulated. To solve this problem, using reinforcement learning, we design a human-like hybrid caching algorithm. Extensive experiments show that compared with popular caching schemes, human-like hybrid caching schemes can improve the cache hit ratio by 20%.

Supernovae (SNe) drive multiphase galactic outflows, impacting galaxy formation; however, cosmological simulations mostly use \textit{ad hoc} feedback models for outflows, making outflow-related predictions from first principles problematic. Recent small-box simulations resolve individual SNe remnants in the interstellar medium (ISM), naturally driving outflows and permitting a determination of the wind loading factors of energy $\eta_E$, mass $\eta_m$, and metals $\eta_Z$. In this Letter, we compile small-box results, and find consensus that the hot outflows are much more powerful than the cool outflows: (i) hot outflows generally dominate the energy flux, and (ii) their specific energy $e_{s,h}$ is 10-1000 times higher than cool outflows. Moreover, the properties of hot outflows are remarkably simple: $e_{s,h} \propto \eta_{E,h}/\eta_{m,h}$ is almost invariant over four orders of magnitude of star formation surface density. Also, we find tentatively that $\eta_{E,h}/\eta_{Z,h} \sim$ 0.5. If corroborated by more simulation data, these correlations reduce the three hot phase loading factors into one. Finally, this one parameter is closely related to whether the ISM has a "breakout" condition. The narrow range of \esh\ indicates that hot outflows cannot escape dark matter halos with log $M_{\rm{halo}}\ [M_\odot] \gtrsim 12$. This mass is also where the galaxy mass-metallicity relation reaches its plateau, implying a deep connection between \textit{hot} outflows and galaxy formation. We argue that hot outflows should be included explicitly in cosmological simulations and (semi-)analytic modeling of galaxy formation., Comment: 9 pages, 4 figures, accepted by ApJL. Compilation of loading factors of supernovae-driven outflows is at https://github.com/limiao0611/loading-factors/wiki

Type Ia supernovae (SNe Ia) provide unique and important feedback in quiescent galaxies, but their impact has been underappreciated. In this paper, we analyze a series of high-resolution simulations to examine the energetics and turbulence of the medium under SNe Ia. We find that when SN remnants are resolved, their effects differ distinctly from a volumetric heating term, as is commonly assumed in unresolved simulations. First, the net heating is significantly higher than expected, by 30$\pm$10\% per cooling time. This is because a large fraction of the medium is pushed into lower densities which cool inefficiently. Second, the medium is turbulent; the root-mean-squared (RMS) velocity of the gas to 20-50 km s$^{-1}$ on a driving scale of tens of parsec. The velocity field of the medium is dominated by compressional modes, which are larger than the solenoidal components by a factor of 3-7. Third, the hot gas has a very broad density distribution. The ratio between the density fluctuations and the RMS Mach number, parameterized as $b$, is 2-20. This is in contrast to previous simulations of turbulent media, which have found $b\lesssim$ 1. The reason for the difference is mainly caused by the \textit{localized} heating of SNe Ia, which creates a large density contrast. Last, the typical length scale of a density fluctuation grows with time, forming increasingly larger bubbles and filamentary ridges. These underlying density fluctuations need to be included when X-ray observations are interpreted., Comment: 16 pages, 12 figures, accepted by ApJ. This work is part of the SMAUG project, see more information at https://www.simonsfoundation.org/flatiron/center-for-computational-astrophysics/galaxy-formation/smaug/papersplash1

A cool phase of the interstellar medium has been observed in many giant elliptical galaxies, but its origin remains unclear. We propose that uneven heating from Type Ia supernovae (SNe Ia), together with radiative cooling, can lead to the formation of the cool phase. The basic idea is that since SNe Ia explode randomly, gas parcels which are not directly heated by SN shocks will cool, forming multiphase gas. We run a series of idealized high-resolution numerical simulations, and find that cool gas develops even when the overall SNe heating rate $H$ exceeds the cooling rate $C$ by a factor as large as 1.4. We also find that the time for multiphase gas development depends on the gas temperature. When the medium has a temperature $T = 3\times 10^6$ K, the cool phase forms within one cooling time \tc; however, the cool phase formation is delayed to a few times \tc\ for higher temperatures. The main reason for the delay is turbulent mixing. Cool gas formed this way would naturally have a metallicity lower than that of the hot medium. For constant $H/C$, there is more turbulent mixing for higher temperature gas. We note that this mechanism of producing cool gas cannot be captured in cosmological simulations, which usually fail to resolve individual SN remnants., Comment: 16 pages, 11 figures, published by ApJ. This work is part of the SMAUG project, see more information at https://www.simonsfoundation.org/flatiron/center-for-computational-astrophysics/galaxy-formation/smaug/papersplash1

The tomographic Alcock-Paczynski (AP) method utilizes the redshift evolution of the AP distortion to place constraints on cosmological parameters. It has proved to be a robust method that can separate the AP signature from the redshift space distortion (RSD) effect, and deliver powerful cosmological constraints using the $\lesssim 40h^{-1}\ \rm Mpc$ clustering region. In previous works, the tomographic AP method was performed via the anisotropic 2-point correlation function statistic. In this work we consider the feasibility of conducting the analysis in the Fourier domain and examine the pros and cons of this approach. We use the integrated galaxy power spectrum (PS) as a function of direction, $\hat P_{\Delta k}(\mu)$, to quantify the magnitude of anisotropy in the large-scale structure clustering, and use its redshift variation to do the AP test. The method is tested on the large, high resolution Big-MultiDark Planck (BigMD) simulation at redshifts $z=0-1$, using the underlying true cosmology $\Omega_m=0.3071,\ w=-1$. Testing the redshift evolution of $\hat P_{\Delta k}(\mu)$ in the true cosmology and cosmologies deviating from the truth with $\delta \Omega_m=0.1,\ \delta w=0.3$, we find that the redshift evolution of the AP distortion overwhelms the effects created by the RSD by a factor of $\sim1.7-3.6$. We test the method in the range of $k\in(0.2,1.8)\ h\ \rm Mpc^{-1}$, and find that it works well throughout the entire regime. We tune the halo mass within the range $2\times 10^{13}$ to $10^{14}\ M_{\odot}$, and find that the change of halo bias results in $\lesssim 5 \%$ change in $\hat P_{\Delta k}(\mu)$, which is less significant compared with the cosmological effect. Our work shows that it is feasible to conduct the tomographic AP analysis in the Fourier space., Comment: 7 pages, 7 figures

The circumgalactic medium (CGM) is closely linked to galaxy formation and evolution, but difficult to characterize observationally and typically poorly resolved in cosmological simulations. We use spherically-symmetric, idealized, high-resolution simulations of the CGM in $10^{12}M_\odot$ and $10^{11}M_\odot$ dark matter halos to characterize the gas pressure, turbulent and radial velocities, and degree of thermal and effective dynamic pressure support in the overall CGM as well as in its high- and low-temperature phases. We find that the $10^{12}M_\odot$ halo contains a CGM mostly formed of a hot gas halo in hydrostatic equilibrium out of which cold gas condenses and falls onto the central galaxy, while the $10^{11}M_\odot$ halo's CGM is not in hydrostatic equilibrium, has a wider spread of properties at a given galactocentric radius, does not have a clear separation of hot and cold phases, and is dominated by bulk motions. We also find that the degree of pressure support in the $10^{11}M_\odot$ halo is strongly dependent on the parameters of the galactic winds of the central galaxy. These results promote the idea that there is no "average" CGM and care must be taken when setting the initial conditions for a small-box simulation of a patch of the CGM., Comment: 20 pages, 17 figures, accepted by MNRAS

We model the ultraviolet spectra of the Seyfert 1 galaxy NGC~5548 obtained with the Hubble Space Telescope during the 6-month reverberation-mapping campaign in 2014. Our model of the emission from NGC 5548 corrects for overlying absorption and deblends the individual emission lines. Using the modeled spectra, we measure the response to continuum variations for the deblended and absorption-corrected individual broad emission lines, the velocity-dependent profiles of Ly$\alpha$ and C IV, and the narrow and broad intrinsic absorption features. We find that the time lags for the corrected emission lines are comparable to those for the original data. The velocity-binned lag profiles of Ly$\alpha$ and C IV have a double-peaked structure indicative of a truncated Keplerian disk. The narrow absorption lines show delayed response to continuum variations corresponding to recombination in gas with a density of $\sim 10^5~\rm cm^{-3}$. The high-ionization narrow absorption lines decorrelate from continuum variations during the same period as the broad emission lines. Analyzing the response of these absorption lines during this period shows that the ionizing flux is diminished in strength relative to the far-ultraviolet continuum. The broad absorption lines associated with the X-ray obscurer decrease in strength during this same time interval. The appearance of X-ray obscuration in $\sim\,2012$ corresponds with an increase in the luminosity of NGC 5548 following an extended low state. We suggest that the obscurer is a disk wind triggered by the brightening of NGC 5548 following the decrease in size of the broad-line region during the preceding low-luminosity state., Comment: 50 pages, 30 figures, uses aastex62.cls. Accepted for publication in ApJ, 07/06/2019. High-level products page in MAST will go live after 7/15/2019. Replaced Figure 4 on 7/12/2019 to be more red/green color-blind friendly

In response to the need for the Astro2020 Decadal Survey to explicitly engage early career astronomers, the National Academies of Sciences, Engineering, and Medicine hosted the Early Career Astronomer and Astrophysicist Focus Session (ECFS) on October 8-9, 2018 under the auspices of Committee of Astronomy and Astrophysics. The meeting was attended by fifty six pre-tenure faculty, research scientists, postdoctoral scholars, and senior graduate students, as well as eight former decadal survey committee members, who acted as facilitators. The event was designed to educate early career astronomers about the decadal survey process, to solicit their feedback on the role that early career astronomers should play in Astro2020, and to provide a forum for the discussion of a wide range of topics regarding the astrophysics career path. This white paper presents highlights and themes that emerged during two days of discussion. In Section 1, we discuss concerns that emerged regarding the coming decade and the astrophysics career path, as well as specific recommendations from participants regarding how to address them. We have organized these concerns and suggestions into five broad themes. These include (sequentially): (1) adequately training astronomers in the statistical and computational techniques necessary in an era of "big data", (2) responses to the growth of collaborations and telescopes, (3) concerns about the adequacy of graduate and postdoctoral training, (4) the need for improvements in equity and inclusion in astronomy, and (5) smoothing and facilitating transitions between early career stages. Section 2 is focused on ideas regarding the decadal survey itself, including: incorporating early career voices, ensuring diverse input from a variety of stakeholders, and successfully and broadly disseminating the results of the survey., Comment: 9 pages; Astro2020 APC White Paper: State of the Profession Consideration

The tomographic Alcock-Paczynski (AP) test is a robust large-scale structure (LSS) measurement that receives little contamination from the redshift space distortion (RSD). It has placed tight cosmological constraints by using small and intermediate clustering scales of the LSS data. However, previous works have neglected the cross-correlation among different redshift bins, which could cause the statistical uncertainty being underestimated by $\sim$20\%. In this work, we further improve this method by including this multi-redshifts full correlation. We apply it to the SDSS DR12 galaxies sample and find out that, for $\Lambda$CDM, the combination of AP with the Planck+BAO dataset slightly reduces (within 1-$\sigma$) $\Omega_m$ to $0.304\pm0.007$ (68.3\% CL). This then leads to a larger $H_0$ and also mildly affects $\Omega_b h^2$, $n_s$ and the derived parameters $z_*$, $r_*$, $z_{re}$ but not $\tau$, $A_s$ and $\sigma_8$. For the flat $w$CDM model, our measurement gives $\Omega_m=0.301\pm 0.010$ and $w=-1.090\pm 0.047$, where the additional AP measurement reduces the error budget by $\sim 25\%$. When including more parameters into the analysis, the AP method also improves the constraints on $\Omega_k$, $\sum m_\mu$, $N_{\rm eff}$ by $20-30\%$. Early universe parameters such as $dn_s/d{\rm ln}k$ and $r$, however, are unaffected. Assuming the dark energy equation of state $w=w_0+w_a \frac{z}{1+z}$, the Planck+BAO+SNIa+$H_0$+AP datasets prefer a dynamical dark energy at $\approx1.5 \sigma$ CL. Finally, we forecast the cosmological constraints expected from the DESI galaxy survey and find that combining AP with CMB+BAO method would improve the $w_0$-$w_a$ constraint by a factor of $\sim 10$., Comment: 13 pages, 12 figures, 3 tables. Accepted for publication in ApJ

We find an interesting fact that fictitious retrograde co-orbitals of Saturn, or small bodies inside the retrograde 1:1 resonance with Saturn, are highly unstable in our numerical simulations. It is shown that in the presence of Jupiter, the retrograde co-orbitals will get ejected from Saturn's co-orbital space within a timescale of 10 Myr. This scenario reminds us of the instability of Saturn Trojans caused by both the Great Inequality and the secular resonances. Therefore, we carry out in-depth inspections on both mechanisms and prove that the retrograde resonance overlap, raised by Great Inequality, cannot serve as an explanation for the instability of retrograde co-orbitals, due to the weakness of the retrograde 2:5 resonance with Jupiter at a low eccentricity. However, we discover that both $\nu_5$ and $\nu_6$ secular resonances contribute to the slow growth of the eccentricity, therefore, are possibly the primary causes of the instability inside Saturn's retrograde co-orbital space., Comment: 6 pages, 6 figures, accepted for publication in MNRAS

The interaction between inflowing gas clouds and galactic outflows at the interface where the galactic disk transitions into the circumgalactic medium is an important process in galaxy fueling, yet remains poorly understood. Using a series of tall-box hydrodynamic ENZO simulations, we have studied the interaction between smooth gas inflow and supernovae-driven outflow at the disk-halo interface with pc-scale resolution. A realistic wind of outflowing material is generated by supernovae explosions in the disk, while inflowing gas is injected at the top boundary of the simulation box with an injection velocity ranging from $10-100 \rm \ km \ s^{-1}$. We find that cooling and hydrodynamic instabilities drive the injected gas to fragment into cold ($\sim 10^{3}$ K) cloud clumps with typical densities of $\sim 1 \rm \ cm^{-3}$. These clumps initially accelerate before interacting and partially mixing with the outflow and decelerating to velocities in the 50-100 $\rm km \ s^{-1}$ range. When the gas clumps hit the disk, $10\%-50 \%$ of the injected material is able to accrete (depending on the injection velocity). Clumps originating from gas injected with a higher initial velocity approach the disk with greater ram pressure, allowing them to penetrate through the disk in low density regions. We use (equilibrium) CLOUDY photoionization models to generate absorption and emission signatures of gas accretion, finding that our mock HI and H$\alpha$ observables are prominent and generally consistent with measurements in the Milky Way. We do not predict enhanced emission/absorption for higher ionization states such as OVI., Comment: 26 pages, 14 figures, 1 table, Accepted for publication in ApJ, Simulation movies available at https://doi.org/10.7916/d8-vgnn-7n73

In this paper, we consider the Cauchy problem for the fractional Schr\"odinger equation $i D_t^\alpha u + (-\Delta)^{\frac{\beta}{2}} u =0$ with $0<\alpha<1$, $\beta>0$. We establish the dispersive estimates for the solutions. In particular, we prove that the decay rates are sharp., Comment: 11 pages

The semilinear space-time fractional Schr\"odinger equation is considered. First, we give the explicit form for the fundamental solutions by using the Fox $H$-functions in order to to establish some $L^s$ decay estimates. After that, we give some space-time estimates for the mild solutions from which the local well-posedness is derived on some proper Banach space., Comment: 23pages