# Posts Tagged numerical simulation

## Recent Postings from numerical simulation

### Formation of sub-horizon black holes from preheating [Cross-Listing]

We study the production of primordial black holes (PBHs) during the preheating stage that follows a chaotic inflationary phase. The scalar fields present in the process are evolved numerically using a modified version of the HLATTICE code. From the output of the numerical simulation we compute the probability distribution of curvature fluctuations paying particular attention to sub-horizon scales. We find that in some specific models these modes grow to large amplitudes developing highly non-Gaussian probability distributions. We then calculate PBH abundances using the standard Press-Schechter criterion and find that overproduction of PBHs is likely in some regions of the chaotic preheating parameter-space.

### Formation of sub-horizon black holes from preheating [Cross-Listing]

We study the production of primordial black holes (PBHs) during the preheating stage that follows a chaotic inflationary phase. The scalar fields present in the process are evolved numerically using a modified version of the HLATTICE code. From the output of the numerical simulation we compute the probability distribution of curvature fluctuations paying particular attention to sub-horizon scales. We find that in some specific models these modes grow to large amplitudes developing highly non-Gaussian probability distributions. We then calculate PBH abundances using the standard Press-Schechter criterion and find that overproduction of PBHs is likely in some regions of the chaotic preheating parameter-space.

### Formation of sub-horizon black holes from preheating

We study the production of primordial black holes (PBHs) during the preheating stage that follows a chaotic inflationary phase. The scalar fields present in the process are evolved numerically using a modified version of the HLATTICE code. From the output of the numerical simulation we compute the probability distribution of curvature fluctuations paying particular attention to sub-horizon scales. We find that in some specific models these modes grow to large amplitudes developing highly non-Gaussian probability distributions. We then calculate PBH abundances using the standard Press-Schechter criterion and find that overproduction of PBHs is likely in some regions of the chaotic preheating parameter-space.

### Numerical simulation of oscillating magnetospheres with resistive electrodynamics

We present a model of the magnetosphere around an oscillating neutron star. The electromagnetic fields are numerically solved by modeling electric charge and current induced by the stellar torsional mode, with particular emphasis on outgoing radiation passing through the magnetosphere. The current is modeled using Ohm’s law, whereby an increase in conductivity results in an increase in the induced current. As a result, the fields are drastically modified, and energy flux is thereby enhanced. This behavior is however localized in the vicinity of the surface since the induced current disappears outwardly in our model, in which the exterior is assumed to gradually approach a vacuum.

### Synthetic Observations of the Evolution of Starless Cores in a Molecular Cloud Simulation: Comparisons with JCMT Data and Predictions for ALMA

Interpreting the nature of starless cores has been a prominent goal in star formation for many years. In order to characterise the evolutionary stages of these objects, we perform synthetic observations of a numerical simulation of a turbulent molecular cloud. We find that nearly all cores that we detect are associated with filaments and eventually form protostars. We conclude that observed starless cores which appear Jeans unstable are only marginally larger than their respective Jeans masses (within a factor of 3). We note single dish observations such as those performed with the JCMT appear to miss significant core structure on small scales due to beam averaging. Finally, we predict that interferometric observations with ALMA Cycle 1 will resolve the important small scale structure, which has so far been missed by mm-wavelength observations.

### Synthetic Observations of the Evolution of Starless Cores in a Molecular Cloud Simulation: Comparisons with JCMT Data and Predictions for ALMA [Replacement]

Interpreting the nature of starless cores has been a prominent goal in star formation for many years. In order to characterise the evolutionary stages of these objects, we perform synthetic observations of a numerical simulation of a turbulent molecular cloud. We find that nearly all cores that we detect are associated with filaments and eventually form protostars. We conclude that observed starless cores which appear Jeans unstable are only marginally larger than their respective Jeans masses (within a factor of 3). We note single dish observations such as those performed with the JCMT appear to miss significant core structure on small scales due to beam averaging. Finally, we predict that interferometric observations with ALMA Cycle 1 will resolve the important small scale structure, which has so far been missed by mm-wavelength observations.

### Chaos and Turbulent Nucleosynthesis Prior to a Supernova Explosion [Replacement]

Three-dimensional (3D), time dependent numerical simulations, of flow of matter in stars, now have sufficient resolution to be fully turbulent. The late stages of the evolution of massive stars, leading up to core collapse to a neutron star (or black hole), and often to supernova explosion and nucleosynthesis, are strongly convective because of vigorous neutrino cooling and nuclear heating. Unlike models based on current stellar evolutionary practice, these simulations show a chaotic dynamics characteristic of highly turbulent flow. Theoretical analysis of this flow, both in the Reynolds-averaged Navier-Stokes (RANS) framework and by simple dynamic models, show an encouraging consistency with the numerical results. It may now be possible to develop physically realistic and robust procedures for convection and mixing which (unlike 3D numerical simulation) may be applied throughout the long life times of stars. In addition, a new picture of the presupernova stages is emerging which is more dynamic and interesting (i.e., predictive of new and newly observed phenomena) than our previous one.

### Chaos and Turbulent Nucleosynthesis Prior to a Supernova Explosion

Three-dimensional (3D), time dependent numerical simulations, of flow of matter in stars, now have sufficient resolution to be fully turbulent. The late stages of the evolution of massive stars, leading up to core collapse to a neutron star (or black hole), and often to supernova explosion and nucleosynthesis, are strongly convective because of vigorous neutrino cooling and nuclear heating. Unlike models based on current stellar evolutionary practice, these simulations show a chaotic dynamics characteristic of highly turbulent flow. Theoretical analysis of this flow, both in the Reynolds-averaged Navier-Stokes (RANS) framework and by simple dynamic models, show an encouraging consistency with the numerical results. It may now be possible to develop physically realistic and robust procedures for convection and mixing which (unlike 3D numerical simulation) may be applied throughout the long life times of stars. In addition, a new picture of the presupernova stages is emerging which is more dynamic and interesting (i.e., predictive of new and newly observed phenomena) than our previous one.

### Numerical simulation of the electron capture process in a magnetar interior

