Researcher profile

Philip Mocz

Philip Mocz contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
15works
0followers
10topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

15 published item(s)

preprint2023arXiv

On the Cosmic Web Elongation in Fuzzy Dark Matter Cosmologies: Effects on Density Profiles, Shapes and Alignments of Halos

The fuzzy dark matter (FDM) scenario has received increased attention in recent years due to the small-scale challenges of the vanilla Lambda cold dark matter ($Λ$CDM) cosmological model and the lack of any experimental evidence for any candidate particle. In this study, we use cosmological $N$-body simulations to investigate high-redshift dark matter halos and their responsiveness to an FDM-like power spectrum cutoff on small scales in the primordial density perturbations. We study halo density profiles, shapes and alignments in FDM-like cosmologies (the latter two for the first time) by providing fits and quantifying departures from $Λ$CDM as a function of the particle mass $m$. Compared to $Λ$CDM, the concentrations of FDM-like halos are lower, peaking at an $m$-dependent halo mass and thus breaking the approximate universality of density profiles in $Λ$CDM. The intermediate-to-major and minor-to-major shape parameter profiles are monotonically increasing with ellipsoidal radius in $N$-body simulations of $Λ$CDM. In FDM-like cosmologies, the monotonicity is broken, halos are more elongated around the virial radius than their $Λ$CDM counterparts and less elongated closer to the center. Finally, intrinsic alignment correlations, stemming from the deformation of initially spherically collapsing halos in an ambient gravitational tidal field, become stronger with decreasing $m$. At $z\sim 4$, we find a $6.4 σ$-significance in the fractional differences between the isotropised linear alignment magnitudes $D_{\text{iso}}$ in the $m=10^{-22}$ eV model and $Λ$CDM. Such FDM-like imprints on the internal properties of virialised halos are expected to be strikingly visible in the high-$z$ Universe.

preprint2022arXiv

Energy balance and Alfvén Mach numbers in compressible magnetohydrodynamic turbulence with a large-scale magnetic field

Energy equipartition is a powerful theoretical tool for understanding astrophysical plasmas. It is invoked, for example, to measure magnetic fields in the interstellar medium (ISM), as evidence for small-scale turbulent dynamo action, and, in general, to estimate the energy budget of star-forming molecular clouds. In this study we motivate and explore the role of the volume-averaged root-mean-squared (rms) magnetic coupling term between the turbulent, $δ\mathbf{B}$ and large-scale, $\mathbf{B}_0$ fields, $\left< (δ\mathbf{B}\cdot\mathbf{B}_0)^{2} \right>^{1/2}_{\mathcal{V}}$. By considering the second moments of the energy balance equations we show that the rms coupling term is in energy equipartition with the volume-averaged turbulent kinetic energy for turbulence with a sub-Alfvénic large-scale field. Under the assumption of exact energy equipartition between these terms, we derive relations for the magnetic and coupling term fluctuations, which provide excellent, parameter-free agreement with time-averaged data from 280 numerical simulations of compressible MHD turbulence. Furthermore, we explore the relation between the turbulent, mean-field and total Alfvén Mach numbers, and demonstrate that sub-Alfvénic turbulence can only be developed through a strong, large-scale magnetic field, which supports an extremely super-Alfvénic turbulent magnetic field. This means that the magnetic field fluctuations are significantly subdominant to the velocity fluctuations in the sub-Alfvénic large-scale field regime. Throughout our study, we broadly discuss the implications for observations of magnetic fields and understanding the dynamics in the magnetised ISM.

preprint2022arXiv

Snowmass2021 Cosmic Frontier White Paper: Cosmological Simulations for Dark Matter Physics

Over the past several decades, unexpected astronomical discoveries have been fueling a new wave of particle model building and are inspiring the next generation of ever-more-sophisticated simulations to reveal the nature of Dark Matter (DM). This coincides with the advent of new observing facilities coming online, including JWST, the Rubin Observatory, the Nancy Grace Roman Space Telescope, and CMB-S4. The time is now to build a novel simulation program to interpret observations so that we can identify novel signatures of DM microphysics across a large dynamic range of length scales and cosmic time. This white paper identifies the key elements that are needed for such a simulation program. We identify areas of growth on both the particle theory side as well as the simulation algorithm and implementation side, so that we can robustly simulate the cosmic evolution of DM for well-motivated models. We recommend that simulations include a fully calibrated and well-tested treatment of baryonic physics, and that outputs should connect with observations in the space of observables. We identify the tools and methods currently available to make predictions and the path forward for building more of these tools. A strong cosmic DM simulation program is key to translating cosmological observations to robust constraints on DM fundamental physics, and provides a connection to lab-based probes of DM physics.

