Hybrid computer techniques for solving partial differential equations

Techniques overcome equipment limitations that restrict other computer techniques in solving trivial cases. The use of curve fitting by quadratic interpolation greatly reduces required digital storage space

Effects of a torsion field on Big Bang nucleosynthesis

In this paper it is investigated whether torsion, which arises naturally in most theories of quantum gravity, has observable implications for the Big Bang nucleosynthesis. Torsion can lead to spin flips amongst neutrinos thus turning them into sterile neutrinos. In the early Universe they can alter the helium abundance which is tightly constrained by observations. Here I calculate to what extent torsion of the string theory type leads to a disagreement with the Big Bang nucleosynthesis predictions.Comment: accepted by General Relativity and Gravitatio

Superconducting niobium thin film slow-wave structures

A superconducting comb structure as a slow-wave element in a traveling-wave maser will significantly improve maser noise temperature and gain by reducing the insertion loss. The results of the insertion loss measurements of superconducting niobium slow-wave structures subjected to maser operating conditions at X-Band frequencies are presented

The antibody loci of the domestic goat (Capra hircus)

The domestic goat (Capra hircus) is an important ruminant species both as a source of antibody-based reagents for research and biomedical applications and as an economically important animal for agriculture, particularly for developing nations that maintain most of the global goat population. Characterization of the loci encoding the goat immune repertoire would be highly beneficial for both vaccine and immune reagent development. However, in goat and other species whose reference genomes were generated using short-read sequencing technologies, the immune loci are poorly assembled as a result of their repetitive nature. Our recent construction of a long-read goat genome assembly (ARS1) has facilitated characterization of all three antibody loci with high confidence and comparative analysis to cattle. We observed broad similarity of goat and cattle antibody-encoding loci but with notable differences that likely influence formation of the functional antibody repertoire. The goat heavy-chain locus is restricted to only four functional and nearly identical IGHV genes, in contrast to the ten observed in cattle. Repertoire analysis indicates that light-chain usage is more balanced in goats, with greater representation of kappa light chains (~ 20-30%) compared to that in cattle (~ 5%). The present study represents the first characterization of the goat antibody loci and will help inform future investigations of their antibody responses to disease and vaccination

Pressure-induced diamond to beta-tin transition in bulk silicon: a near-exact quantum Monte Carlo study

The pressure-induced structural phase transition from diamond to beta-tin in silicon is an excellent test for theoretical total energy methods. The transition pressure provides a sensitive measure of small relative energy changes between the two phases (one a semiconductor and the other a semimetal). Experimentally, the transition pressure is well characterized. Density-functional results have been unsatisfactory. Even the generally much more accurate diffusion Monte Carlo method has shown a noticeable fixed-node error. We use the recently developed phaseless auxiliary-field quantum Monte Carlo (AFQMC) method to calculate the relative energy differences in the two phases. In this method, all but the error due to the phaseless constraint can be controlled systematically and driven to zero. In both structural phases we were able to benchmark the error of the phaseless constraint by carrying out exact unconstrained AFQMC calculations for small supercells. Comparison between the two shows that the systematic error in the absolute total energies due to the phaseless constraint is well within 0.5 mHa/atom. Consistent with these internal benchmarks, the transition pressure obtained by the phaseless AFQMC from large supercells is in very good agreement with experiment.Comment: 9 pages, 5 figure

Two-dimensional molecular para-hydrogen and ortho-deuterium at zero temperature

We study molecular para-hydrogen (p-${\rm H_{2}}$) and ortho-deuterium (o-${\rm D_{2}}$) in two dimensions and in the limit of zero temperature by means of the diffusion Monte Carlo method. We report energetic and structural properties of both systems like the total and kinetic energy per particle, radial pair distribution function, and Lindemann's ratio in the low pressure regime. By comparing the total energy per particle as a function of the density in liquid and solid p-${\rm H_{2}}$, we show that molecular para-hydrogen, and also ortho-deuterium, remain solid at zero temperature. Interestingly, we assess the quality of three different symmetrized trial wave functions, based on the Nosanow-Jastrow model, in the p-${\rm H_{2}}$ solid film at the variational level. In particular, we analyze a new type of symmetrized trial wave function which has been used very recently to describe solid $^{4}$He and found that also characterizes hydrogen satisfactorily. With this wave function, we show that the one-body density matrix $\varrho_{1} (r)$ of solid p-${\rm H_{2}}$ possesses off-diagonal long range order, with a condensate fraction that increases sizably in the negative pressure regime.Comment: 11 pages, 9 figure

The Fermion Monte Carlo revisited

In this work we present a detailed study of the Fermion Monte Carlo algorithm (FMC), a recently proposed stochastic method for calculating fermionic ground-state energies [M.H. Kalos and F. Pederiva, Phys. Rev. Lett. vol. 85, 3547 (2000)]. A proof that the FMC method is an exact method is given. In this work the stability of the method is related to the difference between the lowest (bosonic-type) eigenvalue of the FMC diffusion operator and the exact fermi energy. It is shown that within a FMC framework the lowest eigenvalue of the new diffusion operator is no longer the bosonic ground-state eigenvalue as in standard exact Diffusion Monte Carlo (DMC) schemes but a modified value which is strictly greater. Accordingly, FMC can be viewed as an exact DMC method built from a correlated diffusion process having a reduced Bose-Fermi gap. As a consequence, the FMC method is more stable than any transient method (or nodal release-type approaches). We illustrate the various ideas presented in this work with calculations performed on a very simple model having only nine states but a full sign problem. Already for this toy model it is clearly seen that FMC calculations are inherently uncontrolled.Comment: 49 pages with 4 postscript figure

Generation and measurement of nonstationary random processes technical note no. 3

Generation and measurement of nonstationary stochastic processes related to Monte Carlo studies with analog compute

Dilute Bose gases interacting via power-law potentials

Neutral atoms interact through a van der Waals potential which asymptotically falls off as r^{-6}. In ultracold gases, this interaction can be described to a good approximation by the atom-atom scattering length. However, corrections arise that depend on the characteristic length of the van der Waals potential. We parameterize these corrections by analyzing the energies of two- and few-atom systems under external harmonic confinement, obtained by numerically and analytically solving the Schrodinger equation. We generalize our results to particles interacting through a longer-ranged potential which asymptotically falls off as r^{-4}.Comment: 7 pages, 4 figure