In a superhigh magnetic field, direct Urca reactions can proceed for an arbitrary proton concentration. Since only the electrons with high energy $E$ ($E > Q$, $Q$ is the threshold energy of inverse $\beta-$decay) at large Landau levels can be captured, we introduce the Landau level effect coefficient $q$ and the effective electron capture rate $\Gamma_{\rm eff}$. By using $\Gamma_{\rm eff}$, the values of $L_{\rm X}$ and $L_{\rm \nu}$ are calculated, where and $L_{\rm \nu}$, $L_{\rm X}$ are the average neutrino luminosity of Anomalous X-ray Pulsars (AXPs) and the average X-ray luminosity of AXPs, respectively. The complete process of electron capture inside a magnetar is simulated numerically.

### GLAMER Part II: Multiple Plane Gravitational Lensing

We present an extension to multiple planes of the gravitational lensing code {\small GLAMER}. The method entails projecting the mass in the observed light-cone onto a discrete number of lens planes and inverse ray-shooting from the image to the source plane. The mass on each plane can be represented as halos, simulation particles, a projected mass map extracted form a numerical simulation or any combination of these. The image finding is done in a source oriented fashion, where only regions of interest are iteratively refined on an initially coarse image plane grid. The calculations are performed in parallel on shared memory machines. The code is able to handle different types of analytic halos (NFW, NSIE, power-law, etc.), haloes extracted from numerical simulations and clusters constructed from semi-analytic models ({\small MOKA}). Likewise, there are several different options for modeling the source(s) which can be distributed throughout the light-cone. The distribution of matter in the light-cone can be either taken from a pre-existing N-body numerical simulations, from halo catalogs, or are generated from an analytic mass function. We present several tests of the code and demonstrate some of its applications such as generating mock images of galaxy and galaxy cluster lenses.

### Nonlinear transport of Cosmic Rays in turbulent magnetic field

Recent advances in both the MHD turbulence theory and cosmic ray observations call for revisions in the paradigm of cosmic ray transport. We use the models of magnetohydrodynamic turbulence that were tested in numerical simulation, in which turbulence is injected at large scale and cascades to to small scales. We shall present the nonlinear results for cosmic ray transport, in particular, the cross field transport of CRs and demonstrate that the concept of cosmic ray subdiffusion in general does not apply and the perpendicular motion is well described by normal diffusion with M_A^4 dependence. Moreover, on scales less than injection scale of turbulence, CRs’ transport becomes super-diffusive. Quantitative predictions for both the normal diffusion on large scale and super diffusion are confronted with recent numerical simulations. Implication for shock acceleration is briefly discussed.

### Nonlinear transport of Cosmic Rays in turbulent magnetic field [Replacement]

Recent advances in both the MHD turbulence theory and cosmic ray observations call for revisions in the paradigm of cosmic ray transport. We use the models of magnetohydrodynamic turbulence that were tested in numerical simulation, in which turbulence is injected at large scale and cascades to to small scales. We shall present the nonlinear results for cosmic ray transport, in particular, the cross field transport of CRs and demonstrate that the concept of cosmic ray subdiffusion in general does not apply and the perpendicular motion is well described by normal diffusion with M_A^4 dependence. Moreover, on scales less than injection scale of turbulence, CRs’ transport becomes super-diffusive. Quantitative predictions for both the normal diffusion on large scale and super diffusion are confronted with recent numerical simulations. Implication for shock acceleration is briefly discussed.

### Numerical Simulation of Spectral and Timing Properties of a Two Component Advective Flow around a Black Hole

We study the spectral and timing properties of a two component advective flow (TCAF) around a black hole by numerical simulation. Several cases have been simulated by varying the Keplerian disk rate and the resulting spectra and lightcurves have been produced for all the cases. The dependence of the spectral states and quasi-periodic oscillation (QPO) frequencies on the flow parameters is discussed. We also find the earlier explanation of arising of QPOs as the resonance between infall time scale and cooling time scale remain valid even for Compton cooling.

### A new GPU-accelerated hydrodynamical code for numerical simulation of interacting galaxies

In this paper a new scalable hydrodynamic code GPUPEGAS (GPU-accelerated PErformance Gas Astrophysic Simulation) for simulation of interacting galaxies is proposed. The code is based on combination of Godunov method as well as on the original implementation of FlIC method, specially adapted for GPU-implementation. Fast Fourier Transform is used for Poisson equation solution in GPUPEGAS. Software implementation of the above methods was tested on classical gas dynamics problems, new Aksenov’s test and classical gravitational gas dynamics problems. Collisionless hydrodynamic approach was used for modelling of stars and dark matter. The scalability of GPUPEGAS computational accelerators is shown.

### Probing the structure of local magnetic field of solar features with helioseismology

Motivated by the problem of local solar subsurface magnetic structure, we have used numerical simulation to investigate the propagation of waves through monolithic magnetic flux tubes of different size. A cluster model can be a good approximation to simulate sunspots as well as solar plage regions which are composed of an ensemble of compactly packed thin flux tubes. Simulations of this type is a powerful tool to probe the structure and the dynamic of various solar features which are related directly to solar magnetic field activity.

### Numerical Simulation of Establishment of thermodynamic equilibrium in cosmological model with an arbitrary acceleration

Results of numerical simulation constructed before strict mathematical model of an establishment of thermodynamic equilibrium in originally nonequilibrium cosmological ultrarelativistic plasma for the Universe with any acceleration in the assumption of restoration of a scalling of interactions of elementary particles are presented at energies above a unitary limit. Limiting parametres of residual nonequilibrium distribution of nonequilibrium relic particles of ultrahigh energies are found.

### Using 3D Voronoi grids in radiative transfer simulations

Probing the structure of complex astrophysical objects requires effective three-dimensional (3D) numerical simulation of the relevant radiative transfer (RT) processes. As with any numerical simulation code, the choice of an appropriate discretization is crucial. Adaptive grids with cuboidal cells such as octrees have proven very popular, however several recently introduced hydrodynamical and RT codes are based on a Voronoi tessellation of the spatial domain. Such an unstructured grid poses new challenges in laying down the rays (straight paths) needed in RT codes. We show that it is straightforward to implement accurate and efficient RT on 3D Voronoi grids. We present a method for computing straight paths between two arbitrary points through a 3D Voronoi grid in the context of a RT code. We implement such a grid in our RT code SKIRT, using the open source library Voro++ to obtain the relevant properties of the Voronoi grid cells based solely on the generating points. We compare the results obtained through the Voronoi grid with those generated by an octree grid for two synthetic models, and we perform the well-known Pascucci RT benchmark using the Voronoi grid. The presented algorithm produces correct results for our test models. Shooting photon packages through the geometrically much more complex 3D Voronoi grid is only about three times slower than the equivalent process in an octree grid with the same number of cells, while in fact the total number of Voronoi grid cells may be lower for an equally good representation of the density field. We conclude that the benefits of using a Voronoi grid in RT simulation codes will often outweigh the somewhat slower performance.

