Source author record

J. Osborn

J. Osborn appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

42works
12topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

42 published item(s)

preprint2025arXiv

All-sky neutrino point-source search with IceCube combined track and cascade data

Despite extensive efforts, discovery of high-energy astrophysical neutrino sources remains elusive. We present an event-level simultaneous maximum likelihood analysis of tracks and cascades using IceCube data collected from 04/06/2008 to 05/23/2022 to search the whole sky for neutrino sources and, using a source catalog, for coincidence of neutrino emission with gamma-ray emission. This is the first time a simultaneous fit of different detection channels is used to conduct a time-integrated all-sky scan with IceCube. Combining all-sky tracks, with superior pointing-power and sensitivity in the northern sky, with all-sky cascades, with good energy-resolution and sensitivity in the southern sky, we have developed the most sensitive point-source search to date by IceCube which targets the entire sky. The most significant point in the northern sky aligns with NGC 1068, a Seyfert II galaxy, which, from the catalog search, shows a 3.5$σ$ excess over background after accounting for trials. The most significant point in the southern sky does not align with any source in the catalog and is not significant after accounting for trials. A search for the single most significant Gaussian flare at the locations of NGC 1068, PKS 1424+240, and the southern highest significance point shows results consistent with expectations for steady emission. Notably, this is the first time that a flare shorter than four years has been excluded as being responsible for NGC 1068's emergence as a neutrino source. Our results show that combining tracks and cascades when conducting neutrino source searches improves sensitivity and can lead to new discoveries.

preprint2022arXiv

A low-threshold ultrahigh-energy neutrino search with the Askaryan Radio Array

In the pursuit of the measurement of the still-elusive ultrahigh-energy (UHE) neutrino flux at energies of order EeV, detectors using the in-ice Askaryan radio technique have increasingly targeted lower trigger thresholds. This has led to improved trigger-level sensitivity to UHE neutrinos. Working with data collected by the Askaryan Radio Array (ARA), we search for neutrino candidates at the lowest threshold achieved to date, leading to improved analysis-level sensitivities. A neutrino search on a data set with 208.7~days of livetime from the reduced-threshold fifth ARA station is performed, achieving a 68\% analysis efficiency over all energies on a simulated mixed-composition neutrino flux with an expected background of $0.10_{-0.04}^{+0.06}$ events passing the analysis. We observe one event passing our analysis and proceed to set a neutrino flux limit using a Feldman-Cousins construction. We show that the improved trigger-level sensitivity can be carried through an analysis, motivating the Phased Array triggering technique for use in future radio-detection experiments. We also include a projection using all available data from this detector. Finally, we find that future analyses will benefit from studies of events near the surface to fully understand the background expected for a large-scale detector.

preprint2022arXiv

Triboelectric Backgrounds to radio-based UHE Neutrino Exeperiments

The proposed IceCube-Gen2 (ICG2) seeks to instrument ~500 sq. km of Antarctic ice near the geographic South Pole with radio antennas, in order to observe the highest energy (E>1 EeV) neutrinos in the Universe. To this end, ICG2 will use the impulsive radio-frequency (RF) signal produced by neutrino interactions in polar ice caps. In such experiments, rare single event candidates must be unambiguously separated from background; to date, signal identification strategies primarily reject thermal noise and anthropogenic backgrounds. Here, we consider the possibility that fake neutrino signals may also be naturally generated via the 'triboelectric effect'. This broadly includes any process in which force applied at a boundary layer results in displacement of surface charge, generating a potential difference ΔV. Wind blowing over granular surfaces such as snow can induce such a ΔV, with subsequent discharge. Discharges over nanosecond-timescales can then lead to RF emissions at characteristic MHz-GHz frequencies. We find that such backgrounds are evident in the several neutrino experiments considered, and are generally characterized by: a) a threshold wind velocity which likely depends on the experimental signal trigger threshold and layout; for the experiments considered herein, this value is typically O(10 m/s), b) frequency spectra generally shifted to the low-end of the frequency regime to which current radio experiments are typically sensitive (100-200 MHz), c) for the strongest background signals, an apparent preference for discharges from above-surface structures, although the presence of more isotropic, lower amplitude triboelectric discharges cannot be excluded.

