Kinetic Inductance of Josephson Junction Arrays: Dynamic and Equilibrium Calculations

We show analytically that the inverse kinetic inductance $L^{-1}$ of an overdamped junction array at low frequencies is proportional to the admittance of an inhomogeneous equivalent impedance network. The $ij^{th}$ bond in this equivalent network has an inverse inductance $J_{ij}\cos(\theta_i^0-\theta_j^0-A_{ij})$, where $J_{ij}$ is the Josephson coupling energy of the $ij^{th}$ bond, $\theta_i^0$ is the ground-state phase of the grain $i$, and $A_{ij}$ is the usual magnetic phase factor. We use this theorem to calculate $L^{-1}$ for square arrays as large as $180\times 180$. The calculated $L^{-1}$ is in very good agreement with the low-temperature limit of the helicity modulus $\gamma$ calculated by conventional equilibrium Monte Carlo techniques. However, the finite temperature structure of $\gamma$, as a function of magnetic field, is \underline{sharper} than the zero-temperature $L^{-1}$, which shows surprisingly weak structure. In triangular arrays, the equilibrium calculation of $\gamma$ yields a series of peaks at frustrations $f = \frac{1}{2}(1-1/N)$, where $N$ is an integer $\geq 2$, consistent with experiment.Comment: 14 pages + 6 postscript figures, 3.0 REVTe

Stratification of malaria incidence in Papua New Guinea (2011-2019): contribution towards a sub-national control policy

Malaria risk in Papua New Guinea (PNG) is highly heterogeneous, between and within geographical regions, which is operationally challenging for control. To enhance targeting of malaria interventions in PNG, we investigated risk factors and stratified malaria incidence at the level of health facility catchment areas. Catchment areas and populations of 808 health facilities were delineated using a travel-time accessibility approach and linked to reported malaria cases (2011-2019). Zonal statistics tools were used to calculate average altitude and air temperature in catchment areas before they were spatially joined with incidence rates. In addition, empirical Bayesian kriging (EBK) was employed to interpolate incidence risk strata across PNG. Malaria annual incidence rates are, on average, 186.3 per 1000 population in catchment areas up to 600 m, dropped to 98.8 at (800-1400) m, and to 24.1 cases above 1400 m altitude. In areas above the two altitudinal thresholds 600m and 1400m, the average annual temperature drops below 22°C and 17°C, respectively. EBK models show very low- to low-risk strata ( 200 per 1000) strata are modelled mainly in Momase and Islands Regions. Besides, strata with moderate risk (100-200) predominate throughout the coastal areas. While 35.7% of the PNG population (estimated 3.33 million in 2019) lives in places at high or moderate risk of malaria, 52.2% (estimated 4.88 million) resides in very low-risk areas. In five provinces, relatively large proportions of populations (> 50%) inhabit high-risk areas: New Ireland, East and West New Britain, Sandaun and Milne Bay. Incidence maps show a contrast in malaria risk between coastal and inland areas influenced by altitude. However, the risk is highly variable in low-lying areas. Malaria interventions should be guided by sub-national risk levels in PNG

First-Order Vortex Lattice Melting and Magnetization of YBa2_2Cu3_3O$_{7-\delta}

We present the first non-mean-field calculation of the magnetization $M(T)$ of YBa$_2$Cu$_3$O$_{7-\delta}$ both above and below the flux-lattice melting temperature $T_m(H)$. The results are in good agreement with experiment as a function of transverse applied field $H$. The effects of fluctuations in both order parameter $\psi({\bf r})$ and magnetic induction $B$ are included in the Ginzburg-Landau free energy functional: $\psi({\bf r})$ fluctuates within the lowest Landau level in each layer, while $B$ fluctuates uniformly according to the appropriate Boltzmann factor. The second derivative $(\partial^2 M/\partial T^2)_H$ is predicted to be negative throughout the vortex liquid state and positive in the solid state. The discontinuities in entropy and magnetization at melting are calculated to be $\sim 0.034\, k_B$ per flux line per layer and $\sim 0.0014$~emu~cm$^{-3}$ at a field of 50 kOe.Comment: 11 pages, 4 PostScript figures in one uuencoded fil

Vector composition, abundance, biting patterns and malaria transmission intensity in Madang, Papua New Guinea: assessment after 7 years of an LLIN-based malaria control programme

Background: A malaria control programme based on distribution of long-lasting insecticidal bed nets (LLINs) and artemisinin combination therapy began in Papua New Guinea in 2009. After implementation of the programme, substantial reductions in vector abundance and malaria transmission intensity occurred. The research reported here investigated whether these reductions remained after seven years of sustained effort. Methods: All-night (18:00 to 06:00) mosquito collections were conducted using human landing catches and barrier screen methods in four villages of Madang Province between September 2016 and March 2017. Anopheles species identification and sporozoite infection with Plasmodium vivax and Plasmodium falciparum were determined with molecular methods. Vector composition was expressed as the relative proportion of different species in villages, and vector abundance was quantified as the number of mosquitoes per barrier screen-night and per person-night. Transmission intensity was quantified as the number of sporozoite-infective vector bites per person-night. Results: Five Anopheles species were present, but vector composition varied greatly among villages. Anopheles koliensis, a strongly anthropophilic species was the most prevalent in Bulal, Matukar and Wasab villages, constituting 63.7–73.8% of all Anopheles, but in Megiar Anopheles farauti was the most prevalent species (97.6%). Vector abundance varied among villages (ranging from 2.8 to 72.3 Anopheles per screen-night and 2.2–31.1 Anopheles per person-night), and spatially within villages. Malaria transmission intensity varied among the villages, with values ranging from 0.03 to 0.5 infective Anopheles bites per person-night. Most (54.1–75.1%) of the Anopheles bites occurred outdoors, with a substantial proportion (25.5–50.8%) occurring before 22:00. Conclusion: The estimates of vector abundance and transmission intensity in the current study were comparable to or higher than estimates in the same villages in 2010–2012, indicating impeded programme effectiveness. Outdoor and early biting behaviours of vectors are some of the likely explanatory factors. Heterogeneity in vector composition, abundance and distribution among and within villages challenge malaria control programmes and must be considered when planning them

