The surface densities of molecular gas, $\Sigma_{\rm H_2}$, and the star formation rate (SFR), $\dot\Sigma_\star$, correlate almost linearly on kiloparsec scales in observed star-forming (non-starburst) galaxies. We explore the origin of the linear slope of this correlation using a suite of isolated $L_\star$ galaxy simulations. We show that in simulations with efficient feedback, the slope of the $\dot\Sigma_\star$-$\Sigma_{\rm H_2}$ relation on kiloparsec scales is insensitive to the slope of the $\dot\rho_\star$-$\rho$ relation assumed at the resolution scale. We also find that the slope on kiloparsec scales depends on the criteria used to identify star-forming gas, with a linear slope arising in simulations that identify star-forming gas using a virial parameter threshold. This behavior can be understood using a simple theoretical model based on conservation of interstellar gas mass as the gas cycles between atomic, molecular, and star-forming states under the influence of feedback and dynamical processes. In particular, we show that the linear slope emerges when feedback efficiently regulates and stirs the evolution of dense, molecular gas. We show that the model also provides insights into the likely origin of the relation between the SFR and molecular gas in real galaxies on different scales.

Using a suite of isolated $L_\star$ galaxy simulations, we show that global depletion times and star-forming gas mass fractions in simulated galaxies exhibit systematic and well-defined trends as a function of the local star formation efficiency per freefall time, $\epsilon_{\rm ff}$, strength of stellar feedback, and star formation threshold. We demonstrate that these trends can be reproduced and explained by a simple physical model of global star formation in galaxies. Our model is based on mass conservation and the idea of gas cycling between star-forming and non-star-forming states on certain characteristic time scales under the influence of dynamical and feedback processes. Both the simulation results and our model predictions exhibit two limiting regimes with rather different dependencies of global galactic properties on the local parameters. When $\epsilon_{\rm ff}$ is small and feedback is inefficient, the total star-forming mass fraction, $f_{\rm sf}$, is independent of $\epsilon_{\rm ff}$ and the global depletion time, $\tau_{\rm dep}$, scales inversely with $\epsilon_{\rm ff}$. When $\epsilon_{\rm ff}$ is large or feedback is very efficient, these trends are reversed: $f_{\rm sf} \propto \epsilon_{\rm ff}^{-1}$ and $\tau_{\rm dep}$ is independent of $\epsilon_{\rm ff}$ but scales linearly with the feedback strength. We also compare our results with the observed depletion times and mass fractions of star-forming and molecular gas and show that they provide complementary constraints on $\epsilon_{\rm ff}$ and the feedback strength. We show that useful constraints on $\epsilon_{\rm ff}$ can also be obtained using measurements of the depletion time and its scatter on different spatial scales.

Using an isolated Milky Way-mass galaxy simulation, we compare results from 9 state-of-the-art gravito-hydrodynamics codes widely used in the numerical community. We utilize the infrastructure we have built for the AGORA High-resolution Galaxy Simulations Comparison Project. This includes the common disk initial conditions, common physics models (e.g., radiative cooling and UV background by the standardized package Grackle) and common analysis toolkit yt, all of which are publicly available. Subgrid physics models such as Jeans pressure floor, star formation, supernova feedback energy, and metal production are carefully constrained across code platforms. With numerical accuracy that resolves the disk scale height, we find that the codes overall agree well with one another in many dimensions including: gas and stellar surface densities, rotation curves, velocity dispersions, density and temperature distribution functions, disk vertical heights, stellar clumps, star formation rates, and Kennicutt-Schmidt relations. Quantities such as velocity dispersions are very robust (agreement within a few tens of percent at all radii) while measures like newly-formed stellar clump mass functions show more significant variation (difference by up to a factor of ~3). Systematic differences exist, for example, between mesh-based and particle-based codes in the low density region, and between more diffusive and less diffusive schemes in the high density tail of the density distribution. Yet intrinsic code differences are generally small compared to the variations in numerical implementations of the common subgrid physics such as supernova feedback. Our experiment reassures that, if adequately designed in accordance with our proposed common parameters, results of a modern high-resolution galaxy formation simulation are more sensitive to input physics than to intrinsic differences in numerical schemes.