preprint2020arXiv

Constraints on the Diffuse Flux of Ultra-High Energy Neutrinos from Four Years of Askaryan Radio Array Data in Two Stations

The Askaryan Radio Array (ARA) is an ultra-high energy (UHE, $>10^{17}$ eV) neutrino detector designed to observe neutrinos by searching for the radio waves emitted by the relativistic products of neutrino-nucleon interactions in Antarctic ice. In this paper, we present constraints on the diffuse flux of ultra-high energy neutrinos between $10^{16}-10^{21}$ eV resulting from a search for neutrinos in two complementary analyses, both analyzing four years of data (2013-2016) from the two deep stations (A2, A3) operating at that time. We place a 90 % CL upper limit on the diffuse all flavor neutrino flux at $10^{18}$ eV of $EF(E)=5.6\times10^{-16}$ $\textrm{cm}^{-2}$$\textrm{s}^{-1}$$\textrm{sr}^{-1}$. This analysis includes four times the exposure of the previous ARA result, and represents approximately 1/5 the exposure expected from operating ARA until the end of 2022.

preprint2020arXiv

Limitations imposed by optical turbulence profile structure and evolution on tomographic reconstruction for the ELT

The performance of tomographic adaptive optics systems is intrinsically linked to the vertical profile of optical turbulence. Firstly, a sufficient number of discrete turbulent layers must be reconstructed to model the true continuous turbulence profile. Secondly over the course of an observation, the profile as seen by the telescope changes and the tomographic reconstructor must be updated. These changes can be due to the unpredictable evolution of turbulent layers on meteorological timescales as short as minutes. Here we investigate the effect of changing atmospheric conditions on the quality of tomographic reconstruction by coupling fast analytical adaptive optics simulation to a large database of 10 691 high resolution turbulence profiles measured over two years by the Stereo-SCIDAR instrument at ESO Paranal, Chile. This work represents the first investigation of these effects with a large, statistically significant sample of turbulence profiles. The statistical nature of the study allows us to assess not only the degradation and variability in tomographic error with a set of system parameters (e.g. number of layers, temporal update period) but also the required parameters to meet some error threshold. In the most challenging conditions where the profile is rapidly changing, these parameters must be far more tightly constrained in order to meet this threshold. By providing estimates of these constraints for a wide range of system geometries as well as the impact of different temporal optimisation strategies we may assist the designers of tomographic AO for the ELT to dimension their systems.

preprint2019arXiv

AOtools -- a Python package for adaptive optics modelling and analysis

AOtools is a Python package which is open-source and aimed at providing tools for adaptive optics users and researchers. We present version 1.0 which contains tools for adaptive optics processing, including analysing data in the pupil plane, images and point spread functions in the focal plane, wavefront sensors, modelling of atmospheric turbulence, physical optical propagation of wavefronts, and conversion between frequently used adaptive optics and astronomical units. The main drivers behind AOtools is that it should be easy to install and use. To achieve this the project features extensive documentation, automated unit testing and is registered on the Python Package Index. AOtools is under continuous active development to expand the features available and we encourage everyone involved in adaptive optics to become involved and contribute to the project.

preprint2019arXiv

Characterisation of the ground layer of turbulence at Paranal using a robotic SLODAR system

We describe the implementation of a robotic SLODAR instrument at the Cerro Paranal observatory. The instrument measures the vertical profile of the optical atmospheric turbulence strength, in 8 resolution elements, to a maximum altitude ranging between 100 m and 500 m. We present statistical results of measurements of the turbulence profile on a total of 875 nights between 2014 and 2018. The vertical profile of the ground layer of turbulence is very varied, but in the median case most of the turbulence strength in the ground layer is concentrated within the first 50 m altitude, with relatively weak turbulence at higher altitudes up to 500 m. We find good agreement between measurements of the seeing angle from the SLODAR and from the Paranal DIMM seeing monitor, and also for seeing values extracted from the Shack-Hartmann active optics sensor of VLT UT1, adjusting for the height of each instrument above ground level. The SLODAR data suggest that a median improvement in the seeing angle from 0.689 arcsec to 0.481 arcsec at wavelength 500 nm would be obtained by fully correcting the ground-layer turbulence between the height of the UTs (taken as 10 m) and altitude 500 m.

