Showing 1–38 of 38 results for author: Thomas, D

Multi-axis inertial sensing with 2D arrays of matter waves

Authors: K. Stolzenberg, C. Struckmann, S. Bode, R. Li, A. Herbst, V. Vollenkemper, D. Thomas, E. M. Rasel, N. Gaaloul, D. Schlippert

Abstract: Atom interferometers are an exquisite measurement tool for inertial forces. However, they are commonly limited to one single sensitive axis, allowing high-precision multi-dimensional sensing only through subsequent or postcorrected measurements. Here, we introduce a novel 2D-array-arrangement of Bose-Einstein Condensates (BEC) initialized utilizing time-averaged optical potentials for simultaneous… ▽ More Atom interferometers are an exquisite measurement tool for inertial forces. However, they are commonly limited to one single sensitive axis, allowing high-precision multi-dimensional sensing only through subsequent or postcorrected measurements. Here, we introduce a novel 2D-array-arrangement of Bose-Einstein Condensates (BEC) initialized utilizing time-averaged optical potentials for simultaneous multi-axis inertial sensing. Deploying a 3 x 3 BEC array covering 1.6 mm^2, we perform measurements of angular velocity and acceleration of a rotating reference mirror, as well as a linear acceleration, e.g., induced by gravity, gradients, and higher order derivatives. We anticipate increased sensitivity of our method in interferometers with large scale factors in long-baseline or satellite atom interferometry. Our work paves the way for simple high-precision multi-axis inertial sensing and we envision further applications, e.g., for three-dimensional wave front characterization. △ Less

Submitted 19 July, 2024; v1 submitted 13 March, 2024; originally announced March 2024.

A 2D antiscatter grid and scatter sampling based CBCT method for online dose calculations during CBCT guided radiation therapy of pelvis

Authors: Farhang Bayat, Brian Miller, Yeonok Park, Zhelin Yu, Timur Alexeev, David Thomas, Kelly Stuhr, Brian Kavanagh, Moyed Miften, Cem Altunbas

Abstract: Online dose calculations before radiation treatment have applications in dose delivery verification, plan adaptation, and treatment planning. We propose a novel CBCT imaging pipeline to enhance accuracy. Our approach aims to improve HU accuracy in CBCT images for more precise dose calculations. A quantitative CBCT pipeline was implemented, combining data correction strategies and scatter rejection… ▽ More Online dose calculations before radiation treatment have applications in dose delivery verification, plan adaptation, and treatment planning. We propose a novel CBCT imaging pipeline to enhance accuracy. Our approach aims to improve HU accuracy in CBCT images for more precise dose calculations. A quantitative CBCT pipeline was implemented, combining data correction strategies and scatter rejection, achieving high CT number accuracy. We evaluated the pipeline's effect using pelvis anatomy phantoms and found that dosimetric errors in quantitative CBCT-based dose calculations were minimal. In contrast, clinical CBCT and high-performance ASG CBCT-based plans showed significant errors. The proposed quantitative CBCT pipeline offers comparable dose calculation accuracy to the gold-standard planning CT, eliminating the need for density overrides and enabling precise dose delivery monitoring or online plan adaptations in radiation therapy. △ Less

Submitted 10 October, 2023; originally announced October 2023.

ATHENA Detector Proposal -- A Totally Hermetic Electron Nucleus Apparatus proposed for IP6 at the Electron-Ion Collider

Authors: ATHENA Collaboration, J. Adam, L. Adamczyk, N. Agrawal, C. Aidala, W. Akers, M. Alekseev, M. M. Allen, F. Ameli, A. Angerami, P. Antonioli, N. J. Apadula, A. Aprahamian, W. Armstrong, M. Arratia, J. R. Arrington, A. Asaturyan, E. C. Aschenauer, K. Augsten, S. Aune, K. Bailey, C. Baldanza, M. Bansal, F. Barbosa, L. Barion , et al. (415 additional authors not shown)

Abstract: ATHENA has been designed as a general purpose detector capable of delivering the full scientific scope of the Electron-Ion Collider. Careful technology choices provide fine tracking and momentum resolution, high performance electromagnetic and hadronic calorimetry, hadron identification over a wide kinematic range, and near-complete hermeticity. This article describes the detector design and its e… ▽ More ATHENA has been designed as a general purpose detector capable of delivering the full scientific scope of the Electron-Ion Collider. Careful technology choices provide fine tracking and momentum resolution, high performance electromagnetic and hadronic calorimetry, hadron identification over a wide kinematic range, and near-complete hermeticity. This article describes the detector design and its expected performance in the most relevant physics channels. It includes an evaluation of detector technology choices, the technical challenges to realizing the detector and the R&D required to meet those challenges. △ Less

Submitted 13 October, 2022; originally announced October 2022.

Journal ref: JINST 17 (2022) 10, P10019

Lessons Learned from the Two Largest Galaxy Morphological Classification Catalogues built by Convolutional Neural Networks

Authors: Ting-Yun Cheng, H. Domínguez Sánchez, J. Vega-Ferrero, C. J. Conselice, M. Siudek, A. Aragón-Salamanca, M. Bernardi, R. Cooke, L. Ferreira, M. Huertas-Company, J. Krywult, A. Palmese, A. Pieres, A. A. Plazas Malagón, A. Carnero Rosell, D. Gruen, D. Thomas, D. Bacon, D. Brooks, D. J. James, D. L. Hollowood, D. Friedel, E. Suchyta, E. Sanchez, F. Menanteau , et al. (32 additional authors not shown)

