# Posts Tagged numerical simulation

## Recent Postings from numerical simulation

### The sparkling Universe: a scenario for cosmic void motions

We perform a statistical study of the global motion of cosmic voids using both a numerical simulation and observational data. We analyse their relation to large--scale mass flows and the physical effects that drive those motions. We analyse the bulk motions of voids, defined by the mean velocity of haloes in the surrounding shells in the numerical simulation, and by galaxies in the Sloan Digital Sky Survey Data Release 7. We find void mean bulk velocities close to 400 km/s, comparable to those of haloes (~ 500-600 km/s), depending on void size and the large--scale environment. Statistically, small voids move faster than large ones, and voids in relatively higher density environments have higher bulk velocities than those placed in large underdense regions. Also, we analyze the mean mass density around voids finding, as expected, large--scale overdensities (underdensities) along (opposite to) the void motion direction, suggesting that void motions respond to a pull--push mechanism. This contrasts with massive cluster motions who are mainly governed by the pull of the large-scale overdense regions. Our analysis of void pairwise velocities shows how their relative motions are generated by large--scale density fluctuations. In agreement with linear theory, voids embedded in low (high) density regions mutually recede (attract) each other, providing the general mechanism to understand the bimodal behavior of void motions. In order to compare the theoretical results and the observations we have inferred void motions in the SDSS using linear theory, finding that the estimated observational void motions are in statisticalagreement with the results of the simulation. Regarding large--scale flows, our results suggest a scenario of galaxies and galaxy systems flowing away from void centers with the additional, and morerelevant, contribution of the void bulk motion to the total velocity.

### Estimating SI violation in CMB due to non-circular beam and complex scan in minutes

Mild, unavoidable deviations from circular-symmetry of instrumental beams along with scan strategy can give rise to measurable Statistical Isotropy (SI) violation in Cosmic Microwave Background (CMB) experiments. If not accounted properly, this spurious signal can complicate the extraction of other SI violation signals (if any) in the data. However, estimation of this effect through exact numerical simulation is computationally intensive and time consuming. A generalized analytical formalism not only provides a quick way of estimating this signal, but also gives a detailed understanding connecting the leading beam anisotropy components to a measurable BipoSH characterisation of SI violation. In this paper, we provide an approximate generic analytical method for estimating the SI violation generated due to a non-circular (NC) beam and arbitrary scan strategy, in terms of the Bipolar Spherical Harmonic (BipoSH) spectra. Our analytical method can predict almost all the features introduced by a NC beam in a complex scan and thus reduces the need for extensive numerical simulation worth tens of thousands of CPU hours into minutes long calculations. As an illustrative example, we use WMAP beams and scanning strategy to demonstrate the easability, usability and efficiency of our method. We test all our analytical results against that from exact numerical simulations.

### Light Bridge in a Developing Active Region. II. Numerical Simulation of Flux Emergence and Light Bridge Formation

Light bridges, the bright structure dividing umbrae in sunspot regions, show various activity events. In Paper I, we reported on analysis of multi-wavelength observations of a light bridge in a developing active region (AR) and concluded that the activity events are caused by magnetic reconnection driven by magnetconvective evolution. The aim of this second paper is to investigate the detailed magnetic and velocity structures and the formation mechanism of light bridges. For this purpose, we analyze numerical simulation data from a radiative magnetohydrodynamics model of an emerging AR. We find that a weakly-magnetized plasma upflow in the near-surface layers of the convection zone is entrained between the emerging magnetic bundles that appear as pores at the solar surface. This convective upflow continuously transports horizontal fields to the surface layer and creates a light bridge structure. Due to the magnetic shear between the horizontal fields of the bridge and the vertical fields of the ambient pores, an elongated cusp-shaped current layer is formed above the bridge, which may be favorable for magnetic reconnection. The striking correspondence between the observational results of Paper I and the numerical results of this paper provides a consistent physical picture of light bridges. The dynamic activity phenomena occur as a natural result of the bridge formation and its convective nature, which has much in common with those of umbral dots and penumbral filaments.

### Formation of a condensate during charged collapse [Cross-Listing]

We observe a condensate forming in the interior of a black hole (BH) during numerical simulations of gravitational collapse of a massless charged (complex) scalar field. The magnitude of the scalar field in the interior tends to a non-zero constant; spontaneous breaking of gauge symmetry occurs and a condensate forms. This phenomena occurs in the presence of a BH without the standard symmetry breaking quartic potential; the breaking occurs via the dynamics of the system itself. We also observe that the scalar field in the interior rotates in the complex plane and show that it matches numerically the electric potential to within $1\%$. That a charged scalar condensate can form near the horizon of a black hole in the Abelian Higgs model without the standard symmetry breaking potential had previously been shown analytically in an explicit model involving a massive scalar field in an $AdS_4$ background. Our numerical simulation lends strong support to this finding, although in our case the scalar field is massless and the spacetime is asymptotically flat.

### Formation of a condensate during charged collapse [Replacement]

We observe a condensate forming in the interior of a black hole (BH) during numerical simulations of gravitational collapse of a massless charged (complex) scalar field. The magnitude of the scalar field in the interior tends to a non-zero constant; spontaneous breaking of gauge symmetry occurs and a condensate forms. This phenomena occurs in the presence of a BH without the standard symmetry breaking quartic potential; the breaking occurs via the dynamics of the system itself. We also observe that the scalar field in the interior rotates in the complex plane and show that it matches numerically the electric potential to within $1\%$. That a charged scalar condensate can form near the horizon of a black hole in the Abelian Higgs model without the standard symmetry breaking potential had previously been shown analytically in an explicit model involving a massive scalar field in an $AdS_4$ background. Our numerical simulation lends strong support to this finding, although in our case the scalar field is massless and the spacetime is asymptotically flat.