preprint2016arXiv

A tomographic algorithm to determine tip-tilt information from laser guide stars

Laser Guide Stars (LGS) have greatly increased the sky-coverage of Adaptive Optics (AO) systems. Due to the up-link turbulence experienced by LGSs, a Natural Guide Star (NGS) is still required, preventing full sky-coverage. We present a method of obtaining partial tip-tilt information from LGSs alone in multi-LGS tomographic LGS AO systems. The method of LGS up-link tip-tilt determination is derived using a geometric approach, then an alteration to the Learn and Apply algorithm for tomographic AO is made to accommodate up-link tip-tilt. Simulation results are presented, verifying that the technique shows good performance in correcting high altitude tip-tilt, but not that from low altitudes. We suggest that the method is combined with multiple far off-axis tip-tilt NGSs to provide gains in performance and sky-coverage over current tomographic AO systems.

preprint2016arXiv

Electromagnetic effects on the light pseudoscalar mesons and determination of $m_u/m_d$

The MILC Collaboration has completed production running of electromagnetic effects on light mesons using asqtad improved staggered quarks. In these calculations, we use quenched photons in the noncompact formalism. We study four lattice spacings from $\approx\!0.12\:$fm to $\approx\!0.045\:$fm. To study finite-volume effects, we used six spatial lattice sizes $L/a=12$, 16, 20, 28, 40, and 48, at $a\!\approx\!0.12\:$fm. We update our preliminary values for the correction to Dashen's theorem ($ε$) and the quark-mass ratio $m_u/m_d$.

preprint2016arXiv

Experience with wavefront sensor and deformable mirror interfaces for wide-field adaptive optics systems

Recent advances in adaptive optics (AO) have led to the implementation of wide field-of-view AO systems. A number of wide-field AO systems are also planned for the forthcoming Extremely Large Telescopes. Such systems have multiple wavefront sensors of different types, and usually multiple deformable mirrors (DMs). Here, we report on our experience integrating cameras and DMs with the real-time control systems of two wide-field AO systems. These are CANARY, which has been operating on-sky since 2010, and DRAGON, which is a laboratory adaptive optics real-time demonstrator instrument. We detail the issues and difficulties that arose, along with the solutions we developed. We also provide recommendations for consideration when developing future wide-field AO systems.

preprint2016arXiv

HiPERCAM: A high-speed, quintuple-beam CCD camera for the study of rapid variability in the Universe

HiPERCAM is a high-speed camera for the study of rapid variability in the Universe. The project is funded by a 3.5MEuro European Research Council Advanced Grant. HiPERCAM builds on the success of our previous instrument, ULTRACAM, with very significant improvements in performance thanks to the use of the latest technologies. HiPERCAM will use 4 dichroic beamsplitters to image simultaneously in 5 optical channels covering the u'g'r'i'z' bands. Frame rates of over 1000 per second will be achievable using an ESO CCD controller (NGC), with every frame GPS timestamped. The detectors are custom-made, frame-transfer CCDs from e2v, with 4 low-noise (2.5e-) outputs, mounted in small thermoelectrically-cooled heads operated at 180 K, resulting in virtually no dark current. The two reddest CCDs will be deep-depletion devices with anti-etaloning, providing high quantum efficiencies across the red part of the spectrum with no fringing. The instrument will also incorporate scintillation noise correction via the conjugate-plane photometry technique. The opto-mechanical chassis will make use of additive manufacturing techniques in metal to make a light-weight, rigid and temperature-invariant structure. First light is expected on the 4.2m William Herschel Telescope on La Palma in 2017 (on which the field of view will be 10' with a 0.3"/pixel scale), with subsequent use planned on the 10.4m Gran Telescopio Canarias on La Palma (on which the field of view will be 4' with a 0.11"/pixel scale) and the 3.5m New Technology Telescope in Chile.

preprint2015arXiv

Atmospheric Scintillation in Astronomical Photometry

