Generation of initial molecular dynamics configurations in arbitrary geometries and in parallel

A computational pre-processing tool for generating initial configurations of molecules for molecular dynamics simulations in geometries described by a mesh of unstructured arbitrary polyhedra is described. The mesh is divided into separate zones and each can be filled with a single crystal lattice of atoms. Each zone is filled by creating an expanding cube of crystal unit cells, initiated from an anchor point for the lattice. Each unit cell places the appropriate atoms for the user-specified crystal structure and orientation. The cube expands until the entire zone is filled with the lattice; zones with concave and disconnected volumes may be filled. When the mesh is spatially decomposed into portions for distributed parallel processing, each portion may be filled independently, meaning that the entire molecular system never needs to fit onto a single processor, allowing very large systems to be created. The computational time required to fill a zone with molecules scales linearly with the number of cells in the zone for a fixed number of molecules, and better than linearly with the number of molecules for a fixed number of mesh cells. Our tool, molConfig, has been implemented in the open source C++ code OpenFOAM

mdFoam+: Advanced molecular dynamics in OpenFOAM

This paper introduces mdFoam+, which is an MPI parallelised molecular dynamics (MD) solver implemented entirely within the OpenFOAM software framework. It is open-source and released under the same GNU General Public License (GPL) as OpenFOAM. The source code is released as a publicly open software repository that includes detailed documentation and tutorial cases. Since mdFoam+ is designed entirely within the OpenFOAM C++ object-oriented framework, it inherits a number of key features. The code is designed for extensibility and flexibility, so it is aimed first and foremost as an MD research tool, in which new models and test cases can be developed and tested rapidly. Implementing mdFoam+ in OpenFOAM also enables easier development of hybrid methods that couple MD with continuum-based solvers. Setting up MD cases follows the standard OpenFOAM format, as mdFoam+ also relies upon the OpenFOAM dictionary-based directory structure. This ensures that useful pre- and post-processing capabilities provided by OpenFOAM remain available even though the fully Lagrangian nature of an MD simulation is not typical of most OpenFOAM applications. Results show that mdFoam+ compares well to another well-known MD code (e.g. LAMMPS) in terms of benchmark problems, although it also has additional functionality that does not exist in other open-source MD codes

dsmcFoam+: An OpenFOAM based direct simulation Monte Carlo solver

dsmcFoam+ is a direct simulation Monte Carlo (DSMC) solver for rarefied gas dynamics, implemented within the OpenFOAM software framework, and parallelised with MPI. It is open-source and released under the GNU General Public License in a publicly available software repository that includes detailed documentation and tutorial DSMC gas flow cases. This release of the code includes many features not found in standard dsmcFoam, such as molecular vibrational and electronic energy modes, chemical reactions, and subsonic pressure boundary conditions. Since dsmcFoam+ is designed entirely within OpenFOAM’s C++ object-oriented framework, it benefits from a number of key features: the code emphasises extensibility and flexibility so it is aimed first and foremost as a research tool for DSMC, allowing new models and test cases to be developed and tested rapidly. All DSMC cases are as straightforward as setting up any standard OpenFOAM case, as dsmcFoam+ relies upon the standard OpenFOAM dictionary based directory structure. This ensures that useful pre- and post-processing capabilities provided by OpenFOAM remain available even though the fully Lagrangian nature of a DSMC simulation is not typical of most OpenFOAM applications. We show that dsmcFoam+ compares well to other well-known DSMC codes and to analytical solutions in terms of benchmark results

Factors impacting antiretroviral therapy adherence among human immunodeficiency virus-positive adolescents in Sub-Saharan Africa: a systematic review

© 2018 The Royal Society for Public Health. Published by Elsevier Ltd. All rights reserved.Objectives: Eighty-two percent of human immunodeficiency virus (HIV)–positive adolescents live in Sub-Saharan Africa (SSA). Despite the availability of antiretroviral therapy (ART), adherence levels are suboptimal, leading to poor outcomes. This systematic review investigated factors impacting ART adherence among adolescents in SSA, including religious beliefs and intimate relationships. Methods: A systematic review was conducted between June and August 2016 using eight electronic databases, including Cochrane and PubMed. Published, ongoing and unpublished research, conducted in SSA from 2004 to 2016, was identified and thematic analysis was used to summarise findings. Results: Eleven studies from eight SSA countries, published in English between 2011 and 2016, reported on factors impacting ART adherence among adolescents living with HIV (ALHIV). Forty-four barriers and 29 facilitators to adherence were identified, representing a complex web of factors. The main barriers were stigma, ART side-effects, lack of assistance and forgetfulness. Facilitators included caregiver support, peer support groups and knowledge of HIV status. Conclusions: Stigma reflects difficult relations between ALHIV and their HIV-negative peers and adults. Most interventions target only those with HIV, suggesting a policy shift towards the wider community could be beneficial. Recommendations include engaging religious leaders and schools to change negative societal attitudes. Limitations of the review include the urban settings and recruitment of predominantly vertically infected participants in most included studies. Therefore, the findings cannot be extrapolated to ALHIV residing in rural locations or horizontally infected ALHIV, highlighting the need for further research in those areas.Peer reviewedFinal Accepted Versio

