Recent advances in numerical algorithms and computational power have enabled first-principles simulations of pulsar magnetospheres using Particle-in-Cell (PIC) techniques. These ab-initio simulations seem to indicate that pair creation through photon-photon collision at the light cylinder is required to sustain the pulsar engine. However, for many rotation-powered pulsars pair creation operates effectively only near the stellar surface where magnetic field is high. How these "weak pulsars" fill their magnetospheres without efficient photon-photon pair conversion in the outer magnetosphere is still an open question. In this paper, we present a range of self-consistent solutions to the pulsar magnetosphere that do not require pair production near the light cylinder. When pair production is very efficient near the star, the pulsar magnetosphere converges to previously-reported solutions. However, in the intermediate regime where pair supply is barely enough to sustain the magnetospheric current, we observe a time-dependent solution with quasi-period about half of the rotation period. This new quasi-periodic solution may explain the observed pulsar death line without invoking multipolar components near the star, and can potentially explain the core vs. conal emission patterns observed in pulsar radio signals.

In this paper, we present the classical Bondi accretion theory for the case of non-isothermal accretion processes onto a supermassive black hole (SMBH), including the effects of X-ray heating and the radiation force due to electron scattering and spectral lines. The radiation field is calculated by considering an optically thick, geometrically thin, standard accretion disk as the emitter of UV photons and a spherical central object as a source of X-ray emission. In the present analysis, the UV emission from the accretion disk is assumed to have an angular dependence, while the X-ray/central object radiation is assumed to be isotropic. This allows us to build streamlines in any angular direction we need to. The influence of both types of radiation is evaluated for different flux fractions of the X-ray and UV emissions with and without the effects of spectral line driving. We find that the radiation emitted near the SMBH interacts with the infalling matter and modifies the accretion dynamics. In the presence of line driving, a transition resembles from pure type 1 & 2 to type 5 solutions (see Fig2.1 of Frank etal. 2002), which takes place regardless of whether or not the UV emission dominates over the X-ray emission. We compute the radiative factors at which this transition occurs, and discard type 5 solution from all our models. Estimated values of the accretion radius and accretion rate in terms of the classical Bondi values are also given. The results are useful for the construction of proper initial conditions for time-dependent hydrodynamical simulations of accretion flows onto SMBH at the centre of galaxies.

We present three-dimensional calculations of spherically symmetric Bondi accretion onto a stationary supermassive black hole (SMBH) of mass $10^{8}$ $M_{\odot}$ within a radial range of $0.02-10$ pc, using a modified version of the smoothed particle hydrodynamics (SPH) \gad \sp code, which ensures approximate first-order consistency (i.e., second-order accuracy) for the particle approximation. First-order consistency is restored by allowing the number of neighbours, $n_{\rm neigh}$, and the smoothing length, $h$, to vary with the total number of particles, $N$, such that the asymptotic limits $n_{\rm neigh}\to\infty$ and $h\to 0$ hold as $N\to\infty$. The ability of the method to reproduce the isothermal ($\gamma =1$) and adiabatic ($\gamma =5/3$) Bondi accretion is investigated with increased spatial resolution. In particular, for the isothermal models the numerical radial profiles closely match the Bondi solution, except near the accretor, where the density and radial velocity are slightly underestimated. However, as $n_{\rm neigh}$ is increased and $h$ is decreased, the calculations approach first-order consistency and the deviations from the Bondi solution decrease. The density and radial velocity profiles for the adiabatic models are qualitatively similar to those for the isothermal Bondi accretion. Steady-state Bondi accretion is reproduced by the highly resolved consistent models with a percent relative error of $\lesssim 1$\% for $\gamma =1$ and $\sim 9$\% for $\gamma =5/3$, with the adiabatic accretion taking longer than the isothermal case to reach steady flow. The performance of the method is assessed by comparing the results with those obtained using the standard Gadget and the Gizmo codes.

Recent observations of young stellar systems with the Atacama Large Millimeter/submillimeter Array (ALMA) and the Karl G. Jansky Very Large Array (VLA) are helping to cement the idea that close companion stars form via fragmentation of a gravitationally unstable disk around a protostar early in the star formation process. As the disk grows in mass, it eventually becomes gravitationally unstable and fragments, forming one or more new protostars in orbit with the first at mean separations of 100 astronomical units (AU) or even less. Here we report direct numerical calculations down to scales as small as $\sim 0.1$ AU, using a consistent Smoothed Particle Hydrodynamics (SPH) code, that show the large-scale fragmentation of a cloud core into two protostars accompanied by small-scale fragmentation of their circumstellar disks. Our results demonstrate the two dominant mechanisms of star formation, where the disk forming around a protostar, which in turn results from the large-scale fragmentation of the cloud core, undergoes eccentric ($m=1$) fragmentation to produce a close binary. We generate two-dimensional emission maps and simulated ALMA 1.3 mm continuum images of the structure and fragmentation of the disks that can help explain the dynamical processes occurring within collapsing cloud cores.