### Realistic Simulations of Stellar Surface Convection with ANTARES: I. Boundary Conditions and Model Relaxation [Replacement]

We have implemented open boundary conditions into the ANTARES code to increase the realism of our simulations of stellar surface convection. Even though we greatly benefit from the high accuracy of our fifth order numerical scheme (WENO5), the broader stencils needed for the numerical scheme complicate the implementation of boundary conditions. We show that the effective temperature of a numerical simulation cannot be changed by corrections at the lower boundary since the thermal stratification does only change on the Kelvin-Helmholtz time scale. Except for very shallow models, this time scale cannot be covered by multidimensional simulations due to the enormous computational requirements. We demonstrate to what extent numerical simulations of stellar surface convection are sensitive to the initial conditions and the boundary conditions. An ill-conceived choice of parameters for the boundary conditions can have a severe impact. Numerical simulations of stellar surface convection will only be (physically) meaningful and realistic if the initial model, the extent and position of the simulation box, and the parameters from the boundary conditions are chosen adequately.

### Nonuniformity effects in the negative effective magnetic pressure instability

Using direct numerical simulations (DNS) and mean-field simulations (MFS), the effects of non-uniformity of the magnetic field on the suppression of the turbulent pressure is investigated. This suppression of turbulent pressure can lead to an instability which, in turn, makes the mean magnetic field even more non-uniform. This large-scale instability is caused by a resulting negative contribution of small-scale turbulence to the effective (mean-field) magnetic pressure. We show that enhanced mean current density increases the suppression of the turbulent pressure. The instability leads to magnetic flux concentrations in which the normalized effective mean-field pressure is reduced to a certain value at all depths within a structure.

### The Burst Mode of Accretion in Primordial Star Formation

We present simulation results for the formation and long-term evolution of a primordial protostellar disk harbored by a first star. Using a 2+1D nonaxisymmetric thin disk numerical simulation, together with a barotropic relation for the gas, we are able to probe ~20 kyr of the disk’s evolution. During this time period we observe fragmentation leading to loosely bound gaseous clumps within the disk. These are then torqued inward and accreted onto the growing protostar, giving rise to a burst phenomenon. The luminous feedback produced by this mechanism may have important consequences for the subsequent growth of the protostar.

### Flavor stability analysis of dense supernova neutrinos with flavor-dependent angular distributions [Cross-Listing]

Numerical simulations of the supernova (SN) neutrino self-induced flavor conversions, associated with the neutrino-neutrino interactions in the deepest stellar regions, have been typically carried out assuming the "bulb-model". In this approximation, neutrinos are taken to be emitted half-isotropically by a common neutrinosphere. In the recent Ref. \cite{Mirizzi:2011tu} we have removed this assumption by introducing flavor-dependent angular distributions for SN neutrinos, as suggested by core-collapse simulations. We have found that in this case a novel multi-angle instability in the self-induced flavor transitions can arise. In this work we perform an extensive study of this effect, carrying out a linearized flavor stability analysis for different SN neutrino energy fluxes and angular distributions, in both normal and inverted neutrino mass hierarchy. We confirm that spectra of different nu species which cross in angular space (where F_{\nu_e}=F_{\nu_x} and F_{\bar\nu_e}=F_{\bar\nu_x}) present a significant enhancement of the flavor instability, and a shift of the onset of the flavor conversions at smaller radii with respect to the case of an isotropic neutrino emission. We also illustrate how a qualitative (and sometimes quantitative) understanding of the dynamics of these systems follows from a stability analysis.

### Data assimilation for stratified convection

We show how the 3DVAR data assimilation methodology can be used in the astrophysical context of a two-dimensional convection flow. We study the way this variational approach finds best estimates of the current state of the flow from a weighted average of model states and observations. We use numerical simulations to generate synthetic observations of a vertical two-dimensional slice of the outer part of the solar convection zone for varying noise levels and implement 3DVAR when the covariance matrices are scalar. Our simulation results demonstrate the capability of 3DVAR to produce error estimates of system states between up to tree orders of magnitude below the original noise level present in the observations. This work exemplifies the importance of applying data assimilation techniques in simulations of the stratified convection.

### Synthetic X-ray spectra for simulations of the dynamics of an accretion flow irradiated by a quasar

Ultraviolet and X-ray observations show evidence of outflowing gas around many active galactic nuclei. Some of these outflows may be driven off gas infalling towards the central black hole. We perform radiative transfer calculations to compute the gas ionization state and X-ray spectra for two- and three-dimensional (3D) hydrodynamical simulations of this outflow-from-inflow scenario. By comparison with observations, our results can be used to test the theoretical models and guide future numerical simulations. We predict both absorption and emission features, most of which are formed in a polar funnel of outflowing gas. This outflow causes strong absorption for observer orientation angles of < 35 degrees. Particularly in 3D, the strength of this absorption varies significantly for different lines-of-sight owing to the fragmentary structure of the gas flow. Although infalling material occupies a large fraction of the simulation volume, we do not find that it imprints strong absorption features since the ionization state is very high. Thus, an absence of observed inflow absorption features does not exclude the models. The main spectroscopic consequence of the infalling gas is a scattered continuum component that partially re-fills the absorption features caused by the outflowing polar funnel. Fluorescence and scattering in the outflow is predicted to give rise to several emission features for all observer orientations. For the hydrodynamical simulations considered we find both ionization states and column densities for the outflowing gas that are too high to be quantitatively consistent with well-observed X-ray absorption systems. Nevertheless, our results are qualitatively encouraging and further exploration of the model parameter space is warranted. (Abridged.)

### The Rossiter-McLaughlin effect for exomoons or binary planets