### Formation of a condensate during charged collapse [Replacement]

We observe a condensate forming in the interior of a black hole (BH) during numerical simulations of gravitational collapse of a massless charged (complex) scalar field. The magnitude of the scalar field in the interior tends to a non-zero constant; spontaneous breaking of gauge symmetry occurs and a condensate forms. This phenomena occurs in the presence of a BH without the standard symmetry breaking quartic potential; the breaking occurs via the dynamics of the system itself. We also observe that the scalar field in the interior rotates in the complex plane and show that it matches numerically the electric potential to within $1\%$. That a charged scalar condensate can form near the horizon of a black hole in the Abelian Higgs model without the standard symmetry breaking potential had previously been shown analytically in an explicit model involving a massive scalar field in an $AdS_4$ background. Our numerical simulation lends strong support to this finding, although in our case the scalar field is massless and the spacetime is asymptotically flat.

### Formation of a condensate during charged collapse

We observe a condensate forming in the interior of a black hole (BH) during numerical simulations of gravitational collapse of a massless charged (complex) scalar field. The magnitude of the scalar field in the interior tends to a non-zero constant; spontaneous breaking of gauge symmetry occurs and a condensate forms. This phenomena occurs in the presence of a BH without the standard symmetry breaking quartic potential; the breaking occurs via the dynamics of the system itself. We also observe that the scalar field in the interior rotates in the complex plane and show that it matches numerically the electric potential to within $1\%$. That a charged scalar condensate can form near the horizon of a black hole in the Abelian Higgs model without the standard symmetry breaking potential had previously been shown analytically in an explicit model involving a massive scalar field in an $AdS_4$ background. Our numerical simulation lends strong support to this finding, although in our case the scalar field is massless and the spacetime is asymptotically flat.

### Spin flips in generic black hole binaries [Cross-Listing]

We study the spin dynamics of individual black holes in a binary system. In particular we focus on the polar precession of spins and the possibility of a complete flip of spins with respect to the orbital plane. We perform a full numerical simulation that displays these characteristics. We evolve equal mass binary spinning black holes for $t=20,000M$ from an initial proper separation of $d=25M$ down to merger after 48.5 orbits. We compute the gravitational radiation from this system and compare it to 3.5 post-Newtonian generated waveforms finding close agreement. We then further use 3.5 post-Newtonian evolutions to show the extension of this spin {flip-flop} phenomenon to unequal mass binaries. We also provide analytic expressions to approximate the maximum {flip-flop} angle and frequency in terms of the binary spins and mass ratio parameters at a given orbital radius. Finally we discuss the effect this spin {flip-flop} would have on accreting matter and other potential observational effects.

### Spin flips in generic black hole binaries

We study the spin dynamics of individual black holes in a binary system. In particular we focus on the polar precession of spins and the possibility of a complete flip of spins with respect to the orbital plane. We perform a full numerical simulation that displays these characteristics. We evolve equal mass binary spinning black holes for $t=20,000M$ from an initial proper separation of $d=25M$ down to merger after 48.5 orbits. We compute the gravitational radiation from this system and compare it to 3.5 post-Newtonian generated waveforms finding close agreement. We then further use 3.5 post-Newtonian evolutions to show the extension of this spin {flip-flop} phenomenon to unequal mass binaries. We also provide analytic expressions to approximate the maximum {flip-flop} angle and frequency in terms of the binary spins and mass ratio parameters at a given orbital radius. Finally we discuss the effect this spin {flip-flop} would have on accreting matter and other potential observational effects.

### Visibility moments and power spectrum of turbulence velocity

Here we introduce moments of visibility function and discuss how those can be used to estimate the power spectrum of the turbulent velocity of external spiral galaxies. We perform numerical simulation to confirm the credibility of this method and found that for galaxies with lower inclination angles it works fine. This is the only method to estimate the power spectrum of the turbulent velocity fluctuation in the ISM of the external galaxies.

### Visibility moments and power spectrum of turbulence velocity [Replacement]

Here we introduce moments of visibility function and discuss how those can be used to estimate the power spectrum of the turbulent velocity of external spiral galaxies. We perform numerical simulation to confirm the credibility of this method and found that for galaxies with lower inclination angles it works fine. This is the only method to estimate the power spectrum of the turbulent velocity fluctuation in the ISM of the external galaxies.

### Solar wind turbulence from MHD to sub-ion scales: high-resolution hybrid simulations [Replacement]

We present results from a high-resolution and large-scale hybrid (fluid electrons and particle-in-cell protons) two-dimensional numerical simulation of decaying turbulence. Two distinct spectral regions (separated by a smooth break at proton scales) develop with clear power-law scaling, each one occupying about a decade in wave numbers. The simulation results exhibit simultaneously several properties of the observed solar wind fluctuations: spectral indices of the magnetic, kinetic, and residual energy spectra in the magneto-hydrodynamic (MHD) inertial range along with a flattening of the electric field spectrum, an increase in magnetic compressibility, and a strong coupling of the cascade with the density and the parallel component of the magnetic fluctuations at sub-proton scales. Our findings support the interpretation that in the solar wind large-scale MHD fluctuations naturally evolve beyond proton scales into a turbulent regime that is governed by the generalized Ohm's law.