Scintillation noise due to the Earth's turbulent atmosphere can be a dominant noise source in high-precision astronomical photometry when observing bright targets from the ground. Here we describe the phenomenon of scintillation from its physical origins to its effect on photometry. We show that Young's (1967) scintillation-noise approximation used by many astronomers tends to underestimate the median scintillation noise at several major observatories around the world. We show that using median atmospheric optical turbulence profiles, which are now available for most sites, provides a better estimate of the expected scintillation noise and that real-time turbulence profiles can be used to precisely characterise the scintillation noise component of contemporaneous photometric measurements. This will enable a better understanding and calibration of photometric noise sources and the effectiveness of scintillation correction techniques. We also provide new equations for calculating scintillation noise, including for extremely large telescopes where the scintillation noise will actually be lower than previously thought. These equations highlight the fact that scintillation noise and shot noise have the same dependence on exposure time and so if an observation is scintillation limited, it will be scintillation limited for all exposure times. The ratio of scintillation noise to shot noise is also only weakly dependent on telescope diameter and so a bigger telescope may not yield a reduction in fractional scintillation noise.

preprint2015arXiv

Electromagnetic effects on the light hadron spectrum

For some time, the MILC Collaboration has been studying electromagnetic effects on light mesons. These calculations use fully dynamical QCD, but only quenched photons, which suffices to NLO in XPT. That is, the sea quarks are electrically neutral, while the valence quarks carry charge. For the photons we use the non-compact formalism. We have new results with lattice spacing as small as 0.045 fm and a large range of volumes. We consider how well chiral perturbation theory describes these results and the implications for light quark masses

preprint2014arXiv

Open-loop tomography with artificial neural networks on CANARY: on-sky results

We present recent results from the initial testing of an Artificial Neural Network (ANN) based tomographic reconstructor Complex Atmospheric Reconstructor based on Machine lEarNing (CARMEN) on Canary, an Adaptive Optics demonstrator operated on the 4.2m William Herschel Telescope, La Palma. The reconstructor was compared with contemporaneous data using the Learn and Apply (L&A) tomographic reconstructor. We find that the fully optimised L&A tomographic reconstructor outperforms CARMEN by approximately 5% in Strehl ratio or 15nm rms in wavefront error. We also present results for Canary in Ground Layer Adaptive Optics mode to show that the reconstructors are tomographic. The results are comparable and this small deficit is attributed to limitations in the training data used to build the ANN. Laboratory bench tests show that the ANN can out perform L&A under certain conditions, e.g. if the higher layer of a model two layer atmosphere was to change in altitude by ~300~m (equivalent to a shift of approximately one tenth of a subaperture).

preprint2013arXiv

Electromagnetic contributions to pseudoscalar masses

We report on the calculation by the MILC Collaboration of the electromagnetic effects on kaon and pion masses. These masses are computed in QCD with dynamical (asqtad staggered) quarks plus quenched photons at three lattice spacings varying from 0.12 to 0.06 fm. The masses are fit to staggered chiral perturbation theory with NLO electromagnetic terms, as well as analytic terms at higher order. We extrapolate the results to physical light-quark masses and to the continuum limit. At the current stage of the analysis, most, but not all, of the systematic errors have been estimated. The main goal is the comparison of kaon electromagnetic splittings to those of the pion, i.e., an evaluation of the corrections to "Dashen's theorem." This in turn will allow us to significantly reduce the systematic errors in our determination of m_u/m_d.

preprint2013arXiv

Lattice QCD ensembles with four flavors of highly improved staggered quarks

We present results from our simulations of quantum chromodynamics (QCD) with four flavors of quarks: u, d, s, and c. These simulations are performed with a one-loop Symanzik improved gauge action, and the highly improved staggered quark (HISQ) action. We are generating gauge configurations with four values of the lattice spacing ranging from 0.06 fm to 0.15 fm, and three values of the light quark mass, including the value for which the Goldstone pion mass is equal to the physical pion mass. We discuss simulation algorithms, scale setting, taste symmetry breaking, and the autocorrelations of various quantities. We also present results for the topological susceptibility which demonstrate the improvement of the HISQ configurations relative to those generated earlier with the asqtad improved staggered action.

preprint2013arXiv

Leptonic decay-constant ratio f_{K^+}/f_{pi^+} from lattice QCD with physical light quarks

