We have developed a symbolic algebra approach to automatically produce, verify, and optimize computer code for the Fast Multipole Method (FMM) operators. This approach allows for flexibility in choosing a basis set and kernel, and can generate computer code for any expansion order in multiple languages. The procedure is implemented in the publicly available Python program Mosaic. Optimizations performed at the symbolic level through algebraic manipulations significantly reduce the number of mathematical operations compared with a straightforward implementation of the equations. We find that the optimizer is able to eliminate 20-80% of the floating-point operations and for the expansion orders $p \le 10$ it changes the observed scaling properties. We present our approach using three variants of the operators with the Cartesian basis set for the harmonic potential kernel $1/r$, including the use of totally symmetric and traceless multipole tensors.

We explore to what extent simple algebraic models can be used to describe H II regions when winds, radiation pressure, gravity and photon breakout are included. We a) develop algebraic models to describe the expansion of photoionised H II regions under the influence of gravity and accretion in power-law density fields with $\rho \propto r^{-w}$, b) determine when terms describing winds, radiation pressure, gravity and photon breakout become significant enough to affect the dynamics of the H II region where $w=2$, and c) solve these expressions for a set of physically-motivated conditions. We find that photoionisation feedback from massive stars is the principal mode of feedback on molecular cloud scales, driving accelerating outflows from molecular clouds in cases where the peaked density structure around young massive stars is considered at radii between $\sim$0.1 and 10-100 pc. Under a large range of conditions the effect of winds and radiation on the dynamics of H II regions is around 10% of the contribution from photoionisation. The effect of winds and radiation pressure are most important at high densities, either close to the star or in very dense clouds such as those in the Central Molecular Zone of the Milky Way. Out to $\sim$0.1 pc they are the principal drivers of the H II region. Lower metallicities make the relative effect of photoionisation even stronger as the ionised gas temperature is higher.

The circumgalactic medium (CGM) is one of the frontiers of galaxy formation and intimately connected to the galaxy via accretion of gas on to the galaxy and gaseous outflows from the galaxy. Here we analyse the magnetic field in the CGM of the Milky Way-like galaxies simulated as part of the \textsc{Auriga} project that constitutes a set of high resolution cosmological magnetohydrodynamical zoom simulations. We show that at high redshift the CGM becomes magnetised via galactic outflows that transport magnetised gas from the disk into the halo. At this time the magnetisation of the CGM closely follows its metal enrichment. We then show that at low redshift an in-situ turbulent dynamo that operates on a timescale of Gigayears further amplifies the magnetic field in the CGM and sets its strength down $z=0$. The magnetic field strength reaches a typical value of $0.1\,\mu G$ at the virial radius at $z=0$ and becomes mostly uniform within the virial radius. Its Faraday rotation signal is in excellent agreement with recent observations. For most of its evolution the magnetic field in the CGM is an unordered small scale field. Only strong coherent outflows at low redshift are able to order the magnetic field in parts of the CGM that are directly displaced by these outflows.

Molecular clouds are turbulent structures whose star formation efficiency (SFE) is strongly affected by internal stellar feedback processes. In this paper we determine how sensitive the SFE of molecular clouds is to randomised inputs in the star formation feedback loop, and to what extent relationships between emergent cloud properties and the SFE can be recovered. We introduce the yule suite of 26 radiative magnetohydrodynamic (RMHD) simulations of a 10,000 solar mass cloud similar to those in the solar neighbourhood. We use the same initial global properties in every simulation but vary the initial mass function (IMF) sampling and initial cloud velocity structure. The final SFE lies between 6 and 23 percent when either of these parameters are changed. We use Bayesian mixed-effects models to uncover trends in the SFE. The number of photons emitted early in the cluster's life and the length of the cloud provide are the strongest predictors of the SFE. The HII regions evolve following an analytic model of expansion into a roughly isothermal density field. The more efficient feedback is at evaporating the cloud, the less the star cluster is dispersed. We argue that this is because if the gas is evaporated slowly, the stars are dragged outwards towards surviving gas clumps due to the gravitational attraction between the stars and gas. While star formation and feedback efficiencies are dependent on nonlinear processes, statistical models describing cloud-scale processes can be constructed.