### Solar wind turbulence from MHD to sub-ion scales: high-resolution hybrid simulations

We present results from a high-resolution and large-scale hybrid (fluid electrons and particle-in-cell protons) two-dimensional numerical simulation of decaying turbulence. Two distinct spectral regions (separated by a smooth break at proton scales) develop with clear power-law scaling, each one occupying about a decade in wave numbers. The simulation results exhibit simultaneously several properties of the observed solar wind fluctuations: spectral indices of the magnetic, kinetic, and residual energy spectra in the magneto-hydrodynamic (MHD) inertial range along with a flattening of the electric field spectrum, an increase in magnetic compressibility, and a strong coupling of the cascade with the density and the parallel component of the magnetic fluctuations at sub-proton scales. Our findings support the interpretation that in the solar wind large-scale MHD fluctuations naturally evolve beyond proton scales into a turbulent regime that is governed by the generalized Ohm's law.

### Vacuum high harmonic generation in the shock regime [Cross-Listing]

Electrodynamics becomes nonlinear and permits the self-interaction of fields when the quantised nature of vacuum states is taken into account. The effect on a plane probe pulse propagating through a stronger constant crossed background is calculated using numerical simulation and by analytically solving the corresponding wave equation. The electromagnetic shock resulting from vacuum high harmonic generation is investigated and a nonlinear shock parameter identified.

### Transport by meridional circulations in solar-type stars

Transport by meridional flows has significant consequences for stellar evolution, but is difficult to capture in global-scale numerical simulations because of the wide range of timescales involved. Stellar evolution models therefore usually adopt parameterizations for such transport based on idealized laminar or mean-field models. Unfortunately, recent attempts to model this transport in global simulations have produced results that are not consistent with any of these idealized models. In an effort to explain the discrepancies between global simulations and idealized models, we here use three-dimensional local Cartesian simulations of compressible convection to study the efficiency of transport by meridional flows below a convection zone in several parameter regimes of relevance to the Sun and solar-type stars. In these local simulations we are able to establish the correct ordering of dynamical timescales, although the separation of the timescales remains unrealistic. We find that, even though the generation of internal waves by convective overshoot produces a high degree of time dependence in the meridional flow field, the mean flow has the qualitative behavior predicted by laminar, "balanced" models. In particular, we observe a progressive deepening, or "burrowing", of the mean circulation if the local Eddington-Sweet timescale is shorter than the viscous diffusion timescale. Such burrowing is a robust prediction of laminar models in this parameter regime, but has never been observed in any previous numerical simulation. We argue that previous simulations therefore underestimate the transport by meridional flows.

### On the correction of conserved variables for numerical RMHD with staggered constrained transport

Despite the success of the combination of conservative schemes and staggered constrained transport algorithms in the last fifteen years, the accurate description of highly magnetized, relativistic flows with strong shocks represents still a challenge in numerical RMHD. The present paper focusses in the accuracy and robustness of several correction algorithms for the conserved variables, which has become a crucial ingredient in the numerical simulation of problems where the magnetic pressure dominates over the thermal pressure by more than two orders of magnitude. Two versions of non-relativistic and fully relativistic corrections have been tested and compared using a magnetized cylindrical explosion with high magnetization ($\ge 10^4$) as test. In the non-relativistic corrections, the total energy is corrected for the difference in the classical magnetic energy term between the average of the staggered fields and the conservative ones, before (CA1) and after (CA1') recovering the primitive variables. These corrections are unable to pass the test at any numerical resolution. The two relativistic approaches (CA2 and CA2'), correcting also the magnetic terms depending on the flow speed in both the momentum and the total energy, reveal as much more robust. These algorithms pass the test succesfully and with very small deviations of the energy conservation ($\le 10^{-4}$), and very low values of the total momentum ($\le 10^{-8}$). In particular, the algorithm CA2' (that corrects the conserved variables after recovering the primitive variables) passes the test at all resolutions. The numerical code used to run all the test cases is briefly described.

### Effects of Turbulent Viscosity on A Rotating Gas Ring Around A Black Hole: The Density Profile of Numerical Simulation

In this paper, we present the time evolution of a rotationally axisymmetric gas ring around a non rotating black hole using two dimensional grid-based hydrodynamic simulation. We show the way in which angular momentum transport is included in simulations of non-self-gravitating accretion of matter towards a black hole. We use the Shakura-Sunyaev {\alpha} viscosity prescription to estimate the turbulent viscosity. We investigate how a gas ring which is initially assumed to rotate with Keplerian angular velocity is accreted on to a back hole and hence forms accretion disc in the presence of turbulent viscosity. Furthermore, we also show that increase of the {\alpha} coefficient increases the rate of advection of matter towards the black hole. The density profile we obtain is in good quantitative agreement with that obtained from the analytical results. The dynamics of resulting angular momentum depends strongly on {\alpha}.

### Effects of Turbulent Viscosity on A Rotating Gas Ring Around A Black Hole: Results in Numerical Simulation [Replacement]