preprint2020arXiv

Fuzzy Dark Matter Soliton Cores around Supermassive Black Holes

We explore the effect of a supermassive black hole (SMBH) on the density profile of a fuzzy dark matter (FDM) soliton core at the centre of a dark matter halo. We numerically solve the Schrödinger-Poisson equations, treating the black hole as a gravitational point mass, and demonstrate that this additional perturbing term has a `squeezing&#39; effect on the soliton density profile, decreasing the core radius and increasing the central density. In the limit of large black hole mass, the solution approaches one akin to the hydrogen atom, with radius inversely proportional to the black hole mass. By applying our analysis to two specific galaxies (M87 and the Milky Way) and pairing it with known observational limits on the amount of centrally concentrated dark matter, we obtain a constraint on the FDM particle mass, finding that the range $10^{-22.12}$ eV $\lesssim m \lesssim 10^{-22.06}$ eV should be forbidden (taking into account additional factors concerning the life-time of the soliton in the vicinity of a black hole). Improved observational mass measurements of the black hole and total enclosed masses will significantly extend the lower-bound on the excluded FDM mass region, while self-consistent theoretical modeling of the soliton--black hole system can extend the upper-bound.

preprint2019arXiv

Dynamical Friction in a Fuzzy Dark Matter Universe

We present an in-depth exploration of the phenomenon of dynamical friction in a universe where the dark matter is composed entirely of so-called Fuzzy Dark Matter (FDM), ultralight bosons of mass $m\sim\mathcal{O}(10^{-22})\,$eV. We review the classical treatment of dynamical friction before presenting analytic results in the case of FDM for point masses, extended mass distributions, and FDM backgrounds with finite velocity dispersion. We then test these results against a large suite of fully non-linear simulations that allow us to assess the regime of applicability of the analytic results. We apply these results to a variety of astrophysical problems of interest, including infalling satellites in a galactic dark matter background, and determine that \emph{(1)}~for FDM masses $m\gtrsim 10^{-21}\, {\rm eV}\, c^{-2}$, the timing problem of the Fornax dwarf spheroidal&#39;s globular clusters is no longer solved and \emph{(2)}~the effects of FDM on the process of dynamical friction for satellites of total mass $M$ and relative velocity $v_{\rm rel}$ should require detailed numerical simulations for $\left(M/10^9~M_{\odot}\right) \left(m/10^{-22}~{\rm eV}\right)\left(100~{\rm km}~{\rm s}^{-1}/v_{\rm rel}\right) \sim 1$, parameters which would lie outside the validated range of applicability of any currently developed analytic theory, due to transient wave structures in the time-dependent regime.

preprint2019arXiv

Galaxy Formation with BECDM -- II. Cosmic Filaments and First Galaxies

Bose-Einstein Condensate Dark Matter (BECDM; also known as Fuzzy Dark Matter) is motivated by fundamental physics and has recently received significant attention as a serious alternative to the established Cold Dark Matter (CDM) model. We perform cosmological simulations of BECDM gravitationally coupled to baryons and investigate structure formation at high redshifts ($z \gtrsim 5$) for a boson mass $m=2.5\cdot 10^{-22}~{\rm eV}$, exploring the dynamical effects of its wavelike nature on the cosmic web and the formation of first galaxies. Our BECDM simulations are directly compared to CDM as well as to simulations where the dynamical quantum potential is ignored and only the initial suppression of the power spectrum is considered -- a Warm Dark Matter-like (&#34;WDM&#34;) model often used as a proxy for BECDM. Our simulations confirm that &#34;WDM&#34; is a good approximation to BECDM on large cosmological scales even in the presence of the baryonic feedback. Similarities also exist on small scales, with primordial star formation happening both in isolated haloes and continuously along cosmic filaments; the latter effect is not present in CDM. Global star formation and metal enrichment in these first galaxies are delayed in BECDM/&#34;WDM&#34; compared to the CDM case: in BECDM/&#34;WDM&#34; first stars form at $z\sim 13$/$13.5$ while in CDM star formation starts at $z\sim 35$. The signature of BECDM interference, not present in &#34;WDM&#34;, is seen in the evolved dark matter power spectrum: although the small scale structure is initially suppressed, power on kpc scales is added at lower redshifts. Our simulations lay the groundwork for realistic simulations of galaxy formation in BECDM.

preprint2015arXiv

Improving the convergence properties of the moving-mesh code AREPO

Accurate numerical solutions of the equations of hydrodynamics play an ever more important role in many fields of astrophysics. In this work, we reinvestigate the accuracy of the moving-mesh code \textsc{Arepo} and show how its convergence order can be improved for general problems. In particular, we clarify that for certain problems \textsc{Arepo} only reaches first-order convergence for its original formulation. This can be rectified by simple modifications we propose to the time integration scheme and the spatial gradient estimates of the code, both improving the accuracy of the code. We demonstrate that the new implementation is indeed second-order accurate under the $L^1$ norm, and in particular substantially improves conservation of angular momentum. Interestingly, whereas these improvements can significantly change the results of smooth test problems, we also find that cosmological simulations of galaxy formation are unaffected, demonstrating that the numerical errors eliminated by the new formulation do not impact these simulations. In contrast, simulations of binary stars followed over a large number of orbital times are strongly affected, as here it is particularly crucial to avoid a long-term build up of errors in angular momentum conservation.

preprint2015arXiv

Reducing noise in moving-grid codes with strongly-centroidal Lloyd mesh regularization

A method for improving the accuracy of hydrodynamical codes that use a moving Voronoi mesh is described. Our scheme is based on a new regularization scheme that constrains the mesh to be centroidal to high precision while still allowing the cells to move approximately with the local fluid velocity, thereby retaining the quasi-Lagrangian nature of the approach. Our regularization technique significantly reduces mesh noise that is attributed to changes in mesh topology and deviations from mesh regularity. We demonstrate the advantages of our method on various test problems, and note in particular improvements obtained in handling shear instabilities, mixing, and in angular momentum conservation. Calculations of adiabatic jets in which shear excites Kelvin Helmholtz instability show reduction of mesh noise and entropy generation. In contrast, simulations of the collapse and formation of an isolated disc galaxy are nearly unaffected, showing that numerical errors due to the choice of regularization do not impact the outcome in this case.

preprint2015arXiv

The large-scale properties of simulated cosmological magnetic fields

We perform uniformly sampled large-scale cosmological simulations including magnetic fields with the moving mesh code AREPO. We run two sets of MHD simulations: one including adiabatic gas physics only; the other featuring the fiducial feedback model of the Illustris simulation. In the adiabatic case, the magnetic field amplification follows the $B \propto ρ^{2/3}$ scaling derived from `flux-freezing&#39; arguments, with the seed field strength providing an overall normalization factor. At high baryon overdensities the amplification is enhanced by shear flows and turbulence. Feedback physics and the inclusion of radiative cooling change this picture dramatically. In haloes, gas collapses to much larger densities and the magnetic field is amplified strongly and to the same maximum intensity irrespective of the initial seed field of which any memory is lost. At lower densities a dependence on the seed field strength and orientation, which in principle can be used to constrain models of cosmic magnetogenesis, is still present. Inside the most massive haloes magnetic fields reach values of $\sim 10-100\,\,{\rm μG}$, in agreement with galaxy cluster observations. The topology of the field is tangled and gives rise to rotation measure signals in reasonable agreement with the observations. However, the rotation measure signal declines too rapidly towards larger radii as compared to observational data.

