We examine the properties of the circumgalactic medium (CGM) at low redshift in a range of simulated Milky Way mass halos. The sample is comprised of seven idealized simulations, an adaptive mesh refinement cosmological zoom-in simulation, and two groups of 50 halos with star forming or quiescent galaxies taken from the IllustrisTNG100 simulation. The simulations have very different setups, resolution, and feedback models, but are analyzed in a uniform manner. By comparing median radial profiles and mass distributions of CGM properties, we isolate key similarities and differences. In doing so, we advance the efforts of the SMAUG (Simulating Multiscale Astrophysics to Understand Galaxies) project that aims to understand the inherently multiscale galaxy formation process. In the cosmological simulations, the CGM exhibits nearly flat temperature distributions, and broad pressure and radial velocity distributions. In the idealized simulations, similar distributions are found in the inner CGM ($\lesssim 0.5 \, r_{\rm 200c}$) when strong galactic feedback models are employed, but the outer CGM ($\gtrsim 0.5 \, r_{\rm 200c}$) has a much less prominent cold phase, and narrower pressure and velocity distributions even in models with strong feedback. This comparative analysis demonstrates the dominant role feedback plays in shaping the inner CGM and the increased importance of cosmological effects, such as nonspherical accretion and satellite galaxies, in the outer CGM. Furthermore, our findings highlight that while cosmological simulations are required to capture the multiphase structure of the CGM at large radii, idealized simulations provide a robust framework to study how galactic feedback interacts with the inner CGM and thereby provide a reliable avenue to constrain feedback prescriptions.

We present Powderday, a flexible, fast, open-source dust radiative transfer package designed to interface with galaxy formation simulations. Powderday builds on FSPS population synthesis models, Hyperion dust radiative transfer, and employs yt to interface between different software packages. We include our stellar population synthesis modeling on the fly, which allows for significant run-time flexibility in the assumed stellar physics. We include a model for nebular line emission that can employ either precomputed Cloudy lookup tables (for efficiency), or direct photoionization calculations for all young stars (for flexibility). The dust content follows either observationally-motivated prescriptions, direct modeling from galaxy formation simulations, or a novel approach that includes the dust content via learning-based algorithms from the SIMBA cosmological galaxy formation simulation. AGN can additionally be included via a range of prescriptions. The output of these models are broadband SEDs, as well as filter-convolved images. Powderday is designed to eliminate last-mile efforts by researchers that employ different hydrodynamic galaxy formation models, and seamlessly interfaces with GIZMO, AREPO, GASOLINE, CHANGA, and ENZO. We demonstrate the capabilities of the code via three applications: a model for the star formation rate (SFR) - infrared luminosity relation in galaxies (including the impact of AGN); the impact of circumstellar dust around AGB stars on the mid-infrared emission from galaxy SEDs; and the impact of galaxy inclination angle on dust attenuation laws.

(Abridged) We present an investigation of kinematical imprints of AGN feedback on the Warm Ionized gas Medium (WIM) of massive early-type galaxies (ETGs). To this end, we take a two-fold approach that involves a comparative analysis of Halpha velocity fields in 123 local ETGs from the CALIFA integral field spectroscopy survey with 20 simulated galaxies from high-resolution hydrodynamic cosmological SPHgal simulations. The latter were re-simulated for two modeling setups, one with and another without AGN feedback. In order to quantify the effects of AGN feedback on gas kinematics we measure three parameters that probe deviations from simple regular rotation using the kinemetry package. These indicators trace the possible presence of distinct kinematic components in Fourier space (k3,5/k1), variations in the radial profile of the kinematic major axis (sigma_PA), and offsets between the stellar and gas velocity fields (Delta Phi). These quantities are monitored in the simulations from a redshift 3 to 0.2 to assess the connection between black hole accretion history, stellar mass growth and kinematical perturbation of the WIM. Observed local massive galaxies show a broad range of irregularities, indicating disturbed warm gas motions, irrespective of being classified via diagnostic lines as AGN or not. Simulations of massive galaxies with AGN feedback generally exhibit higher irregularity parameters than without AGN feedback, more consistent with observations. Besides AGN feedback, other processes like major merger events or infalling gas clouds can lead to elevated irregularity parameters, but they are typically of shorter duration. More specifically, k3,5/k1 is most sensitive to AGN feedback, whereas Delta Phi is most strongly affected by gas infall.

The faintness of satellite systems in galaxy groups has contributed to the widely discussed "missing satellite" and "too big to fail" issues. Using techniques based on Tremaine & Richstone (1977), we show that there is no problem with the luminosity function computed from modern codes per se, but that the gap between first and second brightest systems is too big {\it given} the luminosity function, that the same large gap is found in modern, large scale baryonic $\Lambda$CDM simulations such as EAGLE and IllustrisTNG, is even greater in dark matter only simulations, and finally, that this is most likely due to gravitationally induced merging caused by classical dynamical friction. Quantitatively the gap is larger in the computed simulations than in the randomized ones by $1.79 \pm 1.04$, $1.51 \pm 0.93$, $3.43 \pm 1.44$ and $3.33 \pm 1.35$ magnitudes in the EAGLE, IllustrisTNG, and dark matter only simulations of EAGLE and IllustrisTNG respectively. Furthermore the anomalous gaps in the simulated systems are even larger than in the real data by over half a magnitude and are still larger in the dark matter only simulations. Briefly stated, $\Lambda$CDM does not have a problem with an absence of "too big to fail" galaxies. Statistically significant large gaps between first and second brightest galaxies are to be expected.

We explore the effect of active galactic nucleus (AGN) feedback from central galaxies on their satellites by comparing two sets of cosmological zoom-in runs of 27 halos with masses ranging from $10^{12}$ to $10^{13.4}$ solar masses at z=0, with (wAGN) and without (noAGN) AGN feedback. Both simulations include stellar feedback from multiple processes, including powerful winds from supernovae, stellar winds from young massive stars, AGB stars, radiative heating within Str\"omgren spheres and photoelectric heating. Our wAGN model is identical to the noAGN model except that it also includes a model for black hole seeding and accretion, as well as AGN feedback via high-velocity broad absorption line winds and Compton/photoionization heating. We show that the inclusion of AGN feedback from the central galaxy significantly affects the star formation history and the gas content of the satellite galaxies. AGN feedback starts to affect the gas content and the star formation of the satellites as early as z=2. The mean gas rich fraction of satellites at z=0 decreases from 15% in the noAGN simulation to 5% in the wAGN simulation. The difference between the two sets extends as far out as five times the virial radius of the central galaxy at z=1. We investigate the quenching mechanism by studying the physical conditions in the surroundings of pairs of satellites matched across the wAGN and noAGN simulations and find an increase in the temperature and relative velocity of the intergalactic gas.