In this paper, we present the time evolution of a rotationally axisymmetric gas ring around a non rotating black hole using two dimensional grid-based hydrodynamic simulation. We show the way in which angular momentum transport is included in simulations of non-self-gravitating accretion of matter towards a black hole. We use the Shakura-Sunyaev {\alpha} viscosity prescription to estimate the turbulent viscosity for all major viscous stress tensors. We investigate how a gas ring which is initially assumed to rotate with Keplerian angular velocity is accreted on to a black hole and hence forms accretion disc in the presence of turbulent viscosity. We show that the centrifugal pressure supported sub-Keplerian flow with shocks forms when the ring starts to disperse with inclusion of relatively small amount of viscosity. But, if the viscosity is above the critical value, the shock disappears altogether and the whole disc becomes Kepleiran which is subsonic everywhere except in a region close to the horizon, where it supersonically enters to the black hole. We discovered a multiple valued Mach number solution and the corresponding density distributions that connects matter (a) from the initial Keplerian gas ring to a sub-Keplerian disc with shocks in presence of small amount of viscosity and (b) from the sub-Keplerian flow to a Keplerian disc in presence of huge amount of viscosity. We calculate the temporal variations of the magnitude of various time scales which ensure us about the stability of the flow.

### Numerical Simulation of Hot Accretion Flows (III): Revisiting wind properties using trajectory approach

Previous MHD simulations have shown that wind must exist in black hole hot accretion flows. In this paper, we continue our study by investigating the detailed properties of wind, such as mass flux and poloidal speed, and the mechanism of wind production. For this aim, we make use of a three dimensional GRMHD simulation of hot accretion flows around a Schwarzschild black hole. The simulation is designed so that the magnetic flux is not accumulated significantly around the black hole. To distinguish real wind from turbulent outflows, we track the trajectories of the virtual Largrangian particles from simulation data. We find two types of real outflows, i.e., a quasi-relativistic jet close to the axis and a sub-relativistic wind subtending a much larger solid angle. Most of the wind originates from the surface layer of the accretion flow. The poloidal wind speed almost remains constant once they are produced, but the flux-weighted wind speed roughly follows $v_{\rm p, wind}(r)\approx 0.25 v_k(r)$. The mass flux of jet is much lower but the speed is much higher, $v_{\rm p,jet}\sim (0.3-0.4) c$. Consequently, both the energy and momentum fluxes of the wind are much larger than those of the jet. We find that the wind is produced and accelerated primarily by the combination of centrifugal force and magnetic pressure gradient, while the jet is mainly accelerated by magnetic pressure gradient. Finally, we find that the wind production efficiency $\epsilon_{\rm wind}\equiv\dot{E}_{\rm wind}/\dot{M}_{\rm BH}c^2\sim 1/1000$, in good agreement with the value required from large-scale galaxy simulations with AGN feedback.

### Hemispheric Coupling: Comparing Dynamo Simulations and Observations

Numerical simulations that reproduce solar-like magnetic cycles can be used to generate long-term statistics. The variations in N-S hemispheric cycle synchronicity and amplitude produced in simulations has not been widely compared to observations. The observed limits on asymmetry show that hemispheric sunspot area production is no more than 20% asymmetric for cycles 12-23 and phase lags do not exceed 20% (2 yrs) of the total cycle period. Independent studies have found a long-term trend in phase values as one hemisphere leads the other for ~four cycles. Such persistence in phase is not indicative of a stochastic phenomenon. We compare the findings to results from a numerical simulation of solar convection recently produced with the EULAG-MHD model. This simulation spans 1600 yrs and generated 40 regular, sunspot-like cycles. While the simulated cycle length is too long and the toroidal bands remain at too high of latitudes, some solar-like aspects of hemispheric asymmetry are reproduced. The model reproduces the synchrony of polarity inversions and onset of cycle as the simulated phase lags do not exceed 20% of the cycle period. Simulated amplitude variations between the N and S hemispheres are larger than observed in the Sun. The simulations show one hemisphere persistently leads the other for several successive cycles, placing an upper bound on the efficiency of transequatorial magnetic coupling mechanisms. These include magnetic diffusion, cross-equatorial mixing within elongated convective rolls and transequatorial meridional flow cells. One or more of these processes may lead to magnetic flux cancellation whereby the oppositely directed fields come in close proximity and cancel each other across the magnetic equator late in the solar cycle. We discuss the discrepancies between model and observations and the constraints they pose on possible mechanisms of hemispheric coupling.

### Collective coordinate approximation to the scattering of solitons in modified NLS and sine-Gordon models

We use a collective coordinate approximation to model the scattering of two solitons in modified nonlinear Schr\"odinger and sine-Gordon systems. We find that the anomalies of the conservation laws of the charges as calculated using the collective coordinate approximation demonstrate the same dependence on the symmetry of the field configuration as that previously found analytically and using a full numerical simulation. This suggests that the collective coordinate approximation is a suitable method to investigate quasi-integrability in perturbed integrable models. We also discuss the general accuracy of this approximation by comparing our results with those of the full numerical simulations and find that the approximation is often remarkably accurate though less so when the models are a long way from the integrable case.

### The 3D MHD code GOEMHD3 for large-Reynolds-number astrophysical plasmas [Replacement]

The numerical simulation of turbulence and flows in almost ideal, large-Reynolds-number astrophysical plasmas motivates the implementation of almost conservative MHD computer codes. They should efficiently calculate, use highly parallelized schemes scaling well with large numbers of CPU cores, allows to obtain a high grid resolution over large simulation domains and which can easily be adapted to new computer architectures as well as to new initial and boundary conditions, allow modular extensions. The new massively parallel simulation code GOEMHD3 enables efficient and fast simulations of almost ideal, large-Reynolds-number astrophysical plasma flows, well resolved and on huge grids covering large domains. Its abilities are validated by major tests of ideal and weakly dissipative plasma phenomena. The high resolution ($2048^3$ grid points) simulation of a large part of the solar corona above an observed active region proved the excellent parallel scalability of the code using more than 30.000 processor cores.