preprint2014arXiv

A Constrained Transport Scheme for MHD on Unstructured Static and Moving Meshes

Magnetic fields play an important role in many astrophysical systems and a detailed understanding of their impact on the gas dynamics requires robust numerical simulations. Here we present a new method to evolve the ideal magnetohydrodynamic (MHD) equations on unstructured static and moving meshes that preserves the magnetic field divergence-free constraint to machine precision. The method overcomes the major problems of using a cleaning scheme on the magnetic fields instead, which is non-conservative, not fully Galilean invariant, does not eliminate divergence errors completely, and may produce incorrect jumps across shocks. Our new method is a generalization of the constrained transport (CT) algorithm used to enforce the $\nabla\cdot \mathbf{B}=0$ condition on fixed Cartesian grids. Preserving $\nabla\cdot \mathbf{B}=0$ at the discretized level is necessary to maintain the orthogonality between the Lorentz force and $\mathbf{B}$. The possibility of performing CT on a moving mesh provides several advantages over static mesh methods due to the quasi-Lagrangian nature of the former (i.e., the mesh generating points move with the flow), such as making the simulation automatically adaptive and significantly reducing advection errors. Our method preserves magnetic fields and fluid quantities in pure advection exactly.

preprint2014arXiv

Interpreting MAD within multiple accretion regimes