An optimally efficient explicit numerical scheme for solving fluid dynamics equations, or any other parabolic or hyperbolic system of partial differential equations, should allow local regions to advance in time with their own, locally constrained time steps. However, such a scheme can result in violation of the Courant-Friedrichs-Lewy (CFL) condition, which is manifestly non-local. Although the violations can be considered to be "weak" in a certain sense and the corresponding numerical solution may be stable, such calculation does not guarantee the correct propagation speed for arbitrary waves. We use an experimental fluid dynamics code that allows cubic "patches" of grid cells to step with independent, locally constrained time steps to demonstrate how the CFL condition can be enforced by imposing a condition on the time steps of neighboring patches. We perform several numerical tests that illustrate errors introduced in the numerical solutions by weak CFL condition violations and show how strict enforcement of the CFL condition eliminates these errors. In all our tests the strict enforcement of the CFL condition does not impose a significant performance penalty.

We present a model that explains why galaxies form stars on a time scale significantly longer than the time scales of processes governing the evolution of interstellar gas. We show that gas evolves from a non-star-forming to a star-forming state on a relatively short time scale and thus the rate of this evolution does not limit the star formation rate. Instead, the star formation rate is limited because only a small fraction of star-forming gas is converted into stars before star-forming regions are dispersed by feedback and dynamical processes. Thus, gas cycles into and out of star-forming state multiple times, which results in a long time scale on which galaxies convert gas into stars. Our model does not rely on the assumption of equilibrium and can be used to interpret trends of depletion times with the properties of observed galaxies and the parameters of star formation and feedback recipes in simulations. In particular, the model explains how feedback self-regulates the star formation rate in simulations and makes it insensitive to the local star formation efficiency. We illustrate our model using the results of an isolated $L_*$-sized galaxy simulation that reproduces the observed Kennicutt-Schmidt relation for both molecular and atomic gas. Interestingly, the relation for molecular gas is almost linear on kiloparsec scales, although a nonlinear relation is adopted in simulation cells. We discuss how a linear relation emerges from non-self-similar scaling of the gas density PDF with the average gas surface density.

We present a new implementation of star formation in cosmological simulations, by considering star clusters as a unit of star formation. Cluster particles grow in mass over several million years at the rate determined by local gas properties, with high time resolution. The particle growth is terminated by its own energy and momentum feedback on the interstellar medium. We test this implementation for Milky Way-sized galaxies at high redshift, by comparing the properties of model clusters with observations of young star clusters. We find that the cluster initial mass function is best described by a Schechter function rather than a single power law. In agreement with observations, at low masses the logarithmic slope is $\alpha\approx 1.8-2$, while the cutoff at high mass scales with the star formation rate. A related trend is a positive correlation between the surface density of star formation rate and fraction of stars contained in massive clusters. Both trends indicate that the formation of massive star clusters is preferred during bursts of star formation. These bursts are often associated with major merger events. We also find that the median timescale for cluster formation ranges from 0.5 to 4 Myr and decreases systematically with increasing star formation efficiency. Local variations in the gas density and cluster accretion rate naturally lead to the scatter of the overall formation efficiency by an order of magnitude, even when the instantaneous efficiency is kept constant. Comparison of the formation timescale with the observed age spread of young star clusters provides an additional important constraint on the modeling of star formation and feedback schemes.

We present a study of a star formation prescription in which star formation efficiency depends on local gas density and turbulent velocity dispersion, as suggested by direct simulations of SF in turbulent giant molecular clouds (GMCs). We test the model using a simulation of an isolated Milky Way-sized galaxy with a self-consistent treatment of turbulence on unresolved scales. We show that this prescription predicts a wide variation of local star formation efficiency per free-fall time, $\epsilon_{\rm ff} \sim 0.1 - 10\%$, and gas depletion time, $t_{\rm dep} \sim 0.1 - 10$ Gyr. In addition, it predicts an effective density threshold for star formation due to suppression of $\epsilon_{\rm ff}$ in warm diffuse gas stabilized by thermal pressure. We show that the model predicts star formation rates in agreement with observations from the scales of individual star-forming regions to the kiloparsec scales. This agreement is non-trivial, as the model was not tuned in any way and the predicted star formation rates on all scales are determined by the distribution of the GMC-scale densities and turbulent velocities $\sigma$ in the cold gas within the galaxy, which is shaped by galactic dynamics. The broad agreement of the star formation prescription calibrated in the GMC-scale simulations with observations, both gives credence to such simulations and promises to put star formation modeling in galaxy formation simulations on a much firmer theoretical footing.