### The 3D MHD code GOEMHD3 for large-Reynolds-number astrophysical plasmas

The numerical simulation of turbulence and flows in almost ideal, large-Reynolds-number astrophysical plasmas motivates the implementation of almost conservative MHD computer codes. They should efficiently calculate, use highly parallelized schemes scaling well with large numbers of CPU cores, allows to obtain a high grid resolution over large simulation domains and which can easily be adapted to new computer architectures as well as to new initial and boundary conditions, allow modular extensions. The new massively parallel simulation code GOEMHD3 enables efficient and fast simulations of almost ideal, large-Reynolds-number astrophysical plasma flows, well resolved and on huge grids covering large domains. Its abilities are validated by major tests of ideal and weakly dissipative plasma phenomena. The high resolution ($2048^3$ grid points) simulation of a large part of the solar corona above an observed active region proved the excellent parallel scalability of the code using more than 30.000 processor cores.

### Determination of the mass anomalous dimension for $N_f=12$ and $N_f=9$ SU($3$) gauge theories [Cross-Listing]

We show the numerical simulation result for the mass anomalous dimension of the SU($3$) gauge theory coupled to $N_f = 12$ fundamental fermions. We use two independent methods, namely the step scaling method and the hyperscaling method of the Dirac mode number, to determine the anomalous dimension in the vicinity of the infrared fixed point of the theory. We show the continuum extrapolations keeping the renormalized coupling constant as a reference in both analyses. Furthermore, some recent works seems to suggest the lower boundary of the conformal window of the SU($3$) gauge theory exists between $N_f=8$ and $10$. We also briefly report our new project, in which the numerical simulation of the SU($3$) gauge theory coupled to $N_f=9$ fundamental fermions has been performed.

### Three Einstein rings: explicit solution and numerical simulation

We investigated the effects of gravitational lensing for a system in which a lens is a point mass and a homogeneous disc with a central hole. In such system there is a variety of cases resulting in formation of one, two and three Einstein rings. We found an explicit solution and considered conditions for existence of the second Einstein ring arising on the disc. Numerical modelling of the images was made for various ratios of the central mass to the disc one and for various values of the disc surface density. We also analysed dependence of the magnification factor on a source position for such system. The result of our work can be used in search of astrophysical objects with a toroidal (ring) structure.

### Flip-flopping binary black holes

We perform a full numerical simulation of binary spinning black holes to display the long term spin dynamics. We start the simulation at an initial proper separation between holes of d~25M and evolve them down to merger for nearly 48 orbits, 3 precession cycles and half of a flip-flop cycle. The simulation lasts for t=20000M and displays a change in the orientation of the spin of the secondary black hole from initially aligned with the orbital angular momentum to a complete anti-alignment after half of a flip-flop cycle. This process continuously flip-flops the spin during the lifetime of the binary. We discuss the consequences of this oscillation mode for accreting binaries, in particular for the spin growth and binary dynamics as well as the observational consequences for galactic and supermassive black holes.

### Flip-flopping binary black holes [Cross-Listing]

We perform a full numerical simulation of binary spinning black holes to display the long term spin dynamics. We start the simulation at an initial proper separation between holes of d~25M and evolve them down to merger for nearly 48 orbits, 3 precession cycles and half of a flip-flop cycle. The simulation lasts for t=20000M and displays a change in the orientation of the spin of the secondary black hole from initially aligned with the orbital angular momentum to a complete anti-alignment after half of a flip-flop cycle. This process continuously flip-flops the spin during the lifetime of the binary. We discuss the consequences of this oscillation mode for accreting binaries, in particular for the spin growth and binary dynamics as well as the observational consequences for galactic and supermassive black holes.

### Flip-flopping binary black holes [Replacement]

We study binary spinning black holes to display the long term individual spin dynamics. We perform a full numerical simulation starting at an initial proper separation of $d\approx25M$ between equal mass holes and evolve them down to merger for nearly 48 orbits, 3 precession cycles, and half of a flip-flop cycle. The simulation lasts for $t=20000M$ and displays a total change in the orientation of the spin of one of the black holes from initially aligned with the orbital angular momentum to a complete anti-alignment after half of a flip-flop cycle. We compare this evolution with an integration of the 3.5 Post-Newtonian equations of motion and spin evolution to show that this process continuously flip-flops the spin during the lifetime of the binary until merger. We also provide lower order analytic expressions for the maximum flip-flop angle and frequency. We discuss the effects this dynamics may have on spin growth in accreting binaries and on the observational consequences for galactic and supermassive binary black holes.

### Flip-flopping binary black holes [Replacement]

We study binary spinning black holes to display the long term individual spin dynamics. We perform a full numerical simulation starting at an initial proper separation of $d\approx25M$ between equal mass holes and evolve them down to merger for nearly 48 orbits, 3 precession cycles, and half of a flip-flop cycle. The simulation lasts for $t=20000M$ and displays a total change in the orientation of the spin of one of the black holes from initially aligned with the orbital angular momentum to a complete anti-alignment after half of a flip-flop cycle. We compare this evolution with an integration of the 3.5 Post-Newtonian equations of motion and spin evolution to show that this process continuously flip-flops the spin during the lifetime of the binary until merger. We also provide lower order analytic expressions for the maximum flip-flop angle and frequency. We discuss the effects this dynamics may have on spin growth in accreting binaries and on the observational consequences for galactic and supermassive binary black holes.