General relativistic magnetohydrodynamic (GRMHD) simulations of accreting black holes in the radiatively inefficient regime show that systems with sufficient magnetic poloidal flux become magnetically arrested disc (MAD) systems, with a well-defined relationship between the magnetic flux and the mass accretion rate. Recently, Zamaninasab (2014) report that the jet magnetic flux and accretion disc luminosity are tightly correlated over 7 orders of magnitude for a sample of 76 radio-loud active galaxies, concluding that the data are explained by the MAD mode of accretion. Their analysis assumes radiatively efficient accretion, and their sample consists primarily of radiatively efficient sources, while GRMHD simulations of MAD thus far have been carried out in the radiatively inefficient regime. We propose a model to interpret MAD systems in the context of multiple accretion regimes, and apply it to the sample in Zamaninasab (2014), along with additional radiatively inefficient sources from archival data. We show that most of the radiatively inefficient radio-loud galaxies are consistent with being MAD systems. Assuming the MAD relationship found in radiatively inefficient simulations holds at other accretion regimes, a significant fraction of our sample can be candidates for MAD systems. Future GRMHD simulations have yet to verify the validity of this assumption.

preprint2013arXiv

A discontinuous Galerkin method for solving the fluid and MHD equations in astrophysical simulations

A discontinuous Galerkin (DG) method suitable for large-scale astrophysical simulations on Cartesian meshes as well as arbitrary static and moving Voronoi meshes is presented. Most major astrophysical fluid dynamics codes use a finite volume (FV) approach. We demonstrate that the DG technique offers distinct advantages over FV formulations on both static and moving meshes. The DG method is also easily generalized to higher than second-order accuracy without requiring the use of extended stencils to estimate derivatives (thereby making the scheme highly parallelizable). We implement the technique in the AREPO code for solving the fluid and the magnetohydrodynamic (MHD) equations. By examining various test problems, we show that our new formulation provides improved accuracy over FV approaches of the same order, and reduces post-shock oscillations and artificial diffusion of angular momentum. In addition, the DG method makes it possible to represent magnetic fields in a locally divergence-free way, improving the stability of MHD simulations and moderating global divergence errors, and is a viable alternative for solving the MHD equations on meshes where Constrained-Transport (CT) cannot be applied. We find that the DG procedure on a moving mesh is more sensitive to the choice of slope limiter than is its FV method counterpart. Therefore, future work to improve the performance of the DG scheme even further will likely involve the design of optimal slope limiters. As presently constructed, our technique offers the potential of improved accuracy in astrophysical simulations using the moving mesh AREPO code as well as those employing adaptive mesh refinement (AMR).

preprint2012arXiv

The Tully-Fisher Relation for 25,000 SDSS Galaxies as Function of Environment

We construct Tully-Fisher relationships (TFRs) in the $u$, $g$, $r$, $i$ and $z$ bands and stellar mass TFRs (smTFRs) for a sample of $25,698$ late spiral type galaxies (with $0.045<z<0.085$) from the Sloan Digital Sky Survey (SDSS) and study the effects of environment on the relation. We use SDSS-measured Balmer emission line widths, $v_{\rm FWHM}$, as a proxy for disc circular velocity, $v_{\rm circ}$. A priori it is not clear whether we can construct accurate TFRs given the small $3&#34;$ diameter of the fibres used for SDSS spectroscopic measurements. However, we show by modelling the H$α$ emission profile as observed through a $3&#34;$ aperture that for galaxies at appropriate redshifts ($z>0.045$) the fibres sample enough of the disc to obtain a linear relationship between $v_{\rm FWHM}$ and $v_{\rm circ}$, allowing us to obtain a TFR and to investigate dependence on other variables. We also develop a methodology for distinguishing between astrophysical and sample bias in the fibre TFR trends. We observe the well-known steepening of the TFR in redder bands in our sample. We divide the sample of galaxies into four equal groups using projected neighbour density ($Σ$) quartiles and find no significant dependence on environment, extending previous work to a wider range of environments and a much larger sample. Having demonstrated that we can construct SDSS-based TFRs is very useful for future applications because of the large sample size available.

