We investigate the large-scale influence of outflows from AGNs in enriching the IGM with metals in a cosmological context. We combine cosmological simulations of large scale structure formation with a detailed model of metal enrichment, in which outflows expand anisotropically along the direction of least resistance, distributing metals into the IGM. The metals carried by the outflows are generated by two separate stellar populations: stars located near the central AGN, and stars located in the greater galaxy. Using this algorithm, we performed a series of 5 simulations of the propagation of AGN-driven outflows in a cosmological volume of size (128/h Mpc)^3 in a Lambda-CDM universe, and analyze the resulting metal enrichment of the IGM. We found that the metallicity induced in the IGM is greatly dominated by AGNs having bolometric luminosity L > 10^9 L_sun, sources with 10^8 < L / L_sun < 10^9 having a negligible contribution. Our simulations produced an average IGM metallicity of [O/H] = -5 at z = 5.5, which then rises gradually, and remains relatively flat at a value [O/H] = -2.8 between z = 2 and z = 0. The ejection of metals from AGN host galaxies by AGN-driven outflows is found to enrich the IGM to > 10 - 20% of the observed values, the number dependent on redshift. The enriched IGM volume fractions are small at z > 3, then rise rapidly to the following values at z = 0: 6 - 10% of the volume enriched to [O/H] > -2.5, 14 - 24% volume to [O/H] > -3, and 34 - 45% volume to [O/H] > -4. At z > 2, there is a gradient of the induced enrichment, the metallicity decreasing with increasing IGM density, enriching the underdense IGM to higher metallicities, a trend more prominent with increasing anisotropy of the outflows. This can explain observations of metal-enriched low-density IGM at z = 3 - 4.

We investigate the cosmological-scale influence of outflows driven by AGNs on metal enrichment of the intergalactic medium. AGNs are located in dense cosmological structures which tend to be anisotropic. We designed a semi-analytical model for anisotropic AGN outflows which expand away along the direction of least resistance. This model was implemented into a cosmological numerical simulation algorithm for simulating the growth of large-scale structure in the universe. Using this modified algorithm, we perform a series of 9 simulations inside cosmological volumes of size $(128 h^{-1}{\rm Mpc})^3$, in a concordance $\Lambda$CDM universe, varying the opening angle of the outflows, the lifetimes of the AGNs, their kinetic fractions, and their level of clustering. For each simulation, we compute the volume fraction of the IGM enriched in metals by the outflows. The resulting enriched volume fractions are relatively small at $z \gtrsim 2.5$, and then grow rapidly afterward up to $z = 0$. We find that AGN outflows enrich from 65% to 100% of the entire universe at the present epoch, for different values of the model parameters. The enriched volume fraction depends weakly on the opening angle of the outflows. However, increasingly anisotropic outflows preferentially enrich underdense regions, a trend found more prominent at higher redshifts and decreasing at lower redshifts. The enriched volume fraction increases with increasing kinetic fraction and decreasing AGN lifetime and level of clustering.

We simulate anisotropic outflows of AGN, and investigate the large-scale impact of the cosmological population of AGN outflows over the Hubble time by performing N-body LambdaCDM simulations. Using the observed quasar luminosity function to get the redshift and luminosity distribution, and analytical models for the outflow expansion, AGNs are allowed to evolve in a cosmological volume. By the present epoch, 13 - 25% of the total volume is found to be pervaded by AGN outflows, with 10^{-9} G magnetic field.