Pair-distribution functions of the two-dimensional electron gas

Based on its known exact properties and a new set of extensive fixed-node reptation quantum Monte Carlo simulations (both with and without backflow correlations, which in this case turn out to yield negligible improvements), we propose a new analytical representation of (i) the spin-summed pair-distribution function and (ii) the spin-resolved potential energy of the ideal two-dimensional interacting electron gas for a wide range of electron densities and spin polarization, plus (iii) the spin-resolved pair-distribution function of the unpolarized gas. These formulae provide an accurate reference for quantities previously not available in analytic form, and may be relevant to semiconductor heterostructures, metal-insulator transitions and quantum dots both directly, in terms of phase diagram and spin susceptibility, and indirectly, as key ingredients for the construction of new two-dimensional spin density functionals, beyond the local approximation.Comment: 12 pages, 10 figures; misprints correcte

Time-dependent Stochastic Modeling of Solar Active Region Energy

A time-dependent model for the energy of a flaring solar active region is presented based on a stochastic jump-transition model (Wheatland and Glukhov 1998; Wheatland 2008; Wheatland 2009). The magnetic free energy of the model active region varies in time due to a prescribed (deterministic) rate of energy input and prescribed (random) flare jumps downwards in energy. The model has been shown to reproduce observed flare statistics, for specific time-independent choices for the energy input and flare transition rates. However, many solar active regions exhibit time variation in flare productivity, as exemplified by NOAA active region AR 11029 (Wheatland 2010). In this case a time-dependent model is needed. Time variation is incorporated for two cases: 1. a step change in the rates of flare jumps; and 2. a step change in the rate of energy supply to the system. Analytic arguments are presented describing the qualitative behavior of the system in the two cases. In each case the system adjusts by shifting to a new stationary state over a relaxation time which is estimated analytically. The new model retains flare-like event statistics. In each case the frequency-energy distribution is a power law for flare energies less than a time-dependent rollover set by the largest energy the system is likely to attain at a given time. For Case 1, the model exhibits a double exponential waiting-time distribution, corresponding to flaring at a constant mean rate during two intervals (before and after the step change), if the average energy of the system is large. For Case 2 the waiting-time distribution is a simple exponential, again provided the average energy of the system is large. Monte Carlo simulations of Case~1 are presented which confirm the analytic estimates. The simulation results provide a qualitative model for observed flare statistics in active region AR 11029.Comment: 25 pages, 9 figure

Non-linear Evolution of Baryon Acoustic Oscillations from Improved Perturbation Theory in Real and Redshift Spaces

We study the non-linear evolution of baryon acoustic oscillations in the matter power spectrum and correlation function from the improved perturbation theory (PT). Based on the framework of renormalized PT, we apply the {\it closure approximation} that truncates the infinite series of loop contributions at one-loop order, and obtain a closed set of integral equations for power spectrum and non-linear propagator. The resultant integral expressions keep important non-perturbative properties which can dramatically improve the prediction of non-linear power spectrum. Employing the Born approximation, we then derive the analytic expressions for non-linear power spectrum and the predictions are made for non-linear evolution of baryon acoustic oscillations in power spectrum and correlation function. A detailed comparison between improved PT results and N-body simulations shows that a percent-level agreement is achieved in a certain range in power spectrum and in a rather wider range in correlation function. Combining a model of non-linear redshift-space distortion, we also evaluate the power spectrum and correlation function in correlation function. In contrast to the results in real space, the agreement between N-body simulations and improved PT predictions tends to be worse, and a more elaborate modeling for redshift-space distortion needs to be developed. Nevertheless, with currently existing model, we find that the prediction of correlation function has a sufficient accuracy compared with the cosmic-variance errors for future galaxy surveys with volume of a few (Gpc/h)^3 at z>=0.5.Comment: 25 pages, 15 figures, accepted for publication in Phys.Rev.

Primary beam effects of radio astronomy antennas -- II. Modelling the MeerKAT L-band beam