preprint2011arXiv

A Detection of an X-ray Wind and an Ionized Disk in the Chandra HETGS Observation of the Seyfert 2 Galaxy IRAS 18325-5926

We analyze the \textit{Chandra} High Energy Transmission Grating Spectrometer (HETGS) observation of the Seyfert $2$ AGN IRAS\,18325-5926. We detect a $v=-360^{+41}_{-66}$~km~s$^{-1}$ blueshifted ionized absorber in the X-ray spectrum, with photoionization parameter $\logξ=2.0^{+0.1}_{-0.1}$ and hydrogen column density $N_{\rm H}=1.55^{+0.75}_{-0.38}\times10^{21}$~cm$^{-2}$. The absorber may be a photoionized wind originating in the obscuring torus/global covering around the black hole or outer edge of the accretion disk. The estimated mass outflow rate suggests that the supermassive black hole in IRAS\,18325-5926 may significantly affect the large-scale environment of the host galaxy, unless the solid angle subtended by the outflow or the gas filling factor is small. A second warm absorber may be needed to explain the absorption features in the vicinity of the iron K edge, although insufficient counts in the data beyond $7.0$~keV make it difficult to assess the nature of the second absorber. Most plausible is a high ionization ($\logξ\sim 2.3$ to $2.6$), high column density ($N_{\rm H}\sim10^{23}$~cm$^{-2}$) absorber with $v\sim -3000$ to $0$~km~s$^{-1}$, although these parameters are not well constrained. We also examine the broad Fe K emission line in the spectrum, which is likely due to Fe XXV in a highly ionized accretion disk inclined at $25^\circ$, consistent with the \textit{XMM-Newton} EPIC observations of the emission line. Because we are able to view both the obscuring gas and the accretion disk of IRAS\,18325-5926, the surrounding gas of IRAS\,18325-5926 may be patchy or we are viewing the system at an angle just grazing the obscuring torus.

preprint2011arXiv

The inverse-Compton ghost HDF 130 and the giant radio galaxy 6C 0905+3955: matching an analytic model for double radio source evolution

We present new GMRT observations of HDF 130, an inverse-Compton (IC) ghost of a giant radio source that is no longer being powered by jets. We compare the properties of HDF 130 with the new and important constraint of the upper limit of the radio flux density at 240 MHz to an analytic model. We learn what values of physical parameters in the model for the dynamics and evolution of the radio luminosity and X-ray luminosity (due to IC scattering of the cosmic microwave background (CMB)) of a Fanaroff-Riley II (FR II) source are able to describe a source with features (lobe length, axial ratio, X-ray luminosity, photon index and upper limit of radio luminosity) similar to the observations. HDF 130 is found to agree with the interpretation that it is an IC ghost of a powerful double-lobed radio source, and we are observing it at least a few Myr after jet activity (which lasted 5--100 Myr) has ceased. The minimum Lorentz factor of injected particles into the lobes from the hotspot is preferred to be $γ\sim10^3$ for the model to describe the observed quantities well, assuming that the magnetic energy density, electron energy density, and lobe pressure at time of injection into the lobe are linked by constant factors according to a minimum energy argument, so that the minimum Lorentz factor is constrained by the lobe pressure. We also apply the model to match the features of 6C 0905+3955, a classical double FR II galaxy thought to have a low-energy cutoff of $γ\sim10^4$ in the hotspot due to a lack of hotspot inverse-Compton X-ray emission. The models suggest that the low-energy cutoff in the hotspots of 6C 0905+3955 is $γ\gtrsim 10^3$, just slightly above the particles required for X-ray emission.