In this paper we study possible signatures of binary planets or exomoons on the Rossiter-McLaughlin (R-M) effect. Our analyses show that the R-M effect for a binary planet or exomoon during its complete transit phase can be divided into two parts. The first is the conventional one similar to the R-M effect from the transit of a single planet, of which the mass and the projected area are the combinations of the binary components; and the second is caused by the orbital rotation of the binary components, which may add a sine- or linear-mode deviation to the stellar radial velocity curve. We find that the latter effect can be up to several or several ten m/s. By doing numerical simulations as well as analytical analyses, we illustrate that the distribution and dispersion of the latter effects obtained from multiple transit events can be used to constrain the dynamical configuration of the binary planet, such as, how the inner orbit of the binary planet is inclined to its orbit rotating around the central star. We find that the signatures caused by the orbital rotation of the binary components are more likely to be revealed if the two components of binary planet have different masses and mass densities, especially if the heavy one has a high mass density and the light one has a low density. Similar signature on the R-M effect may also be revealed in a hierarchical triple star system containing a dark compact binary and a tertiary star.

### 3D simulations of globules and pillars formation around HII regions: turbulence and shock curvature

We investigate the interplay between the ionization radiation from massive stars and the turbulence inside the surrounding molecular gas thanks to 3D numerical simulations. We used the 3D hydrodynamical code HERACLES to model an initial turbulent medium that is ionized and heated by an ionizing source. Three different simulations are performed with different mean Mach numbers (1, 2 and 4). A non-equilibrium model for the ionization and the associated thermal processes was used. This revealed to be crucial when turbulent ram pressure is of the same order as the ionized-gas pressure. The density structures initiated by the turbulence cause local curvatures of the dense shell formed by the ionization compression. When the curvature of the shell is sufficient, the shell collapse on itself to form a pillar while a smaller curvature leads to the formation of dense clumps that are accelerated with the shell and therefore remain in the shell during the simulation. When the turbulent ram pressure of the cold gas is sufficient to balance the ionized-gas pressure, some dense-gas bubbles have enough kinetic energy to penetrate inside the ionized medium, forming cometary globules. This suggests a direct relation in the observations between the presence of globules and the relative importance of the turbulence compared to the ionized-gas pressure. The probability density functions present a double peak structure when the turbulence is low relative to the ionized-gas pressure. This could be used in observations as an indication of the turbulence inside molecular clouds.

### Black Hole-Neutron Star Mergers: Disk Mass Predictions

Determining the final result of black hole-neutron star mergers, and in particular the amount of matter remaining outside the black hole at late times, has been one of the main motivations behind the numerical simulation of these systems. Black hole-neutron star binaries are amongst the most likely progenitors of short gamma-ray bursts — as long as they result in the formation of massive (at least ~0.1 solar mass) accretion disks around the black hole. Whether this actually happens strongly depends on the physical characteristics of the system, and in particular on the mass ratio, the spin of the black hole, and the radius of the neutron star. We present here a simple two-parameter model, fitted to existing numerical results, for the determination of the mass remaining outside the black hole a few milliseconds after a black hole-neutron star merger. This model predicts the remnant mass within a few percents of the mass of the neutron star, at least for remnant masses up to 20% of the neutron star mass. Results across the range of parameters deemed to be the most likely astrophysically are presented here. We find that, for 10 solar mass black holes, massive disks are only possible for fairly large neutron stars (R>12km), or quasi-extremal black hole spins (a/M>0.9). We also use our model to discuss how the equation of state of the neutron star affects the final remnant, and the strong influence that this can have on the rate of short gamma-ray bursts produced by black hole-neutron star mergers.

### The thick disk rotation-metallicity correlation as a fossil of an "inverse chemical gradient" in the early Galaxy

The thick disk rotation–metallicity correlation, \partial V_\phi/\partial[Fe/H] =40\div 50 km s^{-1}dex^{-1} represents an important signature of the formation processes of the galactic disk. We use nondissipative numerical simulations to follow the evolution of a Milky Way (MW)-like disk to verify if secular dynamical processes can account for this correlation in the old thick disk stellar population. We followed the evolution of an ancient disk population represented by 10 million particles whose chemical abundances were assigned by assuming a cosmologically plausible radial metallicity gradient with lower metallicity in the inner regions, as expected for the 10-Gyr-old MW. Essentially, inner disk stars move towards the outer regions and populate layers located at higher |z|. A rotation–metallicity correlation appears, which well resembles the behaviour observed in our Galaxy at a galactocentric distance between 8 kpc and 10 kpc. In particular,we measure a correlation of \partial V_\phi/\partial[Fe/H]\simeq 60 km s^{-1}dex^{-1} for particles at 1.5 kpc < |z| < 2.0 kpc that persists up to 6 Gyr. Our pure N-body models can account for the V_\phi vs. [Fe/H] correlation observed in the thick disk of our Galaxy, suggesting that processes internal to the disk such as heating and radial migration play a role in the formation of this old stellar component. In this scenario, the positive rotation-metallicity correlation of the old thick disk population would represent the relic signature of an ancient "inverse" chemical (radial) gradient in the inner Galaxy, which resulted from accretion of primordial gas.

### Temporal Evolution of Velocity and Magnetic Field in and around Umbral Dots

We study the temporal evolution of umbral dots (UDs) using measurements from the CRISP imaging spectropolarimeter at the Swedish 1-m Solar Telescope. Scans of the magnetically sensitive 630 nm iron lines were performed under stable atmospheric conditions for 71 min with a cadence of 63 s. These observations allow us to investigate the magnetic field and velocity in and around UDs at a resolution approaching 0.13". From the analysis of 339 UDs, we draw the following conclusions: (1)UDs show clear hints of upflows, as predicted by magnetohydrodynamic (MHD) simulations. By contrast, we could not find systematic downflow signals. Only in very deep layers we detect localized downflows around UDs, but they do not persist in time. (2) We confirm that UDs exhibit weaker and more inclined fields than their surroundings, as reported previously. However, UDs that have strong fields above 2000 G or are in the decay phase show enhanced and more vertical fields. (3)There are enhanced fields at the migration front of UDs detached from penumbral grains, as if their motion were impeded by the ambient field. (4) Long-lived UDs travel longer distances with slower proper motions. Our results appear to confirm some aspects of recent numerical simulations of magnetoconvection in the umbra (e.g., the existence of upflows in UDs), but not others (e.g., the systematic weakening of the magnetic field at the position of UDs.)

### A Mexican Hat with holes: calculating low resolution power spectra from data with gaps

