LF-PPL: A Low-Level First Order Probabilistic Programming Language for Non-Differentiable Models

We develop a new Low-level, First-order Probabilistic Programming Language (LF-PPL) suited for models containing a mix of continuous, discrete, and/or piecewise-continuous variables. The key success of this language and its compilation scheme is in its ability to automatically distinguish parameters the density function is discontinuous with respect to, while further providing runtime checks for boundary crossings. This enables the introduction of new inference engines that are able to exploit gradient information, while remaining efficient for models which are not everywhere differentiable. We demonstrate this ability by incorporating a discontinuous Hamiltonian Monte Carlo (DHMC) inference engine that is able to deliver automated and efficient inference for non-differentiable models. Our system is backed up by a mathematical formalism that ensures that any model expressed in this language has a density with measure zero discontinuities to maintain the validity of the inference engine.Comment: Published in the proceedings of the 22nd International Conference on Artificial Intelligence and Statistics (AISTATS

Density functional theory using an optimized exchange-correlation potential

We have performed self-consistent calculations for first and second row atoms using a variant of density-functional theory, the optimized effective potential method, with an approximation due to Krieger, Li and Iafrate and a correlation-energy functional developed by Colle and Salvetti. The mean absolute deviation of first-row atomic ground-state energies from the exact non-relativistic values is 4.7 mH in our scheme, as compared to 4.5 mH in a recent configuration-interaction calculation. The proposed scheme is significantly more accurate than the conventional Kohn-Sham method while the numerical effort involved is about the same as for an ordinary Hartree-Fock calculation.Comment: To be published in Chemical Physics Letters (1995), latex, 15 pages, no figure

Extrusion-Printing of Multi-Channeled Two-Component Hydrogel Constructs from Gelatinous Peptides and Anhydride-Containing Oligomers

The performance of artificial nerve guidance conduits (NGC) in peripheral nerve regeneration can be improved by providing structures with multiple small channels instead of a single wide lumen. 3D-printing is a strategy to access such multi-channeled structures in a defined and reproducible way. This study explores extrusion-based 3D-printing of two-component hydrogels from a single cartridge printhead into multi-channeled structures under aseptic conditions. The gels are based on a platform of synthetic, anhydride-containing oligomers for cross-linking of gelatinous peptides. Stable constructs with continuous small channels and a variety of footprints and sizes were successfully generated from formulations containing either an organic or inorganic gelation base. The adjustability of the system was investigated by varying the cross-linking oligomer and substituting the gelation bases controlling the cross-linking kinetics. Formulations with organic N-methyl-piperidin-3-ol and inorganic K2HPO4 yielded hydrogels with comparable properties after manual processing and extrusion-based 3D-printing. The slower reaction kinetics of formulations with K2HPO4 can be beneficial for extending the time frame for printing. The two-component hydrogels displayed both slow hydrolytic and activity-dependent enzymatic degradability. Together with satisfying in vitro cell proliferation data, these results indicate the suitability of our cross-linked hydrogels as multi-channeled NGC for enhanced peripheral nerve regeneration

Inferring transmission fitness advantage of SARS-CoV-2 variants of concern from wastewater samples using digital PCR, Switzerland, December 2020 through March 2021

BackgroundThroughout the COVID-19 pandemic, SARS-CoV-2 genetic variants of concern (VOCs) have repeatedly and independently arisen. VOCs are characterised by increased transmissibility, increased virulence or reduced neutralisation by antibodies obtained from prior infection or vaccination. Tracking the introduction and transmission of VOCs relies on sequencing, typically whole genome sequencing of clinical samples. Wastewater surveillance is increasingly used to track the introduction and spread of SARS-CoV-2 variants through sequencing approaches.AimHere, we adapt and apply a rapid, high-throughput method for detection and quantification of the relative frequency of two deletions characteristic of the Alpha, Beta, and Gamma VOCs in wastewater.MethodsWe developed drop-off RT-dPCR assays and an associated statistical approach implemented in the R package WWdPCR to analyse temporal dynamics of SARS-CoV-2 signature mutations (spike Δ69-70 and ORF1a Δ3675-3677) in wastewater and quantify transmission fitness advantage of the Alpha VOC.ResultsBased on analysis of Zurich wastewater samples, the estimated transmission fitness advantage of SARS-CoV-2 Alpha based on the spike Δ69-70 was 0.34 (95% confidence interval (CI): 0.30-0.39) and based on ORF1a Δ3675-3677 was 0.53 (95% CI: 0.49-0.57), aligning with the transmission fitness advantage of Alpha estimated by clinical sample sequencing in the surrounding canton of 0.49 (95% CI: 0.38-0.61).ConclusionDigital PCR assays targeting signature mutations in wastewater offer near real-time monitoring of SARS-CoV-2 VOCs and potentially earlier detection and inference on transmission fitness advantage than clinical sequencing. Keywords: B.1.1.7; SARS-CoV-2; digital PCR; drop-off assays; transmission fitness

Libxc: a library of exchange and correlation functionals for density functional theory

The central quantity of density functional theory is the so-called exchange-correlation functional. This quantity encompasses all non-trivial many-body effects of the ground-state and has to be approximated in any practical application of the theory. For the past 50 years, hundreds of such approximations have appeared, with many successfully persisting in the electronic structure community and literature. Here, we present a library that contains routines to evaluate many of these functionals (around 180) and their derivatives.Comment: 15 page

Measurement of the inclusive and dijet cross-sections of b-jets in pp collisions at sqrt(s) = 7 TeV with the ATLAS detector

The inclusive and dijet production cross-sections have been measured for jets containing b-hadrons (b-jets) in proton-proton collisions at a centre-of-mass energy of sqrt(s) = 7 TeV, using the ATLAS detector at the LHC. The measurements use data corresponding to an integrated luminosity of 34 pb^-1. The b-jets are identified using either a lifetime-based method, where secondary decay vertices of b-hadrons in jets are reconstructed using information from the tracking detectors, or a muon-based method where the presence of a muon is used to identify semileptonic decays of b-hadrons inside jets. The inclusive b-jet cross-section is measured as a function of transverse momentum in the range 20 < pT < 400 GeV and rapidity in the range |y| < 2.1. The bbbar-dijet cross-section is measured as a function of the dijet invariant mass in the range 110 < m_jj < 760 GeV, the azimuthal angle difference between the two jets and the angular variable chi in two dijet mass regions. The results are compared with next-to-leading-order QCD predictions. Good agreement is observed between the measured cross-sections and the predictions obtained using POWHEG + Pythia. MC@NLO + Herwig shows good agreement with the measured bbbar-dijet cross-section. However, it does not reproduce the measured inclusive cross-section well, particularly for central b-jets with large transverse momenta.Comment: 10 pages plus author list (21 pages total), 8 figures, 1 table, final version published in European Physical Journal

Jet energy measurement with the ATLAS detector in proton-proton collisions at root s=7 TeV

The jet energy scale and its systematic uncertainty are determined for jets measured with the ATLAS detector at the LHC in proton-proton collision data at a centre-of-mass energy of √s = 7TeV corresponding to an integrated luminosity of 38 pb-1. Jets are reconstructed with the anti-kt algorithm with distance parameters R=0. 4 or R=0. 6. Jet energy and angle corrections are determined from Monte Carlo simulations to calibrate jets with transverse momenta pT≥20 GeV and pseudorapidities {pipe}η{pipe}<4. 5. The jet energy systematic uncertainty is estimated using the single isolated hadron response measured in situ and in test-beams, exploiting the transverse momentum balance between central and forward jets in events with dijet topologies and studying systematic variations in Monte Carlo simulations. The jet energy uncertainty is less than 2. 5 % in the central calorimeter region ({pipe}η{pipe}<0. 8) for jets with 60≤pT<800 GeV, and is maximally 14 % for pT<30 GeV in the most forward region 3. 2≤{pipe}η{pipe}<4. 5. The jet energy is validated for jet transverse momenta up to 1 TeV to the level of a few percent using several in situ techniques by comparing a well-known reference such as the recoiling photon pT, the sum of the transverse momenta of tracks associated to the jet, or a system of low-pT jets recoiling against a high-pT jet. More sophisticated jet calibration schemes are presented based on calorimeter cell energy density weighting or hadronic properties of jets, aiming for an improved jet energy resolution and a reduced flavour dependence of the jet response. The systematic uncertainty of the jet energy determined from a combination of in situ techniques is consistent with the one derived from single hadron response measurements over a wide kinematic range. The nominal corrections and uncertainties are derived for isolated jets in an inclusive sample of high-pT jets. Special cases such as event topologies with close-by jets, or selections of samples with an enhanced content of jets originating from light quarks, heavy quarks or gluons are also discussed and the corresponding uncertainties are determined. © 2013 CERN for the benefit of the ATLAS collaboration

Single hadron response measurement and calorimeter jet energy scale uncertainty with the ATLAS detector at the LHC

The uncertainty on the calorimeter energy response to jets of particles is derived for the ATLAS experiment at the Large Hadron Collider (LHC). First, the calorimeter response to single isolated charged hadrons is measured and compared to the Monte Carlo simulation using proton-proton collisions at centre-of-mass energies of sqrt(s) = 900 GeV and 7 TeV collected during 2009 and 2010. Then, using the decay of K_s and Lambda particles, the calorimeter response to specific types of particles (positively and negatively charged pions, protons, and anti-protons) is measured and compared to the Monte Carlo predictions. Finally, the jet energy scale uncertainty is determined by propagating the response uncertainty for single charged and neutral particles to jets. The response uncertainty is 2-5% for central isolated hadrons and 1-3% for the final calorimeter jet energy scale.Comment: 24 pages plus author list (36 pages total), 23 figures, 1 table, submitted to European Physical Journal

Measurement of the polarisation of W bosons produced with large transverse momentum in pp collisions at sqrt(s) = 7 TeV with the ATLAS experiment

This paper describes an analysis of the angular distribution of W->enu and W->munu decays, using data from pp collisions at sqrt(s) = 7 TeV recorded with the ATLAS detector at the LHC in 2010, corresponding to an integrated luminosity of about 35 pb^-1. Using the decay lepton transverse momentum and the missing transverse energy, the W decay angular distribution projected onto the transverse plane is obtained and analysed in terms of helicity fractions f0, fL and fR over two ranges of W transverse momentum (ptw): 35 < ptw < 50 GeV and ptw > 50 GeV. Good agreement is found with theoretical predictions. For ptw > 50 GeV, the values of f0 and fL-fR, averaged over charge and lepton flavour, are measured to be : f0 = 0.127 +/- 0.030 +/- 0.108 and fL-fR = 0.252 +/- 0.017 +/- 0.030, where the first uncertainties are statistical, and the second include all systematic effects.Comment: 19 pages plus author list (34 pages total), 9 figures, 11 tables, revised author list, matches European Journal of Physics C versio