We add a suite of chemical abundances to the MACER (Massive AGN Controlled Ellipticals Resolved) 2D code, by solving 12 additional continuity equations for H, He, C, N, O, Ne, Mg, Si, S, Ca, Fe and Ni respectively with sources from AGB stars and supernovae of type Ia and II with metal yields based on standard stellar physics. New stars, formed in Toomre unstable circumnuclear disks (of a size $\lesssim150$ parsec), are assumed to have a top-heavy initial mass function with a power index of 1.65. The metal dilution effects due to cosmic accretion are also included. With the high resolution of few parsecs in central regions, resolved black hole accretion and AGN feedback, we can track the metal enrichment, transportation and dilution throughout the modeled massive elliptical galaxy of velocity dispersion $\sim280$ km/s. We retrieve the chemical composition of the BAL winds launched by the central AGN, synthesize the X-ray features of the hot ISM, and find that (1) the simulated metallicity in the BAL winds could be up to $\sim 8 Z_\odot$, while that of the hot ISM in the host galaxy is $\sim 2.3 Z_\odot$, matching well with SDSS observations of BLR gas; (2) the X-ray emitting hot gas is metal enriched with a typical value $\sim2.5 Z_\odot$; (3) the circumunuclear cold gas disk, where the metals are condensed, further enriched and recycled, plays a critical role in the metal enrichment; (4) the black hole accretion rate $\dot{M}_{\rm BH}$ linearly correlates with the star formation rate $\dot{M}_\star^+$ in the circumnuclear disk, i.e, $\dot{M}_\star^+ \sim 7.7 \dot{M}_{\rm BH}$, but lagged in time by roughly $10^6$ years.

We study the connection between active galactic nuclei (AGN) and their host galaxies through cosmic time in the large-scale cosmological IllustrisTNG simulations. We first compare BH properties, i.e. the hard X-ray BH luminosity function, AGN galaxy occupation fraction, and distribution of Eddington ratios, to available observational constraints. The simulations produce a population of BHs in good agreement with observations, but we note an excess of faint AGN in hard X-ray (L_x ~ 10^{43-44} erg/s), and a lower number of bright AGN (L_x>10^{44} erg/s), a conclusion that varies quantitatively but not qualitatively with BH luminosity estimation method. The lower Eddington ratios of the 10^{9} Msun BHs compared to observations suggest that AGN feedback may be too efficient in this regime. We study galaxy star formation activity and structural properties, and design sample-dependent criteria to identify different galaxy types (star-forming/quiescent, extended/compact) that we apply both to the simulations and observations from the candels fields. We analyze how the simulated and observed galaxies populate the specific star formation rate - stellar mass surface density diagram. A large fraction of the z=0 M_{star}>10^{11} Msun quiescent galaxies first experienced a compaction phase (i.e. reduction of galaxy size) while still forming stars, and then a quenching event. We measure the dependence of AGN fraction on galaxies' locations in this diagram. After correcting the simulations with a redshift and AGN luminosity-dependent model for AGN obscuration, we find good qualitative and quantitative agreement with observations. The AGN fraction is the highest among compact star-forming galaxies (16-20% at z~1.5-2), and the lowest among compact quiescent galaxies (6-10% at z~1.5-2).

We present a series of 20 cosmological zoom simulations of the formation of massive galaxies with and without a model for AGN feedback. Differences in stellar population and kinematic properties are evaluated by constructing mock integral field unit (IFU) maps. The impact of the AGN is weak at high redshift when all systems are mostly fast-rotating and disc-like. After $z \sim 1$ the AGN simulations result in lower mass, older, less metal rich and slower rotating systems with less disky isophotes - in general agreement with observations. Two-dimensional kinematic maps of in-situ and accreted stars show that these differences result from reduced in-situ star formation due to AGN feedback. A full analysis of stellar orbits indicates that galaxies simulated with AGN are typically more triaxial and have higher fractions of x-tubes and box orbits and lower fractions of z-tubes. This trend can also be explained by reduced late in-situ star formation. We introduce a global parameter, $\xi_3$ , to characterise the anti-correlation between the third-order kinematic moment $h_3$ and the line-of-sight velocity ($v_{los}/{\sigma}$), and compare to ATLAS$^{3D}$ observations. The kinematic asymmetry parameter $\xi_3$ might be a useful diagnostic for large integral field surveys as it is a kinematic indicator for intrinsic shape and orbital content.

We compute synthetic optical and ultraviolet (UV) emission-line properties of galaxies in a full cosmological framework by coupling, in post-processing, new-generation nebular-emission models with high-resolution, cosmological zoom-in simulations of massive galaxies. Our self-consistent modelling accounts for nebular emission from young stars and accreting black holes (BHs). We investigate which optical- and UV-line diagnostic diagrams can best help to discern between the main ionizing sources, as traced by the ratio of BH accretion to star formation rates in model galaxies, over a wide range of redshifts. At low redshift, simulated star-forming galaxies, galaxies dominated by active galactic nuclei and composite galaxies are appropriately differentiated by standard selection criteria in the classical [OIII]$\lambda$5007/H$\beta$ versus [NII]$\lambda$6584/H$\alpha$ diagram. At redshifts $z \gt 1$, however, this optical diagram fails to discriminate between active and inactive galaxies at metallicities below $0.5\ Z_\odot$. To robustly classify the ionizing radiation of such metal-poor galaxies, which dominate in the early Universe, we confirm 3 previous, and propose 11 novel diagnostic diagrams based on equivalent widths and luminosity ratios of UV emission lines, such as EW(OIII]$\lambda$1663) versus OIII]$\lambda$1663/HeII$\lambda$1640, CIII]$\lambda$1908/HeII$\lambda$1640 versus OIII]$\lambda$1663/HeII$\lambda$1640, and CIV$\lambda$1550/CIII]$\lambda$1908 versus CIII]$\lambda$1908/CII$\lambda$2326. We formulate associated UV selection criteria and discuss some caveats of our results (e.g., uncertainties in the modelling of the HeII$\lambda$1640 line). These UV diagnostic diagrams are potentially important for the interpretation of high-quality spectra of very distant galaxies to be gathered by next-generation telescopes, such as the James Webb Space Telescope.