A simple method for calculating a low-resolution power spectrum from data with gaps is described. The method is a modification of the $\Delta$-variance method previously described by Stutzki and Ossenkopf. A Mexican Hat filter is used to single out fluctuations at a given spatial scale and the variance of the convolved image is calculated. The gaps in the image, defined by the mask, are corrected for by representing the Mexican Hat filter as a difference between two Gaussian filters with slightly different widths, convolving the image and mask with these filters and dividing the results before calculating the final filtered image. This method cleanly compensates for data gaps even if these have complicated shapes and cover a significant fraction of the data. The method was developed to deal with problematic 2D images, where irregular detector edges and masking of contaminating sources compromise the power spectrum estimates, but it can also be straightforwardly applied to 1D timing analysis or 3D data cubes from numerical simulations.

### Dark Matter Subhalos in the Ursa Minor Dwarf Galaxy

Through numerical simulations, we study the dissolution timescale of the Ursa Minor cold stellar clump, due to the combination of phase-mixing and gravitational encounters with compact dark substructures in the halo of Ursa Minor. We compare two scenarios; one where the dark halo is made up by a smooth mass distribution of light particles and one where the halo contains 10% of its mass in the form of substructures (subhalos). In a smooth halo, the stellar clump survives for a Hubble time provided that the dark matter halo has a big core. In contrast, when the point-mass dark substructures are added, the clump survives barely for \sim 1.5 Gyr. These results suggest a strong test to the \Lambda-cold dark matter scenario at dwarf galaxy scale.

### The impact of magnetic fields on the IMF in star-forming clouds near a supermassive black hole

Star formation in the centers of galaxies is thought to yield massive stars with a possibly top-heavy stellar mass distribution. It is likely that magnetic fields play a crucial role in the distribution of stellar masses inside star-forming molecular clouds. In this context, we explore the effects of magnetic fields, with a typical field strength of 38 {\mu}G, such as in RCW 38, and a field strength of 135 {\mu}G, similar to NGC 2024 and the infrared dark cloud G28.34+0.06, on the initial mass function (IMF) near (\leq 10 pc) a 10^7 solar mass black hole. Using these conditions, we perform a series of numerical simulations with the hydrodynamical code FLASH to elucidate the impact of magnetic fields on the IMF and the star-formation efficiency (SFE) emerging from an 800 solar mass cloud. We find that the collapse of a gravitationally unstable molecular cloud is slowed down with increasing magnetic field strength and that stars form along the field lines. The total number of stars formed during the simulations increases by a factor of 1.5-2 with magnetic fields. The main component of the IMF has a lognormal shape, with its peak shifted to sub-solar (\leq 0.3 Mo) masses in the presence of magnetic fields, due to a decrease in the accretion rates from the gas reservoir. In addition, we see a top-heavy, nearly flat IMF above \sim 2 solar masses, from regions that were supported by magnetic pressure until high masses are reached. We also consider the effects of X-ray irradiation if the central black hole is active. X-ray feedback inhibits the formation of sub-solar masses and decreases the SFEs even further. Thus, the second contribution is no longer visible. We conclude that magnetic fields potentially change the SFE and the IMF both in active and inactive galaxies, and need to be taken into account in such calculations.

### The impact of magnetic fields on the IMF in star-forming clouds near a supermassive black hole [Replacement]

Star formation in the centers of galaxies is thought to yield massive stars with a possibly top-heavy stellar mass distribution. It is likely that magnetic fields play a crucial role in the distribution of stellar masses inside star-forming molecular clouds. In this context, we explore the effects of magnetic fields, with a typical field strength of 38 microG, such as in RCW 38, and a field strength of 135 microG, similar to NGC 2024 and the infrared dark cloud G28.34+0.06, on the initial mass function (IMF) near (< 10 pc) a 10^7 solar mass black hole. Using these conditions, we perform a series of numerical simulations with the hydrodynamical code FLASH to elucidate the impact of magnetic fields on the IMF and the star-formation efficiency (SFE) emerging from an 800 solar mass cloud. We find that the collapse of a gravitationally unstable molecular cloud is slowed down with increasing magnetic field strength and that stars form along the field lines. The total number of stars formed during the simulations increases by a factor of 1.5-2 with magnetic fields. The main component of the IMF has a lognormal shape, with its peak shifted to sub-solar (< 0.3 M_sun) masses in the presence of magnetic fields, due to a decrease in the accretion rates from the gas reservoir. In addition, we see a top-heavy, nearly flat IMF above ~2 solar masses, from regions that were supported by magnetic pressure until high masses are reached. We also consider the effects of X-ray irradiation if the central black hole is active. X-ray feedback inhibits the formation of sub-solar masses and decreases the SFEs even further. Thus, the second contribution is no longer visible. We conclude that magnetic fields potentially change the SFE and the IMF both in active and inactive galaxies, and need to be taken into account in such calculations.

### Rotational effects on the negative magnetic pressure instability

The surface layers of the Sun are strongly stratified. In the presence of turbulence with a weak mean magnetic field, a large-scale instability resulting in the formation of non-uniform magnetic structures, can be excited over the scale of many turbulent eddies or convection cells. This instability is caused by a negative contribution of turbulence to the effective (mean-field) magnetic pressure and has previously been discussed in connection with the formation of active regions and perhaps sunspots. We want to understand the effects of rotation on this instability in both two and three dimensions. We use mean-field magnetohydrodynamics in a parameter regime in which the properties of the negative effective magnetic pressure instability have previously been found to be in agreement with those of direct numerical simulations. We find that the instability is suppressed already for relatively slow rotation with Coriolis numbers (i.e. inverse Rossby numbers) around 0.2. The suppression is strongest at the equator. In the nonlinear regime, we find traveling wave solutions with propagation in the prograde direction at the equator with additional poleward migration away from the equator. The prograde rotation of the magnetic pattern near the equator is argued to be a possible explanation for the faster rotation speed of magnetic tracers found on the Sun. In the bulk of the domain, kinetic and current helicities are negative in the northern hemisphere and positive in the southern.

### Numerical simulation code for self-gravitating Bose-Einstein condensates [Replacement]