A calculation of the ratio of leptonic decay constants f_{K^+}/f_{π^+} makes possible a precise determination of the ratio of CKM matrix elements |V_{us}|/|V_{ud}| in the Standard Model, and places a stringent constraint on the scale of new physics that would lead to deviations from unitarity in the first row of the CKM matrix. We compute f_{K^+}/f_{π^+} numerically in unquenched lattice QCD using gauge-field ensembles recently generated that include four flavors of dynamical quarks: up, down, strange, and charm. We analyze data at four lattice spacings a ~ 0.06, 0.09, 0.12, and 0.15 fm with simulated pion masses down to the physical value 135 MeV. We obtain f_{K^+}/f_{π^+} = 1.1947(26)(37), where the errors are statistical and total systematic, respectively. This is our first physics result from our N_f = 2+1+1 ensembles, and the first calculation of f_{K^+}/f_{π^+} from lattice-QCD simulations at the physical point. Our result is the most precise lattice-QCD determination of f_{K^+}/f_{π^+}, with an error comparable to the current world average. When combined with experimental measurements of the leptonic branching fractions, it leads to a precise determination of |V_{us}|/|V_{ud}| = 0.2309(9)(4) where the errors are theoretical and experimental, respectively.

preprint2013arXiv

Stereo-SCIDAR: Optical turbulence profiling with high sensitivity using a modified SCIDAR instrument

The next generation of adaptive optics (AO) systems will require tomographic reconstruction techniques to map the optical refractive index fluctuations, generated by the atmospheric turbulence, along the line of sight to the astronomical target. These systems can be enhanced with data from an external atmospheric profiler. This is important for Extremely Large Telescope scale tomography. Here we propose a new instrument which utilises the generalised SCIntillation Detection And Ranging (SCIDAR) technique to allow high sensitivity vertical profiles of the atmospheric optical turbulence and wind velocity profile above astronomical observatories. The new approach, which we refer to as 'Stereo-SCIDAR', uses a stereoscopic system with the scintillation pattern from each star of a double-star target incident on a separate detector. Separating the pupil images for each star has several advantages including: increased magnitude difference tolerance for the target stars; negating the need for re-calibration due to the normalisation errors usually associated with SCIDAR; an increase of at least a factor of two in the signal-to-noise ratio of the cross-covariance function and hence the profile for equal magnitude target stars and up to a factor of 16 improvement for targets of 3 magnitudes difference; and easier real-time reconstruction of the wind-velocity profile. Theoretical response functions are calculated for the instrument, and the performance is investigated using a Monte-Carlo simulation. The technique is demonstrated using data recorded at the 2.5 m Nordic Optical Telescope and the 1.0 m Jacobus Kapteyn Telescope, both on La Palma.

preprint2013arXiv

Update on the 2+1+1 flavor QCD equation of state with HISQ

We present recent results on the QCD equation of state with 2+1+1 flavors of highly improved staggered quarks (HISQ). We focus on three sets of ensembles with temporal extent 6, 8 and 10, that reach up to temperatures of 967, 725 and 580 MeV, respectively. The strange and charm quark masses are tuned to the physical values and the light quarks mass is set to one fifth of the strange. This corresponds to a Goldstone pion of about 300 MeV.

preprint2012arXiv

Status of the MILC calculation of electromagnetic contributions to pseudoscalar masses

We calculate pseudoscalar masses on gauge configurations containing the effects of 2+1 flavors of dynamical asqtad quarks and quenched electromagnetism. The lattice spacings vary from 0.12 to 0.06 fm. The masses are fit with staggered chiral perturbation theory including NLO electromagnetic terms. We attempt to extract the fit parameters for the electromagnetic contributions, while taking into account the finite volume effects, and extrapolate them to the physical limit.

preprint2011arXiv

Properties of light pseudoscalars from lattice QCD with HISQ ensembles

We fit lattice-QCD data for light-pseudoscalar masses and decay constants, from HISQ configurations generated by MILC, to SU(3) staggered chiral perturbation theory. At present such fits have rather high values of chi^2/d.o.f., possibly due to the lack of ensembles with lighter-than-physical sea strange-quark masses. We propose solutions to this problem for future work. We also perform simple linear interpolations near the physical point on two ensembles with different lattice spacings, and obtain the preliminary result (f_K / f_pi)^phys = 1.1872(41) in the continuum limit.