Abstract: We compare the two largest galaxy morphology catalogues, which separate early and late type galaxies at intermediate redshift. The two catalogues were built by applying supervised deep learning (convolutional neural networks, CNNs) to the Dark Energy Survey data down to a magnitude limit of $\sim$21 mag. The methodologies used for the construction of the catalogues include differences such as the… ▽ More We compare the two largest galaxy morphology catalogues, which separate early and late type galaxies at intermediate redshift. The two catalogues were built by applying supervised deep learning (convolutional neural networks, CNNs) to the Dark Energy Survey data down to a magnitude limit of $\sim$21 mag. The methodologies used for the construction of the catalogues include differences such as the cutout sizes, the labels used for training, and the input to the CNN - monochromatic images versus $gri$-band normalized images. In addition, one catalogue is trained using bright galaxies observed with DES ($i<18$), while the other is trained with bright galaxies ($r<17.5$) and `emulated' galaxies up to $r$-band magnitude $22.5$. Despite the different approaches, the agreement between the two catalogues is excellent up to $i<19$, demonstrating that CNN predictions are reliable for samples at least one magnitude fainter than the training sample limit. It also shows that morphological classifications based on monochromatic images are comparable to those based on $gri$-band images, at least in the bright regime. At fainter magnitudes, $i>19$, the overall agreement is good ($\sim$95\%), but is mostly driven by the large spiral fraction in the two catalogues. In contrast, the agreement within the elliptical population is not as good, especially at faint magnitudes. By studying the mismatched cases we are able to identify lenticular galaxies (at least up to $i<19$), which are difficult to distinguish using standard classification approaches. The synergy of both catalogues provides an unique opportunity to select a population of unusual galaxies. △ Less

Submitted 14 September, 2022; originally announced September 2022.

Comments: 17 pages, 14 figures (1 appendix for galaxy examples including 3 figures)

Design of the ECCE Detector for the Electron Ion Collider

Authors: J. K. Adkins, Y. Akiba, A. Albataineh, M. Amaryan, I. C. Arsene, C. Ayerbe Gayoso, J. Bae, X. Bai, M. D. Baker, M. Bashkanov, R. Bellwied, F. Benmokhtar, V. Berdnikov, J. C. Bernauer, F. Bock, W. Boeglin, M. Borysova, E. Brash, P. Brindza, W. J. Briscoe, M. Brooks, S. Bueltmann, M. H. S. Bukhari, A. Bylinkin, R. Capobianco , et al. (259 additional authors not shown)

Abstract: The EIC Comprehensive Chromodynamics Experiment (ECCE) detector has been designed to address the full scope of the proposed Electron Ion Collider (EIC) physics program as presented by the National Academy of Science and provide a deeper understanding of the quark-gluon structure of matter. To accomplish this, the ECCE detector offers nearly acceptance and energy coverage along with excellent track… ▽ More The EIC Comprehensive Chromodynamics Experiment (ECCE) detector has been designed to address the full scope of the proposed Electron Ion Collider (EIC) physics program as presented by the National Academy of Science and provide a deeper understanding of the quark-gluon structure of matter. To accomplish this, the ECCE detector offers nearly acceptance and energy coverage along with excellent tracking and particle identification. The ECCE detector was designed to be built within the budget envelope set out by the EIC project while simultaneously managing cost and schedule risks. This detector concept has been selected to be the basis for the EIC project detector. △ Less

Submitted 20 July, 2024; v1 submitted 6 September, 2022; originally announced September 2022.

Comments: 34 pages, 30 figures, 9 tables

Report number: JLAB-PHY-24-4124

Detector Requirements and Simulation Results for the EIC Exclusive, Diffractive and Tagging Physics Program using the ECCE Detector Concept

Authors: A. Bylinkin, C. T. Dean, S. Fegan, D. Gangadharan, K. Gates, S. J. D. Kay, I. Korover, W. B. Li, X. Li, R. Montgomery, D. Nguyen, G. Penman, J. R. Pybus, N. Santiesteban, R. Trotta, A. Usman, M. D. Baker, J. Frantz, D. I. Glazier, D. W. Higinbotham, T. Horn, J. Huang, G. Huber, R. Reed, J. Roche , et al. (258 additional authors not shown)

Abstract: This article presents a collection of simulation studies using the ECCE detector concept in the context of the EIC's exclusive, diffractive, and tagging physics program, which aims to further explore the rich quark-gluon structure of nucleons and nuclei. To successfully execute the program, ECCE proposed to utilize the detecter system close to the beamline to ensure exclusivity and tag ion beam/fr… ▽ More This article presents a collection of simulation studies using the ECCE detector concept in the context of the EIC's exclusive, diffractive, and tagging physics program, which aims to further explore the rich quark-gluon structure of nucleons and nuclei. To successfully execute the program, ECCE proposed to utilize the detecter system close to the beamline to ensure exclusivity and tag ion beam/fragments for a particular reaction of interest. Preliminary studies confirmed the proposed technology and design satisfy the requirements. The projected physics impact results are based on the projected detector performance from the simulation at 10 or 100 fb^-1 of integrated luminosity. Additionally, a few insights on the potential 2nd Interaction Region can (IR) were also documented which could serve as a guidepost for the future development of a second EIC detector. △ Less

Submitted 6 March, 2023; v1 submitted 30 August, 2022; originally announced August 2022.

Open Heavy Flavor Studies for the ECCE Detector at the Electron Ion Collider

Authors: X. Li, J. K. Adkins, Y. Akiba, A. Albataineh, M. Amaryan, I. C. Arsene, C. Ayerbe Gayoso, J. Bae, X. Bai, M. D. Baker, M. Bashkanov, R. Bellwied, F. Benmokhtar, V. Berdnikov, J. C. Bernauer, F. Bock, W. Boeglin, M. Borysova, E. Brash, P. Brindza, W. J. Briscoe, M. Brooks, S. Bueltmann, M. H. S. Bukhari, A. Bylinkin , et al. (262 additional authors not shown)

Abstract: The ECCE detector has been recommended as the selected reference detector for the future Electron-Ion Collider (EIC). A series of simulation studies have been carried out to validate the physics feasibility of the ECCE detector. In this paper, detailed studies of heavy flavor hadron and jet reconstruction and physics projections with the ECCE detector performance and different magnet options will… ▽ More The ECCE detector has been recommended as the selected reference detector for the future Electron-Ion Collider (EIC). A series of simulation studies have been carried out to validate the physics feasibility of the ECCE detector. In this paper, detailed studies of heavy flavor hadron and jet reconstruction and physics projections with the ECCE detector performance and different magnet options will be presented. The ECCE detector has enabled precise EIC heavy flavor hadron and jet measurements with a broad kinematic coverage. These proposed heavy flavor measurements will help systematically study the hadronization process in vacuum and nuclear medium especially in the underexplored kinematic region. △ Less

Submitted 23 July, 2022; v1 submitted 21 July, 2022; originally announced July 2022.

Comments: Open heavy flavor studies with the EIC reference detector design by the ECCE consortium. 11 pages, 11 figures, to be submitted to the Nuclear Instruments and Methods A

Report number: LANL report number: LA-UR-22-27181

Exclusive J/$ψ$ Detection and Physics with ECCE

Authors: X. Li, J. K. Adkins, Y. Akiba, A. Albataineh, M. Amaryan, I. C. Arsene, C. Ayerbe Gayoso, J. Bae, X. Bai, M. D. Baker, M. Bashkanov, R. Bellwied, F. Benmokhtar, V. Berdnikov, J. C. Bernauer, F. Bock, W. Boeglin, M. Borysova, E. Brash, P. Brindza, W. J. Briscoe, M. Brooks, S. Bueltmann, M. H. S. Bukhari, A. Bylinkin , et al. (262 additional authors not shown)

Abstract: Exclusive heavy quarkonium photoproduction is one of the most popular processes in EIC, which has a large cross section and a simple final state. Due to the gluonic nature of the exchange Pomeron, this process can be related to the gluon distributions in the nucleus. The momentum transfer dependence of this process is sensitive to the interaction sites, which provides a powerful tool to probe the… ▽ More Exclusive heavy quarkonium photoproduction is one of the most popular processes in EIC, which has a large cross section and a simple final state. Due to the gluonic nature of the exchange Pomeron, this process can be related to the gluon distributions in the nucleus. The momentum transfer dependence of this process is sensitive to the interaction sites, which provides a powerful tool to probe the spatial distribution of gluons in the nucleus. Recently the problem of the origin of hadron mass has received lots of attention in determining the anomaly contribution $M_{a}$. The trace anomaly is sensitive to the gluon condensate, and exclusive production of quarkonia such as J/$ψ$ and $Υ$ can serve as a sensitive probe to constrain it. In this paper, we present the performance of the ECCE detector for exclusive J/$ψ$ detection and the capability of this process to investigate the above physics opportunities with ECCE. △ Less

Submitted 21 July, 2022; originally announced July 2022.

Comments: 11 pages, 14 figures, 1 table

Design and Simulated Performance of Calorimetry Systems for the ECCE Detector at the Electron Ion Collider

Authors: F. Bock, N. Schmidt, P. K. Wang, N. Santiesteban, T. Horn, J. Huang, J. Lajoie, C. Munoz Camacho, J. K. Adkins, Y. Akiba, A. Albataineh, M. Amaryan, I. C. Arsene, C. Ayerbe Gayoso, J. Bae, X. Bai, M. D. Baker, M. Bashkanov, R. Bellwied, F. Benmokhtar, V. Berdnikov, J. C. Bernauer, W. Boeglin, M. Borysova, E. Brash , et al. (263 additional authors not shown)

Abstract: We describe the design and performance the calorimeter systems used in the ECCE detector design to achieve the overall performance specifications cost-effectively with careful consideration of appropriate technical and schedule risks. The calorimeter systems consist of three electromagnetic calorimeters, covering the combined pseudorapdity range from -3.7 to 3.8 and two hadronic calorimeters. Key… ▽ More We describe the design and performance the calorimeter systems used in the ECCE detector design to achieve the overall performance specifications cost-effectively with careful consideration of appropriate technical and schedule risks. The calorimeter systems consist of three electromagnetic calorimeters, covering the combined pseudorapdity range from -3.7 to 3.8 and two hadronic calorimeters. Key calorimeter performances which include energy and position resolutions, reconstruction efficiency, and particle identification will be presented. △ Less

Submitted 19 July, 2022; originally announced July 2022.

Comments: 19 pages, 22 figures, 5 tables

Modeling in higher dimensions to improve diagnostic testing accuracy: theory and examples for multiplex saliva-based SARS-CoV-2 antibody assays

Authors: Rayanne A. Luke, Anthony J. Kearsley, Nora Pisanic, Yukari C. Manabe, David L. Thomas, Christopher D. Heaney, Paul N. Patrone

Abstract: The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) pandemic has emphasized the importance and challenges of correctly interpreting antibody test results. Identification of positive and negative samples requires a classification strategy with low error rates, which is hard to achieve when the corresponding measurement values overlap. Additional uncertainty arises when classification s… ▽ More The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) pandemic has emphasized the importance and challenges of correctly interpreting antibody test results. Identification of positive and negative samples requires a classification strategy with low error rates, which is hard to achieve when the corresponding measurement values overlap. Additional uncertainty arises when classification schemes fail to account for complicated structure in data. We address these problems through a mathematical framework that combines high dimensional data modeling and optimal decision theory. Specifically, we show that appropriately increasing the dimension of data better separates positive and negative populations and reveals nuanced structure that can be described in terms of mathematical models. We combine these models with optimal decision theory to yield a classification scheme that better separates positive and negative samples relative to traditional methods such as confidence intervals (CIs) and receiver operating characteristics. We validate the usefulness of this approach in the context of a multiplex salivary SARS-CoV-2 immunoglobulin G assay dataset. This example illustrates how our analysis: (i) improves the assay accuracy (e.g. lowers classification errors by up to 42 % compared to CI methods); (ii) reduces the number of indeterminate samples when an inconclusive class is permissible (e.g. by 40 % compared to the original analysis of the example multiplex dataset); and (iii) decreases the number of antigens needed to classify samples. Our work showcases the power of mathematical modeling in diagnostic classification and highlights a method that can be adopted broadly in public health and clinical settings. △ Less

Submitted 9 November, 2022; v1 submitted 28 June, 2022; originally announced June 2022.

Comments: 14 pages, 4 figures, 1 table, 6 supplemental figures; fixed typos/updated table, corrected PDF in Appendix

Optimal Decision Theory for Diagnostic Testing: Minimizing Indeterminate Classes with Applications to Saliva-Based SARS-CoV-2 Antibody Assays

Authors: Paul N. Patrone, Prajakta Bedekar, Nora Pisanic, Yukari C. Manabe, David L. Thomas, Christopher D. Heaney, Anthony J. Kearsley

Abstract: In diagnostic testing, establishing an indeterminate class is an effective way to identify samples that cannot be accurately classified. However, such approaches also make testing less efficient and must be balanced against overall assay performance. We address this problem by reformulating data classification in terms of a constrained optimization problem that (i) minimizes the probability of lab… ▽ More In diagnostic testing, establishing an indeterminate class is an effective way to identify samples that cannot be accurately classified. However, such approaches also make testing less efficient and must be balanced against overall assay performance. We address this problem by reformulating data classification in terms of a constrained optimization problem that (i) minimizes the probability of labeling samples as indeterminate while (ii) ensuring that the remaining ones are classified with an average target accuracy X. We show that the solution to this problem is expressed in terms of a bathtub principle that holds out those samples with the lowest local accuracy up to an X-dependent threshold. To illustrate the usefulness of this analysis, we apply it to a multiplex, saliva-based SARS-CoV-2 antibody assay and demonstrate up to a 30 % reduction in the number of indeterminate samples relative to more traditional approaches. △ Less

Submitted 31 January, 2022; originally announced February 2022.

Acquisition-invariant brain MRI segmentation with informative uncertainties

Authors: Pedro Borges, Richard Shaw, Thomas Varsavsky, Kerstin Klaser, David Thomas, Ivana Drobnjak, Sebastien Ourselin, M Jorge Cardoso

Abstract: Combining multi-site data can strengthen and uncover trends, but is a task that is marred by the influence of site-specific covariates that can bias the data and therefore any downstream analyses. Post-hoc multi-site correction methods exist but have strong assumptions that often do not hold in real-world scenarios. Algorithms should be designed in a way that can account for site-specific effects,… ▽ More Combining multi-site data can strengthen and uncover trends, but is a task that is marred by the influence of site-specific covariates that can bias the data and therefore any downstream analyses. Post-hoc multi-site correction methods exist but have strong assumptions that often do not hold in real-world scenarios. Algorithms should be designed in a way that can account for site-specific effects, such as those that arise from sequence parameter choices, and in instances where generalisation fails, should be able to identify such a failure by means of explicit uncertainty modelling. This body of work showcases such an algorithm, that can become robust to the physics of acquisition in the context of segmentation tasks, while simultaneously modelling uncertainty. We demonstrate that our method not only generalises to complete holdout datasets, preserving segmentation quality, but does so while also accounting for site-specific sequence choices, which also allows it to perform as a harmonisation tool. △ Less

Submitted 7 November, 2021; originally announced November 2021.

Comments: 25 pages, 8 figures

The role of MRI physics in brain segmentation CNNs: achieving acquisition invariance and instructive uncertainties

Authors: Pedro Borges, Richard Shaw, Thomas Varsavsky, Kerstin Klaser, David Thomas, Ivana Drobnjak, Sebastien Ourselin, M Jorge Cardoso

Abstract: Being able to adequately process and combine data arising from different sites is crucial in neuroimaging, but is difficult, owing to site, sequence and acquisition-parameter dependent biases. It is important therefore to design algorithms that are not only robust to images of differing contrasts, but also be able to generalise well to unseen ones, with a quantifiable measure of uncertainty. In th… ▽ More Being able to adequately process and combine data arising from different sites is crucial in neuroimaging, but is difficult, owing to site, sequence and acquisition-parameter dependent biases. It is important therefore to design algorithms that are not only robust to images of differing contrasts, but also be able to generalise well to unseen ones, with a quantifiable measure of uncertainty. In this paper we demonstrate the efficacy of a physics-informed, uncertainty-aware, segmentation network that employs augmentation-time MR simulations and homogeneous batch feature stratification to achieve acquisition invariance. We show that the proposed approach also accurately extrapolates to out-of-distribution sequence samples, providing well calibrated volumetric bounds on these. We demonstrate a significant improvement in terms of coefficients of variation, backed by uncertainty based volumetric validation. △ Less

Submitted 4 November, 2021; originally announced November 2021.

Comments: 10 pages, 3 figures, published in: Simulation and Synthesis in Medical Imaging 6th International Workshop, SASHIMI 2021, Held in Conjunction with MICCAI 2021, Strasbourg, France, September 27, 2021, Proceedings

Evaluation of silicon carbide as a divertor armor material in DIII-D H-mode discharges

Authors: T. Abrams, S. Bringuier, D. M. Thomas, G. Sinclair, S. Gonderman, L. Holland, D. L. Rudakov, R. S. Wilcox, E. A. Unterberg, F. Scotti

Abstract: Silicon carbide (SiC) represents a promising but largely untested plasma-facing material (PFM) for next-step fusion devices. In this work, an analytic mixed-material erosion model is developed by calculating the physical (via SDTrimSP) and chemical (via empirical scalings) sputtering yield from SiC, Si, and C. The Si content in the near-surface SiC layer is predicted to increase during D plasma bo… ▽ More Silicon carbide (SiC) represents a promising but largely untested plasma-facing material (PFM) for next-step fusion devices. In this work, an analytic mixed-material erosion model is developed by calculating the physical (via SDTrimSP) and chemical (via empirical scalings) sputtering yield from SiC, Si, and C. The Si content in the near-surface SiC layer is predicted to increase during D plasma bombardment due to more efficient physical and chemical sputtering of C relative to Si. Silicon erosion from SiC thereby occurs primarily from sputtering of the enriched Si layer, rather than directly from the SiC itself. SiC coatings on ATJ graphite, manufactured via chemical vapor deposition, were exposed to repeated H-mode plasma discharges in the DIII-D tokamak to test this model. The qualitative trends from analytic modeling are reproduced by the experimental measurements, obtained via spectroscopic inference using the S/XB method. Quantitatively the model slightly under-predicts measured erosion rates, which is attributed to uncertainties in the ion impact angle distribution, as well as the effect of edge-localized modes. After exposure, minimal changes to the macroscopic or microscopic surface morphology of the SiC coatings were observed. Compositional analysis reveals Si enrichment of about 10%, in line with expectations from the erosion model. Extrapolating to a DEMO-type device, an order-of-magnitude decrease in impurity sourcing, and up to a factor of 2 decrease in impurity radiation, is expected with SiC walls, relative to graphite, if low C plasma impurity content can be achieved. These favorable erosion properties motivate further investigations of SiC as a low-Z, non-metallic PFM. △ Less

Submitted 8 April, 2021; originally announced April 2021.

Comments: 17 pages, 15 figures

Exploring the impacts of conformer selection methods on ion mobility collision cross section predictions

Authors: Felicity F. Nielson, Sean M. Colby, Dennis G. Thomas, Ryan S. Renslow, Thomas O. Metz

Abstract: The prediction of structure dependent molecular properties, such as collision cross sections as measured using ion mobility spectrometry, are crucially dependent on the selection of the correct population of molecular conformers. Here, we report an in-depth evaluation of multiple conformation selection techniques, including simple averaging, Boltzmann weighting, lowest energy selection, low energy… ▽ More The prediction of structure dependent molecular properties, such as collision cross sections as measured using ion mobility spectrometry, are crucially dependent on the selection of the correct population of molecular conformers. Here, we report an in-depth evaluation of multiple conformation selection techniques, including simple averaging, Boltzmann weighting, lowest energy selection, low energy threshold reductions, and similarity reduction. Generating 50,000 conformers each for 18 molecules, we used the In Silico Chemical Library Engine (ISiCLE) to calculate the collision cross sections for the entire dataset. First, we employed Monte Carlo simulations to understand the variability between conformer structures as generated using simulated annealing. Then we employed Monte Carlo simulations to the aforementioned conformer selection techniques applied on the simulated molecular property - the ion mobility collision cross section. Based on our analyses, we found Boltzmann weighting to be a good tradeoff between precision and theoretical accuracy. Combining multiple techniques revealed that energy thresholds and root-mean-squared deviation-based similarity reductions can save considerable computational expense while maintaining property prediction accuracy. Molecular dynamic conformer generation tools like AMBER can continue to generate new lowest energy conformers even after tens of thousands of generations, decreasing precision between runs. This reduced precision can be ameliorated and theoretical accuracy increased by running density functional theory geometry optimization on carefully selected conformers. △ Less

Submitted 14 October, 2020; originally announced October 2020.

Physics-informed brain MRI segmentation

Authors: Pedro Borges, Carole Sudre, Thomas Varsavsky, David Thomas, Ivana Drobnjak, Sebastien Ourselin, M. Jorge Cardoso

Abstract: Magnetic Resonance Imaging (MRI) is one of the most flexible and powerful medical imaging modalities. This flexibility does however come at a cost; MRI images acquired at different sites and with different parameters exhibit significant differences in contrast and tissue appearance, resulting in downstream issues when quantifying brain anatomy or the presence of pathology. In this work, we propose… ▽ More Magnetic Resonance Imaging (MRI) is one of the most flexible and powerful medical imaging modalities. This flexibility does however come at a cost; MRI images acquired at different sites and with different parameters exhibit significant differences in contrast and tissue appearance, resulting in downstream issues when quantifying brain anatomy or the presence of pathology. In this work, we propose to combine multiparametric MRI-based static-equation sequence simulations with segmentation convolutional neural networks (CNN), to make these networks robust to variations in acquisition parameters. Results demonstrate that, when given both the image and their associated physics acquisition parameters, CNNs can produce segmentations that exhibit robustness to acquisition variations. We also show that the proposed physics-informed methods can be used to bridge multi-centre and longitudinal imaging studies where imaging acquisition varies across a site or in time. △ Less

Submitted 29 January, 2020; originally announced January 2020.

Comments: Accepted at SASHIMI2019

Impact of ELM control techniques on tungsten sputtering in the DIII-D divertor and extrapolations to ITER

Authors: T. Abrams, E. A. Unterberg, D. L. Rudakov, A. W. Leonard, O. Schmitz, D. Shiraki, L. R. Baylor, P. C. Stangeby, D. M. Thomas, H. Q. Wang

Abstract: The free-streaming plus recycling model (FSRM) has recently been developed to understand and predict tungsten gross erosion rates from the divertor during edge localized modes (ELMs). In this work, the FSRM was tested against experimental measurements of W sputtering during ELMs, conducted via fast WI spectroscopy. Good agreement is observed using a variety of controlling techniques, including gas… ▽ More The free-streaming plus recycling model (FSRM) has recently been developed to understand and predict tungsten gross erosion rates from the divertor during edge localized modes (ELMs). In this work, the FSRM was tested against experimental measurements of W sputtering during ELMs, conducted via fast WI spectroscopy. Good agreement is observed using a variety of controlling techniques, including gas puffing, neutral beam heating, and plasma shaping to modify the pedestal stability boundary and thus the ELM behavior. ELM mitigation by pellet pacing was observed to strongly reduce W sputtering by flushing C impurities from the pedestal and reducing the divertor target electron temperature. No reduction of W sputtering was observed during the application of resonant magnetic perturbations (RMPs), in contrast to the prediction of the FSRM. Potential sources of this discrepancy are discussed. Finally, the framework of the FSRM is utilized to predict intra-ELM W sputtering rates in ITER. It is concluded that W erosion during ELMs in ITER will be caused mainly by free-streaming fuel ions, but free-streaming seeded impurities (N or Ne) may increase the erosion rate significantly if present in the pedestal at even the 1% level. Impurity recycling is not expected to cause significant W erosion in ITER due to the very low target electron temperature. △ Less

Submitted 6 November, 2019; originally announced November 2019.

Journal ref: Phys. Plasmas 26, 062504 (2019)

ISiCLE: A molecular collision cross section calculation pipeline for establishing large in silico reference libraries for compound identification

Authors: Sean M. Colby, Dennis G. Thomas, Jamie R. Nunez, Douglas J. Baxter, Kurt R. Glaesemann, Joseph M. Brown, Meg A Pirrung, Niranjan Govind, Justin G. Teeguarden, Thomas O. Metz, Ryan S. Renslow

Abstract: Comprehensive and confident identifications of metabolites and other chemicals in complex samples will revolutionize our understanding of the role these chemically diverse molecules play in biological systems. Despite recent advances, metabolomics studies still result in the detection of a disproportionate number of features than cannot be confidently assigned to a chemical structure. This inadequ… ▽ More Comprehensive and confident identifications of metabolites and other chemicals in complex samples will revolutionize our understanding of the role these chemically diverse molecules play in biological systems. Despite recent advances, metabolomics studies still result in the detection of a disproportionate number of features than cannot be confidently assigned to a chemical structure. This inadequacy is driven by the single most significant limitation in metabolomics: the reliance on reference libraries constructed by analysis of authentic reference chemicals. To this end, we have developed the in silico chemical library engine (ISiCLE), a high-performance computing-friendly cheminformatics workflow for generating libraries of chemical properties. In the instantiation described here, we predict probable three-dimensional molecular conformers using chemical identifiers as input, from which collision cross sections (CCS) are derived. The approach employs state-of-the-art first-principles simulation, distinguished by use of molecular dynamics, quantum chemistry, and ion mobility calculations to generate structures and libraries, all without training data. Importantly, optimization of ISiCLE included a refactoring of the popular MOBCAL code for trajectory-based mobility calculations, improving its computational efficiency by over two orders of magnitude. Calculated CCS values were validated against 1,983 experimentally-measured CCS values and compared to previously reported CCS calculation approaches. An online database is introduced for sharing both calculated and experimental CCS values (metabolomics.pnnl.gov), initially including a CCS library with over 1 million entries. Finally, three successful applications of molecule characterization using calculated CCS are described. This work represents a promising method to address the limitations of small molecule identification. △ Less

Submitted 21 September, 2018; originally announced September 2018.

Radiative lifetimes of the bound excited states of $\text{Pt}^{-}$

Authors: K. C. Chartkunchand, M. Kamińska, E. K. Anderson, M. K. Kristiansson, G. Eklund, O. M. Hole, R. F. Nascimento, M. Blom, M. Björkhage, A. Källberg, P. Löfgren, P. Reinhed, S. Rosén, A. Simonsson, R. D. Thomas, S. Mannervik, V. T. Davis, P. A. Neill, J. S. Thompson, D. Hanstorp, H. Zettergren, H. Cederquist, H. T. Schmidt

Abstract: The intrinsic radiative lifetimes of the $5d^{10}6s$ $^{2}\text{S}_{1/2}$ and $5d^{9}6s^{2}$ $^{2}\text{D}_{3/2}$ bound excited states in the platinum anion $\text{Pt}^{-}$ have been studied at cryogenic temperatures at the Double ElectroStatic Ion Ring Experiment (DESIREE) facility at Stockholm University. The intrinsic lifetime of the higher-lying $5d^{10}6s$ $^{2}\text{S}_{1/2}$ state was measu… ▽ More The intrinsic radiative lifetimes of the $5d^{10}6s$ $^{2}\text{S}_{1/2}$ and $5d^{9}6s^{2}$ $^{2}\text{D}_{3/2}$ bound excited states in the platinum anion $\text{Pt}^{-}$ have been studied at cryogenic temperatures at the Double ElectroStatic Ion Ring Experiment (DESIREE) facility at Stockholm University. The intrinsic lifetime of the higher-lying $5d^{10}6s$ $^{2}\text{S}_{1/2}$ state was measurement to be 2.54$\pm$0.10 s, while only a lifetime in the range of 50 - 200 ms could be estimated for the $5d^{9}6s^{2}$ $^{2}\text{D}_{3/2}$ fine-structure level. The storage lifetime of the $\text{Pt}^{-}$ ion beam was measured to be a little over 15 minutes at a ring temperature of 13 K. The present study reports the lifetime of an atomic negative ion in an excited bound state with an electron configuration different from that of the ground state. △ Less

Submitted 31 March, 2018; originally announced April 2018.

Comments: 6 pages, 5 figures

Journal ref: Physical Review A 94, 032501 (2016)

First storage of ion beams in the Double Electrostatic Ion-Ring Experiment - DESIREE

Authors: H. T. Schmidt, R. D. Thomas, M. Gatchell, S. Rosén, P. Reinhed, P. Löfgren, L. Brännholm, M. Blom, M. Björkhage, E. Bäckström, J. D. Alexander, S. Leontein, D. Hanstorp, H. Zettergren, L. Liljeby, A. Källberg, A. Simonsson, F. Hellberg, S. Mannervik, M. Larsson, W. D. Geppert, K. G. Rensfelt, H. Danared, A. Paál, M. Masuda , et al. (9 additional authors not shown)

Abstract: We report on the first storage of ion beams in the Double ElectroStatic Ion Ring ExpEriment; DESIREE, at Stockholm University. We have produced beams of atomic carbon anions and small carbon anion molecules (C$_n^-$, $n=1,2,3,4$) in a sputter ion source. The ion beams were accelerated to 10 keV kinetic energy and stored in an electrostatic ion storage ring enclosed in a vacuum chamber at 13 K. For… ▽ More We report on the first storage of ion beams in the Double ElectroStatic Ion Ring ExpEriment; DESIREE, at Stockholm University. We have produced beams of atomic carbon anions and small carbon anion molecules (C$_n^-$, $n=1,2,3,4$) in a sputter ion source. The ion beams were accelerated to 10 keV kinetic energy and stored in an electrostatic ion storage ring enclosed in a vacuum chamber at 13 K. For 10 keV C$_2^-$ molecular anions we measure the residual-gas limited beam storage lifetime to be 448 s $\pm$ 18 s with two independent detector systems. Using the measured storage lifetimes we estimate that the residual gas pressure is in the 10$^{-14}$ mbar range. When high current ion beams are injected, the number of stored particles does not follow a single exponential decay law as would be expected for stored particles lost solely due to electron detachment in collision with the residual-gas. Instead, we observe a faster initial decay rate, which we ascribe to the effect of the space charge of the ion beam on the storage capacity. %The latter effect becomes insignificant after longer storage times of typically 100-150 seconds and we then observe a constant decay rate due to residual-gas collisions. △ Less

Submitted 31 March, 2018; originally announced April 2018.

Journal ref: Rev. Sci. Instrum. 84, 055115 (2013)

Rotationally cold OH$^-$ ions in the cryogenic electrostatic ion-beam storage ring DESIREE

Authors: H. T. Schmidt, G. Eklund, K. C. Chartkunchand, E. K. Anderson, M. Kamińska, N. de Ruette, R. D. Thomas, M. K. Kristiansson, M. Gatchell, P. Reinhed, S. Rosén, A. Simonsson, A. Källberg, P. Löfgren, S. Mannervik, H. Zettergren, H. Cederquist

Abstract: We apply near-threshold laser photodetachment to characterize the rotational quantum level distribution of OH$^-$ ions stored in the cryogenic ion-beam storage ring, DESIREE, at Stockholm University. We find that the stored ions relax to a rotational temperature of 13.4$\pm$0.2 K with 94.9$\pm$0.3 % of the ions in the rotational ground state. This is consistent with the storage ring temperature of… ▽ More We apply near-threshold laser photodetachment to characterize the rotational quantum level distribution of OH$^-$ ions stored in the cryogenic ion-beam storage ring, DESIREE, at Stockholm University. We find that the stored ions relax to a rotational temperature of 13.4$\pm$0.2 K with 94.9$\pm$0.3 % of the ions in the rotational ground state. This is consistent with the storage ring temperature of 13.5$\pm$0.5 K as measured with eight silicon diodes, but in contrast to all earlier studies in cryogenic traps and rings where the rotational temperatures were always much higher than those of the storage devices at their lowest temperatures. Furthermore, we actively modify the rotational distribution through selective photodetachment to produce an OH$^-$ beam where 99.1$\pm$0.1 % of approximately one million stored ions are in the $J$=0 rotational ground state. △ Less

Submitted 31 March, 2018; v1 submitted 28 March, 2018; originally announced March 2018.

Journal ref: Phys. Rev. Lett. 119, 073001 (2017)

Prefer Nested Segmentation to Compound Segmentation

Authors: Haley D. Clark, Stefan A. Reinsberg, Vitali Moiseenko, Jonn Wu, Steven D. Thomas

Abstract: Introduction: Intra-organ radiation dose sensitivity is becoming increasingly relevant in clinical radiotherapy. One method for assessment involves partitioning delineated regions of interest and comparing the relative contributions or importance to clinical outcomes. We show that an intuitive method for dividing organ contours, compound (sub-)segmentation, can unintentionally lead to sub-segments… ▽ More Introduction: Intra-organ radiation dose sensitivity is becoming increasingly relevant in clinical radiotherapy. One method for assessment involves partitioning delineated regions of interest and comparing the relative contributions or importance to clinical outcomes. We show that an intuitive method for dividing organ contours, compound (sub-)segmentation, can unintentionally lead to sub-segments with inconsistent volumes, which will bias relative importance assessment. An improved technique, nested segmentation, is introduced and compared. Methods: Clinical radiotherapy planning parotid contours from 510 patients were segmented. Counts of radiotherapy dose matrix voxels interior to sub-segments were used to determine the equivalency of sub-segment volumes. The distribution of voxel counts within sub-segments were compared using Kolmogorov-Smirnov tests and characterized by their dispersion. Analytical solutions for 2D/3D analogues were derived and sub-segment area/volume were compared directly. Results: Both parotid and 2D/3D region of interest analogue segmentation confirmed compound segmentation intrinsically produces sub-segments with volumes that depend on the region of interest shape and selection location. Significant volume differences were observed when sub-segmenting parotid contours into 18ths, and vanishingly small sub-segments were observed when sub-segmenting into 96ths. Central sub-segments were considerably smaller than sub-segments on the periphery. Nested segmentation did not exhibit these shortcomings and produced sub-segments with equivalent volumes when dose grid and contour collinearity was addressed, even when dividing the parotid into 96ths. Nested segmentation was always faster or equivalent in runtime to compound segmentation. Conclusions: Nested segmentation is more suited than compound segmentation for analyses requiring equal weighting of sub-segments. △ Less

Submitted 3 May, 2017; originally announced May 2017.

Comments: 7 figures

Overview of recent physics results from MAST

Authors: A Kirk, J Adamek, RJ Akers, S Allan, L Appel, F Arese Lucini, M Barnes, T Barrett, N Ben Ayed, W Boeglin, J Bradley, P K Browning, J Brunner, P Cahyna, M Carr, F Casson, M Cecconello, C Challis, IT Chapman, S Chapman, S Conroy, N Conway, WA Cooper, M Cox, N Crocker , et al. (138 additional authors not shown)

Abstract: New results from MAST are presented that focus on validating models in order to extrapolate to future devices. Measurements during start-up experiments have shown how the bulk ion temperature rise scales with the square of the reconnecting field. During the current ramp up models are not able to correctly predict the current diffusion. Experiments have been performed looking at edge and core turbu… ▽ More New results from MAST are presented that focus on validating models in order to extrapolate to future devices. Measurements during start-up experiments have shown how the bulk ion temperature rise scales with the square of the reconnecting field. During the current ramp up models are not able to correctly predict the current diffusion. Experiments have been performed looking at edge and core turbulence. At the edge detailed studies have revealed how filament characteristic are responsible for determining the near and far SOL density profiles. In the core the intrinsic rotation and electron scale turbulence have been measured. The role that the fast ion gradient has on redistributing fast ions through fishbone modes has led to a redesign of the neutral beam injector on MAST Upgrade. In H-mode the turbulence at the pedestal top has been shown to be consistent with being due to electron temperature gradient modes. A reconnection process appears to occur during ELMs and the number of filaments released determines the power profile at the divertor. Resonant magnetic perturbations can mitigate ELMs provided the edge peeling response is maximised and the core kink response minimised. The mitigation of intrinsic error fields with toroidal mode number n>1 has been shown to be important for plasma performance. △ Less

Submitted 18 November, 2016; originally announced November 2016.

Comments: 34 pages, 10 figures. This is an author-created, un-copyedited version of an article submitted for publication in Nuclear Fusion. IoP Publishing Ltd is not responsible for any errors or omissions in this version of the manuscript or any version derived from it

A treecode to simulate dust-plasma interactions

Authors: D. M. Thomas, J. T. Holgate

Abstract: The interaction of a small object with surrounding plasma is an area of plasma-physics research with a multitude of applications. This paper introduces the plasma octree code pot, a microscopic simulator of a spheroidal dust grain in a plasma. pot uses the Barnes-Hut treecode algorithm to perform $N$-body simulations of electrons and ions in the vicinity of a chargeable spheroid, employing also th… ▽ More The interaction of a small object with surrounding plasma is an area of plasma-physics research with a multitude of applications. This paper introduces the plasma octree code pot, a microscopic simulator of a spheroidal dust grain in a plasma. pot uses the Barnes-Hut treecode algorithm to perform $N$-body simulations of electrons and ions in the vicinity of a chargeable spheroid, employing also the Boris particle-motion integrator and Hutchinson's reinjection algorithm from SCEPTIC; a description of the implementation of all three algorithms is provided. We present results from pot simulations of the charging of spheres in magnetized plasmas, and of spheroids in unmagnetized plasmas. The results call into question the validity of using the Boltzmann relation in hybrid PIC codes. △ Less

Submitted 1 November, 2016; v1 submitted 19 May, 2016; originally announced May 2016.

Comments: 20 pages, 7 figures, revised for PPCF

The lifetime of the bound excited level in Ni$^-$

Authors: M. Kamińska, V. T. Davis, O. M. Hole, R. F. Nascimento, K. C. Chartakunchand, M. Blom, M. Björkhage, A. Källberg, P. Löfgren, P. Reinhed, S. Rosén, A. Simonsson, R. D. Thomas, S. Mannervik, P. A. Neill, J. S. Thompson, H. T. Schmidt, H. Cederquist, D. Hanstorp

Abstract: The intrinsic lifetime of the upper level in the bound-bound 3d$^9$4s$^2$ $^2$D$_{3/2}$ $\rightarrow$ 3d$^9$4s$^2$ $^2$D$_{5/2}$ radiative transition in Ni$^-$ was measured to be 15.1 $\pm$ 0.4 s. The experiment was performed at cryogenic temperatures in one of the ion-beam storage rings of the DESIREE (Double ElectroStatic Ion Ring Experiment) facility at Stockholm University. The storage lifetim… ▽ More The intrinsic lifetime of the upper level in the bound-bound 3d$^9$4s$^2$ $^2$D$_{3/2}$ $\rightarrow$ 3d$^9$4s$^2$ $^2$D$_{5/2}$ radiative transition in Ni$^-$ was measured to be 15.1 $\pm$ 0.4 s. The experiment was performed at cryogenic temperatures in one of the ion-beam storage rings of the DESIREE (Double ElectroStatic Ion Ring Experiment) facility at Stockholm University. The storage lifetime of the Ni$^-$ ion-beam was measured to be close to five minutes at a ring temperature of 13 K. △ Less

Submitted 2 December, 2015; v1 submitted 13 November, 2015; originally announced November 2015.

Comments: 12 pages, 5 figures

Journal ref: Phys. Rev. A 93, 012512 (2016)

Cylindrical plasmas generated by an annular beam of ultraviolet light

Authors: D. M. Thomas, J. E. Allen

Abstract: We investigate a cylindrical plasma system with ionization, by an annular beam of ultraviolet light, taking place only in the cylinder's outer region. In the steady state, both the outer and inner regions contain a plasma, with that in the inner region being uniform and field-free. At the interface between the two regions, there is an infinitesimal jump in ion density, the magnitude approaching ze… ▽ More We investigate a cylindrical plasma system with ionization, by an annular beam of ultraviolet light, taking place only in the cylinder's outer region. In the steady state, both the outer and inner regions contain a plasma, with that in the inner region being uniform and field-free. At the interface between the two regions, there is an infinitesimal jump in ion density, the magnitude approaching zero in the quasi-neutral $(λ_D \rightarrow 0)$ limit. The system offers the possibility of producing a uniform stationary plasma in the laboratory, hitherto obtained only with thermally produced alkali plasmas. △ Less

Submitted 29 July, 2015; originally announced July 2015.

Comments: 5 figures

Journal ref: Phys. Plasmas 22 (2015) 073507

Observations of 2D Doppler backscattering on MAST

Authors: D. A. Thomas, K. J. Brunner, S. J. Freethy, B. K. Huang, V. F. Shevchenko, R. G. L. Vann

Abstract: The Synthetic Aperture Microwave Imaging (SAMI) diagnostic has conducted proof-of-principle 2D Doppler backscattering (DBS) experiments on MAST. SAMI actively probes the plasma edge using a wide (+-40 degrees vertical and horizontal) and tuneable (10-35.5 GHz) beam. The Doppler backscattered signal is digitised in vector form using an array of eight Vivaldi PCB antennas. This allows the receiving… ▽ More The Synthetic Aperture Microwave Imaging (SAMI) diagnostic has conducted proof-of-principle 2D Doppler backscattering (DBS) experiments on MAST. SAMI actively probes the plasma edge using a wide (+-40 degrees vertical and horizontal) and tuneable (10-35.5 GHz) beam. The Doppler backscattered signal is digitised in vector form using an array of eight Vivaldi PCB antennas. This allows the receiving array to be focused in any direction within the field of view simultaneously to an angular range of 6-24 degrees FWHM at 10-34.5 GHz. This capability is unique to SAMI and is an entirely novel way of conducting DBS experiments. In this paper the feasibility of conducting 2D DBS experiments is explored. Initial measurements of phenomena observed on conventional DBS experiments are presented; such as momentum injection from neutral beams and an abrupt change in power and turbulence velocity coinciding with the onset of H-mode. In addition, being able to carry out 2D DBS imaging allows a measurement of magnetic pitch angle to be made; preliminary results are presented. Capabilities gained through steering a beam using a phased array and the limitations of this technique are discussed. △ Less

Submitted 23 July, 2015; originally announced July 2015.

The role of correlation and solvation in ion interactions with B-DNA

Authors: Maria L. Sushko, Dennis G. Thomas, Suzette A. Pabit, Lois Pollack, Alexey V. Onufriev, Nathan A. Baker

Abstract: The ionic atmospheres around nucleic acids play important roles in biological function. Large-scale explicit solvent simulations coupled to experimental assays such as anomalous small-angle X-ray scattering (ASAXS) can provide important insights into the structure and energetics of such atmospheres but are time- and resource-intensive. In this paper, we use classical density functional theory (cDF… ▽ More The ionic atmospheres around nucleic acids play important roles in biological function. Large-scale explicit solvent simulations coupled to experimental assays such as anomalous small-angle X-ray scattering (ASAXS) can provide important insights into the structure and energetics of such atmospheres but are time- and resource-intensive. In this paper, we use classical density functional theory (cDFT) to explore the balance between ion-DNA, ion-water, and ion-ion interactions in ionic atmospheres of RbCl, SrCl$_2$, and CoHexCl$_3$ (cobalt hexammine chloride) around a B-form DNA molecule. The accuracy of the cDFT calculations was assessed by comparison between simulated and experimental ASAXS curves, demonstrating that an accurate model should take into account ion-ion correlation and ion hydration forces, DNA topology, and the discrete distribution of charges on DNA strands. As expected, these calculations revealed significant differences between monovalent, divalent, and trivalent cation distributions around DNA. About half of the DNA-bound Rb$^+$ ions penetrate into the minor groove of the DNA and half adsorb on the DNA strands. The fraction of cations in the minor groove decreases for the larger Sr$^{2+}$ ions and becomes zero for CoHex$^{3+}$ ions, which all adsorb on the DNA strands. The distribution of CoHex$^{3+}$ ions is mainly determined by Coulomb and steric interactions, while ion-correlation forces play a central role in the monovalent Rb$^+$ distribution and a combination of ion-correlation and hydration forces affect the Sr$^{2+}$ distribution around DNA. △ Less

Submitted 6 December, 2015; v1 submitted 25 June, 2015; originally announced June 2015.

Asymptotic theory of double layer and shielding of electric field at the edge of illuminated plasma

Authors: M. S. Benilov, D. M. Thomas

Abstract: The method of matched asymptotic expansions is applied to the problem of a collisionless plasma generated by UV illumination localized in a central part of the plasma in the limiting case of small Debye length $λ_D$. A second-approximation asymptotic solution is found for the double layer positioned at the boundary of the illuminated region and for the un-illuminated plasma for the plane geometry.… ▽ More The method of matched asymptotic expansions is applied to the problem of a collisionless plasma generated by UV illumination localized in a central part of the plasma in the limiting case of small Debye length $λ_D$. A second-approximation asymptotic solution is found for the double layer positioned at the boundary of the illuminated region and for the un-illuminated plasma for the plane geometry. Numerical calculations for different values of $λ_D$ are reported and found to confirm the asymptotic results. The net integral space charge of the double layer is asymptotically small, although in the plane geometry it is just sufficient to shield the ambipolar electric field existing in the illuminated region and thus to prevent it from penetrating into the un-illuminated region. The double layer has the same mathematical nature as the intermediate transition layer separating an active plasma and a collisionless sheath, and the underlying physics is also the same. In essence, the two layers represent the same physical object: a transonic layer. △ Less

Submitted 27 February, 2014; originally announced February 2014.

Comments: 5 figures, submitted to Physics of Plasmas

Journal ref: Phys. Plasmas 21 (2014) 043501

Plasmas generated by ultra-violet light rather than electron impact

Authors: R. N. Franklin, J. E. Allen, D. M. Thomas, M. S. Benilov

Abstract: We analyze, in both plane and cylindrical geometries, a collisionless plasma consisting of an inner region where generation occurs by UV illumination, and an un-illuminated outer region with no generation. Ions generated in the inner region flow outwards through the outer region and into a wall. We solve for this system's steady state, first in the quasi-neutral regime (where the Debye length… ▽ More We analyze, in both plane and cylindrical geometries, a collisionless plasma consisting of an inner region where generation occurs by UV illumination, and an un-illuminated outer region with no generation. Ions generated in the inner region flow outwards through the outer region and into a wall. We solve for this system's steady state, first in the quasi-neutral regime (where the Debye length $λ_D$ vanishes and analytic solutions exist) and then in the general case, which we solve numerically. In the general case a double layer forms where the illuminated and un-illuminated regions meet, and an approximately quasi-neutral plasma connects the double layer to the wall sheath; in plane geometry the ions coast through the quasi-neutral section at slightly more than the Bohm speed $c_s$. The system, although simple, therefore has two novel features: a double layer that does not require counter-streaming ions and electrons, and a quasi-neutral plasma where ions travel in straight lines with at least the Bohm speed. We close with a précis of our asymptotic solutions of this system, and suggest how our theoretical conclusions might be extended and tested in the laboratory. △ Less

Submitted 2 December, 2013; originally announced December 2013.

Comments: 10 pages, 3 figures, accepted by Physics of Plasmas

Journal ref: Phys. Plasmas 20 (2013) 123508

The PHENIX Forward Silicon Vertex Detector

Authors: C. Aidala, L. Anaya, E. Anderssen, A. Bambaugh, A. Barron, J. G. Boissevain, J. Bok, S. Boose, M. L. Brooks, S. Butsyk, M. Cepeda, P. Chacon, S. Chacon, L. Chavez, T. Cote, C. D'Agostino, A. Datta, K. DeBlasio, L. DelMonte, E. J. Desmond, J. M. Durham, D. Fields, M. Finger, C. Gingu, B. Gonzales , et al. (60 additional authors not shown)

Abstract: A new silicon detector has been developed to provide the PHENIX experiment with precise charged particle tracking at forward and backward rapidity. The Forward Silicon Vertex Tracker (FVTX) was installed in PHENIX prior to the 2012 run period of the Relativistic Heavy Ion Collider (RHIC). The FVTX is composed of two annular endcaps, each with four stations of silicon mini-strip sensors, covering a… ▽ More A new silicon detector has been developed to provide the PHENIX experiment with precise charged particle tracking at forward and backward rapidity. The Forward Silicon Vertex Tracker (FVTX) was installed in PHENIX prior to the 2012 run period of the Relativistic Heavy Ion Collider (RHIC). The FVTX is composed of two annular endcaps, each with four stations of silicon mini-strip sensors, covering a rapidity range of $1.2<|η|<2.2$ that closely matches the two existing PHENIX muon arms. Each station consists of 48 individual silicon sensors, each of which contains two columns of mini-strips with 75 $μ$m pitch in the radial direction and lengths in the $φ$ direction varying from 3.4 mm at the inner radius to 11.5 mm at the outer radius. The FVTX has approximately 0.54 million strips in each endcap. These are read out with FPHX chips, developed in collaboration with Fermilab, which are wire bonded directly to the mini-strips. The maximum strip occupancy reached in central Au-Au collisions is approximately 2.8%. The precision tracking provided by this device makes the identification of muons from secondary vertices away from the primary event vertex possible. The expected distance of closest approach (DCA) resolution of 200 $μ$m or better for particles with a transverse momentum of 5 GeV/$c$ will allow identification of muons from relatively long-lived particles, such as $D$ and $B$ mesons, through their broader DCA distributions. △ Less

Submitted 14 February, 2014; v1 submitted 14 November, 2013; originally announced November 2013.

Comments: version submitted to NIM

Delay differential equation models for single and coupled bubble dynamics in a compressible liquid

Authors: Derek C. Thomas, Yurii A. Ilinskii, Mark F. Hamilton

Abstract: Various models for interacting spherical bubbles in a compressible liquid based on delay differential equations are considered. It is shown that most previously proposed models for interacting spherical bubbles in a compressible liquid based on the Keller-Miksis and Gilmore-Akulichev models are unstable for closely spaced bubbles. A new model for a single spherical bubble in a compressible liquid… ▽ More Various models for interacting spherical bubbles in a compressible liquid based on delay differential equations are considered. It is shown that most previously proposed models for interacting spherical bubbles in a compressible liquid based on the Keller-Miksis and Gilmore-Akulichev models are unstable for closely spaced bubbles. A new model for a single spherical bubble in a compressible liquid is proposed and used to derive a stable model for interacting bubbles. A qualitative comparison to the results of direct numerical integration of the fluid equations of motion suggests that the new model provides more accurate results than the standard Keller-Miksis or Gilmore-Akulichev models for single bubble dynamics. △ Less

Submitted 10 October, 2013; originally announced October 2013.

Comments: 13 pages, 4 figures

Equilibrium probability distribution of a conductive sphere's floating charge in a collisionless, drifting Maxwellian plasma

Authors: Drew M. Thomas, Michael Coppins

Abstract: A dust grain in a plasma has a fluctuating electric charge, and past work concludes that spherical grains in a stationary, collisionless plasma have an essentially Gaussian charge probability distribution. This paper extends that work to flowing plasmas and arbitrarily large spheres, deriving analytic charge probability distributions up to normalizing constants. We find that these distributions al… ▽ More A dust grain in a plasma has a fluctuating electric charge, and past work concludes that spherical grains in a stationary, collisionless plasma have an essentially Gaussian charge probability distribution. This paper extends that work to flowing plasmas and arbitrarily large spheres, deriving analytic charge probability distributions up to normalizing constants. We find that these distributions also have good Gaussian approximations, with analytic expressions for their mean and variance. △ Less

Submitted 27 January, 2014; v1 submitted 24 May, 2013; originally announced May 2013.

Comments: 10 pages & 6 figures. v3 corrects two typos in published version (missing slash added in denom. of eq. 66 and, on the same page, an order-of estimate's fixed to O(1/α))

Journal ref: Phys. Rev. E 88 (2013) 023110

Endocortical bone loss in osteoporosis: The role of bone surface availability

Authors: Pascal R. Buenzli, C. David L. Thomas, John G. Clement, Peter Pivonka

Abstract: Age-related bone loss and postmenopausal osteoporosis are disorders of bone remodelling, in which less bone is reformed than resorbed. Yet, this dysregulation of bone remodelling does not occur equally in all bone regions. Loss of bone is more pronounced near and at the endocortex, leading to cortical wall thinning and medullary cavity expansion, a process sometimes referred to as "trabecularisati… ▽ More Age-related bone loss and postmenopausal osteoporosis are disorders of bone remodelling, in which less bone is reformed than resorbed. Yet, this dysregulation of bone remodelling does not occur equally in all bone regions. Loss of bone is more pronounced near and at the endocortex, leading to cortical wall thinning and medullary cavity expansion, a process sometimes referred to as "trabecularisation" or "cancellisation". Cortical wall thinning is of primary concern in osteoporosis due to the strong deterioration of bone mechanical properties that it is associated with. In this paper, we examine the possibility that the non-uniformity of microscopic bone surface availability could explain the non-uniformity of bone loss in osteoporosis. We use a computational model of bone remodelling in which microscopic bone surface availability influences bone turnover rate and simulate the evolution of the bone volume fraction profile across the midshaft of a long bone. We find that bone loss is accelerated near the endocortical wall where the specific surface is highest. Over time, this leads to a substantial reduction of cortical wall thickness from the endosteum. The associated expansion of the medullary cavity can be made to match experimentally observed cross-sectional data from the Melbourne Femur Collection. Finally, we calculate the redistribution of the mechanical stresses in this evolving bone structure and show that mechanical load becomes critically transferred to the periosteal cortical bone. △ Less

Submitted 6 August, 2012; v1 submitted 26 June, 2012; originally announced June 2012.

Comments: 13 pages, 3 figures. V2: minor stylistic improvements in text/figures; more accurately referenced subsection "Internal mechanical stress distribution"; some improved remarks in the Discussion section

Journal ref: Int J Numer Meth Biomed Engng (2013) 29:1307-1322

Experimental Verification of the Chemical Sensitivity of Two-Site Double Core-Hole States Formed by an X-ray FEL

Authors: P. Salen, P. van der Meulen, H. T. Schmidt, R. D. Thomas, M. Larsson, R. Feifel, M. N. Piancastelli, L. Fang, B. Murphy, T. Osipov, N. Berrah, E. Kukk, K. Ueda, J. D. Bozek, C. Bostedt, S. Wada, R. Richter, V. Feyer, K. C. Prince

Abstract: We have performed X-ray two-photon photoelectron spectroscopy (XTPPS) using the Linac Coherent Light Source (LCLS) X-ray free-electron laser (FEL) in order to study double core-hole (DCH) states of CO2, N2O and N2. The experiment verifies the theory behind the chemical sensitivity of two-site (ts) DCH states by comparing a set of small molecules with respect to the energy shift of the tsDCH state… ▽ More We have performed X-ray two-photon photoelectron spectroscopy (XTPPS) using the Linac Coherent Light Source (LCLS) X-ray free-electron laser (FEL) in order to study double core-hole (DCH) states of CO2, N2O and N2. The experiment verifies the theory behind the chemical sensitivity of two-site (ts) DCH states by comparing a set of small molecules with respect to the energy shift of the tsDCH state and by extracting the relevant parameters from this shift. △ Less

Submitted 2 May, 2012; originally announced May 2012.

Comments: 11 pages, 2 figures

Journal ref: Phys. Rev. Lett. 108, 153003 (2012)

A Not-So-Fundamental Limitation on Studying Complex Systems with Statistics: Comment on Rabin (2011)

Authors: Drew M. Thomas

Abstract: Although living organisms are affected by many interrelated and unidentified variables, this complexity does not automatically impose a fundamental limitation on statistical inference. Nor need one invoke such complexity as an explanation of the "Truth Wears Off" or "decline" effect; similar "decline" effects occur with far simpler systems studied in physics. Selective reporting and publication bi… ▽ More Although living organisms are affected by many interrelated and unidentified variables, this complexity does not automatically impose a fundamental limitation on statistical inference. Nor need one invoke such complexity as an explanation of the "Truth Wears Off" or "decline" effect; similar "decline" effects occur with far simpler systems studied in physics. Selective reporting and publication bias, and scientists' biases in favour of reporting eye-catching results (in general) or conforming to others' results (in physics) better explain this feature of the "Truth Wears Off" effect than Rabin's suggested limitation on statistical inference. △ Less

Submitted 26 August, 2012; v1 submitted 26 April, 2012; originally announced April 2012.

Comments: 4 pages, no figures, 11 references. Commentary on Y. Rabin, Journal of Statistical Physics, 144, 213-216 (2011)

Journal ref: Journal of Statistical Physics, 149 (2012), 1168-1171

The Pierre Auger Observatory V: Enhancements

Authors: The Pierre Auger Collaboration, P. Abreu, M. Aglietta, E. J. Ahn, I. F. M. Albuquerque, D. Allard, I. Allekotte, J. Allen, P. Allison, J. Alvarez Castillo, J. Alvarez-Muñiz, M. Ambrosio, A. Aminaei, L. Anchordoqui, S. Andringa, T. Antičić, A. Anzalone, C. Aramo, E. Arganda, F. Arqueros, H. Asorey, P. Assis, J. Aublin, M. Ave, M. Avenier , et al. (471 additional authors not shown)

Abstract: Ongoing and planned enhancements of the Pierre Auger Observatory Ongoing and planned enhancements of the Pierre Auger Observatory △ Less

Submitted 24 July, 2011; originally announced July 2011.

Comments: Contributions to 32nd International Cosmic Ray Conference, Beijing, China, August 2011

The Pierre Auger Observatory IV: Operation and Monitoring

Authors: The Pierre Auger Collaboration, P. Abreu, M. Aglietta, E. J. Ahn, I. F. M. Albuquerque, D. Allard, I. Allekotte, J. Allen, P. Allison, J. Alvarez Castillo, J. Alvarez-Muñiz, M. Ambrosio, A. Aminaei, L. Anchordoqui, S. Andringa, T. Antičić, A. Anzalone, C. Aramo, E. Arganda, F. Arqueros, H. Asorey, P. Assis, J. Aublin, M. Ave, M. Avenier , et al. (471 additional authors not shown)

Abstract: Technical reports on operations and monitoring of the Pierre Auger Observatory Technical reports on operations and monitoring of the Pierre Auger Observatory △ Less

Submitted 24 July, 2011; originally announced July 2011.

Comments: Constributions to 32nd International Cosmic Ray Conference, Beijing, China, August 2011