We completed the development of simulation code that is designed to study the behavior of a conjectured dark matter galactic halo that is in the form of a Bose-Einstein Condensate (BEC). The BEC is described by the Gross-Pitaevskii equation, which can be solved numerically using the Crank-Nicholson method. The gravitational potential, in turn, is described by Poisson’s equation, that can be solved using the relaxation method. Our code combines these two methods to study the time evolution of a self-gravitating BEC. The inefficiency of the relaxation method is balanced by the fact that in subsequent time iterations, previously computed values of the gravitational field serve as very good initial estimates. The code is robust (as evidenced by its stability on coarse grids) and efficient enough to simulate the evolution of a system over the course of 1E9 years using a finer (100x100x100) spatial grid, in less than a day of processor time on a contemporary desktop computer.

### Numerical simulation code for self-gravitating Bose-Einstein condensates

We completed the development of simulation code that is designed to study the behavior of a conjectured dark matter galactic halo that is in the form of a Bose-Einstein Condensate (BEC). The BEC is described by the Gross-Pitaevskii equation, which can be solved numerically using the Crank-Nicholson method. The gravitational potential, in turn, is described by Poisson’s equation, that can be solved using the relaxation method. Our code combines these two methods to study the time evolution of a self-gravitating BEC. The inefficiency of the relaxation method is balanced by the fact that in subsequent time iterations, previously computed values of the gravitational field serve as very good initial estimates. The code is robust (as evidenced by its stability on coarse grids) and efficient enough to simulate the evolution of a system over the course of 1E9 years using a finer (100\times100\times100) spatial grid, in less than a day of processor time on a contemporary desktop computer.

### Hysteresis and thermal limit cycles in MRI simulations of accretion discs