The near-Earth asteroid 3200 Phaethon (1983 TB) is an attractive object not only from a scientific viewpoint but also because of JAXA's DESTINY+ target. The rotational lightcurve and spin properties were investigated based on the data obtained in the ground-based observation campaign of Phaethon. We aim to refine the lightcurves and shape model of Phaethon using all available lightcurve datasets obtained via optical observation, as well as our time-series observation data from the 2017 apparition. Using eight 1-2-m telescopes and an optical imager, we acquired the optical lightcurves and derived the spin parameters of Phaethon. We applied the lightcurve inversion method and SAGE algorithm to deduce the convex and non-convex shape model and pole orientations. We analysed the optical lightcurve of Phaethon and derived a synodic and a sidereal rotational period of 3.6039 h, with an axis ratio of a/b = 1.07. The ecliptic longitude (lambda) and latitude (beta) of the pole orientation were determined as (308, -52) and (322, -40) via two independent methods. A non-convex model from the SAGE method, which exhibits a concavity feature, is also presented.

We use cosmological hydrodynamical simulations to investigate the role of feedback from accreting black holes on the evolution of sizes, compactness, stellar core density and specific star-formation of massive galaxies with stellar masses of $ M_{star} > 10^{10.9} M_{\odot}$. We perform two sets of cosmological zoom-in simulations of 30 halos to z=0: (1) without black holes and Active Galactic Nucleus (AGN) feedback and (2) with AGN feedback arising from winds and X-ray radiation. We find that AGN feedback can alter the stellar density distribution, reduce the core density within the central 1 kpc by 0.3 dex from z=1, and enhance the size growth of massive galaxies. We also find that galaxies simulated with AGN feedback evolve along similar tracks to those characterized by observations in specific star formation versus compactness. We confirm that AGN feedback plays an important role in transforming galaxies from blue compact galaxies into red extended galaxies in two ways: (1) it effectively quenches the star formation, transforming blue compact galaxies into compact quiescent galaxies and (2) it also removes and prevents new accretion of cold gas, shutting down in-situ star formation and causing subsequent mergers to be gas-poor or mixed. Gas poor minor mergers then build up an extended stellar envelope. AGN feedback also puffs up the central region through the fast AGN driven winds as well as the slow expulsion of gas while the black hole is quiescent. Without AGN feedback, large amounts of gas accumulate in the central region, triggering star formation and leading to overly massive blue galaxies with dense stellar cores.

A tightly correlated star formation rate-stellar mass relation of star forming galaxies, or star-forming sequence (SFS), is a key feature in galaxy property-space that is predicted by modern galaxy formation models. We present a flexible data-driven approach for identifying this SFS over a wide range of star formation rates and stellar masses using Gaussian mixture modeling (GMM). Using this method, we present a consistent comparison of the $z=0$ SFSs of central galaxies in the Illustris, EAGLE, and Mufasa hydrodynamic simulations and the Santa Cruz semi-analytic model (SC-SAM), alongside data from the Sloan Digital Sky Survey. We find, surprisingly, that the amplitude of the SFS varies by up to ${\sim} 0.7\,\mathrm{dex}$ (factor of ${\sim} 5$) among the simulations with power-law slopes range from $0.7$ to $1.2$. In addition to the SFS, our GMM method also identifies sub-components in the star formation rate-stellar mass relation corresponding to star-burst, transitioning, and quiescent sub-populations. The hydrodynamic simulations are similarly dominated by SFS and quiescent sub-populations unlike the SC-SAM, which predicts substantial fractions of transitioning and star-burst galaxies at stellar masses above and below $10^{10} M_\odot$, respectively. All of the simulations also produce an abundance of low-mass quiescent central galaxies in apparent tension with observations. These results illustrate that, even among models that well reproduce many observables of the galaxy population, the $z=0$ SFS and other sub-populations still show marked differences that can provide strong constraints on galaxy formation models.