We study the consistency and convergence of smoothed particle hydrodynamics (SPH), as a function of the interpolation parameters, namely the number of particles $N$, the number of neighbors $n$, and the smoothing length $h$, using simulations of the collapse and fragmentation of protostellar rotating cores. The calculations are made using a modified version of the GADGET-2 code that employs an improved scheme for the artificial viscosity and power-law dependences of $n$ and $h$ on $N$, as was recently proposed by Zhu et al., which comply with the combined limit $N\to\infty$, $h\to 0$, and $n\to\infty$ with $n/N\to 0$ for full SPH consistency, as the domain resolution is increased. We apply this realization to the "standard isothermal test case" in the variant calculated by Burkert & Bodenheimer and the Gaussian cloud model of Boss to investigate the response of the method to adaptive smoothing lengths in the presence of large density and pressure gradients. The degree of consistency is measured by tracking how well the estimates of the consistency integral relations reproduce their continuous counterparts. In particular, $C^{0}$ and $C^{1}$ particle consistency is demonstrated, meaning that the calculations are close to second-order accuracy. As long as $n$ is increased with $N$, mass resolution also improves as the minimum resolvable mass $M_{\rm min}\sim n^{-1}$. This aspect allows proper calculation of small-scale structures in the flow associated with the formation and instability of protostellar disks around the growing fragments, which are seen to develop a spiral structure and fragment into close binary/multiple systems as supported by recent observations.

We present digital tables for the radiative terms that appear in the energy and momentum equations used to simulate the accretion onto supermassive black holes (SMBHs) in the center of galaxies. Cooling and heating rates and radiative accelerations are calculated with two different Spectral Energy Distributions (SEDs). One SED is composed of an accretion disk + [X-ray]-powerlaw, while the other is made of an accretion disk + [Corona]-bremsstrahlung with $T_X=1.16 \times 10^8$ K, where precomputed conditions of adiabatic expansion are included. Quantification of different physical mechanisms at operation are presented, showing discrepancies and similarities between both SEDs in different ranges of fundamental physical parameters (i.e., ionization parameter, density, and temperature). With the recent discovery of outflows originating at sub-parsec scales, these tables may provide a useful tool to model gas accretion processes onto a SMBH.

In this paper we show a plausible mechanism that could lead to the formation of the Dark Lanes in Lunar Swirls, and the electromagnetic shielding of the lunar surface that results in the preservation of the white colour of the lunar regolith. We present the results of a fully self-consistent 2 and 3 dimensional particle-in-cell simulations of mini-magnetospheres that form above the lunar surface and show that they are consistent with the formation of `lunar swirls' such as the archetypal formation Reiner Gamma. The simulations show how the microphysics of the deflection/shielding of plasma operates from a kinetic-scale cavity, and show that this interaction leads to a footprint with sharp features that could be the mechanism behind the generation of `dark lanes'. The physics of mini-magnetospheres is described and shown to be controlled by space-charge fields arising due to the magnetized electrons and unmagnetized ions. A comparison between model and observation is shown for a number of key plasma parameters.

We determine both representations of the Fundamental Plane [Re Sig^a Ie^b and Re (Sig^2/Ie)^lambda] and the luminosity-effective phase space density [L fe^(-gamma)] scaling relation for N-body remnants of binary mergers of spiral-like galaxies. The main set of merger simulations involves a mass-ratio of the progenitors in the range of about 1:1 to 1:5, harboring or not a bulge-like component, and are constructed using a cosmological motivated model. Equal-mass mergers are also considered. Remnants lead to average values for the scaling indices of <a>~1.6, <b>~0.6, <lambda>~0.7, and <gamma>~0.65. These values are consistent with those of K--band observations (Mobasher et al. 1999) of ellipticals: <a>~1.5, <b>~0.8, <lambda>~0.7, and <gamma>~0.60. The b index is, however, not well reproduced. This study does not allow us to establish a conclusive preference for models with or without a bulge as progenitors. Our results indicate that the L--fe and FP scalings might be determined to a large extent by dissipationless processes, a result that appears to be in contradiction to other dissipationless results.