Petschek-type time-dependent reconnection (TDR) and quasi-stationary reconnection (QSR) models are considered to understand reconnection outflow structures and the features of the associated locally generated turbulence in the solar wind. We show that the outflow structures, such as discontinuites, Kelvin-Helmholtz (KH) unstable flux tubes or continuous space filling flows cannot be distinguished from one-point WIND measurements. In both models the reconnection outflows can generate more or less spatially extended turbulent boundary layers (TBDs). The structure of an unique extended reconnection outflow is investigated in detail. The analysis of spectral scalings and break locations show that reconnection outflows can control the local field and plasma conditions which may play in favor of one or another turbulent dissipation mechanisms with their characteristic scales and wavenumbers.

The exact mechanism by which astrophysical jets are formed is still unknown. It is believed that necessary elements are a rotating (Kerr) black hole and a magnetised accreting plasma. We model the accreting plasma as a collection of magnetic flux tubes/strings. If such a tube falls into a Kerr black hole, then the leading portion loses angular momentum and energy as the string brakes, and to compensate for this loss, momentum and energy is redistributed to the trailing portion of the tube.} {We found that buoyancy creates a pronounced helical magnetic field structure aligned with the spin axis. Along the field lines, the plasma is centrifugally accelerated close to the speed of light. This process leads to unlimited stretching of the flux tube since one part of the tube continues to fall into the black hole and simultaneously the other part of the string is pushed outward. Eventually, reconnection cuts the tube, the inner part is filled with new material and the outer part forms a collimated bubble-structured relativistic jet. Each plasmoid can be considered like an outgoing particle in the Penrose mechanism: it carries extracted rotational energy away from the black hole while the falling part with the corresponding negative energy is left inside the ergosphere.

This work investigates a set of cosmological collisionless N-body simulations with featured power spectra of initial perturbations in the context of the core-cusp and satellites problems. On the studied power spectra some scales of fluctuations were suppressed. Such spectral features can be caused by multicomponent dark matter. The density profiles innermost resolved slopes $\alpha\equiv d \log(\rho) /d \log(r) $ of the five largest haloes were measured and its dependence on the parameters of the suppression was traced. In a certain range of the parameters the slopes flatten from initial value of about -1.2 to -0.6 or even to -0.2 in one of the cases. This qualitatively demonstrates that (i) profiles shape depends on initial power spectrum and (ii) this effect may be responsible for the solution of the core-cusp problem. The suppression of some part of the initial power spectrum also leads to the decrease of the number of massive subhaloes.

The origin of jets emitted from black holes is not well understood, however there are two possible energy sources, the accretion disk or the rotating black hole. Magnetohydrodynamic simulations show a well-defined jet that extracts energy from a black hole. If plasma near the black hole is threaded by large-scale magnetic flux, it will rotate with respect to asymptotic infinity creating large magnetic stresses. These stresses are released as a relativistic jet at the expense of black hole rotational energy. The physics of the jet initiation in the simulations is described by the theory of black hole gravitohydromagnetics.

Relativistic magnetohydrodynamics can be reformulated in terms of magnetic flux tubes which turn out to obey equations of non-linear strings. The string approach is applied to study how a test flux tube falls into a Kerr black hole. Analytical treatment and numerical simulations show that the leading portion of the falling tube loses angular momentum and energy as the string breaks, and to compensate for this loss, momentum and energy has to be generated to conserve energy and momentum for the tube. Inside the ergosphere the energy of the leading part can be negative, and the rest of the tube then extracts energy from the hole. Increasing centrifugal forces eject the part of the tube with extra positive energy from the ergosphere after a time producing a relativistic jet.

In the MHD description of plasma phenomena the concept of magnetic field lines frozen into the plasma turns out to be very useful. We present here a method of introducing Lagrangian coordinates into relativistic MHD equations in general relativity, which enables a convenient mathematical formulation for the behaviour of flux tubes. With the introduction of these Lagrangian, so--called ``frozen--in'' coordinates, the relativistic MHD equations reduce to a set of nonlinear 1D string equations, and the plasma may therefore be regarded as a gas of nonlinear strings corresponding to flux tubes. Numerical simulation shows that if such a tube/string falls into a Kerr black hole, then the leading portion loses angular momentum and energy as the string brakes, and to compensate for this loss, momentum and energy is radiated to infinity to conserve energy and momentum for the tube. Inside the ergosphere the energy of the leading part turns out to be negative after some time, and the rest of the tube then gets energy from the hole. In our simulations most of the compensated positive energy is also localized inside the ergosphere because the inward speed of the plasma is approximately equal to the velocity of the MHD wave which transports energy outside. Therefore, an additional physical process has to be included which can remove energy from the ergophere. Magnetic reconnection seems fills this role releasing Maxwellian stresses and producing a relativistic jet.