preprint2010arXiv

Results for light pseudoscalar mesons

We present the current status of the MILC collaboration's calculations of the properties of the light pseudoscalar meson sector. We use asqtad staggered ensembles with 2+1 dynamical flavors down to $a \approx 0.045$ fm and light quark mass down to 0.05 $m_s$. Here we describe fits to the data using chiral forms from SU(3) chiral perturbation theory, including all staggered taste violations at NLO and the continuum NNLO chiral logarithms. We emphasize issues of convergence of the chiral expansion.

preprint2010arXiv

Scaling studies of QCD with the dynamical HISQ action

We study the lattice spacing dependence, or scaling, of physical quantities using the highly improved staggered quark (HISQ) action introduced by the HPQCD/UKQCD collaboration, comparing our results to similar simulations with the asqtad fermion action. Results are based on calculations with lattice spacings approximately 0.15, 0.12 and 0.09 fm, using four flavors of dynamical HISQ quarks. The strange and charm quark masses are near their physical values, and the light-quark mass is set to 0.2 times the strange-quark mass. We look at the lattice spacing dependence of hadron masses, pseudoscalar meson decay constants, and the topological susceptibility. In addition to the commonly used determination of the lattice spacing through the static quark potential, we examine a determination proposed by the HPQCD collaboration that uses the decay constant of a fictitious "unmixed s bar s" pseudoscalar meson. We find that the lattice artifacts in the HISQ simulations are much smaller than those in the asqtad simulations at the same lattice spacings and quark masses.

preprint2010arXiv

Simulations with dynamical HISQ quarks

We report on the status of a program of generating and using configurations with four flavors of dynamical quarks, using the HISQ action. We study the lattice spacing dependence of physical quantities in these simulations, using runs at several lattice spacings, but with the light quark mass held fixed at two tenths of the strange quark mass. We find that the lattice artifacts in the HISQ simulations are much smaller than those in the asqtad simulations at the same lattice spacings and quark masses. We also discuss methods for setting the scale, or assigning a lattice spacing to ensembles run at unphysical parameters.

preprint2010arXiv

Topological susceptibility with the asqtad action

Chiral perturbation theory predicts that in quantum chromodynamics (QCD), light dynamical quarks suppress the gauge-field topological susceptibility of the vacuum. The degree of suppression depends on quark multiplicity and masses. It provides a strong consistency test for fermion formulations in lattice QCD. Such tests are especially important for staggered fermion formulations that lack a full chiral symmetry and use the "fourth-root" procedure to achieve the desired number of sea quarks. Over the past few years we have measured the topological susceptibility on a large database of 18 gauge field ensembles, generated in the presence of 2+1 flavors of dynamical asqtad quarks with up and down quark masses ranging from 0.05 to 1 in units of the strange quark mass and lattice spacings ranging from 0.045 fm to 0.12 fm. Our study also includes three quenched ensembles with lattice spacings ranging from 0.06 to 0.12 fm. We construct the topological susceptibility from the integrated point-to-point correlator of the discretized topological charge density F-Fdual. To reduce its variance, we model the asymptotic tail of the correlator. The continuum extrapolation of our results for the topological susceptibility agrees nicely at small quark mass with the predictions of lowest-order SU(3) chiral perturbation theory, thus lending support to the validity of the fourth-root procedure.

preprint2009arXiv

Chebyshev type lattice path weight polynomials by a constant term method

We prove a constant term theorem which is useful for finding weight polynomials for Ballot/Motzkin paths in a strip with a fixed number of arbitrary `decorated' weights as well as an arbitrary `background' weight. Our CT theorem, like Viennot's lattice path theorem from which it is derived primarily by a change of variable lemma, is expressed in terms of orthogonal polynomials which in our applications of interest often turn out to be non-classical. Hence we also present an efficient method for finding explicit closed form polynomial expressions for these non-classical orthogonal polynomials. Our method for finding the closed form polynomial expressions relies on simple combinatorial manipulations of Viennot's diagrammatic representation for orthogonal polynomials. In the course of the paper we also provide a new proof of Viennot's original orthogonal polynomial lattice path theorem. The new proof is of interest because it uses diagonalization of the transfer matrix, but gets around difficulties that have arisen in past attempts to use this approach. In particular we show how to sum over a set of implicitly defined zeros of a given orthogonal polynomial, either by using properties of residues or by using partial fractions. We conclude by applying the method to two lattice path problems important in the study of polymer physics as models of steric stabilization and sensitized flocculation.