Sersic parameters characterising the density profiles of remnants formed in collision-less disc galaxy mergers are obtained; no bulge is included in our simulations. For the luminous component we find that the Sersic index is n ~ (1.5,5.3) with <n> ~ 3 +/- 1 and an effective radius of R_e ~ (1.6,12.9) kpc with <R_e> ~ 5 +/- 3 kpc. A strong correlation of n with the central projected density I_0 is found [n ~ I_0^(-0.14)] which is consistent with observations. No positive linear correlation between the size (R_e) and structure (n) of our remnants is found; we do not advocate the existence of this. The photometric plane (PHP) of the luminous component [n ~ R_e^(0.05) I_0^(0.15)] agrees well, within the uncertainties and the assumption of a constant mass-to-light ratio, with those observationally determined for ellipticals. We found that the surface defined by Sersic parameters (n, R_e, mu_0) in log-space is not a true plane, but a pseudo-plane with a small curvature at low values of n owed to intrinsic properties of the Sersic model. The dark haloes of the remnants have a 3-dimensional Sersic index of <n> ~ 4 +/- 0.5 that are smaller than the ones obtained for dark haloes in LCDM cosmologies <n> ~ 6 +/- 1. A tight dark Sersic ``plane'' (DSP) is also defined by the parameters of the remnants haloes with n \~ r_e^(0.07) rho_0^(0.10). We conclude that collision-less merger remnants of pure disc galaxies have Sersic properties and correlations consistent with those of observed in early-type galaxies and local remnants. It seems that a ``primordial'' bulge in spirals is not a necessary condition to form bona fide ellipticals on grounds of the Sersic properties of remnants.

We consider our universe as a 3d domain wall embedded in a 5d dimensional Minkowski space-time. We address the problem of inflation and late time acceleration driven by bulk particles colliding with the 3d domain wall. The expansion of our universe is mainly related to these bulk particles. Since our universe tends to be permeated by a large number of isolated structures, as temperature diminishes with the expansion, we model our universe with a 3d domain wall with increasing internal structures. These structures could be unstable 2d domain walls evolving to fermi-balls which are candidates to cold dark matter. The momentum transfer of bulk particles colliding with the 3d domain wall is related to the reflection coefficient. We show a nontrivial dependence of the reflection coefficient with the number of internal dark matter structures inside the 3d domain wall. As the population of such structures increases the velocity of the domain wall expansion also increases. The expansion is exponential at early times and polynomial at late times. We connect this picture with string/M-theory by considering BPS 3d domain walls with structures which can appear through the bosonic sector of a five-dimensional supergravity theory.

We present the results of N-body simulations of the accretion of high-density dwarf galaxies by low-density giant galaxies. Both galaxies contain power-law central density cusps and point masses representing supermassive black holes; the ratio of galaxy masses is 3:1. The cusp of the dwarf galaxy is always disrupted during the merger, leading to a remnant with a weak power law in the intrinsic density and a ``core'' in the projected density. Removing both black holes from the giant and dwarf galaxies allows the dwarf galaxy to remain intact and leads to a remnant with a high central density, contrary to what is observed. Our results support the hypothesis than the persistence of low-density cores in giant galaxies following mergers is a consequence of the existence of supermassive central black holes.

We present 0.55 x 10^6 particle simulations of the accretion of high-density dwarf galaxies by low-density giant galaxies, using models that contain both power-law central density cusps and point masses representing supermassive black holes. The cusp of the dwarf galaxy is disrupted during the merger, producing a remnant with a central density that is only slightly higher than that of the giant galaxy initially. Removing the black hole from the giant galaxy allows the dwarf galaxy to remain intact and leads to a remnant with a high central density, contrary to what is observed. Our results support the hypothesis that the persistence of low-density cores in giant galaxies is a consequence of supermassive black holes.

Space-based observations are beginning to yield detailed information about the stellar kinematics at the very centers of galaxies, within the sphere of gravitational influence of the black hole. The structure and dynamics of these regions is probably determined in part by the infall and coalescence of black holes during galaxy mergers. A goal of N-body simulations is to reproduce the kinematics near the black holes as well as the relations that exist between the nuclear and global properties of galaxies. However, the problem is computationally difficult due to the wide range of length and time scales, and no single N-body code can efficiently follow the evolution from kiloparsec to sub-parsec scales. We review existing N-body work on this problem and present the first, fully self-consistent merger simulations of galaxies containing dense stellar cusps and black holes.