Feedback by Active Galactic Nuclei is often divided into quasar and radio mode, powered by radiation or radio jets, respectively. Both are fundamental in galaxy evolution, especially in late-type galaxies, as shown by cosmological simulations and observations of jet-ISM interactions in these systems. We compare AGN feedback by radiation and by collimated jets through a suite of simulations, in which a central AGN interacts with a clumpy, fractal galactic disc. We test AGN of $10^{43}$ and $10^{46}$ erg/s, considering jets perpendicular or parallel to the disc. Mechanical jets drive the more powerful outflows, exhibiting stronger mass and momentum coupling with the dense gas, while radiation heats and rarifies the gas more. Radiation and perpendicular jets evolve to be quite similar in outflow properties and effect on the cold ISM, while inclined jets interact more efficiently with all the disc gas, removing the densest $20\%$ in $20$ Myr, and thereby reducing the amount of cold gas available for star formation. All simulations show small-scale inflows of $0.01-0.1$ M$_\odot$/yr, which can easily reach down to the Bondi radius of the central supermassive black hole (especially for radiation and perpendicular jets), implying that AGN modulate their own duty cycle in a feedback/feeding cycle.

Feedback from Active Galactic Nuclei (AGN) and subsequent jet cocoons and outflow bubbles can have a significant impact on star formation in the host galaxy. To investigate feedback physics on small scales, we perform hydrodynamic simulations of realistically fast AGN winds striking Bonnor-Ebert (BE) spheres and examine gravitational collapse and ablation. We test AGN wind velocities ranging from 300--3,000 km s$^{-1}$ and wind densities ranging from 0.5--10 $m_\mathrm{p}\,\mathrm{cm}^{-3}$. We include heating and cooling of low- and high-temperature gas, self-gravity, and spatially correlated perturbations in the shock, with a maximum resolution of 0.01 pc. We find that the ram pressure is the most important factor that determines the fate of the cloud. High ram pressure winds increase fragmentation and decrease the star formation rate, but also cause star formation to occur on a much shorter time scale and with increased velocities of the newly formed stars. We find a threshold ram pressure of $\sim 2\times10^{-8}$ dyne cm$^{-2}$ above which stars are not formed because the resulting clumps have internal velocities large enough to prevent collapse. Our results indicate that simultaneous positive and negative feedback will be possible in a single galaxy as AGN wind parameters will vary with location within a galaxy.

The quasar mode of Active Galactic Nuclei (AGN) in the high-redshift Universe is routinely observed in gas-rich galaxies together with large-scale AGN-driven winds. It is crucial to understand how photons emitted by the central AGN source couple to the ambient interstellar-medium to trigger large-scale outflows. By means of radiation-hydrodynamical simulations of idealised galactic discs, we study the coupling of photons with the multiphase galactic gas, and how it varies with gas cloud sizes, and the radiation bands included in the simulations, which are ultraviolet (UV), optical, and infrared (IR). We show how a quasar with a luminosity of $10^{46}$ erg/s can drive large-scale winds with velocities of $10^2-10^3$ km/s and mass outflow rates around $10^3$ M$_\odot$/yr for times of order a few million years. Infrared radiation is necessary to efficiently transfer momentum to the gas via multi-scattering on dust in dense clouds. However, IR multi-scattering, despite being extremely important at early times, quickly declines as the central gas cloud expands and breaks up, allowing the radiation to escape through low gas density channels. The typical number of multi-scattering events for an IR photon is only about a quarter of the mean optical depth from the center of the cloud. Our models account for the observed outflow rates of $\sim$500-1000 M$_\odot$/yr and high velocities of $\sim 10^3$ km/s, favouring winds that are energy-driven via extremely fast nuclear outflows, interpreted here as being IR-radiatively-driven winds.

Feedback from active galactic nuclei (AGN) has often been invoked both in simulations and in interpreting observations for regulating star formation and quenching cooling flows in massive galaxies. AGN activity can, however, also over-pressurise the dense star-forming regions of galaxies and thus enhance star formation, leading to a positive feedback effect. To understand this pressurisation better, we investigate the effect of an ambient external pressure on gas fragmentation and triggering of starburst activity by means of hydrodynamical simulations. We find that moderate levels of over-pressurisation of the galaxy boost the global star formation rate of the galaxy by an order of magnitude, turn stable discs unstable, and lead to significant fragmentation of the gas content of the galaxy, similar to what is observed in high redshift galaxies.

Feedback in massive galaxies generally involves quenching of star formation, a favored candidate being outflows from a central supermassive black hole. At high redshifts however, explanation of the huge rates of star formation often found in galaxies containing AGN may require a more vigorous mode of star formation than attainable by simply enriching the gas content of galaxies in the usual gravitationally-driven mode that is associated with the nearby Universe. Using hydrodynamical simulations, we demonstrate that AGN-pressure-driven star formation potentially provides the positive feedback that may be required to generate the accelerated star formation rates observed in the distant Universe.