### Analysis of the spiral structure in a simulated galaxy

We analyze the spiral structure that results in a numerical simulation of a galactic disk with stellar and gaseous components evolving in a potential that includes an axisymmetric halo and bulge. We perform a second simulation without the gas component to observe how it affects the spiral structure in the disk. To quantify this, we use a Fourier analysis and obtain values for the pitch angle and the velocity of the self-excited spiral pattern of the disk. The results show a tighter spiral in the simulation with gaseous component. The spiral structure is consistent with a superposition of waves, each with a constant pattern velocity in given radial ranges.

### Numerical Simulation of Radio Signal from Extended Air Showers

The burst of radio emission by the extensive air shower provides a promising alternative for detecting ultra-high energy cosmic rays.We have developed an independent numerical program to simulate these radio signals. Our code is based on a microscopic treatment, with both the geosynchrotron radiation and charge excess effect included. Here we make a first presentation of our basic program and its results. The time signal for different polarizations are computed, we find that the pulses take on a bipolar pattern, the spectrum is suppressed towards the lower frequencies.We investigate how the shower at different heights in atmosphere contribute to the total signal, and examine the signal strength and distribution at sites of different elevations. We also study the signal from showers of different inclination angles and azimuth directions. In all these cases we find the charge excess effect important.

### Towards Collisions of Inhomogeneous Shockwaves in AdS

We perform a numerical simulation of the evolution of inhomogeneities with transverse profile in a collision of gravitational shockwaves in asymptotically anti-de Sitter spacetime. This constitutes a step closer towards an accurate holographic description of the thermalization of a strongly coupled plasma, which can model the dynamics of heavy ion collisions. The results indicate that previous calculations of the thermalization time may have been underestimates.

### Towards Collisions of Inhomogeneous Shockwaves in AdS [Replacement]

We perform a numerical simulation of the evolution of inhomogeneities with transverse profile in a collision of gravitational shockwaves in asymptotically anti-de Sitter spacetime. This constitutes a step closer towards an accurate holographic description of the thermalization of a strongly coupled plasma, which can model the dynamics of heavy ion collisions. The results indicate that the considered inhomogeneities typically become hydrodynamical earlier or at the same moment when hydrodynamics applies to the background, even though they decay slowly.

### Towards Collisions of Inhomogeneous Shockwaves in AdS [Cross-Listing]

We perform a numerical simulation of the evolution of inhomogeneities with transverse profile in a collision of gravitational shockwaves in asymptotically anti-de Sitter spacetime. This constitutes a step closer towards an accurate holographic description of the thermalization of a strongly coupled plasma, which can model the dynamics of heavy ion collisions. The results indicate that previous calculations of the thermalization time may have been underestimates.

### Towards Collisions of Inhomogeneous Shockwaves in AdS [Replacement]

We perform a numerical simulation of the evolution of inhomogeneities with transverse profile in a collision of gravitational shockwaves in asymptotically anti-de Sitter spacetime. This constitutes a step closer towards an accurate holographic description of the thermalization of a strongly coupled plasma, which can model the dynamics of heavy ion collisions. The results indicate that the considered inhomogeneities typically become hydrodynamical earlier or at the same moment when hydrodynamics applies to the background, even though they decay slowly.

### Towards Collisions of Inhomogeneous Shockwaves in AdS

We perform a numerical simulation of the evolution of inhomogeneities with transverse profile in a collision of gravitational shockwaves in asymptotically anti-de Sitter spacetime. This constitutes a step closer towards an accurate holographic description of the thermalization of a strongly coupled plasma, which can model the dynamics of heavy ion collisions. The results indicate that previous calculations of the thermalization time may have been underestimates.

### Towards Collisions of Inhomogeneous Shockwaves in AdS [Replacement]

We perform a numerical simulation of the evolution of inhomogeneities with transverse profile in a collision of gravitational shockwaves in asymptotically anti-de Sitter spacetime. This constitutes a step closer towards an accurate holographic description of the thermalization of a strongly coupled plasma, which can model the dynamics of heavy ion collisions. The results indicate that the considered inhomogeneities typically become hydrodynamical earlier or at the same moment when hydrodynamics applies to the background, even though they decay slowly.

### Numerical Simulation of Superhalo Electrons Generated by Magnetic Reconnection in the Solar Wind Source Region

Superhalo electrons appear to be continuously present in the interplanetary medium, even at very quiet times, with a power-law spectrum at energies above $\sim$2 keV. Here we numerically investigate the generation of superhalo electrons by magnetic reconnection in the solar wind source region, using the MHD and test particle simulations for both single X-line reconnection and multiple X-line reconnection. We find that the direct current electric field, produced in the magnetic reconnection region, can accelerate electrons from an initial thermal energy of T $\sim10^5$ K up to hundreds of keV. After acceleration, some of the accelerated electrons, together with the nascent solar wind flow driven by the reconnection, propagate upwards along the newly-opened magnetic field lines into the interplanetary space, while the rest move downwards into the lower atmosphere. Similar to the observed superhalo electrons at 1 AU, the flux of the upward-traveling accelerated electrons versus energy displays a power-law distribution at $\sim$ 2 $-$ 100 keV, $f(E) \sim E^{-\delta}$, with a $\delta$ of $\sim$ 1.5 $-$ 2.4. For single (multiple) X-line reconnection, the spectrum becomes harder (softer) as the anomalous resistivity parameter $\alpha$ (uniform resistivity $\eta$) increases. These modeling results suggest that the acceleration in the solar wind source region may contribute to superhalo electrons.