Vortices and 2D bosons: A Path-Integral Monte Carlo Study

The vortex system in a high-T_c superconductor has been studied numerically using the mapping to 2D bosons and the path-integral Monte Carlo method. We find a single first-order transition from an Abrikosov lattice to an entangled vortex liquid. The transition is characterized by an entropy jump dS = 0.4 k_B per vortex and layer (parameters for YBCO) and a Lindemann number c_L = 0.25. The increase in density at melting is given by d\rho = 6.0*10^{-4} / \lambda(T)^2. The vortex liquid corresponds to a bosonic superfluid, with \rho_s = \rho even in the limit \lambda -> \infty.Comment: 9 pages, RevTeX, 4 PostScript figures. The entropy jump at the transition has been recomputed and is now in agreement with experiments on YBCO. Some minor modifications were made in the tex

Flux Pinning and Phase Transitions in Model High-Temperature Superconductors with Columnar Defects

We calculate the degree of flux pinning by defects in model high-temperature superconductors (HTSC's). The HTSC is modeled as a three-dimensional network of resistively-shunted Josephson junctions in an external magnetic field, corresponding to a HTSC in the extreme Type-II limit. Disorder is introduced either by randomizing the coupling between grains (Model A disorder) or by removing grains (Model B disorder). Three types of defects are considered: point disorder, random line disorder, and periodic line disorder; but the emphasis is on random line disorder. Static and dynamic properties of the models are determined by Monte Carlo simulations and by solution of the analogous coupled overdamped Josephson equations in the presence of thermal noise. Random line defects considerably raise the superconducting transition temperature T$_c(B)$, and increase the apparent critical current density J$_c(B,T)$, in comparison to the defect-free crystal. They are more effective in these respects than a comparable volume density of point defects, in agreement with the experiments of Civale {\it et al}. Periodic line defects commensurate with the flux lattice are found to raise T$_c(B)$ even more than do random line defects. Random line defects are most effective when their density approximately equals the flux density. Near T$_c(B)$, our static and dynamic results appear consistent with the anisotropic Bose glass scaling hypotheses of Nelson and Vinokur, but with possibly different critical indices:Comment: 10 pages, LaTeX(REVTeX v3.0, twocolumn), 11 figures (not included

First-Order Melting and Dynamics of Flux Lines in a Model for YBa2_2Cu3_3O7−δ_{7-\delta}

We have studied the statics and dynamics of flux lines in a model for YBCO, using both Monte Carlo simulations and Langevin dynamics. For a clean system, both approaches yield the same melting curve, which is found to be weakly first order with a heat of fusion of about $0.02 k_BT_m$ per vortex pancake at a field of $50 {\rm kG}.$ The time averaged magnetic field distribution experienced by a fixed spin is found to undergo a qualitative change at freezing, in agreement with NMR and $\mu {\rm SR}$ experiments. Melting in the clean system is accompanied by a proliferation of free disclinations which show a clear B-dependent 3D-2D crossover from long disclination lines parallel to the c-axis at low fields, to 2D ``pancake'' disclinations at higher fields. Strong point pins produce a logarithmical $\ln t$ relaxation which results from slow annealing out of disclinations in disordered samples.Comment: 31 pages, latex, revtex, 12 figures available upon request, No major changes to the original text, but some errors in the axes scale for Figures 6 and 7 were corrected(new figures available upon request), to be published in Physical Review B, July 199

Grundstate Properties of the 3D Ising Spin Glass

We study zero--temperature properties of the 3d Edwards--Anderson Ising spin glass on finite lattices up to size $12^3$. Using multicanonical sampling we generate large numbers of groundstate configurations in thermal equilibrium. Finite size scaling with a zero--temperature scaling exponent $y = 0.74 \pm 0.12$ describes the data well. Alternatively, a descriptions in terms of Parisi mean field behaviour is still possible. The two scenarios give significantly different predictions on lattices of size $\ge 12^3$.Comment: LATEX 9pages,figures upon request ,SCRI-9

Metamagnetism in the 2D Hubbard Model with easy axis

Although the Hubbard model is widely investigated, there are surprisingly few attempts to study the behavior of such a model in an external magnetic field. Using the Projector Quantum Monte Carlo technique, we show that the Hubbard model with an easy axis exhibits metamagnetic behavior if an external field is turned on. For the case of intermediate correlations strength $U$, we observe a smooth transition from an antiferromagnetic regime to a paramagnetic phase. While the staggered magnetization will decrease linearly up to a critical field $B_c$, uniform magnetization develops only for fields higher than $B_c$.Comment: RevTeX 5 pages + 2 postscript figures (included), accepted for PRB Rapid Communication