preprint2009arXiv

HISQ action in dynamical simulations

We report on recent progress in employing the Highly Improved Staggered Quark (HISQ) action introduced by the HPQCD/UKQCD collaboration in simulations with dynamical fermions. The HISQ action is an order $a^2$ Symanzik-improved action with further suppressed taste symmetry violations. The improvement in taste symmetry is achieved by introducing Fat7 smearing of the original gauge links and reunitarization (projection to an element of U(3) or SU(3)) followed by Asq-type smearing. Major challenges for calculating the fermion force are related to the reunitarization step. We present a preliminary study of the HISQ action on two 2+1+1 flavor ensembles with the lattice spacing roughly equivalent to the MILC asqtad a=0.125 and 0.09 fm ensembles.

preprint2009arXiv

MILC results for light pseudoscalars

We present the latest preliminary results of the MILC collaboration's analysis of the light pseudoscalar meson sector. The analysis includes data from new ensembles with smaller lattice spacings, smaller light quark masses and lighter-than-physical strange quark masses. Both SU(2) and SU(3) chiral fits, including NNLO chiral logarithms, are shown. We give results for decay constants, quark masses, Gasser-Leutwyler low energy constants, and condensates in the two- and three-flavor chiral limits.

preprint2009arXiv

Progress on four flavor QCD with the HISQ action

We describe recent progress on generation of gauge configurations using the Highly Improved Staggered Quark (HISQ) action that was designed by the HPQCD/UKQCD collaboration. The HISQ action requires two levels of smearing with a reunitarization of the links before the second smearing. We describe how we deal with the occurrence of occasional large forces arising from the reunitarization step. The MILC collaboration is currently generating ensembles with approximate lattice spacings of 0.15, 0.12, 0.09, and 0.06 fm, with the strange and charm quark masses close to their physical values and the mass of the light quarks m_l set to 0.2 m_s. We present recent results for pion taste splittings, light hadron masses, the static potential, the eta_c dispersion relation and the topological susceptibility.

preprint2009arXiv

Results from the MILC collaboration's SU(3) chiral perturbation theory analysis

We present the status of the MILC collaboration's analysis of the light pseudoscalar meson sector with SU(3) chiral fits. The analysis includes data from new ensembles with smaller lattice spacing, smaller light quark masses and lighter than physical strange quark masses. Our fits include the NNLO chiral logarithms. We present results for decay constants, quark masses, Gasser-Leutwyler low energy constants, and condensates in the two- and three-flavor chiral limits.

preprint2009arXiv

SU(2) chiral fits to light pseudoscalar masses and decay constants

We present the results of fits to recent asqtad data in the light pseudoscalar sector using SU(2) partially-quenched staggered chiral perturbation theory. Superfine (a~0.06 fm) and ultrafine (a~0.045 fm) ensembles are used, where light sea quark masses and taste splittings are small compared to the strange quark mass. Our fits include continuum NNLO chiral logarithms and analytic terms. We give preliminary results for the pion decay constant, SU(2) low-energy constants and the chiral condensate in the two-flavor chiral limit.

preprint2008arXiv

Electromagnetic splittings of hadrons from improved staggered quarks in full QCD

We present our initial study of the electromagnetic splittings of charged and neutral mesons, and the violation of Dashen's theorem. Hadron masses are calculated on MILC N_f=2+1 QCD ensembles at lattice spacing \approx 0.15fm, together with quenched non-compact U(1) configurations. The O(a^2) tadpole improved staggered quark (asqtad) action is used both for the sea quarks and for six different valence quark masses. Chiral extrapolations are performed using partially quenched chiral perturbation theory incorporating electromagnetic corrections.

preprint2007arXiv

Baryon masses with improved staggered quarks

The MILC collaboration's simulations with improved staggered quarks are being extended with runs at a lattice spacing of 0.06 fm with quark masses down to one tenth the strange quark mass. We give a brief introduction to these new simulations and the determination of the lattice spacing. Then we combine these new runs with older results to study the masses of the nucleon and the Omega minus in the continuum and chiral limits.