### Numerical Simulation Of Spectral And Timing Properties Of Galactic Black Holes

A black hole accretion may have both the Keplerian and the sub-Keplerian components. We consider the most general accretion flow configuration, namely, two-component advective flow (TCAF) in which the Keplerian disk is immersed inside a low angular momentum, accreting sub-Keplerian halo component around a black hole. Low energy (soft) photons from the Keplerian component and hot electrons in the sub-Keplerian component exchange their energy through Comptonization or inverse-Comptonization processes. In the sub-Keplerian component, a shock is generally formed due to the centrifugal force. The post-shock region is known as the CENtrifugal pressure dominated BOundary Layer (CENBOL). The spectral and the timing properties of TCAF have been extensively studied using mostly analytical and some time dependent numerical simulations since the model was proposed by Chakrabarti & Titarchuk in 1995. The findings are the key inputs of understanding several observed features of black hole candidates. In this thesis, using numerical simulation, we rigorously prove some of the conjectures of the TCAF model. In the work presented in this thesis, we have considered for the first time the presence of both the Keplerian and the sub-Keplerian flow in a single simulation. The Keplerian disk resides on the equatorial plane and is the standard disk from which low energy photons having multi-color blackbody spectrum is emitted. The hydrodynamics as well as the thermal properties of the sub-Keplerian halo are simulated using a finite difference code which uses the principle of total variation diminishing (TVD). The Comptonization between the photons and the hot electrons is simulated using a Monte Carlo code. These two codes are then coupled and the resulting localized heating and cooling are included in the coupled code. Using this code, we study the spectral and timing properties of the TCAF.

### Rotating Solar Jets in Simulations of Flux Emergence with Thermal Conduction

We study the formation of coronal jets through numerical simulation of the emergence of a twisted magnetic flux rope into a pre-existing open magnetic field. Reconnection inside the emerging flux rope in addition to that between the emerging and pre-existing fields give rise to the violent eruption studied. The simulated event closely resembles the coronal jets ubiquitously observed by Hinode/XRT and demonstrates that heated plasma is driven into the extended atmosphere above. Thermal conduction implemented in the model allows us to qualitatively compare simulated and observed emission from such events. We find that untwisting field lines after the reconnection drive spinning outflows of plasma in the jet column. The Poynting flux in the simulated jet is dominated by the untwisting motions of the magnetic fields loaded with high-density plasma. The simulated jet is comprised of spires of untwisting field that are loaded with a mixture of cold and hot plasma and exhibit rotational motion of order 20 km/s and match contemporary observations.

### An Analysis of Chandra Deep Follow-up GRBs: Implications for Off-Axis Jets [Replacement]

We present a sample of 27 GRBs with detailed Swift light curves supplemented by late time Chandra observations. To answer the missing jet-break problem in general, we develop a numerical simulation based model which can be directly fit to the data using Monte Carlo methods. Our numerical model takes into account all the factors that can shape a jet-break: (i) lateral expansion (ii) edge effects and (iii) off-axis effects. Our results provide improved fits to the light curves and constraints on physical parameters. More importantly, our results suggest that off-axis effects are important and must be included in interpretations of GRB jet breaks.

### Restricted phase-space approximation in real-time stochastic quantization

We perform and extend real-time numerical simulation of a scalar quantum field theory using stochastic quantization. After a brief review of the quantization method, we calculate the propagator and the perturbative series and compare with analytical results. This is a first step toward general applications, and we focus only on the vacuum properties of the theory; this enables us to handle the boundary condition by the $i\epsilon$ prescription. Then, we explicitly check the convergence and solve the differential equation in frequency space. For clarity we drop the spatial-derivative terms and make a comparison between our results and the numerically exact results obtained by diagonalization of the Hamiltonian. While we can control stability of the numerical simulation for any coupling strength, our results turn out to flow into an unphysical attractor if the simulation is out of the weak-coupling regime. We propose a simple truncation scheme to incorporate the interaction terms, which we name the "restricted phase-space approximation." With this method, we obtain results with stable simulation at good accuracy. Finally we give a short discussion on the closed-time path formalism.

### Restricted phase-space approximation in real-time stochastic quantization [Cross-Listing]

We perform and extend real-time numerical simulation of a scalar quantum field theory using stochastic quantization. After a brief review of the quantization method, we calculate the propagator and the perturbative series and compare with analytical results. This is a first step toward general applications, and we focus only on the vacuum properties of the theory; this enables us to handle the boundary condition by the $i\epsilon$ prescription. Then, we explicitly check the convergence and solve the differential equation in frequency space. For clarity we drop the spatial-derivative terms and make a comparison between our results and the numerically exact results obtained by diagonalization of the Hamiltonian. While we can control stability of the numerical simulation for any coupling strength, our results turn out to flow into an unphysical attractor if the simulation is out of the weak-coupling regime. We propose a simple truncation scheme to incorporate the interaction terms, which we name the "restricted phase-space approximation." With this method, we obtain results with stable simulation at good accuracy. Finally we give a short discussion on the closed-time path formalism.

### Restricted phase-space approximation in real-time stochastic quantization [Replacement]