Search for R-Parity Violating Decays of Scalar Fermions at LEP

A search for pair-produced scalar fermions under the assumption that R-parity is not conserved has been performed using data collected with the OPAL detector at LEP. The data samples analysed correspond to an integrated luminosity of about 610 pb-1 collected at centre-of-mass energies of sqrt(s) 189-209 GeV. An important consequence of R-parity violation is that the lightest supersymmetric particle is expected to be unstable. Searches of R-parity violating decays of charged sleptons, sneutrinos and squarks have been performed under the assumptions that the lightest supersymmetric particle decays promptly and that only one of the R-parity violating couplings is dominant for each of the decay modes considered. Such processes would yield final states consisting of leptons, jets, or both with or without missing energy. No significant single-like excess of events has been observed with respect to the Standard Model expectations. Limits on the production cross- section of scalar fermions in R-parity violating scenarios are obtained. Constraints on the supersymmetric particle masses are also presented in an R-parity violating framework analogous to the Constrained Minimal Supersymmetric Standard Model.Comment: 51 pages, 24 figures, Submitted to Eur. Phys. J.

Measurement of the Hadronic Photon Structure Function F_2^gamma at LEP2

The hadronic structure function of the photon F_2^gamma is measured as a function of Bjorken x and of the factorisation scale Q^2 using data taken by the OPAL detector at LEP. Previous OPAL measurements of the x dependence of F_2^gamma are extended to an average Q^2 of 767 GeV^2. The Q^2 evolution of F_2^gamma is studied for average Q^2 between 11.9 and 1051 GeV^2. As predicted by QCD, the data show positive scaling violations in F_2^gamma. Several parameterisations of F_2^gamma are in agreement with the measurements whereas the quark-parton model prediction fails to describe the data.Comment: 4 pages, 2 figures, to appear in the proceedings of Photon 2001, Ascona, Switzerlan

A measurement of the tau mass and the first CPT test with tau leptons

We measure the mass of the tau lepton to be 1775.1+-1.6(stat)+-1.0(syst.) MeV using tau pairs from Z0 decays. To test CPT invariance we compare the masses of the positively and negatively charged tau leptons. The relative mass difference is found to be smaller than 3.0 10^-3 at the 90% confidence level.Comment: 10 pages, 4 figures, Submitted to Phys. Letts.

Measurement of the B0 Lifetime and Oscillation Frequency using B0->D*+l-v decays

The lifetime and oscillation frequency of the B0 meson has been measured using B0->D*+l-v decays recorded on the Z0 peak with the OPAL detector at LEP. The D*+ -> D0pi+ decays were reconstructed using an inclusive technique and the production flavour of the B0 mesons was determined using a combination of tags from the rest of the event. The results t_B0 = 1.541 +- 0.028 +- 0.023 ps, Dm_d = 0.497 +- 0.024 +- 0.025 ps-1 were obtained, where in each case the first error is statistical and the second systematic.Comment: 17 pages, 4 figures, submitted to Phys. Lett.

First Measurement of Z/gamma* Production in Compton Scattering of Quasi-real Photons

We report the first observation of Z/gamma* production in Compton scattering of quasi-real photons. This is a subprocess of the reaction e+e- to e+e-Z/gamma*, where one of the final state electrons is undetected. Approximately 55 pb-1 of data collected in the year 1997 at an e+e- centre-of-mass energy of 183 GeV with the OPAL detector at LEP have been analysed. The Z/gamma* from Compton scattering has been detected in the hadronic decay channel. Within well defined kinematic bounds, we measure the product of cross-section and Z/gamma* branching ratio to hadrons to be (0.9+-0.3+-0.1) pb for events with a hadronic mass larger than 60 GeV, dominated by (e)eZ production. In the hadronic mass region between 5 GeV and 60 GeV, dominated by (e)egamma* production, this product is found to be (4.1+-1.6+-0.6) pb. Our results agree with the predictions of two Monte Carlo event generators, grc4f and PYTHIA.Comment: 18 pages, LaTeX, 5 eps figures included, submitted to Physics Letters