We explore the effect of momentum-driven winds representing radiation pressure driven outflows from accretion onto supermassive black holes in a set of numerical hydrodynamical simulations. We explore two matched sets of cosmological zoom-in runs of 24 halos with masses ~$10^{12.0}-10^{13.4}$ M_sun run with two different feedback models. Our `NoAGN' model includes stellar feedback via UV heating, stellar winds and supernovae, photoelectric heating and cosmic X-ray background heating from a meta-galactic background. Our fiducial `MrAGN' model is identical except that it also includes a model for black hole seeding and accretion, as well as heating and momentum injection associated with the radiation from black hole accretion. Our MrAGN model launches galactic outflows which result in both `ejective' feedback - the outflows themselves which drive gas out of galaxies - and `preventative' feedback, which suppresses the inflow of new and recycling gas. As much as 80 % of outflowing galactic gas can be expelled, and accretion can be suppressed by as much as a factor of 30 in the MrAGN runs when compared with the NoAGN runs. The histories of NoAGN galaxies are recycling-dominated, with ~70% of material that leaves the galaxy eventually returning, and the majority of outflowing gas re-accretes on 1 Gyr timescales without AGN feedback. Outflowing gas in the MrAGN runs has higher characteristic velocity (500 - 1,000 km/s versus 100-300 km/s for outflowing NoAGN gas) and travels as far as a few Mpcs. Only ~10% of ejected material is re-accreted in the MrAGN galaxies.

Single-photon detectors in space must retain useful performance characteristics despite being bombarded with sub-atomic particles. Mitigating the effects of this space radiation is vital to enabling new space applications which require high-fidelity single-photon detection. To this end, we conducted proton radiation tests of various models of avalanche photodiodes (APDs) and one model of photomultiplier tube potentially suitable for satellite-based quantum communications. The samples were irradiated with 106 MeV protons at doses approximately equivalent to lifetimes of 0.6 , 6, 12 and 24 months in a low-Earth polar orbit. Although most detection properties were preserved, including efficiency, timing jitter and afterpulsing probability, all APD samples demonstrated significant increases in dark count rate (DCR) due to radiation-induced damage, many orders of magnitude higher than the 200 counts per second (cps) required for ground-to-satellite quantum communications. We then successfully demonstrated the mitigation of this DCR degradation through the use of deep cooling, to as low as -86 degrees C. This achieved DCR below the required 200 cps over the 24 months orbit duration. DCR was further reduced by thermal annealing at temperatures of +50 to +100 degrees C.

Galaxies occupy different regions of the [OIII]$\lambda5007$/H$\beta$-versus-[NII]$\lambda6584$/H$\alpha$ emission-line ratio diagram in the distant and local Universe. We investigate the origin of this intriguing result by modelling self-consistently, for the first time, nebular emission from young stars, accreting black holes (BHs) and older, post-asymptotic-giant-branch (post-AGB) stellar populations in galaxy formation simulations in a full cosmological context. In post-processing, we couple new-generation nebular-emission models with high-resolution, cosmological zoom-in simulations of massive galaxies to explore which galaxy physical properties drive the cosmic evolution of the optical-line ratios [OIII]$\lambda5007$/H$\beta$, [NII]$\lambda6584$/H$\alpha$, [SII]$\lambda\lambda6717,6731$/H$\alpha$ and [OI]$\lambda6300$/H$\alpha$. The line ratios of simulated galaxies agree well with observations of both star-forming and active local SDSS galaxies. Towards higher redshifts, at fixed galaxy stellar mass, the average [OIII]/H$\beta$ increases and [NII]/H$\alpha$, [SII]/H$\alpha$ and [OI]/H$\alpha$ decrease -- widely consistent with observations. At fixed stellar mass, we identify star formation history, which controls nebular emission from young stars via the ionization parameter, as the primary driver of the cosmic evolution of [OIII]/H$\beta$ and [NII]/H$\alpha$. For [SII]/H$\alpha$ and [OI]/H$\alpha$, this applies only to redshifts above $z=1.5$, the evolution at lower redshift being driven in roughly equal parts by nebular emission from AGN and post-AGB stars. Instead, changes in the hardness of ionizing radiation, ionized-gas density, the prevalence of BH accretion relative to star formation and the dust-to-metal mass ratio (whose impact on the gas-phase N/O ratio we model at fixed O/H) play at most a minor role in the cosmic evolution of simulated galaxy line ratios.

We explore observational and theoretical constraints on how galaxies might transition between the "star-forming main sequence" (SFMS) and varying "degrees of quiescence" out to $z=3$. Our analysis is focused on galaxies with stellar mass $M_*>10^{10}M_{\odot}$, and is enabled by GAMA and CANDELS observations, a semi-analytic model (SAM) of galaxy formation, and a cosmological hydrodynamical "zoom in" simulation with momentum-driven AGN feedback. In both the observations and the SAM, transition galaxies tend to have intermediate S\'ersic indices, half-light radii, and surface stellar mass densities compared to star-forming and quiescent galaxies out to $z=3$. We place an observational upper limit on the average population transition timescale as a function of redshift, finding that the average high-redshift galaxy is on a "fast track" for quenching whereas the average low-redshift galaxy is on a "slow track" for quenching. We qualitatively identify four physical origin scenarios for transition galaxies in the SAM: oscillations on the SFMS, slow quenching, fast quenching, and rejuvenation. Quenching timescales in both the SAM and the hydrodynamical simulation are not fast enough to reproduce the quiescent population that we observe at $z\sim3$. In the SAM, we do not find a clear-cut morphological dependence of quenching timescales, but we do predict that the mean stellar ages, cold gas fractions, SMBH masses, and halo masses of transition galaxies tend to be intermediate relative to those of star-forming and quiescent galaxies at $z<3$.

We ask how the inclusion of various physical heating processes due to the metal content of gas affect the evolution of massive galaxies and compute a suite of cosmological hydrodynamical simulations that follow these systems and their supermassive black holes. We use a smoothed particle hydrodynamics code with a pressure-entropy formulation and a more accurate treatment of the metal production, turbulent diffusion and cooling rate based on individual element abundances. The feedback models include (1) AGN feedback via high velocity BAL winds and Compton/photoionization heating, (2) explicit stellar feedback from multiple processes including powerful winds from supernova events, stellar winds from young massive stars and AGB stars as well as radiative heating within Stromgren spheres around massive stars, and (3) additional heating effects due to the presence of metals including grain photoelectric heating, metallicity dependent X-ray heating by nearby accreting black holes and from the cosmic X-ray background, which are the major improvement in our feedback model. With a suite of zoom-in simulations of 30 halos with $M_{vir} \sim 10^{12-13.4}$, we show that energy and momentum budget from all feedback effects generate realistic galaxy properties. We explore the detailed role of each feedback model with three additional sets of simulations with varying input physics. We show that the metal induced heating mechanisms reduce the fraction of accreted stellar material by mainly suppressing the growth of diffuse small stellar systems at high redshift but overall have a relatively minor effect on the final stellar and gas properties of massive galaxies. The inclusion of AGN feedback significantly improves the ability of our cosmological simulations to yield realistic gas and stellar properties of massive galaxies with reasonable fraction of the final stellar mass accreted from other galaxies.

We have performed the first three-dimensional relativistic hydrodynamic simulations of extragalactic jets of pure leptonic and baryonic plasma compositions propagating into a hydrostatic intracluster medium environment. The numerical simulations use a general equation of state for a multi-component relativistic gas which closely reproduces the Synge equation of state for a relativistic perfect gas. We find that morphological and dynamical differences between leptonic and baryonic jets are much less evident than those between hot and cold jets. In all these models, the jets first propagate with essentially constant velocities within the core radius of the intracluster medium and then accelerate progressively so as to increase the jet advance velocity by a factor of between 1.2 and 1.6 at the end of simulations, depending upon the models. The temporal evolution of the average cavity pressure is not consistent with that expected by the extended theoretical model even if the average cavity pressure decreases as a function of time with a power-law. Our simulations produce synthetic radio images that are dominated by bright hot spots and appear similar to observations of the extended radio galaxies with collimated radio jets. These bright radio lobes would be visible as dark regions in X-ray images and are morphologically similar to observed X-ray cavities in the intracluster medium. This supports the expectation that the bow shock surrounding the head of the jet is important mechanism for producing X-ray cavities in the intracluster medium. Although there are quantitative differences among the models the total radio and X-ray intensity curves show qualitatively similar trends in all of them.

We present three-dimensional hydrodynamical simulations showing the effect of kinetic and radiative AGN feedback on a model galaxy representing a massive quiescent low-redshift early-type galaxy of $M_* = 8.41\times 10^{10} M_\odot$, harbouring a $M_\mathrm{BH} = 4\times 10^8 M_\odot $ black hole surrounded by a cooling gaseous halo. We show that, for a total baryon fraction of $\sim 20\%$ of the cosmological value, feedback from the AGN can keep the galaxy quiescent for about 4.35 Gyr and with properties consistent with black hole mass and X-ray luminosity scaling relations. However, this can only be achieved if the AGN feedback model includes both kinetic and radiative feedback modes. The simulation with only kinetic feedback fails to keep the model galaxy fully quiescent, while one with only radiative feedback leads to excessive black-hole growth. For higher baryon fractions (e.g. 50\% of the cosmological value), the X-ray luminosities exceed observed values by at least one order of magnitude, and rapid cooling results in a star-forming galaxy. The AGN plays a major role in keeping the circumgalactic gas at observed metallicities of $Z/Z_\odot \gtrsim 0.3 $ within the central $\sim 30$ kpc by venting nuclear gas enriched with metals from residual star formation activity. As indicated by previous cosmological simulations, our results are consistent with a model for which the black hole mass and the total baryon fraction are set at higher redshifts $z > 1$ and the AGN alone can keep the model galaxy on observed scaling relations. Models without AGN feedback violate both the quiescence criterion as well as CGM metallicity constraints.

Various heuristic approaches to model unresolved supernova (SN) feedback in galaxy formation simulations exist to reproduce the formation of spiral galaxies and the overall inefficient conversion of gas into stars. Some models, however, require resolution dependent scalings. We present a sub-resolution model representing the three major phases of supernova blast wave evolution $-$free expansion, energy conserving Sedov-Taylor, and momentum conserving snowplow$-$ with energy scalings adopted from high-resolution interstellar-medium simulations in both uniform and multiphase media. We allow for the effects of significantly enhanced SN remnant propagation in a multiphase medium with the cooling radius scaling with the hot volume fraction, $f_{\mathrm{hot}}$, as $(1 - f_{\mathrm{hot}})^{-4/5}$. We also include winds from young massive stars and AGB stars, Str\"omgren sphere gas heating by massive stars, and a gas cooling limiting mechanism driven by radiative recombination of dense HII regions. We present initial tests for isolated Milky-Way like systems simulated with the GADGET based code SPHgal with improved SPH prescription. Compared to pure thermal SN input, the model significantly suppresses star formation at early epochs, with star formation extended both in time and space in better accord with observations. Compared to models with pure thermal SN feedback, the age at which half the stellar mass is assembled increases by a factor of 2.4, and the mass loading parameter and gas outflow rate from the galactic disk increase by a factor of 2. Simulation results are converged for a two order of magnitude variation in particle mass in the range (1.3$-$130)$\times 10^4$ solar masses.

We employ cosmological hydrodynamical simulations to investigate the effects of AGN feedback on the formation of massive galaxies with present-day stellar masses of $M_{stel} = 8.8 \times 10^{10} - 6.0 \times 10^{11} M_{sun}$. Using smoothed particle hydrodynamics simulations with a pressure-entropy formulation that allows an improved treatment of contact discontinuities and fluid mixing, we run three sets of simulations of 20 halos with different AGN feedback models: (1) no feedback, (2) thermal feedback, and (3) mechanical and radiation feedback. We assume that seed black holes are present at early cosmic epochs at the centre of emerging dark matter halos and trace their mass growth via gas accretion and mergers with other black holes. Both feedback models successfully recover the observed M_BH - sigma relation and black hole-to-stellar mass ratio for simulated central early-type galaxies. The baryonic conversion efficiencies are reduced by a factor of two compared to models without any AGN feedback at all halo masses. However, massive galaxies simulated with thermal AGN feedback show a factor of ~10-100 higher X-ray luminosities than observed. The mechanical/radiation feedback model reproduces the observed correlation between X-ray luminosities and velocity dispersion, e.g. for galaxies with sigma = 200 km/s, the X-ray luminosity is reduced from $10^{42}$ erg/s to $10^{40}$ erg/s. It also efficiently suppresses late time star formation, reducing the specific star formation rate from $10^{-10.5}$ $yr^{-1}$ to $10^{-14}$ $yr^{-1}$ on average and resulting in quiescent galaxies since z=2, whereas the thermal feedback model shows higher late time in-situ star formation rates than observed.

We have investigated particle acceleration and emission from shocks and shear flows associated with an unmagnetized relativistic jet plasma propagating into an unmagnetized ambient plasma. Strong electro-magnetic fields are generated in the jet shock via the filamentation (Weibel) instability. Shock field strength and structure depend on plasma composition (($e^{\pm}$ or $e^-$- $p^+$ plasmas) and Lorentz factor. In the velocity shear between jet and ambient plasmas, strong AC ($e^{\pm}$ plasmas) or DC ($e^-$- $p^+$ plasmas) magnetic fields are generated via the kinetic Kelvin-Helmholtz instability (kKHI), and the magnetic field structure also depends on the jet Lorentz factor. We have calculated, self-consistently, the radiation from electrons accelerated in shock generated magnetic fields. The spectra depend on the jet's initial Lorentz factor and temperature via the resulting particle acceleration and magnetic field generation. Our ongoing "Global" jet simulations containing shocks and velocity shears will provide us with the ability to calculate and model the complex time evolution and/or spectral structure observed from gamma-ray bursts, AGN jets, and supernova remnants.

We report the results of a 3D particle-in-cell (PIC) simulation carried out to study the early-stage evolution of the shock formed when an unmagnetized relativistic jet interacts with an ambient electron-ion plasma. Full-shock structures associated with the interaction are observed in the ambient frame. When open boundaries are employed in the direction of the jet; the forward shock is seen as a hybrid structure consisting of an electrostatic shock combined with a double layer, while the reverse shock is seen as a double layer. The ambient ions show two distinct features across the forward shock: a population penetrating into the shocked region from the precursor region and an accelerated population escaping from the shocked region into the precursor region. This behavior is a signature of a combination of an electrostatic shock and a double layer. Jet electrons are seen to be electrostatically trapped between the forward and reverse shock structures showing a ring-like distribution in a phase-space plot, while ambient electrons are thermalized and become essentially isotropic in the shocked region. The magnetic energy density grows to a few percent of the jet kinetic energy density at both the forward and the reverse shock transition layers in a rather short time scale. We see little disturbance of the jet ions over this time scale.

We study the effect of AGN mechanical and radiation feedback on the formation of bulge dominated galaxies via mergers of disc galaxies. The merging galaxies have mass-ratios of 1:1 to 6:1 and include pre-existing hot gaseous halos to properly account for the global impact of AGN feedback. Using smoothed particle hydrodynamics simulation code (GADGET-3) we compare three models with different AGN feedback models: (1) no black hole and no AGN feedback; (2) thermal AGN feedback; and (3) mechanical and radiative AGN feedback. The last model is motivated by observations of broad line quasars which show winds with initial velocities of $v_w \ge$ 10,000 km/s and also heating associated with the central AGN X-ray radiation. The primary changes in gas properties due to mechanical AGN feedback are lower thermal X-ray luminosity from the final galaxy - in better agreement with observations - and galactic outflows with higher velocity $\sim 1000$ km/s similar to recent direct observations of nearby merger remnants. The kinetic energy of the outflowing gas is a factor of $\sim$ 20 higher than in the thermal feedback case. All merger remnants with momentum-based AGN feedback with $v_w \sim 10,000$ km/s and $\epsilon_w=2 \times 10^{-3}$, independent of their progenitor mass-ratios, reproduce the observed relations between stellar velocity dispersion and black hole mass ($M_{\rm bh} - \sigma$) as well as X-ray luminosity ($L_X - \sigma$) with $10^{37.5}$ erg/s $\lesssim L_X (0.3-8~{\rm keV}) \lesssim 10^{39.2}$ erg/s for velocity dispersions in the range of 120 km/s $\lesssim \sigma \lesssim$ 190 km/s. In addition, the mechanical feedback produces a much greater AGN variability. We also show that gas is more rapidly and impulsively stripped from the galactic centres driving a moderate increase in galaxy size and decrease in central density with the mechanical AGN feedback model.

We have investigated generation of magnetic fields associated with velocity shear between an unmagnetized relativistic jet and an unmagnetized sheath plasma. We have examined the strong magnetic fields generated by kinetic shear (Kelvin-Helmholtz) instabilities. Compared to the previous studies using counter-streaming performed by Alves et al. (2012), the structure of KKHI of our jet-sheath configuration is slightly different even for the global evolution of the strong transverse magnetic field. In our simulations the major components of growing modes are the electric field $E_{\rm z}$ and the magnetic field $B_{\rm y}$. After the $B_{\rm y}$ component is excited, an induced electric field $E_{\rm x}$ becomes significant. However, other field components remain small. We find that the structure and growth rate of KKHI with mass ratios $m_{\rm i}/m_{\rm e} = 1836$ and $m_{\rm i}/m_{\rm e} = 20$ are similar. In our simulations saturation in the nonlinear stage is not as clear as in counter-streaming cases. The growth rate for a mildly-relativistic jet case ($\gamma_{\rm j} = 1.5$) is larger than for a relativistic jet case ($\gamma_{\rm j} = 15$).

We have investigated particle acceleration and shock structure associated with an unmagnetized relativistic jet propagating into an unmagnetized plasma. Strong magnetic fields generated in the trailing jet shock lead to transverse deflection and acceleration of the electrons. We have self-consistently calculated the radiation from the electrons accelerated in the turbulent magnetic fields. We find that the synthetic spectra depend on the bulk Lorentz factor of the jet, the jet temperature, and the strength of the magnetic fields generated in the shock. We have also begun study of electron acceleration in the strong magnetic fields generated by kinetic shear (Kelvin-Helmholtz) instabilities. Our calculated spectra should lead to a better understanding of the complex time evolution and/or spectral structure from gamma-ray bursts, relativistic jets, and supernova remnants.

We study the growth of black holes (BHs) in galaxies using three-dimensional smoothed particle hydrodynamic simulations with new implementations of the momentum mechanical feedback, and restriction of accreted elements to those that are gravitationally bound to the BH. We also include the feedback from the X-ray radiation emitted by the BH, which heats the surrounding gas in the host galaxies, and adds radial momentum to the fluid. We perform simulations of isolated galaxies and merging galaxies and test various feedback models with the new treatment of the Bondi radius criterion. We find that overall the BH growth is similar to what has been obtained by earlier workers using the Springel, Di Matteo, & Hernquist algorithms. However, the outflowing wind velocities and mechanical energy emitted by winds are considerably higher (v_w ~ 1000-3000 km/s) compared to the standard thermal feedback model (v_w ~ 50-100 km/s). While the thermal feedback model emits only 0.1 % of BH released energy in winds, the momentum feedback model emits more than 30 % of the total energy released by the BH in winds. In the momentum feedback model, the degree of fluctuation in both radiant and wind output is considerably larger than in the standard treatments. We check that the new model of the BH mass accretion agrees with analytic results for the standard Bondi problem.

Thermal instability (TI) can strongly affect the structure and dynamics of the interstellar medium (ISM) in the Milky Way and other disk galaxies. Thermal conduction plays an important role in the TI by stabilizing small scales and limiting the size of the smallest condensates. In the magnetized ISM, however, heat is conducted anisotropically (primarily along magnetic field lines). We investigate the effects of anisotropic thermal conduction on the nonlinear regime of the TI by performing two-dimensional magnetohydrodynamic simulations. We present models with magnetic fields of different initial geometries and strengths, and compare them to hydrodynamic models with isotropic conduction. We find anisotropic conduction does not significantly alter the overall density and temperature statistics in the saturated state of the TI. However, it can strongly affect the shapes and sizes of cold clouds formed by the TI. For example, for uniform initial fields long filaments of cold gas are produced that are reminiscent of some observed HI clouds. For initially tangled fields, such filaments are not produced. We also show that anisotropic conduction suppresses turbulence generated by evaporative flows from the surfaces of cold blobs, which may have implications for mechanisms for driving turbulence in the ISM.

lasma instabilities excited in collisionless shocks are responsible for particle acceleration. We have investigated the particle acceleration and shock structure associated with an unmagnetized relativistic electron-positron jet propagating into an unmagnetized electron-positron plasma. Cold jet electrons are thermalized and slowed while the ambient electrons are swept up to create a partially developed hydrodynamic-like shock structure. In the leading shock, electron density increases by a factor of about 3.5 in the simulation frame. Strong electromagnetic fields are generated in the trailing shock and provide an emission site. These magnetic fields contribute to the electron's transverse deflection behind the shock. Our initial results of a jet-ambient interaction with anti-parallel magnetic fields show pile-up of magnetic fields at the colliding shock, which may lead to reconnection and associated particle acceleration. We will investigate the radiation in transient stage as a possible generation mechanism of precursors of prompt emission. In our simulations we calculate the radiation from electrons in the shock region. The detailed properties of this radiation are important for understanding the complex time evolution and spectral structure in gamma-ray bursts, relativistic jets, and supernova remnants.

The ideal gas equation of state with a constant adiabatic index, although commonly used in relativistic hydrodynamics, is a poor approximation for most relativistic astrophysical flows. Here we propose a new general equation of state for a multi-component relativistic gas which is consistent with the Synge equation of state for a relativistic perfect gas and is suitable for numerical (special) relativistic hydrodynamics. We also present a multidimensional relativistic hydrodynamics code incorporating the proposed general equation of state, based on the HLL scheme, which does not make use of a full characteristic decomposition of the relativistic hydrodynamic equations. The accuracy and robustness of this code is demonstrated in multidimensional calculations through several highly relativistic test problems taking into account nonvanishing tangential velocities. Results from three-dimensional simulations of relativistic jets show that the morphology and dynamics of the relativistic jets are significantly influenced by the different equation of state and by different compositions of relativistic perfect gases. Our new numerical code, combined with our proposed equation of state is very efficient and robust, and unlike previous codes, it gives very accurate results for thermodynamic variables in relativistic astrophysical flows.

The deposition of mechanical feedback from a supermassive black hole (SMBH) in an active galactic nucleus (AGN) into the surrounding galaxy occurs via broad-line winds which must carry mass and radial momentum as well as energy. The effect can be summarized by the dimensionless parameter $\eta=dot{M_outflow}/dot{M_accretion}= (2 \epsilon_w c^2)/v_w^2$ where ($\epslion_w \equiv dot{E}_w/(dot{M_accretion} c^2)$) is the efficiency by which accreted matter is turned into wind energy in the disc surrounding the central SMBH. The outflowing mass and omentum are proportional to $\eta$, and many prior treatments have essentially assumed that $\eta=0$. We perform one- and two-dimensional simulations and find that the growth of the central SMBH is very sensitive to the inclusion of the mass and momentum driving but is insensitive to the assumed mechanical efficiency. For example in representative calculations, the omission of momentum and mass feedback leads to an hundred fold increase in the mass of the SMBH to over $10^{10} \Msun$. When allowance is made for momentum driving, the final SMBH mass is much lower and the wind efficiencies which lead to the most observationally acceptable results are relatively low with $\epsilon_w \lesssim 10^{-4}$.

We study filamentary structure in the galaxy distribution at z ~ 0.8 using data from the Deep Extragalactic Evolutionary Probe 2 (DEEP2) Redshift Survey and its evolution to z ~ 0.1 using data from the Sloan Digital Sky Survey (SDSS). We trace individual filaments for both surveys using the Smoothed Hessian Major Axis Filament Finder, an algorithm which employs the Hessian matrix of the galaxy density field to trace the filamentary structures in the distribution of galaxies. We extract 33 subsamples from the SDSS data with a geometry similar to that of DEEP2. We find that the filament length distribution has not significantly changed since z ~ 0.8, as predicted in a previous study using a $\Lamda$CDM cosmological N-body simulation. However, the filament width distribution, which is sensitive to the non-linear growth of structure, broadens and shifts to smaller widths for smoothing length scales of 5-10 Mpc/h from z ~ 0.8 to z ~ 0.1, in accord with N-body simulations.

We describe a method for incorporating ambipolar diffusion in the strong coupling approximation into a multidimensional magnetohydrodynamics code based on the total variation diminishing scheme. Contributions from ambipolar diffusion terms are included by explicit finite difference operators in a fully unsplit way, maintaining second order accuracy. The divergence-free condition of magnetic fields is exactly ensured at all times by a flux-interpolated constrained transport scheme. The super time stepping method is used to accelerate the timestep in high resolution calculations and/or in strong ambipolar diffusion. We perform two test problems, the steady-state oblique C-type shocks and the decay of Alfv\'en waves, confirming the accuracy and robustness of our numerical approach. Results from the simulations of the compressible MHD turbulence with ambipolar diffusion show the flexibility of our method as well as its ability to follow complex MHD flows in the presence of ambipolar diffusion. These simulations show that the dissipation rate of MHD turbulence is strongly affected by the strength of ambipolar diffusion.

The most massive Galactic globular cluster, Omega centauri, appears to have multiple populations. Its bluest main sequence and extended horizontal branch stars are suggested to have the common origin, that is, an extremely high helium abundance of Y ~ 0.4. The high helium abundance is most often attributed to asymptotic giant branch (AGB) stars. In this study we test the AGB hypothesis. We simulate the maximum-AGB models where the impact of AGB stars is maximised by assuming that supernova explosions do not affect the chemical evolution of the proto cloud. We compare the enrichment history of helium, metals, carbon and nitrogen to the observed values. Even under the most generous condition, the maximum-AGB models fail to reproduce the large values of helium Y ~ 0.4 and helium enrichment parameter Delta Y / Delta Z ~ 70 which were deduced from the colour-magnitude diagram fits. They also fail to reproduce the C and N contents of the blue population spectroscopically determined. We conclude that the AGB scenario with the canonical stellar evolution theory cannot explain the observational constraints and that the self chemical enrichment does not provide a viable solution. Alternative processes are desperately called for.

The most massive Galactic globular cluster omega Cen appears to have two, or perhaps more, distinct main sequences. Its bluest main sequence is at the centre of debate because it has been suggested to have an extremely high helium abundance of Y ~ 0.4. The same helium abundance is claimed to explain the presence of extreme horizontal branch stars of omega Cen as well. This demands a relative helium to metal enrichment of deltaY/deltaZ ~ 70; that is, more than one order of magnitude larger than the generally accepted value. Candidate solutions, namely, AGB stars, massive stars, and supernovae, have been suggested; but in this study, we show that none of them is a viable channel, in terms of reproducing the high value of deltaY/deltaZ for the constrained age difference between the red and blue populations. Essentially no populations with an ordinary initial mass function, including those candidates, can produce such a high deltaY/deltaZ because they all produce metals as well as helium. As an alternative, we investigate the possibility of the stochastic ``first star'' contamination to the gas from which the younger generation of omega Cen formed. This requires the assumption that Population III star formation episode overlaps with that of Population II. While the required condition appears extreme, very massive objects in the first star generation provide a solution that is at least as plausible as any other suggestions made before.

We have studied three-dimensional hydrodynamic interactions of relativistic extragalactic jets with two-phase ambient media. These jets propagate through a denser homogeneous gas and then impact clouds with densities 100 to 1000 times higher than the initial beam density. The deflection angle of the jet is influenced more by the density contrast of the cloud than by the beam Mach number of the jet. A relativistic jet with low relativistic beam Mach number can eventually be slightly bent after it crosses the dense cloud; however, we have not seen permanently bent structures in the interaction of a high relativistic beam Mach number jet with a cloud. The relativistic jet impacts on dense clouds do not necessarily destroy the clouds completely, and much of the cloud body can survive as a coherent blob. This enhancement of cloud durability is partly due to the geometric influence of the off-axis collisions we consider and also arises from the lower rate of cloud fragmentation through the Kelvin-Helmholtz instability for relativistic jets. To compare our simulations with observed extragalactic radio jets, we have computed the approximate surface distributions of synchrotron emission at different viewing angles. These surface intensity maps show relativistic jets interacting with clouds can produce synchrotron emission knots similar to structures observed in many VLBI-scale radio sources. We find that the synchrotron emission increases steeply at the moment of impact and the emission peaks right before the jet passes through the cloud.

Relativistic temperature of gas raises the issue of the equation of state (EoS) in relativistic hydrodynamics. We study the EoS for numerical relativistic hydrodynamics, and propose a new EoS that is simple and yet approximates very closely the EoS of the single-component perfect gas in relativistic regime. We also discuss the calculation of primitive variables from conservative ones for the EoS's considered in the paper, and present the eigenstructure of relativistic hydrodynamics for a general EoS, in a way that they can be used to build numerical codes. Tests with a code based on the Total Variation Diminishing (TVD) scheme are presented to highlight the differences induced by different EoS's.

This paper describes a multidimensional hydrodynamic code which can be used for studies of relativistic astrophysical flows. The code solves the special relativistic hydrodynamic equations as a hyperbolic system of conservation laws based on the total variation diminishing (TVD) scheme. It uses a new set of conserved quantities and employs an analytic formula for transformation from the conserved quantities in the reference frame to the physical quantities in the local rest frame. Several standard tests, including relativistic shock tubes, a relativistic wall shock, and a relativistic blast wave, are presented to demonstrate that the code captures discontinuities correctly and sharply in ultrarelativistic regimes. The robustness and flexibility of the code are demonstrated through test simulations of the relativistic Hawley-Zabusky shock and a relativistic extragalactic jet.

We present an investigation of H-alpha emission line variations observed in the massive Algol binary, RY Per. We give new radial velocity data for the secondary based upon our optical spectra and for the primary based upon high dispersion UV spectra. We present revised orbital elements and an estimate of the primary's projected rotational velocity (which indicates that the primary is rotating 7 times faster than synchronous). We use a Doppler tomography algorithm to reconstruct the individual primary and secondary spectra in the region of H-alpha, and we subtract the latter from each of our observations to obtain profiles of the primary and its disk alone. Our H-alpha observations of RY Per show that the mass gaining primary is surrounded by a persistent but time variable accretion disk. The profile that is observed outside-of-eclipse has weak, double-peaked emission flanking a deep central absorption, and we find that these properties can be reproduced by a disk model that includes the absorption of photospheric light by the band of the disk seen in projection against the face of the star. We developed a new method to reconstruct the disk surface density distribution from the ensemble of H-alpha profiles observed around the orbit, and this method accounts for the effects of disk occultation by the stellar components, the obscuration of the primary by the disk, and flux contributions from optically thick disk elements. The resulting surface density distribution is elongated along the axis joining the stars, in the same way as seen in hydrodynamical simulations of gas flows that strike the mass gainer near trailing edge of the star. This type of gas stream configuration is optimal for the transfer of angular momentum, and we show that rapid rotation is found in other Algols that have passed through a similar stage.

We report a study to constrain the fraction of baryonic matter in the cold plus hot dark matter (CHDM) universe by numerical simulations which include the hydrodynamics of baryonic matter as well as the particle dynamics of dark matter. Spatially flat, COBE-normalized CHDM models with the fraction of hot component $\Omega_h\leq0.2$ are considered. We show that the models with $h/n/\Omega_h=0.5/0.9/0.1$ and $0.5/0.9/0.2$ give a linear power spectrum which agrees well with observations. Here, $h$ is the Hubble constant in unit of $100~km/s/Mpc$ and $n$ is the spectral index of the initial power spectrum. Then, for the models with $h/n/\Omega_h=0.5/0.9/0.2$ and baryonic fraction $\Omega_b=0.05$ and 0.1 we calculate the properties of X-ray clusters, such as luminosity function, temperature distribution function, luminosity-temperature relation, histogram of gas to total mass ratio, and change of average temperature with redshift $z$. Comparison with the observed data of X-ray clusters indicates that the model with $\Omega_b=0.05$ is preferred. The COBE-normalized CHDM model with $\Omega_b>0.1$ may be ruled out by the present work, since it produces too many X-ray bright clusters.