-
Pushing the limits of cell segmentation models for imaging mass cytometry
Authors:
Kimberley M. Bird,
Xujiong Ye,
Alan M. Race,
James M. Brown
Abstract:
Imaging mass cytometry (IMC) is a relatively new technique for imaging biological tissue at subcellular resolution. In recent years, learning-based segmentation methods have enabled precise quantification of cell type and morphology, but typically rely on large datasets with fully annotated ground truth (GT) labels. This paper explores the effects of imperfect labels on learning-based segmentation…
▽ More
Imaging mass cytometry (IMC) is a relatively new technique for imaging biological tissue at subcellular resolution. In recent years, learning-based segmentation methods have enabled precise quantification of cell type and morphology, but typically rely on large datasets with fully annotated ground truth (GT) labels. This paper explores the effects of imperfect labels on learning-based segmentation models and evaluates the generalisability of these models to different tissue types. Our results show that removing 50% of cell annotations from GT masks only reduces the dice similarity coefficient (DSC) score to 0.874 (from 0.889 achieved by a model trained on fully annotated GT masks). This implies that annotation time can in fact be reduced by at least half without detrimentally affecting performance. Furthermore, training our single-tissue model on imperfect labels only decreases DSC by 0.031 on an unseen tissue type compared to its multi-tissue counterpart, with negligible qualitative differences in segmentation. Additionally, bootstrapping the worst-performing model (with 5% of cell annotations) a total of ten times improves its original DSC score of 0.720 to 0.829. These findings imply that less time and work can be put into the process of producing comparable segmentation models; this includes eliminating the need for multiple IMC tissue types during training, whilst also providing the potential for models with very few labels to improve on themselves. Source code is available on GitHub: https://github.com/kimberley/ISBI2024.
△ Less
Submitted 6 February, 2024;
originally announced February 2024.
-
An assessment of measuring local levels of homelessness through proxy social media signals
Authors:
Yoshi Meke Bird,
Sarah E. Grobe,
Michael V. Arnold,
Sean P. Rogers,
Mikaela I. Fudolig,
Julia Witte Zimmerman,
Christopher M. Danforth,
Peter Sheridan Dodds
Abstract:
Recent studies suggest social media activity can function as a proxy for measures of state-level public health, detectable through natural language processing. We present results of our efforts to apply this approach to estimate homelessness at the state level throughout the US during the period 2010-2019 and 2022 using a dataset of roughly 1 million geotagged tweets containing the substring ``hom…
▽ More
Recent studies suggest social media activity can function as a proxy for measures of state-level public health, detectable through natural language processing. We present results of our efforts to apply this approach to estimate homelessness at the state level throughout the US during the period 2010-2019 and 2022 using a dataset of roughly 1 million geotagged tweets containing the substring ``homeless.'' Correlations between homelessness-related tweet counts and ranked per capita homelessness volume, but not general-population densities, suggest a relationship between the likelihood of Twitter users to personally encounter or observe homelessness in their everyday lives and their likelihood to communicate about it online. An increase to the log-odds of ``homeless'' appearing in an English-language tweet, as well as an acceleration in the increase in average tweet sentiment, suggest that tweets about homelessness are also affected by trends at the nation-scale. Additionally, changes to the lexical content of tweets over time suggest that reversals to the polarity of national or state-level trends may be detectable through an increase in political or service-sector language over the semantics of charity or direct appeals. An analysis of user account type also revealed changes to Twitter-use patterns by accounts authored by individuals versus entities that may provide an additional signal to confirm changes to homelessness density in a given jurisdiction. While a computational approach to social media analysis may provide a low-cost, real-time dataset rich with information about nationwide and localized impacts of homelessness and homelessness policy, we find that practical issues abound, limiting the potential of social media as a proxy to complement other measures of homelessness.
△ Less
Submitted 15 May, 2023;
originally announced May 2023.
-
A Near Surface Temperature Model of Arrokoth
Authors:
O. M. Umurhan,
W. M. Grundy,
M. K. Bird,
R. Beyer,
J. T. Keane,
I. R. Linscott,
S. Birch,
C. Bierson,
L. A. Young,
S. A. Stern,
C. M. Lisse,
C. J. A. Howett,
S. Protopapa,
J. R. Spencer,
R. P. Binzel,
W. B. Mckinnon,
T. R. Lauer,
H. A. Weaver,
C. B. Olkin,
K. N. Singer,
A. J. Verbiscer,
A. H. Parker
Abstract:
A near surface thermal model for Arrokoth is developed based on the recently released $10^5$ facet model of the body. This thermal solution takes into account Arrokoth's surface re-radiation back onto itself. The solution method exploits Arrokoth's periodic orbital character to develop a thermal response using a time-asymptotic solution method, which involves a Fourier transform solution of the he…
▽ More
A near surface thermal model for Arrokoth is developed based on the recently released $10^5$ facet model of the body. This thermal solution takes into account Arrokoth's surface re-radiation back onto itself. The solution method exploits Arrokoth's periodic orbital character to develop a thermal response using a time-asymptotic solution method, which involves a Fourier transform solution of the heat equation, an approach recently used by others. We display detailed thermal solutions assuming that Arrokoth's near surface material's thermal inertia ${\cal I} = $ 2.5 W/m$^{-2}$K$^{-1}$s$^{1/2}$. We predict that at New Horizons' encounter with Arrokoth its encounter hemisphere surface temperatures were $\sim$ 57-59 K in its polar regions, 30-40 K on its equatorial zones, and 11-13 K for its winter hemisphere. Arrokoth's orbitally averaged temperatures are around 30-35 K in its polar regions, and closer to 40 K near its equatorial zones. Thermal reradiation from the surrounding surface amounts to less than 5\% of the total energy budget, while the total energy ensconced into and exhumed out Arrokoth's interior via thermal conduction over one orbit is about 0.5\% of the total energy budget. As a generalized application of this thermal modeling together with other KBO origins considerations, we favor the interpretation that New Horizons' REX instrument's $29 \pm 5$K brightness temperature measurement is consistent with Arrokoth's near surface material's being made of sub-to-few mm sized tholin-coated amorphous \water ice grains with 1 W/m$^{-2}$K$^{-1}$s$^{1/2}$ $< {\cal I} < $10-20 W/m$^{-2}$K$^{-1}$s$^{1/2}$, and which are characterized by an X-band emissivity in the range 0.9 and 1.
△ Less
Submitted 21 February, 2022;
originally announced February 2022.
-
Initial results from the New Horizons exploration of 2014 MU69, a small Kuiper Belt Object
Authors:
S. A. Stern,
H. A. Weaver,
J. R. Spencer,
C. B. Olkin,
G. R. Gladstone,
W. M. Grundy,
J. M. Moore,
D. P. Cruikshank,
H. A. Elliott,
W. B. McKinnon,
J. Wm. Parker,
A. J. Verbiscer,
L. A. Young,
D. A. Aguilar,
J. M. Albers,
T. Andert,
J. P. Andrews,
F. Bagenal,
M. E. Banks,
B. A. Bauer,
J. A. Bauman,
K. E. Bechtold,
C. B. Beddingfield,
N. Behrooz,
K. B. Beisser
, et al. (180 additional authors not shown)
Abstract:
The Kuiper Belt is a distant region of the Solar System. On 1 January 2019, the New Horizons spacecraft flew close to (486958) 2014 MU69, a Cold Classical Kuiper Belt Object, a class of objects that have never been heated by the Sun and are therefore well preserved since their formation. Here we describe initial results from these encounter observations. MU69 is a bi-lobed contact binary with a fl…
▽ More
The Kuiper Belt is a distant region of the Solar System. On 1 January 2019, the New Horizons spacecraft flew close to (486958) 2014 MU69, a Cold Classical Kuiper Belt Object, a class of objects that have never been heated by the Sun and are therefore well preserved since their formation. Here we describe initial results from these encounter observations. MU69 is a bi-lobed contact binary with a flattened shape, discrete geological units, and noticeable albedo heterogeneity. However, there is little surface color and compositional heterogeneity. No evidence for satellites, ring or dust structures, gas coma, or solar wind interactions was detected. By origin MU69 appears consistent with pebble cloud collapse followed by a low velocity merger of its two lobes.
△ Less
Submitted 2 April, 2020;
originally announced April 2020.
-
Color, Composition, and Thermal Environment of Kuiper Belt Object (486958) Arrokoth
Authors:
W. M. Grundy,
M. K. Bird,
D. T. Britt,
J. C. Cook,
D. P. Cruikshank,
C. J. A. Howett,
S. Krijt,
I. R. Linscott,
C. B. Olkin,
A. H. Parker,
S. Protopapa,
M. Ruaud,
O. M. Umurhan,
L. A. Young,
C. M. Dalle Ore,
J. J. Kavelaars,
J. T. Keane,
Y. J. Pendleton,
S. B. Porter,
F. Scipioni,
J. R. Spencer,
S. A. Stern,
A. J. Verbiscer,
H. A. Weaver,
R. P. Binzel
, et al. (24 additional authors not shown)
Abstract:
The outer Solar System object (486958) Arrokoth (provisional designation 2014 MU$_{69}$) has been largely undisturbed since its formation. We study its surface composition using data collected by the New Horizons spacecraft. Methanol ice is present along with organic material, which may have formed through radiation of simple molecules. Water ice was not detected. This composition indicates hydrog…
▽ More
The outer Solar System object (486958) Arrokoth (provisional designation 2014 MU$_{69}$) has been largely undisturbed since its formation. We study its surface composition using data collected by the New Horizons spacecraft. Methanol ice is present along with organic material, which may have formed through radiation of simple molecules. Water ice was not detected. This composition indicates hydrogenation of carbon monoxide-rich ice and/ or energetic processing of methane condensed on water ice grains in the cold, outer edge of the early Solar System. There are only small regional variations in color and spectra across the surface, suggesting Arrokoth formed from a homogeneous or well-mixed reservoir of solids. Microwave thermal emission from the winter night side is consistent with a mean brightness temperature of 29$\pm$5 K.
△ Less
Submitted 16 February, 2020;
originally announced February 2020.
-
DeepGestalt - Identifying Rare Genetic Syndromes Using Deep Learning
Authors:
Yaron Gurovich,
Yair Hanani,
Omri Bar,
Nicole Fleischer,
Dekel Gelbman,
Lina Basel-Salmon,
Peter Krawitz,
Susanne B Kamphausen,
Martin Zenker,
Lynne M. Bird,
Karen W. Gripp
Abstract:
Facial analysis technologies have recently measured up to the capabilities of expert clinicians in syndrome identification. To date, these technologies could only identify phenotypes of a few diseases, limiting their role in clinical settings where hundreds of diagnoses must be considered.
We developed a facial analysis framework, DeepGestalt, using computer vision and deep learning algorithms,…
▽ More
Facial analysis technologies have recently measured up to the capabilities of expert clinicians in syndrome identification. To date, these technologies could only identify phenotypes of a few diseases, limiting their role in clinical settings where hundreds of diagnoses must be considered.
We developed a facial analysis framework, DeepGestalt, using computer vision and deep learning algorithms, that quantifies similarities to hundreds of genetic syndromes based on unconstrained 2D images. DeepGestalt is currently trained with over 26,000 patient cases from a rapidly growing phenotype-genotype database, consisting of tens of thousands of validated clinical cases, curated through a community-driven platform. DeepGestalt currently achieves 91% top-10-accuracy in identifying over 215 different genetic syndromes and has outperformed clinical experts in three separate experiments.
We suggest that this form of artificial intelligence is ready to support medical genetics in clinical and laboratory practices and will play a key role in the future of precision medicine.
△ Less
Submitted 23 January, 2018;
originally announced January 2018.
-
US Cosmic Visions: New Ideas in Dark Matter 2017: Community Report
Authors:
Marco Battaglieri,
Alberto Belloni,
Aaron Chou,
Priscilla Cushman,
Bertrand Echenard,
Rouven Essig,
Juan Estrada,
Jonathan L. Feng,
Brenna Flaugher,
Patrick J. Fox,
Peter Graham,
Carter Hall,
Roni Harnik,
JoAnne Hewett,
Joseph Incandela,
Eder Izaguirre,
Daniel McKinsey,
Matthew Pyle,
Natalie Roe,
Gray Rybka,
Pierre Sikivie,
Tim M. P. Tait,
Natalia Toro,
Richard Van De Water,
Neal Weiner
, et al. (226 additional authors not shown)
Abstract:
This white paper summarizes the workshop "U.S. Cosmic Visions: New Ideas in Dark Matter" held at University of Maryland on March 23-25, 2017.
This white paper summarizes the workshop "U.S. Cosmic Visions: New Ideas in Dark Matter" held at University of Maryland on March 23-25, 2017.
△ Less
Submitted 14 July, 2017;
originally announced July 2017.
-
Inverse correlation between quasiparticle mass and Tc in a cuprate high-Tc superconductor
Authors:
C. Putzke,
L. Malone,
S. Badoux,
B. Vignolle,
D. Vignolles,
W. Tabis,
P. Walmsley,
M. Bird,
N. E. Hussey,
C. Proust,
A. Carrington
Abstract:
Close to a zero temperature transition between ordered and disordered electronic phases, quantum fluctuations can lead to a strong enhancement of the electron mass and to the emergence of competing phases such as superconductivity. A correlation between the existence of such a quantum phase transition and superconductivity is quite well established in some heavy fermion and iron-based superconduct…
▽ More
Close to a zero temperature transition between ordered and disordered electronic phases, quantum fluctuations can lead to a strong enhancement of the electron mass and to the emergence of competing phases such as superconductivity. A correlation between the existence of such a quantum phase transition and superconductivity is quite well established in some heavy fermion and iron-based superconductors and there have been suggestions that high temperature superconductivity in the copper oxide materials (cuprates) may also be driven by the same mechanism. Close to optimal doping, where the superconducting transition temperature $T_c$ is maximum in the cuprates, two different phases are known to compete with superconductivity: a poorly understood pseudogap phase and a charge ordered phase. Recent experiments have shown a strong increase in quasiparticle mass $m^*$ in the cuprate YBa$_2$Cu$_3$O$_{7-δ}$ as optimal doping is approached suggesting that quantum fluctuations of the charge ordered phase may be responsible for the high-$T_c$ superconductivity. We have tested the robustness of this correlation between $m^*$ and $T_c$ by performing quantum oscillation studies on the stoichiometric compound YBa$_2$Cu$_4$O$_8$ under hydrostatic pressure. In contrast to the results for YBa$_2$Cu$_3$O$_{7-δ}$, we find that in YBa$_2$Cu$_4$O$_8$ the mass decreases as $T_c$ increases under pressure. This inverse correlation between $m^*$ and $T_c$ suggests that quantum fluctuations of the charge order enhance $m^*$ but do not enhance $T_c$.
△ Less
Submitted 29 March, 2016;
originally announced March 2016.
-
The Pluto system: Initial results from its exploration by New Horizons
Authors:
S. A. Stern,
F. Bagenal,
K. Ennico,
G. R. Gladstone,
W. M. Grundy,
W. B. McKinnon,
J. M. Moore,
C. B. Olkin,
J. R. Spencer,
H. A. Weaver,
L. A. Young,
T. Andert,
J. Andrews,
M. Banks,
B. Bauer,
J. Bauman,
O. S. Barnouin,
P. Bedini,
K. Beisser,
R. A. Beyer,
S. Bhaskaran,
R. P. Binzel,
E. Birath,
M. Bird,
D. J. Bogan
, et al. (126 additional authors not shown)
Abstract:
The Pluto system was recently explored by NASA's New Horizons spacecraft, making closest approach on 14 July 2015. Pluto's surface displays diverse landforms, terrain ages, albedos, colors, and composition gradients. Evidence is found for a water-ice crust, geologically young surface units, surface ice convection, wind streaks, volatile transport, and glacial flow. Pluto's atmosphere is highly ext…
▽ More
The Pluto system was recently explored by NASA's New Horizons spacecraft, making closest approach on 14 July 2015. Pluto's surface displays diverse landforms, terrain ages, albedos, colors, and composition gradients. Evidence is found for a water-ice crust, geologically young surface units, surface ice convection, wind streaks, volatile transport, and glacial flow. Pluto's atmosphere is highly extended, with trace hydrocarbons, a global haze layer, and a surface pressure near 10 microbars. Pluto's diverse surface geology and long-term activity raise fundamental questions about how small planets remain active many billions of years after formation. Pluto's large moon Charon displays tectonics and evidence for a heterogeneous crustal composition, its north pole displays puzzling dark terrain. Small satellites Hydra and Nix have higher albedos than expected.
△ Less
Submitted 26 October, 2015;
originally announced October 2015.
-
The Hera Saturn Entry Probe Mission
Authors:
O. Mousis,
D. H. Atkinson,
T. Spilker,
E. Venkatapathy,
J. Poncy,
R. Frampton,
A. Coustenis,
K. Reh,
J. -P. Lebreton,
L. N. Fletcher,
R. Hueso,
M. J. Amato,
A. Colaprete,
F. Ferri,
D. Stam,
P. Wurz,
S. Atreya,
S. Aslam,
D. J. Banfield,
S. Calcutt,
G. Fischer,
A. Holland,
C. Keller,
E. Kessler,
M. Leese
, et al. (19 additional authors not shown)
Abstract:
The Hera Saturn entry probe mission is proposed as an M--class mission led by ESA with a contribution from NASA. It consists of one atmospheric probe to be sent into the atmosphere of Saturn, and a Carrier-Relay spacecraft. In this concept, the Hera probe is composed of ESA and NASA elements, and the Carrier-Relay Spacecraft is delivered by ESA. The probe is powered by batteries, and the Carrier-R…
▽ More
The Hera Saturn entry probe mission is proposed as an M--class mission led by ESA with a contribution from NASA. It consists of one atmospheric probe to be sent into the atmosphere of Saturn, and a Carrier-Relay spacecraft. In this concept, the Hera probe is composed of ESA and NASA elements, and the Carrier-Relay Spacecraft is delivered by ESA. The probe is powered by batteries, and the Carrier-Relay Spacecraft is powered by solar panels and batteries. We anticipate two major subsystems to be supplied by the United States, either by direct procurement by ESA or by contribution from NASA: the solar electric power system (including solar arrays and the power management and distribution system), and the probe entry system (including the thermal protection shield and aeroshell). Hera is designed to perform in situ measurements of the chemical and isotopic compositions as well as the dynamics of Saturn's atmosphere using a single probe, with the goal of improving our understanding of the origin, formation, and evolution of Saturn, the giant planets and their satellite systems, with extrapolation to extrasolar planets. Hera's aim is to probe well into the cloud-forming region of the troposphere, below the region accessible to remote sensing, to the locations where certain cosmogenically abundant species are expected to be well mixed. By leading to an improved understanding of the processes by which giant planets formed, including the composition and properties of the local solar nebula at the time and location of giant planet formation, Hera will extend the legacy of the Galileo and Cassini missions by further addressing the creation, formation, and chemical, dynamical, and thermal evolution of the giant planets, the entire solar system including Earth and the other terrestrial planets, and formation of other planetary systems.
△ Less
Submitted 26 October, 2015;
originally announced October 2015.
-
Maximum likelihood and pseudo score approaches for parametric time-to-event analysis with informative entry times
Authors:
Brian D. M. Tom,
Vernon T. Farewell,
Sheila M. Bird
Abstract:
We develop a maximum likelihood estimating approach for time-to-event Weibull regression models with outcome-dependent sampling, where sampling of subjects is dependent on the residual fraction of the time left to developing the event of interest. Additionally, we propose a two-stage approach which proceeds by iteratively estimating, through a pseudo score, the Weibull parameters of interest (i.e.…
▽ More
We develop a maximum likelihood estimating approach for time-to-event Weibull regression models with outcome-dependent sampling, where sampling of subjects is dependent on the residual fraction of the time left to developing the event of interest. Additionally, we propose a two-stage approach which proceeds by iteratively estimating, through a pseudo score, the Weibull parameters of interest (i.e., the regression parameters) conditional on the inverse probability of sampling weights; and then re-estimating these weights (given the updated Weibull parameter estimates) through the profiled full likelihood. With these two new methods, both the estimated sampling mechanism parameters and the Weibull parameters are consistently estimated under correct specification of the conditional referral distribution. Standard errors for the regression parameters are obtained directly from inverting the observed information matrix in the full likelihood specification and by either calculating bootstrap or robust standard errors for the hybrid pseudo score/profiled likelihood approach. Loss of efficiency with the latter approach is considered. Robustness of the proposed methods to misspecification of the referral mechanism and the time-to-event distribution is also briefly examined. Further, we show how to extend our methods to the family of parametric time-to-event distributions characterized by the generalized gamma distribution. The motivation for these two approaches came from data on time to cirrhosis from hepatitis C viral infection in patients referred to the Edinburgh liver clinic. We analyze these data here.
△ Less
Submitted 31 July, 2014;
originally announced July 2014.
-
Controlling Tokamak Geometry with 3D Magnetic Perturbations
Authors:
Thomas M Bird,
Chris C Hegna
Abstract:
It is shown that small externally applied magnetic perturbations can significantly alter important geometric properties of magnetic flux surfaces in tokamaks. Through 3D shaping, experimentally relevant perturbation levels are large enough to influence turbulent transport and MHD stability in the pedestal region. It is shown that the dominant pitch-resonant flux surface deformations are primarily…
▽ More
It is shown that small externally applied magnetic perturbations can significantly alter important geometric properties of magnetic flux surfaces in tokamaks. Through 3D shaping, experimentally relevant perturbation levels are large enough to influence turbulent transport and MHD stability in the pedestal region. It is shown that the dominant pitch-resonant flux surface deformations are primarily induced by non-resonant 3D fields, particularly in the presence of significant axisymmetric shaping. The spectral content of the applied 3D field can be used to control these effects.
△ Less
Submitted 23 June, 2014;
originally announced June 2014.
-
Electron-hole pair creation by atoms incident on a metal surface
Authors:
J. R. Trail,
D. M. Bird,
M. Persson,
S. Holloway
Abstract:
Electron-hole pair creation by an adsorbate incident on a metal surface is described using \textit{ab initio} methods. The approach starts with standard first principles electronic structure theory, and proceeds to combine classical, quantum oscillator and time dependent density functional methods to provide a consistent description of the non-adiabatic energy transfer from adsorbate to substrat…
▽ More
Electron-hole pair creation by an adsorbate incident on a metal surface is described using \textit{ab initio} methods. The approach starts with standard first principles electronic structure theory, and proceeds to combine classical, quantum oscillator and time dependent density functional methods to provide a consistent description of the non-adiabatic energy transfer from adsorbate to substrate. Of particular interest is the conservation of the total energy at each level of approximation, and the importance of a spin transition as a function of the adsorbate/surface separation. Results are presented and discussed for H and D atoms incident on the Cu(111) surface.
△ Less
Submitted 30 September, 2009;
originally announced September 2009.
-
Energy loss of atoms at metal surfaces due to electron-hole pair excitations: First-principles theory of "chemicurrents"
Authors:
J. R. Trail,
M. C. Graham,
D. M. Bird,
M. Persson,
S. Holloway
Abstract:
A method is presented for calculating electron-hole pair excitation due to an incident atom or molecule interacting with a metal surface. Energy loss is described using an \textit{ab initio} approach that obtains a position-dependent friction coefficient for an adsorbate moving near a metal surface from a total energy pseudopotential calculation. A semi-classical forced oscillator model is const…
▽ More
A method is presented for calculating electron-hole pair excitation due to an incident atom or molecule interacting with a metal surface. Energy loss is described using an \textit{ab initio} approach that obtains a position-dependent friction coefficient for an adsorbate moving near a metal surface from a total energy pseudopotential calculation. A semi-classical forced oscillator model is constructed, using the same friction coefficient description of the energy loss, to describe excitation of the electron gas due to the incident molecule. This approach is applied to H and D atoms incident on a Cu(111) surface, and we obtain theoretical estimates of the `chemicurrents' measured by Nienhaus et al [Phys. Rev. Lett. \textbf{82}, 446 (1999)] for these atoms incident on the surface of a Schottky diode.
△ Less
Submitted 30 September, 2009;
originally announced September 2009.
-
Electronic damping of molecular motion at metal surfaces
Authors:
J. R. Trail,
M. C. Graham,
D. M. Bird
Abstract:
A method for the calculation of the damping rate due to electron-hole pair excitation for atomic and molecular motion at metal surfaces is presented. The theoretical basis is provided by Time Dependent Density Functional Theory (TDDFT) in the quasi-static limit and calculations are performed within a standard plane-wave, pseudopotential framework. The artificial periodicity introduced by using a…
▽ More
A method for the calculation of the damping rate due to electron-hole pair excitation for atomic and molecular motion at metal surfaces is presented. The theoretical basis is provided by Time Dependent Density Functional Theory (TDDFT) in the quasi-static limit and calculations are performed within a standard plane-wave, pseudopotential framework. The artificial periodicity introduced by using a super-cell geometry is removed to derive results for the motion of an isolated atom or molecule, rather than for the coherent motion of an ordered over-layer. The algorithm is implemented in parallel, distributed across both ${\bf k}$ and ${\bf g}$ space, and in a form compatible with the CASTEP code. Test results for the damping of the motion of hydrogen atoms above the Cu(111) surface are presented.
△ Less
Submitted 30 September, 2009;
originally announced September 2009.
-
Embedding methods for large-scale surface calculations
Authors:
J. R. Trail,
D. M. Bird
Abstract:
One of the goals in the development of large scale electronic structure methods is to perform calculations explicitly for a localised region of a system, while still taking into account the rest of the system outside of this region. An example of this in surface physics would be to embed an adsorbate and a few surface atoms into an extended substrate, hence considerably reducing computational co…
▽ More
One of the goals in the development of large scale electronic structure methods is to perform calculations explicitly for a localised region of a system, while still taking into account the rest of the system outside of this region. An example of this in surface physics would be to embed an adsorbate and a few surface atoms into an extended substrate, hence considerably reducing computational costs. Here we apply the constrained electron density method of embedding a Kohn-Sham system in a substrate system (first described by P. Cortona\cite{1} and T.A. Wesolowski\cite{2}), within a plane-wave basis and pseudopotential framework. This approach divides the charge density of the system into substrate and embedded charge densities, the sum of which is the charge density of the actual system of interest. Two test cases are considered. First we construct fcc bulk aluminium by embedding one cubic lattice of atoms within another. Second, we examine a model surface/adsorbate system of aluminium on aluminium and compare with full Kohn-Sham results.
△ Less
Submitted 30 September, 2009;
originally announced September 2009.
-
Density-functional embedding using a plane-wave basis
Authors:
J. R. Trail,
D. M. Bird
Abstract:
The constrained electron density method of embedding a Kohn-Sham system in a substrate system (first described by P. Cortona, Phys. Rev. B {\bf 44}, 8454 (1991) and T.A. Wesolowski and A. Warshel, J. Phys. Chem {\bf 97}, 8050 (1993)) is applied with a plane-wave basis and both local and non-local pseudopotentials. This method divides the electron density of the system into substrate and embedded…
▽ More
The constrained electron density method of embedding a Kohn-Sham system in a substrate system (first described by P. Cortona, Phys. Rev. B {\bf 44}, 8454 (1991) and T.A. Wesolowski and A. Warshel, J. Phys. Chem {\bf 97}, 8050 (1993)) is applied with a plane-wave basis and both local and non-local pseudopotentials. This method divides the electron density of the system into substrate and embedded electron densities, the sum of which is the electron density of the system of interest. Coupling between the substrate and embedded systems is achieved via approximate kinetic energy functionals. Bulk aluminium is examined as a test case for which there is a strong interaction between the substrate and embedded systems. A number of approximations to the kinetic-energy functional, both semi-local and non-local, are investigated. It is found that Kohn-Sham results can be well reproduced using a non-local kinetic energy functional, with the total energy accurate to better than 0.1 eV per atom and good agreement between the electron densities.
△ Less
Submitted 30 September, 2009;
originally announced September 2009.
-
Core reconstruction in pseudopotential calculations
Authors:
J. R. Trail,
D. M. Bird
Abstract:
A new method is presented for obtaining all-electron results from a pseudopotential calculation. This is achieved by carrying out a localised calculation in the region of an atomic nucleus using the embedding potential method of Inglesfield [J.Phys. C {\bf 14}, 3795 (1981)]. In this method the core region is \emph{reconstructed}, and none of the simplifying approximations (such as spherical symm…
▽ More
A new method is presented for obtaining all-electron results from a pseudopotential calculation. This is achieved by carrying out a localised calculation in the region of an atomic nucleus using the embedding potential method of Inglesfield [J.Phys. C {\bf 14}, 3795 (1981)]. In this method the core region is \emph{reconstructed}, and none of the simplifying approximations (such as spherical symmetry of the charge density/potential or frozen core electrons) that previous solutions to this problem have required are made. The embedding method requires an accurate real space Green function, and an analysis of the errors introduced in constructing this from a set of numerical eigenstates is given. Results are presented for an all-electron reconstruction of bulk aluminium, for both the charge density and the density of states.
△ Less
Submitted 30 September, 2009;
originally announced September 2009.
-
Accurate structure factors from pseudopotential methods
Authors:
J. R. Trail,
D. M. Bird
Abstract:
Highly accurate experimental structure factors of silicon are available in the literature, and these provide the ideal test for any \emph{ab initio} method for the construction of the all-electron charge density. In a recent paper [J. R. Trail and D. M. Bird, Phys. Rev. B {\bf 60}, 7863 (1999)] a method has been developed for obtaining an accurate all-electron charge density from a first princip…
▽ More
Highly accurate experimental structure factors of silicon are available in the literature, and these provide the ideal test for any \emph{ab initio} method for the construction of the all-electron charge density. In a recent paper [J. R. Trail and D. M. Bird, Phys. Rev. B {\bf 60}, 7863 (1999)] a method has been developed for obtaining an accurate all-electron charge density from a first principles pseudopotential calculation by reconstructing the core region of an atom of choice. Here this method is applied to bulk silicon, and structure factors are derived and compared with experimental and Full-potential Linear Augmented Plane Wave results (FLAPW). We also compare with the result of assuming the core region is spherically symmetric, and with the result of constructing a charge density from the pseudo-valence density + frozen core electrons. Neither of these approximations provide accurate charge densities. The aspherical reconstruction is found to be as accurate as FLAPW results, and reproduces the residual error between the FLAPW and experimental results.
△ Less
Submitted 30 September, 2009;
originally announced September 2009.
-
Newns-Anderson model of chemicurrents in H/Cu and H/Ag
Authors:
M. S. Mizielinski,
D. M. Bird,
M. Persson,
S. Holloway
Abstract:
The excitation of the electronic system induced by the adsorption of a hydrogen atom on the (111) surfaces of copper and silver is investigated using the time-dependent, mean-field Newns-Anderson model. Parameters for the model are obtained by fitting to density functional theory calculations, allowing the charge and energy transfer between adsorbate and surface to be calculated, together with t…
▽ More
The excitation of the electronic system induced by the adsorption of a hydrogen atom on the (111) surfaces of copper and silver is investigated using the time-dependent, mean-field Newns-Anderson model. Parameters for the model are obtained by fitting to density functional theory calculations, allowing the charge and energy transfer between adsorbate and surface to be calculated, together with the spectrum of electronic excitations. These results are used to make direct comparisons with experimental measurements of chemicurrents, yielding good agreement for both the magnitude of the current and the ratio of the currents for H and D adsorption.
△ Less
Submitted 16 April, 2008;
originally announced April 2008.
-
Discussion of: Statistical analysis of an archeological find--skeptical counting challenges to an archaeological find
Authors:
Sheila M. Bird
Abstract:
Discussion of ``Statistical analysis of an archeological find'' by Andrey Feuerverger [arXiv:0804.0079]
Discussion of ``Statistical analysis of an archeological find'' by Andrey Feuerverger [arXiv:0804.0079]
△ Less
Submitted 1 April, 2008;
originally announced April 2008.
-
Spectrum of electronic excitations due to the adsorption of atoms on metal surfaces
Authors:
M. S. Mizielinski,
D. M. Bird,
M. Persson,
S. Holloway
Abstract:
The time-dependent, mean-field Newns-Anderson model for a spin-polarised adsorbate approaching a metallic surface is solved in the wide-band limit. Equations for the time-evolution of the electronic structure of the adsorbate-metal system are derived and the spectrum of electronic excitations is found. The behaviour of the model is demonstrated for a set of physically reasonable parameters.
The time-dependent, mean-field Newns-Anderson model for a spin-polarised adsorbate approaching a metallic surface is solved in the wide-band limit. Equations for the time-evolution of the electronic structure of the adsorbate-metal system are derived and the spectrum of electronic excitations is found. The behaviour of the model is demonstrated for a set of physically reasonable parameters.
△ Less
Submitted 9 November, 2006;
originally announced November 2006.
-
International Colloquium "Scattering and Scintillation in Radio Astronomy" was held on June 19-23, 2006 in Pushchino, Moscow region, Russia
Authors:
V. I. Shishov,
W. A. Coles,
B. J. Rickett,
M. K. Bird,
A. I. Efimov,
L. N. Samoznaev,
V. K. Rudash,
I. V. Chashei,
D. Plettemeier,
S. R. Spangler,
Yu. Tokarev,
Yu. Belov,
G. Boiko,
G. Komrakov,
J. Chau,
J. Harmon,
M. Sulzer,
M. Kojima,
M. Tokumaru,
K. Fujiki,
P. Janardhan,
B. V. Jackson,
P. P. Hick,
A. Buffington,
M. R. Olyak
, et al. (32 additional authors not shown)
Abstract:
Topics of the Colloquium: a) Interplanetary scintillation b) Interstellar scintillation c) Modeling and physical origin of the interplanetary and the interstellar plasma turbulence d) Scintillation as a tool for investigation of radio sources e) Seeing through interplanetary and interstellar turbulent media Ppt-presentations are available on the Web-site: http://www.prao.ru/conf/Colloquium/main.…
▽ More
Topics of the Colloquium: a) Interplanetary scintillation b) Interstellar scintillation c) Modeling and physical origin of the interplanetary and the interstellar plasma turbulence d) Scintillation as a tool for investigation of radio sources e) Seeing through interplanetary and interstellar turbulent media Ppt-presentations are available on the Web-site: http://www.prao.ru/conf/Colloquium/main.html
△ Less
Submitted 19 September, 2006;
originally announced September 2006.
-
Recent Searches for the Radio Lines of NH_3 in Comets
Authors:
J. Hatchell,
M. K. Bird,
F. F. S. van der Tak,
W. A. Sherwood
Abstract:
Radio observations in the ammonia inversion lines of four comets, C/2001 A2 (LINEAR), 153P/Ikeya-Zhang, C/2001 Q4 (NEAT) and C/2002 T7 (LINEAR), were performed at the Effelsberg 100-m Radio Telescope during their respective close approaches to Earth. None of the four lowest energy metastable lines (J,K=J),J=1--4, could be detected in these comets. We derive the following 3 sigma upper bounds on…
▽ More
Radio observations in the ammonia inversion lines of four comets, C/2001 A2 (LINEAR), 153P/Ikeya-Zhang, C/2001 Q4 (NEAT) and C/2002 T7 (LINEAR), were performed at the Effelsberg 100-m Radio Telescope during their respective close approaches to Earth. None of the four lowest energy metastable lines (J,K=J),J=1--4, could be detected in these comets. We derive the following 3 sigma upper bounds on the NH_3 production rate, and c omparing to the corresponding water production rates, percentage NH_3 abundances relative to H_2O:
Q(NH_3) < 1.9x10^26 s^-1 (0.63%) for C/2001 A2 (LINEAR),
Q(NH_3) < 2.7x10^26 s^-1 (0.13%) for C/2001 Q4 (NEAT),
Q(NH_3) < 2.3x10^27 s^-1 (0.74%) for C/2002 T7 (LINEAR) and Q(NH_3) <= 6.3x10^26 s^-1 (0.63%) for Comet 153P/Ikeya-Zhang.
At 0.74% or less, the ammonia-to-water ratios are factors of 2 below the value for C/1995 O1 (Hale-Bopp) and 1P/Halley, suggesting chemical diversity between comets. The 18-cm lines of OH were clearly detected in the two comets observed during the 2004 campaign, thereby validating the cometary ephemerides.
△ Less
Submitted 20 May, 2005;
originally announced May 2005.
-
A Mixed Basis Approach for the Efficient Calculation of Potential Energy Surfaces
Authors:
O. Gulseren,
D. M. Bird,
S. E. Humphreys
Abstract:
First principles calculations based on density functional theory are having an incerasing impact on our understanding of molecule-surface interactions. For example, calculations of the multi-dimensional potential energy surface have provided considerable insight into th edynamics of dissociation processes. However, these calculations using a plane-wave basis set are very compute expensive if the…
▽ More
First principles calculations based on density functional theory are having an incerasing impact on our understanding of molecule-surface interactions. For example, calculations of the multi-dimensional potential energy surface have provided considerable insight into th edynamics of dissociation processes. However, these calculations using a plane-wave basis set are very compute expensive if they are to be fully converged with respect to the plane-wave energy cutoff, k-point sampling, supercell size, slab thickness, etc. Because of this, in this study, we have implemented a mixed-basis set approach which uses pseudo-atomic orbitals and a few low-energy plane waves as the basis set within a density functional, pseudopotential calculation. We show that the method offers a computationally cheap but accurate alternative. The energy barrier for hydrogen dissociation on Cu(111) is calculated as an example.
△ Less
Submitted 16 October, 1997;
originally announced October 1997.
-
Radio Science Investigation on a Mercury Orbiter Mission
Authors:
J. D. Anderson,
S. G. Turyshev,
S. W. Asmar,
M. K. Bird,
A. S. Konopliv,
T. P. Krisher,
E. L. Lau,
G. Schubert,
W. L. Sjogren
Abstract:
We review the results from {\it Mariner 10} regarding Mercury's gravity field and the results from radar ranging regarding topography. We discuss the implications of improving these results, including a determination of the polar component, as well as the opportunity to perform relativistic gravity tests with a future {\it Mercury Orbiter}. With a spacecraft placed in orbit with periherm at 400…
▽ More
We review the results from {\it Mariner 10} regarding Mercury's gravity field and the results from radar ranging regarding topography. We discuss the implications of improving these results, including a determination of the polar component, as well as the opportunity to perform relativistic gravity tests with a future {\it Mercury Orbiter}. With a spacecraft placed in orbit with periherm at 400 km altitude, apherm at 16,800 km, period 13.45 hr and latitude of periherm at +30 deg, one can expect a significant improvement in our knowledge of Mercury's gravity field and geophysical properties. The 2000 Plus mission that evolved during the European Space Agency (ESA) {\it Mercury Orbiter} assessment study can provide a global gravity field complete through the 25th degree and order in spherical harmonics. If after completion of the main mission, the periherm could be lowered to 200 km altitude, the gravity field could be extended to 50th degree and order. We discuss the possibility that a search for a Hermean ionosphere could be performed during the mission phases featuring Earth occultations.
Because of its relatively large eccentricity and close proximity to the Sun, Mercury's orbital motion provides one of the best solar-system tests of general relativity. Consequently, we emphasize the number of feasible relativistic gravity tests that can be performed within the context of the parameterized post-Newtonian formalism - a useful framework for testing modern gravitational theories. We pointed out that current results on relativistic precession of Mercury's perihelion are uncertain by 0.5 %, and we discuss the expected improvement using {\it Mercury Orbiter}. We discuss the importance of {\it Mercury Orbiter} for setting limits on a possible time variation in the
△ Less
Submitted 16 October, 1995;
originally announced October 1995.
-
The Velocity Dispersion -- Temperature Correlation from a Limited Cluster Sample
Authors:
Christina M. Bird,
Richard F. Mushotzky,
Christopher A. Metzler
Abstract:
Most studies of correlations between X-ray and optical properties of galaxy clusters have used the largest samples of data available, regardless of the morphological types of clusters included. Given the increasing evidence that morphology is related to a cluster's degree of dynamical evolution, we approach the study of X-ray and optical correlations differently. We evaluate the relationship bet…
▽ More
Most studies of correlations between X-ray and optical properties of galaxy clusters have used the largest samples of data available, regardless of the morphological types of clusters included. Given the increasing evidence that morphology is related to a cluster's degree of dynamical evolution, we approach the study of X-ray and optical correlations differently. We evaluate the relationship between velocity dispersion and temperature for a limited set of galaxy clusters taken from Bird (1994), which all possess dominant central galaxies and which have been explicitly corrected for the presence of substructure. We find that $σ_r \propto T^{0.61 \pm 0.13}$. We use a Monte Carlo computer routine to estimate the significance of this deviation from the $σ_r \propto T^{0.5}$ relationship predicted by the virial theorem. We find that the simulated correlation is steeper than the observed value only 4\% of the time, suggesting that the deviation is significant. The combination of protogalactic winds and dynamical friction reproduces nearly exactly the observed relationship between $σ_r$ and $T$.
△ Less
Submitted 4 May, 1995;
originally announced May 1995.
-
The Effects of Substructure on Galaxy Cluster Mass Determinations
Authors:
Christina M. Bird
Abstract:
Although numerous studies of individual galaxy clusters have demonstrated the presence of significant substructure, previous studies of the distribution of masses of galaxy clusters determined from optical observations have failed to explicitly correct for substructure in those systems. In this {\it Letter} I present the distributions of velocity dispersion, mean separation, and dynamical masses o…
▽ More
Although numerous studies of individual galaxy clusters have demonstrated the presence of significant substructure, previous studies of the distribution of masses of galaxy clusters determined from optical observations have failed to explicitly correct for substructure in those systems. In this {\it Letter} I present the distributions of velocity dispersion, mean separation, and dynamical masses of clusters when substructure is eliminated from the cluster datasets. I also discuss the changes in these distributions because of the substructure correction. Comparing the masses of clusters with central galaxies before and after correction for the presence of substructure reveals a significant change. This change is driven by reductions in the mean separation of galaxies, not by a decrease in the velocity dispersions as has generally been assumed. Correction for substructure reduces most significantly the masses of systems with cool X-ray temperatures, suggesting that the use of a constant linear radius (1.5$h _{100} ^{-1}$ Mpc in this study) to determine cluster membership is inappropriate for clusters spanning a range of temperatures and/or morphologies.
△ Less
Submitted 9 March, 1995; v1 submitted 8 March, 1995;
originally announced March 1995.
-
Subcluster Merger and Galaxy Infall in A2151
Authors:
C. M. Bird,
D. S. Davis,
T. C. Beers
Abstract:
We have obtained a 12.5 ksec image of the Hercules Cluster, A2151, with the {\it ROSAT} PSPC. Comparison of the optical and X-ray data suggest the presence of at least three distinct subclusters in A2151. The brightest X-ray emission coincides with the highest-density peak in the galaxy distribution, and is bimodal. The northern subclump, distinct in position and velocity, has {\it no} detectabl…
▽ More
We have obtained a 12.5 ksec image of the Hercules Cluster, A2151, with the {\it ROSAT} PSPC. Comparison of the optical and X-ray data suggest the presence of at least three distinct subclusters in A2151. The brightest X-ray emission coincides with the highest-density peak in the galaxy distribution, and is bimodal. The northern subclump, distinct in position and velocity, has {\it no} detectable X-ray gas. The eastern subclump, apparent in the optical contour map, is indistinguishable from the main clump in velocity space but is clearly visible in the X-ray image. X-ray spectra derived from the central peak of emission yield a best-fit temperature of 1.6 keV. The emission coincident with the eastern clump of galaxies is cooler, 0.8 keV, and is outside the 90\% confidence intervals of the central peak temperature. We suggest that the eastern and central subclusters have recently undergone a merger event. The lack of X-ray emission to the north suggests that those galaxies do not form a physically-distinct structure (i.e. they are not located within a distinct gravitational potential), but rather that they are falling into the cluster core along the filament defined by the Hercules Supercluster.
△ Less
Submitted 1 December, 1994;
originally announced December 1994.
-
Detecting Bimodality in Astronomical Datasets
Authors:
Keith M. Ashman,
Christina M. Bird,
Steven E. Zepf
Abstract:
We discuss statistical techniques for detecting and quantifying bimodality in astronomical datasets. We concentrate on the KMM algorithm, which estimates the statistical significance of bimodality in such datasets and objectively partitions data into sub-populations. By simulating bimodal distributions with a range of properties we investigate the sensitivity of KMM to datasets with varying char…
▽ More
We discuss statistical techniques for detecting and quantifying bimodality in astronomical datasets. We concentrate on the KMM algorithm, which estimates the statistical significance of bimodality in such datasets and objectively partitions data into sub-populations. By simulating bimodal distributions with a range of properties we investigate the sensitivity of KMM to datasets with varying characteristics. Our results facilitate the planning of optimal observing strategies for systems where bimodality is suspected. Mixture-modeling algorithms similar to the KMM algorithm have been used in previous studies to partition the stellar population of the Milky Way into subsystems. We illustrate the broad applicability of KMM by analysing published data on globular cluster metallicity distributions, velocity distributions of galaxies in clusters, and burst durations of gamma-ray sources. PostScript versions of the tables and figures, as well as FORTRAN code for KMM and instructions for its use, are available by anonymous ftp from kula.phsx.ukans.edu.
△ Less
Submitted 10 August, 1994;
originally announced August 1994.
-
Substructure in Clusters and Central Galaxy Peculiar Velocities
Authors:
Christina M. Bird
Abstract:
Formation theories for central dominant galaxies in clusters require them to be located at the minimum of the cluster gravitational potential. However, 32\% (8 out of 25) of the clusters with more than 50 measured redshifts have central galaxies with significant velocity offsets (with respect to other cluster members). By studying their velocity distributions and correlations between velocity an…
▽ More
Formation theories for central dominant galaxies in clusters require them to be located at the minimum of the cluster gravitational potential. However, 32\% (8 out of 25) of the clusters with more than 50 measured redshifts have central galaxies with significant velocity offsets (with respect to other cluster members). By studying their velocity distributions and correlations between velocity and position, I show that the presence of a large peculiar velocity is strongly correlated with the presence of substructure in these massive systems. About 85\% (21 of 25) of all well-studied clusters show some evidence for substructure, in contrast to the 30-40\% found when using only galaxy or gas distributions. The correlation between substructure and central galaxy location verifies the hypothesis of Merritt (1985) and Tremaine (1990) that high peculiar velocities are indicative of recent merger events between less-massive systems of galaxies. Dynamical friction should act quickly to pull the central galaxy, the most massive discrete object in a cluster, to the minimum of the potential. The less-massive galaxies retain information about their primordial subclusters for a longer period of time. I use an objective partitioning algorithm to assign cluster galaxies to their host subclumps. When galaxies are allocated in this fashion to their subclusters, 75\% of the significant velocity offsets are eliminated. Only 2 out of the 25 clusters have central galaxies which are not centrally-located when substructure is considered in the analysis.
△ Less
Submitted 23 February, 1994;
originally announced February 1994.