After a decade of design and construction, South Africa's SKA-MID precursor MeerKAT has begun its science operations. To make full use of the widefield capability of the array, it is imperative that we have an accurate model of the primary beam of its antennas. We have taken available L-band full-polarization 'astro-holographic' observations of three antennas and a generic electromagnetic simulation and created sparse representations of the beams using principal components and Zernike polynomials. The spectral behaviour of the spatial coefficients has been modelled using discrete cosine transform. We have provided the Zernike-based model over a diameter of 10 deg averaged over the beams of three antennas in an associated software tool (EIDOS) that can be useful in direction-dependent calibration and imaging. The model is more accurate for the diagonal elements of the beam Jones matrix and at lower frequencies. As we get more accurate beam measurements and simulations in the future, especially for the cross-polarization patterns, our pipeline can be used to create more accurate sparse representations of MeerKAT beams.Comment: 16 pages, 18 figures. This is a pre-copyedited, author-produced PDF of an article accepted for publication in MNRAS following peer review. The version of record [K. M. B. Asad et al., 2021] is available online at: https://doi.org/10.1093/mnras/stab10

A 3-mode, Variable Velocity Jet Model for HH 34

Variable ejection velocity jet models can qualitatively explain the appearance of successive working surfaces in Herbig-Haro (HH) jets. This paper presents an attempt to explore which features of the HH-34 jet can indeed be reproduced by such a model. From previously published data on this object, we find evidence for the existence of a 3-mode ejection velocity variability, and then explore the implications of such a variability. From simple, analytic considerations it is possible to show that the longer period modes produce a modulation on the shorter period modes, resulting in the formation of ``trains'' of multiple knots. The knots observed close to the source of HH-34 could correspond to such a structure. Finally, a numerical simulation with the ejection velocity variability deduced from the HH-34 data is computed. This numerical simulation shows a quite remarkable resemblance with the observed properties of the HH-34 jet.Comment: 28 pages LaTex, 10 postscript figure

Full-sky maps for gravitational lensing of the CMB

We use the large cosmological Millennium Simulation (MS) to construct the first all-sky maps of the lensing potential and the deflection angle, aiming at gravitational lensing of the CMB, with the goal of properly including small-scale non-linearities and non-Gaussianity. Exploiting the Born approximation, we implement a map-making procedure based on direct ray-tracing through the gravitational potential of the MS. We stack the simulation box in redshift shells up to $z\sim 11$, producing continuous all-sky maps with arcminute angular resolution. A randomization scheme avoids repetition of structures along the line of sight and structures larger than the MS box size are added to supply the missing contribution of large-scale (LS) structures to the lensing signal. The angular power spectra of the projected lensing potential and the deflection-angle modulus agree quite well with semi-analytic estimates on scales down to a few arcminutes, while we find a slight excess of power on small scales, which we interpret as being due to non-linear clustering in the MS. Our map-making procedure, combined with the LS adding technique, is ideally suited for studying lensing of CMB anisotropies, for analyzing cross-correlations with foreground structures, or other secondary CMB anisotropies such as the Rees-Sciama effect.Comment: LaTeX file, 10 pages, MNRAS in press, scales larger than the Millennium Simulation box size semi-analytically added, maps changed, references added, typos correcte

Magnetic heating across the cosmological recombination era: Results from 3D MHD simulations

The origin of cosmic magnetic fields is an unsolved problem and magnetogenesis could have occurred in the early Universe. We study the evolution of such primordial magnetic fields across the cosmological recombination epoch via 3D magnetohydrodynamic numerical simulations. We compute the effective or net heating rate of baryons due to decaying magnetic fields and its dependence on the magnetic field strength and spectral index. In the drag-dominated regime ($z \gtrsim 1500$), prior to recombination, we find no real heating is produced. Our simulations allow us to smoothly trace a new transition regime ($600 \lesssim z \lesssim 1500$), where magnetic energy decays, at first, into the kinetic energy of baryons. A turbulent velocity field is built up until it saturates, as the net heating rate rises from a low value at recombination to its peak towards the end of the transition regime. This is followed by a turbulent decay regime ($z \lesssim 600$) where magnetic energy dissipates via turbulent decay of both magnetic and velocity fields while net heating remains appreciable and declines slowly. Both the peak of the net heating rate and the onset of turbulent decay are delayed significantly beyond recombination, by up to 0.5 Myr (until $z\simeq 600-700$), for scale-invariant magnetic fields. We provide analytic approximations and present numerical results for a range of field strengths and spectral indices, illustrating the redshift-dependence of dissipation and net heating rates. These can be used to study cosmic microwave background constraints on primordial magnetic fields.Comment: Submitted to MNRAS, comments are welcome; 22 pages, 26 figures, 2 table