We perform and extend real-time numerical simulation of a low-dimensional scalar field theory or a quantum mechanical system using stochastic quantization. After a brief review of the quantization method and the complex Langevin dynamics, we calculate the propagator and make a comparison with analytical results. This is a first step toward general applications, and we focus only on the vacuum properties of the theory; this enables us to handle the boundary condition with the $i\epsilon$ prescription in frequency space. While we can control stability of the numerical simulation for any coupling strength, our results turn out to flow into an unphysical fixed-point, which is qualitatively understood from the corresponding Fokker-Planck equation. We propose a simple truncation scheme, "restricted phase-space approximation," to avoid the unphysical fixed-point. With this method, we obtain stable results at reasonably good accuracy. Finally we give a short discussion on the closed-time path formalism and demonstrate the direct computation of the vacuum expectation value not with the $i\epsilon$ prescription but from an explicit construction of the Feynman kernel.

### Restricted phase-space approximation in real-time stochastic quantization [Replacement]

We perform and extend real-time numerical simulation of a low-dimensional scalar field theory or a quantum mechanical system using stochastic quantization. After a brief review of the quantization method and the complex Langevin dynamics, we calculate the propagator and make a comparison with analytical results. This is a first step toward general applications, and we focus only on the vacuum properties of the theory; this enables us to handle the boundary condition with the $i\epsilon$ prescription in frequency space. While we can control stability of the numerical simulation for any coupling strength, our results turn out to flow into an unphysical fixed-point, which is qualitatively understood from the corresponding Fokker-Planck equation. We propose a simple truncation scheme, "restricted phase-space approximation," to avoid the unphysical fixed-point. With this method, we obtain stable results at reasonably good accuracy. Finally we give a short discussion on the closed-time path formalism and demonstrate the direct computation of the vacuum expectation value not with the $i\epsilon$ prescription but from an explicit construction of the Feynman kernel.

### Spatially Extended 21 cm Signal from Strongly Clustered UV and X-Ray Sources in the Early Universe [Replacement]

We present our prediction for the local 21 cm differential brightness temperature ($\delta T_{b}$) from a set of strongly clustered sources of Population III (Pop III) and II (Pop II) objects in the early Universe, by a numerical simulation of their formation and radiative feedback. These objects are located inside a highly biased environment, which is a rare, high-density peak ("Rarepeak") extending to $\sim7$ comoving Mpc. We study the impact of ultraviolet (UV) and X-ray photons on the intergalactic medium (IGM) and the resulting $\delta T_{b}$, when Pop III stars are assumed to emit X-ray photons by forming X-ray binaries very efficiently. We parameterize the rest-frame spectral energy distribution (SED) of X-ray photons, which regulates X-ray photon-trapping, IGM-heating, secondary Lyman-alpha pumping and the resulting morphology of $\delta T_{b}$. A combination of emission ($\delta T_{b}>0$) and absorption ($\delta T_{b}<0$) regions appears in varying amplitudes and angular scales. The boost of the signal by the high-density environment ($\delta\sim0.64$) and on a relatively large scale combine to make Rarepeak a discernible, spatially-extended ($\theta\sim10'$) object for 21 cm observation at $13\lesssim z\lesssim17$, which is found to be detectable as a single object by SKA with integration time of $\sim1000$ hours. Power spectrum analysis by some of the SKA precursors (LOFAR, MWA, PAPER) of such rare peaks is found difficult due to the rarity of these peaks, and the contribution only by these rare peaks to the total power spectrum remains subdominant compared to that by all astrophysical sources.

### Spinning black hole in the puncture method: Numerical experiments

The strong-field region inside a black hole needs special attention during numerical simulation. One approach for handling the problem is the moving puncture method, which has become an important tool in numerical relativity since it allows long term simulations of binary black holes. An essential component of this method is the choice of the '1+log'-slicing condition. We present an investigation of this slicing condition in rotating black hole spacetimes. We discuss how the results of the stationary Schwarzschild '1+log'-trumpet change when spin is added. This modification enables a simple and cheap algorithm for determining the spin of a non-moving black hole for this particular slicing condition. Applicability of the algorithm is verified in simulations of single black hole, binary neutron star and mixed binary simulations.

### Formation of a Flare-Productive Active Region: Observation and Numerical Simulation of NOAA AR 11158

We present a comparison of the Solar Dynamics Observatory (SDO) analysis of NOAA Active Region (AR) 11158 and numerical simulations of flux-tube emergence, aiming to investigate the formation process of this flare-productive AR. First, we use SDO/Helioseismic and Magnetic Imager (HMI) magnetograms to investigate the photospheric evolution and Atmospheric Imaging Assembly (AIA) data to analyze the relevant coronal structures. Key features of this quadrupolar region are a long sheared polarity inversion line (PIL) in the central delta-sunspots and a coronal arcade above the PIL. We find that these features are responsible for the production of intense flares, including an X2.2-class event. Based on the observations, we then propose two possible models for the creation of AR 11158 and conduct flux-emergence simulations of the two cases to reproduce this AR. Case 1 is the emergence of a single flux tube, which is split into two in the convection zone and emerges at two locations, while Case 2 is the emergence of two isolated but neighboring tubes. We find that, in Case 1, a sheared PIL and a coronal arcade are created in the middle of the region, which agrees with the AR 11158 observation. However, Case 2 never builds a clear PIL, which deviates from the observation. Therefore, we conclude that the flare-productive AR 11158 is, between the two cases, more likely to be created from a single split emerging flux than from two independent flux bundles.

### Formation of subhorizon black holes from preheating [Replacement]

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.