preprint2007arXiv

Status of the MILC light pseudoscalar meson project

We discuss the current status of our calculation of the physics of pi and K mesons using three dynamical flavors of improved staggered quarks. This year, we have a new ensemble with a lattice spacing of 0.06 fm and a light sea mass of 0.2 m_s, as well as significant increases in statistics at several coarser lattice spacings and/or heavier sea masses. Results for decay constants, quark masses, low energy constants, condensates, and V_{us} are presented.

preprint2007arXiv

The 2+1 flavor topological susceptibility from the asqtad action at 0.06 fm

We report new data for the topological susceptibility computed on 2+1 flavor dynamical configurations with lattice spacing 0.06 fm, generated with the asqtad action. The topological susceptibility is computed by HYP smearing and compared with rooted staggered chiral perturbation theory as the pion mass goes to zero. At 0.06 fm, the raw data is already quite close to the continuum extrapolated values obtained from coarser lattices. These results provide a further test of the asqtad action with rooted staggered flavors.

preprint2006arXiv

Update on the physics of light pseudoscalar mesons

We present an update of the MILC investigation of the properties of light pseudoscalar mesons using three flavors of improved staggered quarks. Results are presented for the $π$ and $K$ leptonic decay constants, the CKM matrix element $V_{us}$, the up, down and strange quark masses, and the coefficients of the $O(p^4)$ chiral lagrangian. We have new data for lattice spacing $a \approx 0.15$ fm with several values of the light quark mass down to one-tenth the strange quark mass, higher statistics for $a \approx 0.09$ fm with the light quark mass equal to one-tenth the strange quark mass, and initial results for our smallest lattice spacing, $a \approx 0.06$ fm with light quark mass two-fifths of the strange quark mass.

preprint2005arXiv

Update on pi and K Physics

We present an update of the MILC studies of the physics of light pseudoscalars using improved staggered fermions. New runs at lighter quark mass, as well as increased statistics for older sets, are enabling us to improve the results for decay constants in full QCD. In addition, we have analyzed quenched runs at two different lattice spacings. This makes possible a test of the applicability of staggered chiral perturbation theory in a different context.

preprint2004arXiv

First determination of the strange and light quark masses from full lattice QCD

We compute the strange quark mass $m_s$ and the average of the $u$ and $d$ quark masses $\hat m$ using full lattice QCD with three dynamical quarks combined with experimental values for the pion and kaon masses. The simulations have degenerate $u$ and $d$ quarks with masses $m_u=m_d\equiv \hat m$ as low as $m_s/8$, and two different values of the lattice spacing. The bare lattice quark masses obtained are converted to the $\msbar$ scheme using perturbation theory at $O(alpha_s)$. Our results are: $m_s^\msbar$(2 GeV) = 76(0)(3)(7)(0) MeV, $\hat m^\msbar$(2 GeV) = 2.8(0)(1)(3)(0) MeV and $m_s/\hat m$ = 27.4(1)(4)(0)(1), where the errors are from statistics, simulation, perturbation theory, and electromagnetic effects, respectively.

preprint2004arXiv

Leptonic decay constants f_Ds and f_D in three flavor lattice QCD

We determine the leptonic decay constants in three flavor unquenched lattice QCD. We use O(a^2)-improved staggered light quarks and O(a)-improved charm quarks in the Fermilab heavy quark formalism. Our preliminary results, based upon an analysis at a single lattice spacing, are f_Ds = 263(+5-9)(+/-24) MeV and f_D = 225(+11-13)(+/-21) MeV. In each case, the first reported error is statistical while the is the combined systematic uncertainty.

preprint2004arXiv

The scaling dimension of low lying Dirac eigenmodes and of the topological charge density

As a quantitative measure of localization, the inverse participation ratio of low lying Dirac eigenmodes and topological charge density is calculated on quenched lattices over a wide range of lattice spacings and volumes. Since different topological objects (instantons, vortices, monopoles, and artifacts) have different co-dimension, scaling analysis provides information on the amount of each present and their correlation with the localization of low lying eigenmodes.