The recurrent outbursts that characterise low-mass binary systems reflect thermal state changes in their associated accretion discs. The observed outbursts are connected to the strong variation in disc opacity as hydrogen ionises near 5000 K. This leads to accretion disc models that exhibit bistability and thermal limit cycles, whereby the disc jumps between a family of cool and low accreting states and a family of hot and efficiently accreting states. Previous models have parametrised the turbulence via an alpha (or `eddy’) viscosity. In this paper we treat the turbulence more realistically via a suite of numerical simulations of the magnetorotational instability (MRI) in local geometry. Radiative cooling is included via a simple but physically motivated prescription. We show the existence of bistable equilibria and thus the prospect of thermal limit cycles, and in so doing demonstrate that MRI-induced turbulence is compatible with the classical theory. Our simulations also show that the turbulent stress and pressure perturbations are only weakly dependent on each other; as a consequence, thermal instability connected to variations in turbulent heating (as opposed to radiative cooling) is unlikely to operate, in agreement with previous numerical results. Our work presents a first step towards unifying simulations of full MHD turbulence with the correct thermal and radiative physics of the outbursting discs associated with dwarf novae, low-mass X-ray binaries, and possibly young stellar objects.

### The impact of baryons on the spins and shapes of dark matter haloes

We use numerical simulations to investigate how the statistical properties of dark matter (DM) haloes are affected by the baryonic processes associated with galaxy formation. We focus on how these processes influence the spin and shape of a large number of DM haloes covering a wide range of mass scales, from dwarf galaxies to clusters, at redshifts zero, one and two. The haloes are extracted from the OverWhelmingly Large Simulations, a suite of state-of-the-art high-resolution cosmological simulations run with a range of feedback prescriptions. We find that the median spin parameter in DM-only simulations is independent of mass, redshift and cosmology. Baryons increase the spin of the DM in the central region (< 0.25r_{200}) by up to 50 per cent when feedback is weak or absent. This increase can be attributed to the transfer of angular momentum from baryons to the DM. We also present fits to the mass dependence of the DM halo shape at both low and high redshift. At z=0 the sphericity (triaxiality) is negatively (positively) correlated with halo mass and both results are independent of cosmology. Interestingly, these mass-dependent trends are markedly weaker at z=2. While the cooling of baryons acts to make the overall DM halo more spherical, stronger feedback prescriptions (e.g. from active galactic nuclei) tend to reduce the impact of baryons by reducing the central halo mass concentration. More generally, we demonstrate a strongly positive (negative) correlation between halo sphericity (triaxiality) and galaxy formation efficiency, with the latter measured using the central halo baryon fraction. In conclusion, our results suggest that the effects of baryons on the DM halo spin and shape are minor when the effects of cooling are mitigated, as required by realistic models of galaxy formation, although they remain significant for the inner halo.

### The nature of assembly bias - II. Halo spin

We study an assembly-type bias parametrized by the dimensionless spin parameter that affects massive structures. In numerical simulations higher spin haloes are more strongly clustered than lower spin haloes of equal mass. We detect a difference of over a 30 per cent in the clustering strength for dark matter haloes of 10^13-10^14 Msun, which is similar to the result of Bett et al. We explore whether the dependence of clustering strength on halo spin is removed if we apply the redefinition of overdensity peak height proposed by Lacerna & Padilla (Paper I) obtained using assembly ages. We find that this is not the case due to two reasons. Firstly, only a few objects of low-virial mass are moved into the mass range where the spin introduces an assembly bias after using this redefinition. Secondly, this formalism does not alter the mass of massive objects. We then repeat the process of finding the redefined peak height of Paper I but using the spin. In this case, the new masses show no spin-related assembly bias but they introduce a previously absent assembly bias with respect to relative age. From this result, we conclude that the assembly-type bias with respect to the halo spin has a different origin than with respect to assembly age. The former may be due to the material from filaments, which is accreted by massive haloes, that is enhanced in high-density environments, thus causing more extreme spin values without significantly changing the formation age of the halo. In addition, high-mass objects may correspond, in some cases, to a different peak height than that suggested by their mass in numerical simulations, providing a possible explanation for the assembly bias with respect to spin. (abridged)

### The stability of stratified, rotating systems and the generation of vorticity in the Sun

We examine the linear behavior of three-dimensional Lagrangian displacements in a stratified, shearing background. The isentropic and iso-rotation surfaces of the equilibrium flow are assumed to be axisymmetric, but otherwise fully two-dimensional. Three-dimensional magnetic fields are included in the perturbation equations; however the equilibrium is assumed to be well-described by purely hydrodynamic forces. The model, in principle very general, is used to study the behavior of fluid displacements in an environment resembling the solar convection zone. Some very suggestive results emerge. All but high-latitude displacements align themselves with the observed surfaces of constant angular velocity. The tendency for the angular velocity to remain constant with depth in the bulk of the convective zone, together with other critical features of the rotation profile, emerge from little more than a visual inspection of the governing equation. In the absence of a background axial angular velocity gradient, displacements exhibit no poleward bias, suggesting that solar convection "plays-off" of prexisting shear rather than creates it. We argue that baroclinic vorticity of precisely the right order is generated at the radiative/convective zone boundary due to centrifugal distortion of equipotential surfaces that is not precisely followed by isothermal surfaces. If so, many features of the Sun’s internal rotation become more clear, including: i) the general appearance of the tachocline; ii) the extension of differential rotation well into the radiative zone; iii) the abrupt change of morphology of convective zone isorotation surfaces; and iv) the inability of current numerical simulations to reproduce the solar rotation profile without imposed entropy boundary conditions.

### Thermal and Photophoretic Properties of Dust Mantled Chondrules and Sorting in the Solar Nebula

Many chondrules are found to be surrounded by a fine grained rim. Supposedly, these rims were acquired from dust accreted to the chondrules in a protoplanetary disk. In numerical simulations we study the heat transfer in illuminated bare and dust mantled chondrules. The calculations consider the chondrule size, the dust mantle size, and the thermal conductivities of both components as parameters. We calculate the photophoretic force and compare the numerical results to analytical approximations. We give an expression to quantify the photophoretic force on a spherical particle in the free molecular regime to better than 2%. We describe the influence of a dust mantle on the photophoretic strength by an effective thermal conductivity of the core-mantle particle. The effective thermal conductivity significantly depends on the size ratio between mantle and chondrule but not on the absolute sizes. It also strongly depends on the thermal conductivity of the mantle with minor influence of the thermal conductivity of the chondrule. The size ratio between rim and chondrule in meteorites is found to vary systematically with overall size by other authors. Based on this, our calculations show that a photophoretic size sorting can occur for dust mantled chondrules in optically thin disks or at the evolving inner edge of the solar nebula.

### Infinite impulse response modal filtering in visible adaptive optics

Diffraction limited resolution adaptive optics (AO) correction in visible wavelengths requires a high performance control. In this paper we investigate infinite impulse response filters that optimize the wavefront correction: we tested these algorithms through full numerical simulations of a single-conjugate AO system comprising an adaptive secondary mirror with 1127 actuators and a pyramid wavefront sensor (WFS). The actual practicability of the algorithms depends on both robustness and knowledge of the real system: errors in the system model may even worsen the performance. In particular we checked the robustness of the algorithms in different conditions, proving that the proposed method can reject both disturbance and calibration errors.

### Extremely fast focal-plane wavefront sensing for extreme adaptive optics

We present a promising approach to the extremely fast sensing and correction of small wavefront errors in adaptive optics systems. As our algorithm’s computational complexity is roughly proportional to the number of actuators, it is particularly suitable to systems with 10,000 to 100,000 actuators. Our approach is based on sequential phase diversity and simple relations between the point-spread function and the wavefront error in the case of small aberrations. The particular choice of phase diversity, introduced by the deformable mirror itself, minimizes the wavefront error as well as the computational complexity. The method is well suited for high-contrast astronomical imaging of point sources such as the direct detection and characterization of exoplanets around stars, and it works even in the presence of a coronagraph that suppresses the di?raction pattern. The accompanying paper in these proceedings by Korkiakoski et al. describes the performance of the algorithm using numerical simulations and laboratory tests.

### Vertical density waves in the Milky Way disc induced by the Sagittarius Dwarf Galaxy

Recently, Widrow and collaborators announced the discovery of vertical density waves in the Milky Way disk. Here we investigate a scenario where these waves were induced by the Sagittarius dwarf galaxy as it plunged through the Galaxy. Using numerical simulations, we find that the Sagittarius impact produces North-South asymmetries and vertical wave-like behavior that qualitatively agrees with what is observed. The extent to which vertical modes can radially penetrate into the disc, as well as their amplitudes, depend on the mass of the perturbing satellite. We show that the mean height of the disc is expected to vary more rapidly in the radial than in the azimuthal direction. If the observed vertical density asymmetry is indeed caused by vertical oscillations, we predict radial and azimuthal variations of the mean vertical velocity, correlating with the spatial structure. These variations can have amplitudes as large as 8 km/s.

### Long Duration X-Ray Flash and X-Ray Rich Gamma Ray Burst from Low Mass Population III Star

Recent numerical simulations suggest that Population III (Pop III) stars are born with masses not larger than $\sim 100M_\odot$ but typically $\sim 40M_{\odot}$. We investigate whether such a low mass Pop III star can raise a Gamma Ray Burst (GRB) by considering the propagation of a jet, which is launched from the black hole, in the stellar envelope. It is generally believed that a super giant star is not an appropriate progenitor of a GRB, since the large envelope prevents the successful jet breakout. Especially for Pop III stars, the mass loss is not expected and the large hydrogen envelope is kept due to the low opacity envelope. We find, however, that those Pop III stars who end as blue super giants are compact enough for jets to break out the stellar envelopes successfully. We evaluate observational characters of Pop III GRBs and predict that Pop III GRBs have the duration of $\sim 10^5$ sec in the observer frame and the peak luminosity of $\sim 5 \times 10^{50}{\rm erg/sec}$. Moreover, assuming that the $E_p-L_p$ correlation (or the $E_p-E_{\gamma, \rm iso}$ correlation) holds for Pop III GRBs, we find that the spectrum peak energy falls $\sim$ a few keV (or $\sim 100$ keV) in the observer frame. We discuss the detectability of Pop III GRBs by future satellite missions such as EXIST and Lobster. If the $E_p-E_{\gamma, \rm iso}$ correlation holds for Pop III GRBs, we find that EXIST is more appropriate for GRB detections and that EXIST can detect Pop III GRBs at $z \lesssim 9$. We observe such Pop III GRBs at $z \sim 9$ as long duration X-ray rich GRBs by EXIST. On the other hand, if the $E_p-L_p$ correlation holds, we find that Lobster is more appropriate for detecting GRBs and that Lobster can detect very high z Pop III GRBs up to $z \sim 19$. We observe Pop III GRBs as long duration X-ray flashes by Lobster.

### Three-dimensional magnetohydrodynamic simulations of M-dwarf chromospheres

We present first results from three-dimensional radiation magnetohydrodynamic simulations of M-type dwarf stars with CO5BOLD. The local models include the top of the convection zone, the photosphere, and the chromosphere. The results are illustrated for models with an effective temperature of 3240 K and a gravitational acceleration of log g = 4.5, which represent analogues of AD Leo. The models have different initial magnetic field strengths and field topologies. This first generation of models demonstrates that the atmospheres of M-dwarfs are highly dynamic and intermittent. Magnetic fields and propagating shock waves produce a complicated fine-structure, which is clearly visible in synthetic intensity maps in the core of the Ca II K spectral line and also at millimeter wavelengths. The dynamic small-scale pattern cannot be described by means of one-dimensional models, which has important implications for the construction of semi-empirical model atmospheres and thus for the interpretation of observations in general. Detailed three-dimensional numerical simulations are valuable in this respect. Furthermore, such models facilitate the analysis of small-scale processes, which cannot be observed on stars but nevertheless might be essential for understanding M-dwarf atmospheres and their activity. An example are so-called "magnetic tornadoes", which have recently been found on the Sun and are presented here in M-dwarf models for the first time.

### Modelling non-dust fluids in cosmology

Often in cosmology Newtonian physics is used to model the dynamics of matter inhomogeneities. For pressureless dark matter, or ‘dust’, this approach gives the correct results, but for scenarios in which the fluid has pressure this is no longer the case. In this article, we explicitly highlight the relationship between the variables in Newtonian and cosmological perturbation theory, showing exact equivalence for pressureless matter, and giving the relativistic corrections for matter with pressure. As an example, we focus on the scalar field dark matter model, which has recently gained popularity but which has non-zero pressure perturbations. We discuss some problems which may arise when modelling this theory with numerical simulations, and when using CMB Boltzmann codes.

### Semi-analytic galaxy formation in early dark energy cosmologies

We study the impact of early dark energy (EDE) cosmologies on galaxy properties by coupling high-resolution numerical simulations with semi-analytic modeling (SAM) of galaxy formation and evolution. EDE models are characterized by a non-vanishing high-redshift contribution of dark energy, producing an earlier growth of structures and a modification of large-scale structure evolution. They can be viewed as typical representatives of non-standard dark energy models in which only the expansion history is modified, and hence the impact on galaxy formation is indirect. We show that in EDE cosmologies the predicted space density of galaxies is enhanced at all scales with respect to the standard LCDM scenario, and the corresponding cosmic star formation history and stellar mass density is increased at high-redshift. We compare these results with a set of theoretical predictions obtained with alternative SAMs applied to our reference LCDM simulation, yielding a rough measure of the systematic uncertainty of the models. We find that the modifications in galaxy properties induced by EDE cosmologies are of the same order of magnitude as intra-SAM variations for a standard LCDM realization (unless rather extreme EDE models are considered), suggesting that is difficult to use such predictions alone to disentangle between different cosmological scenarios. However, when independent information on the underlying properties of host dark matter haloes is included, the SAM predictions on galaxy bias may provide important clues on the expansion history and the equation-of-state evolution.

### MPTbreeze: A fast renormalized perturbative scheme

We put forward and test a simple description of multi-point propagators (MP), which serve as building-blocks to calculate the nonlinear matter power spectrum. On large scales these propagators reduce to the well-known kernels in standard perturbation theory, while at smaller scales they are suppresed due to nonlinear couplings. Through extensive testing with numerical simulations we find that this decay is characterized by the same damping scale for both two and three-point propagators. In turn this transition can be well modeled with resummation results that exponentiate one-loop computations. For the first time, we measure the four components of the non-linear (two-point) propagator using dedicated simulations started from two independent random Gaussian fields for positions and velocities, verifying in detail the fundamentals of propagator resummation. We use these results to develop an implementation of the MP-expansion for the nonlinear power spectrum that only requires seconds to evaluate at BAO scales. To test it we construct six suites of large numerical simulations with different cosmologies. From these and LasDamas runs we show that the nonlinear power spectrum can be described at the ~ 2% level at BAO scales for redshifts in the range [0-2.5]. We make a public release of the MPTbreeze code with the hope that it can be useful to the community.

### An improved model for the nonlinear velocity power spectrum

The velocity divergence power spectrum is a key ingredient in modelling redshift space distortion effects on quasi-linear and nonlinear scales. We present an improved model for the z=0 velocity divergence auto and cross power spectrum which was originally suggested by Jennings et al. 2011. Using numerical simulations we measure the velocity fields using a Delaunay tesselation and obtain an accurate prediction of the velocity divergence power spectrum on scales k < 1 hMpc^{-1}. We use this to update the model which is now accurate to 2% for both P_{\theta \theta} and P_{\theta \delta} at z=0 on scales k <0.7 hMpc^{-1} and k <0.5 hMpc^{-1} respectively. We find that the formula for the redshift dependence of the velocity divergence power spectra proposed by Jennings et al. 2011 recovers the measured z>0 P(k) to even greater accuracy with the new model. The nonlinear P_{\theta \theta} and P_{\theta \delta} at z =1 are recovered accurately to better than 2% on scales k<0.2 hMpc^{-1}. Recently it was shown that the velocity field shows larger differences between modified gravity cosmologies and \Lambda CDM compared to the matter field. An accurate model for the velocity divergence power spectrum, such as the one presented here, is a valuable tool for analysing redshift space distortion effects in future galaxy surveys and for constraining deviations from general relativity.

### Magnetic Fields from QCD Phase Transitions

We study the evolution of QCD phase transition-generated magnetic fields in freely decaying MHD turbulence of the expanding Universe. We consider a magnetic field generation model that starts from basic non-perturbative QCD theory and predicts stochastic magnetic fields with an amplitude of the order of 0.02 $\mu$G and small magnetic helicity. We employ direct numerical simulations to model the MHD turbulence decay and identify two different regimes: "weakly helical" turbulence regime, when magnetic helicity increases during decay, and "fully helical" turbulence, when maximal magnetic helicity is reached and an inverse cascade develops. The results of our analysis show that in the most optimistic scenario the magnetic correlation length in the comoving frame can reach 10 kpc with the amplitude of the effective magnetic field being 0.007 nG. We demonstrate that the considered model of magneto-genesis can provide the seed magnetic field for galaxies and clusters.