-
TCP-Diffusion: A Multi-modal Diffusion Model for Global Tropical Cyclone Precipitation Forecasting with Change Awareness
Authors:
Cheng Huang,
Pan Mu,
Cong Bai,
Peter AG Watson
Abstract:
Precipitation from tropical cyclones (TCs) can cause disasters such as flooding, mudslides, and landslides. Predicting such precipitation in advance is crucial, giving people time to prepare and defend against these precipitation-induced disasters. Developing deep learning (DL) rainfall prediction methods offers a new way to predict potential disasters. However, one problem is that most existing m…
▽ More
Precipitation from tropical cyclones (TCs) can cause disasters such as flooding, mudslides, and landslides. Predicting such precipitation in advance is crucial, giving people time to prepare and defend against these precipitation-induced disasters. Developing deep learning (DL) rainfall prediction methods offers a new way to predict potential disasters. However, one problem is that most existing methods suffer from cumulative errors and lack physical consistency. Second, these methods overlook the importance of meteorological factors in TC rainfall and their integration with the numerical weather prediction (NWP) model. Therefore, we propose Tropical Cyclone Precipitation Diffusion (TCP-Diffusion), a multi-modal model for global tropical cyclone precipitation forecasting. It forecasts TC rainfall around the TC center for the next 12 hours at 3 hourly resolution based on past rainfall observations and multi-modal environmental variables. Adjacent residual prediction (ARP) changes the training target from the absolute rainfall value to the rainfall trend and gives our model the ability of rainfall change awareness, reducing cumulative errors and ensuring physical consistency. Considering the influence of TC-related meteorological factors and the useful information from NWP model forecasts, we propose a multi-model framework with specialized encoders to extract richer information from environmental variables and results provided by NWP models. The results of extensive experiments show that our method outperforms other DL methods and the NWP method from the European Centre for Medium-Range Weather Forecasts (ECMWF).
△ Less
Submitted 16 October, 2024;
originally announced October 2024.
-
Unveiling Multiple Physical Processes on a Cluster Galaxy at z=0.3 Using JWST
Authors:
Peter J. Watson,
Benedetta Vulcani,
Ariel Werle,
Bianca Poggianti,
Marco Gullieuszik,
Michele Trenti,
Xin Wang,
Namrata Roy
Abstract:
We present a study of a previously identified candidate ``jellyfish'' galaxy in the Abell 2744 cluster, F0083, which showed faint signs of a tidal interaction in archival imaging. We used publicly available, PSF-matched deep photometric data from the Hubble and James Webb Space Telescopes, to infer the spatially-resolved star-formation history of this galaxy. F0083 shows clear signs of ram-pressur…
▽ More
We present a study of a previously identified candidate ``jellyfish'' galaxy in the Abell 2744 cluster, F0083, which showed faint signs of a tidal interaction in archival imaging. We used publicly available, PSF-matched deep photometric data from the Hubble and James Webb Space Telescopes, to infer the spatially-resolved star-formation history of this galaxy. F0083 shows clear signs of ram-pressure stripping (RPS), with a recently enhanced star-formation rate (SFR) orientated towards the South-West quadrant of the stellar disc. The stellar mass surface density is heavily asymmetric, with a variation of nearly 1 dex between the Western spiral arm and the postulated tidal feature. This feature appears to contain a high proportion of older stars, ruling out RPS as the cause of this ``unwinding''. An investigation of nearby galaxies reveals the most probable interaction candidate to be a compact elliptical galaxy, 28171, of mass $\log_{10}(M_*/\rm{M}_{\odot})=8.53\pm0.04$ embedded in the tidal stream along our line-of-sight. The star-formation history of the tidal feature indicates a steep change in SFR at lookback times $t_L\lesssim1\,$Gyr. The most probable formation scenario of F0083 thus indicates a low-speed interaction with 28171, followed by RPS as the combined system approaches pericentre passage. Our results demonstrate the ability of photometric data to distinguish between these consecutive processes.
△ Less
Submitted 23 September, 2024;
originally announced September 2024.
-
ASTRODEEP-JWST: NIRCam-HST multiband photometry and redshifts for half a million sources in six extragalactic deep fields
Authors:
E. Merlin,
P. Santini,
D. Paris,
M. Castellano,
A. Fontana,
T. Treu,
S. L. Finkelstein,
J. S. Dunlop,
P. Arrabal Haro,
M. Bagley,
K. Boyett,
A. Calabrò,
M. Correnti,
K. Davis,
M. Dickinson,
C. T. Donnan,
H. C. Ferguson,
F. Fortuni,
M. Giavalisco,
K. Glazebrook,
A. Grazian,
N. A. Grogin,
N. Hathi,
M. Hirschmann,
J. S. Kartaltepe
, et al. (30 additional authors not shown)
Abstract:
We present a set of photometric catalogs primarily aimed at providing the community with a comprehensive database for the study of galaxy populations in the high redshift Universe. The set gathers data from eight JWST NIRCam observational programs, targeting the Abell 2744 (GLASS-JWST, UNCOVER, DDT2756 and GO3990), EGS (CEERS), COSMOS and UDS (PRIMER), and GOODS North and South (JADES and NGDEEP)…
▽ More
We present a set of photometric catalogs primarily aimed at providing the community with a comprehensive database for the study of galaxy populations in the high redshift Universe. The set gathers data from eight JWST NIRCam observational programs, targeting the Abell 2744 (GLASS-JWST, UNCOVER, DDT2756 and GO3990), EGS (CEERS), COSMOS and UDS (PRIMER), and GOODS North and South (JADES and NGDEEP) deep fields, for a total area of $\sim$0.2 sq. degrees. Photometric estimates are obtained by means of well-established techniques, including tailored improvements designed to enhance the performance on the specific dataset. We also include new measurements from HST archival data, thus collecting 16 bands spanning from 0.44 to 4.44 $μ$m. A grand total of $\sim$530 thousand sources is detected on stacks of NIRCam 3.56 and 4.44 $μ$m mosaics. We assess the photometric accuracy by comparing fluxes and colors against archival catalogs. We also provide photometric redshift estimates, statistically validated against a large set of robust spectroscopic data. The catalogs are publicly available on the Astrodeep website.
△ Less
Submitted 9 October, 2024; v1 submitted 30 August, 2024;
originally announced September 2024.
-
The GLASS-JWST Early Release Science Program. IV. Data release of 263 spectra from 245 unique sources
Authors:
S. Mascia,
G. Roberts-Borsani,
T. Treu,
L. Pentericci,
W. Chen,
A. Calabrò,
E. Merlin,
D. Paris,
P. Santini,
G. Brammer,
A. Henry,
P. L. Kelly,
C. Mason,
T. Morishita,
T. Nanayakkara,
N. Roy,
X. Wang,
H. Williams,
K. Boyett,
M. Bradač,
M. Castellano,
K. Glazebrook,
T. Jones,
L. Napolitano,
B. Vulcani
, et al. (2 additional authors not shown)
Abstract:
We release fully reduced spectra obtained with NIRSpec onboard JWST as part of the GLASS-JWST Early Release Science Program and a follow-up Director's Discretionary Time program 2756. From these 263 spectra of 245 unique sources, acquired with low ($R =30-300$) and high dispersion ($R\sim2700$) gratings, we derive redshifts for 200 unique sources in the redshift range $z=0-10$. We describe the sam…
▽ More
We release fully reduced spectra obtained with NIRSpec onboard JWST as part of the GLASS-JWST Early Release Science Program and a follow-up Director's Discretionary Time program 2756. From these 263 spectra of 245 unique sources, acquired with low ($R =30-300$) and high dispersion ($R\sim2700$) gratings, we derive redshifts for 200 unique sources in the redshift range $z=0-10$. We describe the sample selection and characterize its high completeness as a function of redshift and apparent magnitude. Comparison with independent estimates based on different methods and instruments shows that the redshifts are accurate, with 80\% differing less than 0.005. We stack the GLASS-JWST spectra to produce the first high-resolution ($R \sim 2700$) JWST spectral template extending in the rest frame wavelength from 2000~Å to 20, 000~Å. Catalogs, reduced spectra, and template are made publicly available to the community.
△ Less
Submitted 29 August, 2024;
originally announced August 2024.
-
Machine learning emulation of precipitation from km-scale regional climate simulations using a diffusion model
Authors:
Henry Addison,
Elizabeth Kendon,
Suman Ravuri,
Laurence Aitchison,
Peter AG Watson
Abstract:
High-resolution climate simulations are very valuable for understanding climate change impacts and planning adaptation measures. This has motivated use of regional climate models at sufficiently fine resolution to capture important small-scale atmospheric processes, such as convective storms. However, these regional models have very high computational costs, limiting their applicability. We presen…
▽ More
High-resolution climate simulations are very valuable for understanding climate change impacts and planning adaptation measures. This has motivated use of regional climate models at sufficiently fine resolution to capture important small-scale atmospheric processes, such as convective storms. However, these regional models have very high computational costs, limiting their applicability. We present CPMGEM, a novel application of a generative machine learning model, a diffusion model, to skilfully emulate precipitation simulations from such a high-resolution model over England and Wales at much lower cost. This emulator enables stochastic generation of high-resolution (8.8km), daily-mean precipitation samples conditioned on coarse-resolution (60km) weather states from a global climate model. The output is fine enough for use in applications such as flood inundation modelling. The emulator produces precipitation predictions with realistic intensities and spatial structures and captures most of the 21st century climate change signal. We show evidence that the emulator has skill for extreme events up to and including 1-in-100 year intensities. Potential applications include producing high-resolution precipitation predictions for large-ensemble climate simulations and downscaling different climate models and climate change scenarios to better sample uncertainty in climate changes at local-scale.
△ Less
Submitted 19 July, 2024;
originally announced July 2024.
-
Machine learning emulation of a local-scale UK climate model
Authors:
Henry Addison,
Elizabeth Kendon,
Suman Ravuri,
Laurence Aitchison,
Peter AG Watson
Abstract:
Climate change is causing the intensification of rainfall extremes. Precipitation projections with high spatial resolution are important for society to prepare for these changes, e.g. to model flooding impacts. Physics-based simulations for creating such projections are very computationally expensive. This work demonstrates the effectiveness of diffusion models, a form of deep generative models, f…
▽ More
Climate change is causing the intensification of rainfall extremes. Precipitation projections with high spatial resolution are important for society to prepare for these changes, e.g. to model flooding impacts. Physics-based simulations for creating such projections are very computationally expensive. This work demonstrates the effectiveness of diffusion models, a form of deep generative models, for generating much more cheaply realistic high resolution rainfall samples for the UK conditioned on data from a low resolution simulation. We show for the first time a machine learning model that is able to produce realistic samples of high-resolution rainfall based on a physical model that resolves atmospheric convection, a key process behind extreme rainfall. By adding self-learnt, location-specific information to low resolution relative vorticity, quantiles and time-mean of the samples match well their counterparts from the high-resolution simulation.
△ Less
Submitted 29 November, 2022;
originally announced November 2022.
-
Padded Helmet Shell Covers in American Football: A Comprehensive Laboratory Evaluation with Preliminary On-Field Findings
Authors:
Nicholas J. Cecchi,
Ashlyn A. Callan,
Landon P. Watson,
Yuzhe Liu,
Xianghao Zhan,
Ramanand V. Vegesna,
Collin Pang,
Enora Le Flao,
Gerald A. Grant,
Michael M. Zeineh,
David B. Camarillo
Abstract:
Protective headgear effects measured in the laboratory may not always translate to the field. In this study, we evaluated the impact attenuation capabilities of a commercially available padded helmet shell cover in the laboratory and field. In the laboratory, we evaluated the efficacy of the padded helmet shell cover in attenuating impact magnitude across six impact locations and three impact velo…
▽ More
Protective headgear effects measured in the laboratory may not always translate to the field. In this study, we evaluated the impact attenuation capabilities of a commercially available padded helmet shell cover in the laboratory and field. In the laboratory, we evaluated the efficacy of the padded helmet shell cover in attenuating impact magnitude across six impact locations and three impact velocities when equipped to three different helmet models. In a preliminary on-field investigation, we used instrumented mouthguards to monitor head impact magnitude in collegiate linebackers during practice sessions while not wearing the padded helmet shell covers (i.e., bare helmets) for one season and whilst wearing the padded helmet shell covers for another season. The addition of the padded helmet shell cover was effective in attenuating the magnitude of angular head accelerations and two brain injury risk metrics (DAMAGE, HARM) across most laboratory impact conditions, but did not significantly attenuate linear head accelerations for all helmets. Overall, HARM values were reduced in laboratory impact tests by an average of 25% at 3.5 m/s (range: 9.7 - 39.6%), 18% at 5.5 m/s (range: -5.5 - 40.5%), and 10% at 7.4 m/s (range: -6.0 - 31.0%). However, on the field, no significant differences in any measure of head impact magnitude were observed between the bare helmet impacts and padded helmet impacts. Further laboratory tests were conducted to evaluate the ability of the padded helmet shell cover to maintain its performance after exposure to repeated, successive impacts and across a range of temperatures. This research provides a detailed assessment of padded helmet shell covers and supports the continuation of in vivo helmet research to validate laboratory testing results.
△ Less
Submitted 27 November, 2022;
originally announced November 2022.
-
Correlative extinction and single fluorophore bleaching microscopy for ligand quantification on gold nanoparticles
Authors:
Nicole Slesiona,
Lukas Payne,
Iestyn Pope,
Paola Borri,
Wolfgang Langbein,
Peter Watson
Abstract:
Nanoparticles (NPs) are promising therapeutic delivery agents, yet it is increasingly apparent that the number and manner of presentation of cell binding ligands on the NP can affect the eventual fate of the therapeutic. Whenever NPs are conjugated with biomolecules, a heterogenous population of decorated NPs will be produced and the details of the subpopulations of particle-ligand structures need…
▽ More
Nanoparticles (NPs) are promising therapeutic delivery agents, yet it is increasingly apparent that the number and manner of presentation of cell binding ligands on the NP can affect the eventual fate of the therapeutic. Whenever NPs are conjugated with biomolecules, a heterogenous population of decorated NPs will be produced and the details of the subpopulations of particle-ligand structures needs to be characterised for a reliable interpretation of NP-based data. We report an optical microscopy method to quantitatively evaluate the conjugation on a single particle basis in samples consisting of gold NPs (GNPs) decorated with human holo-transferrin fluorescently labelled with Alexa647 (Tf). We employed widefield fluorescence and extinction microscopy on NP-ligand constructs sparesly deposited onto a glass surface, alongside a correlative analysis which spatially co-localises diffraction-limited sources of fluorescence with the optical extinction by individual GNPs. A photobleaching step analysis of the fluorescence emission was employed to estimate the number of fluorophores contributing to the detected emission rate. The method quantifies the number of fluorescent biomolecules attached per GNP, the numbers of unconjugated GNPs and unbound Tf present within the mixed population, and the size and intraparticle clustering propensity of conjugated GNPs. We found a high variability in the number of Tf ligands per GNP within the GNP population, when analysed at the single-particle level, unraveling a non-trivial statistical distribution not accessible in ensemble averaged approaches
△ Less
Submitted 27 September, 2022;
originally announced October 2022.
-
Machine learning applications for weather and climate need greater focus on extremes
Authors:
Peter AG Watson
Abstract:
Multiple studies have now demonstrated that machine learning (ML) can give improved skill for predicting or simulating fairly typical weather events, for tasks such as short-term and seasonal weather forecasting, downscaling simulations to higher resolution and emulating and speeding up expensive model parameterisations. Many of these used ML methods with very high numbers of parameters, such as n…
▽ More
Multiple studies have now demonstrated that machine learning (ML) can give improved skill for predicting or simulating fairly typical weather events, for tasks such as short-term and seasonal weather forecasting, downscaling simulations to higher resolution and emulating and speeding up expensive model parameterisations. Many of these used ML methods with very high numbers of parameters, such as neural networks, which are the focus of the discussion here. Not much attention has been given to the performance of these methods for extreme event severities of relevance for many critical weather and climate prediction applications, with return periods of more than a few years. This leaves a lot of uncertainty about the usefulness of these methods, particularly for general purpose prediction systems that must perform reliably in extreme situations. ML models may be expected to struggle to predict extremes due to there usually being few samples of such events. However, there are some studies that do indicate that ML models can have reasonable skill for extreme weather, and that it is not hopeless to use them in situations requiring extrapolation. This article reviews these studies and argues that this is an area that needs researching more. Ways to get a better understanding of how well ML models perform at predicting extreme weather events are discussed.
△ Less
Submitted 20 September, 2022; v1 submitted 15 July, 2022;
originally announced July 2022.
-
The SAMI Galaxy Survey: The Link Between [$α$/Fe] and Kinematic Morphology
Authors:
Peter J. Watson,
Roger L. Davies,
Jesse van de Sande,
Sarah Brough,
Scott M. Croom,
Francesco D'Eugenio,
Karl Glazebrook,
Brent Groves,
Ángel R. López-Sánchez,
Nicholas Scott,
Sam P. Vaughan,
C. Jakob Walcher,
Joss Bland-Hawthorn,
Julia J. Bryant,
Michael Goodwin,
Jon S. Lawrence,
Nuria P. F. Lorente,
Matt S. Owers,
Samuel Richards
Abstract:
We explore a sample of 1492 galaxies with measurements of the mean stellar population properties and the spin parameter proxy, $λ_{R_{\rm{e}}}$, drawn from the SAMI Galaxy Survey. We fit a global $\left[α/\rm{Fe}\right]$-$σ$ relation, finding that $\left[α/\rm{Fe}\right]=(0.395\pm0.010)\rm{log}_{10}\left(σ\right)-(0.627\pm0.002)$. We observe an anti-correlation between the residuals…
▽ More
We explore a sample of 1492 galaxies with measurements of the mean stellar population properties and the spin parameter proxy, $λ_{R_{\rm{e}}}$, drawn from the SAMI Galaxy Survey. We fit a global $\left[α/\rm{Fe}\right]$-$σ$ relation, finding that $\left[α/\rm{Fe}\right]=(0.395\pm0.010)\rm{log}_{10}\left(σ\right)-(0.627\pm0.002)$. We observe an anti-correlation between the residuals $Δ\left[α/\rm{Fe}\right]$ and the inclination-corrected $λ_{\,R_{\rm{e}}}^{\rm{\,eo}}$, which can be expressed as $Δ\left[α/\rm{Fe}\right]=(-0.057\pm0.008)λ_{\,R_{\rm{e}}}^{\rm{\,eo}}+(0.020\pm0.003)$. The anti-correlation appears to be driven by star-forming galaxies, with a gradient of $Δ\left[α/\rm{Fe}\right]\sim(-0.121\pm0.015)λ_{\,R_{\rm{e}}}^{\rm{\,eo}}$, although a weak relationship persists for the subsample of galaxies for which star formation has been quenched. We take this to be confirmation that disk-dominated galaxies have an extended duration of star formation. At a reference velocity dispersion of 200 km s$^{-1}$, we estimate an increase in half-mass formation time from $\sim$0.5 Gyr to $\sim$1.2 Gyr from low- to high-$λ_{\,R_{\rm{e}}}^{\rm{\,eo}}$ galaxies. Slow rotators do not appear to fit these trends. Their residual $α$-enhancement is indistinguishable from other galaxies with $λ_{\,R_{\rm{e}}}^{\rm{\,eo}}\lessapprox0.4$, despite being both larger and more massive. This result shows that galaxies with $λ_{\,R_{\rm{e}}}^{\rm{\,eo}}\lessapprox0.4$ experience a similar range of star formation histories, despite their different physical structure and angular momentum.
△ Less
Submitted 26 April, 2022;
originally announced April 2022.
-
On the Complexity of Object Detection on Real-world Public Transportation Images for Social Distancing Measurement
Authors:
Nik Khadijah Nik Aznan,
John Brennan,
Daniel Bell,
Jennine Jonczyk,
Paul Watson
Abstract:
Social distancing in public spaces has become an essential aspect in helping to reduce the impact of the COVID-19 pandemic. Exploiting recent advances in machine learning, there have been many studies in the literature implementing social distancing via object detection through the use of surveillance cameras in public spaces. However, to date, there has been no study of social distance measuremen…
▽ More
Social distancing in public spaces has become an essential aspect in helping to reduce the impact of the COVID-19 pandemic. Exploiting recent advances in machine learning, there have been many studies in the literature implementing social distancing via object detection through the use of surveillance cameras in public spaces. However, to date, there has been no study of social distance measurement on public transport. The public transport setting has some unique challenges, including some low-resolution images and camera locations that can lead to the partial occlusion of passengers, which make it challenging to perform accurate detection. Thus, in this paper, we investigate the challenges of performing accurate social distance measurement on public transportation. We benchmark several state-of-the-art object detection algorithms using real-world footage taken from the London Underground and bus network. The work highlights the complexity of performing social distancing measurement on images from current public transportation onboard cameras. Further, exploiting domain knowledge of expected passenger behaviour, we attempt to improve the quality of the detections using various strategies and show improvement over using vanilla object detection alone.
△ Less
Submitted 14 February, 2022;
originally announced February 2022.
-
Frame Shift Prediction
Authors:
Zheng-Xin Yong,
Patrick D. Watson,
Tiago Timponi Torrent,
Oliver Czulo,
Collin F. Baker
Abstract:
Frame shift is a cross-linguistic phenomenon in translation which results in corresponding pairs of linguistic material evoking different frames. The ability to predict frame shifts enables automatic creation of multilingual FrameNets through annotation projection. Here, we propose the Frame Shift Prediction task and demonstrate that graph attention networks, combined with auxiliary training, can…
▽ More
Frame shift is a cross-linguistic phenomenon in translation which results in corresponding pairs of linguistic material evoking different frames. The ability to predict frame shifts enables automatic creation of multilingual FrameNets through annotation projection. Here, we propose the Frame Shift Prediction task and demonstrate that graph attention networks, combined with auxiliary training, can learn cross-linguistic frame-to-frame correspondence and predict frame shifts.
△ Less
Submitted 5 January, 2022;
originally announced January 2022.
-
A comparison of sports-related head accelerations with and without direct head impacts
Authors:
Samuel J. Raymond,
Yuzhe Liu,
Nicholas J. Cecchi,
Eli Rice,
Ashlyn A. Callan,
Landon P. Watson,
Sohrab Sami,
Zhou Zhou,
Xiaogai Li,
Svein Kleiven,
Michael Zeineh,
David B. Camarillo
Abstract:
Concussion and repeated exposure to mild traumatic brain injury are risks for athletes in many sports. While direct head impacts are analyzed to improve the detection and awareness of head acceleration events so that an athlete's brain health can be appropriately monitored and treated. However, head accelerations can also be induced by impacts with little or no head involvement. In this work we ev…
▽ More
Concussion and repeated exposure to mild traumatic brain injury are risks for athletes in many sports. While direct head impacts are analyzed to improve the detection and awareness of head acceleration events so that an athlete's brain health can be appropriately monitored and treated. However, head accelerations can also be induced by impacts with little or no head involvement. In this work we evaluated if impacts that do not involve direct head contact, such as being pushed in the torso, can be sufficient in collegiate American football to induce head accelerations comparable to direct head impacts. Datasets of impacts with and without direct head contact were collected and compared. These datasets were gathered using a state-of-the-art impact detection algorithm embedded in an instrumented mouthguard to record head kinematics. Video analysis was used to differentiate between impact types. In total, 15 impacts of each type were used in comparison, with clear video screenshots available to distinguish each impact type. Analysis of the kinematics showed that the impacts without direct head contact achieved similar levels of linear and angular accelerations during impact compared to those from direct head impacts. Finite element analyses using the median and peak kinematic signals were used to calculate maximum principal strain of the brain. Statistical analysis revealed that no significant difference was found between the two datasets based on a Bonferroni-adjusted p-value threshold of 0.017 , with the exception of peak linear acceleration. Impacts without direct head contact showed higher mean values of peak linear acceleration values of 17.6 g compared to the direct-head impact mean value of 6.1g. These results indicated that impacts other than direct head impacts could still produce meaningful kinematic loads in the head and as such should be included in athlete health monitoring.
△ Less
Submitted 12 October, 2021; v1 submitted 19 August, 2021;
originally announced August 2021.
-
The SAMI Galaxy Survey: Trends in [α/Fe] as a Function of Morphology and Environment
Authors:
Peter J. Watson,
Roger L. Davies,
Sarah Brough,
Scott M. Croom,
Francesco D'Eugenio,
Karl Glazebrook,
Brent Groves,
Ángel R. López-Sánchez,
Jesse van de Sande,
Nicholas Scott,
Sam P. Vaughan,
Jakob Walcher,
Joss Bland-Hawthorn,
Julia J. Bryant,
Michael Goodwin,
Jon S. Lawrence,
Nuria P. F. Lorente,
Matt S. Owers,
Samuel Richards
Abstract:
We present a new set of index-based measurements of [$α$/Fe] for a sample of 2093 galaxies in the SAMI Galaxy Survey. Following earlier work, we fit a global relation between [$α$/Fe] and the galaxy velocity dispersion $σ$ for red sequence galaxies, [$α$/Fe]=(0.378$\pm$0.009)log($σ$/100)+(0.155$\pm$0.003). We observe a correlation between the residuals and the local environmental surface density,…
▽ More
We present a new set of index-based measurements of [$α$/Fe] for a sample of 2093 galaxies in the SAMI Galaxy Survey. Following earlier work, we fit a global relation between [$α$/Fe] and the galaxy velocity dispersion $σ$ for red sequence galaxies, [$α$/Fe]=(0.378$\pm$0.009)log($σ$/100)+(0.155$\pm$0.003). We observe a correlation between the residuals and the local environmental surface density, whereas no such relation exists for blue cloud galaxies. In the full sample, we find that elliptical galaxies in high-density environments are $α$-enhanced by up to 0.057$\pm$0.014 dex at velocity dispersions $σ$<100 km/s, compared with those in low-density environments. This $α$-enhancement is morphology-dependent, with the offset decreasing along the Hubble sequence towards spirals, which have an offset of 0.019$\pm$0.014 dex. At low velocity dispersion and controlling for morphology, we estimate that star formation in high-density environments is truncated $\sim1$ Gyr earlier than in low-density environments. For elliptical galaxies only, we find support for a parabolic relationship between [$α$/Fe] and $σ$, with an environmental $α$-enhancement of at least 0.03 dex. This suggests strong contributions from both environment and mass-based quenching mechanisms. However, there is no evidence for this behaviour in later morphological types.
△ Less
Submitted 24 August, 2021; v1 submitted 3 June, 2021;
originally announced June 2021.
-
High-resolution Probabilistic Precipitation Prediction for use in Climate Simulations
Authors:
Sherman Lo,
Peter Watson,
Peter Dueben,
Ritabrata Dutta
Abstract:
The accurate prediction of precipitation is important to allow for reliable warnings of flood or drought risk in a changing climate. However, to make trust-worthy predictions of precipitation, at a local scale, is one of the most difficult challenges for today's weather and climate models. This is because important features, such as individual clouds and high-resolution topography, cannot be resol…
▽ More
The accurate prediction of precipitation is important to allow for reliable warnings of flood or drought risk in a changing climate. However, to make trust-worthy predictions of precipitation, at a local scale, is one of the most difficult challenges for today's weather and climate models. This is because important features, such as individual clouds and high-resolution topography, cannot be resolved explicitly within simulations due to the significant computational cost of high-resolution simulations. Climate models are typically run at $\sim$50-100 km resolution which is insufficient to represent local precipitation events in satisfying detail. Here, we develop a method to make probabilistic precipitation predictions based on features that climate models can resolve well and that is not highly sensitive to the approximations used in individual models. To predict, we will use a temporal compound Poisson distribution dependent on the output of climate models at a location. We use the output of Earth System models at coarse resolution $\sim$50 km as input and train the statistical models towards precipitation observations over Wales at $\sim$10 km resolution. A Bayesian inferential scheme is provided so that the compound-Poisson model can be inferred using a Gibbs-within-Metropolis-Elliptic-Slice sampling scheme which enables us to quantify the uncertainty of our predictions. In addition, we use a Gaussian process regressor on the posterior samples of the model parameters, to infer a spatially coherent model and hence to produce spatially coherent rainfall prediction. We illustrate the prediction performance of our model by training over 5 years of the data up to 31st December 1999 and predicting precipitation for 20 years afterwards for Cardiff and Wales.
△ Less
Submitted 25 February, 2021; v1 submitted 17 December, 2020;
originally announced December 2020.
-
TRU-NET: A Deep Learning Approach to High Resolution Prediction of Rainfall
Authors:
Rilwan Adewoyin,
Peter Dueben,
Peter Watson,
Yulan He,
Ritabrata Dutta
Abstract:
Climate models (CM) are used to evaluate the impact of climate change on the risk of floods and strong precipitation events. However, these numerical simulators have difficulties representing precipitation events accurately, mainly due to limited spatial resolution when simulating multi-scale dynamics in the atmosphere. To improve the prediction of high resolution precipitation we apply a Deep Lea…
▽ More
Climate models (CM) are used to evaluate the impact of climate change on the risk of floods and strong precipitation events. However, these numerical simulators have difficulties representing precipitation events accurately, mainly due to limited spatial resolution when simulating multi-scale dynamics in the atmosphere. To improve the prediction of high resolution precipitation we apply a Deep Learning (DL) approach using an input of CM simulations of the model fields (weather variables) that are more predictable than local precipitation. To this end, we present TRU-NET (Temporal Recurrent U-Net), an encoder-decoder model featuring a novel 2D cross attention mechanism between contiguous convolutional-recurrent layers to effectively model multi-scale spatio-temporal weather processes. We use a conditional-continuous loss function to capture the zero-skewed %extreme event patterns of rainfall. Experiments show that our model consistently attains lower RMSE and MAE scores than a DL model prevalent in short term precipitation prediction and improves upon the rainfall predictions of a state-of-the-art dynamical weather model. Moreover, by evaluating the performance of our model under various, training and testing, data formulation strategies, we show that there is enough data for our deep learning approach to output robust, high-quality results across seasons and varying regions.
△ Less
Submitted 12 February, 2021; v1 submitted 20 August, 2020;
originally announced August 2020.
-
A semi-supervised learning framework for quantitative structure-activity regression modelling
Authors:
Oliver P Watson,
Isidro Cortes-Ciriano,
James A Watson
Abstract:
Supervised learning models, also known as quantitative structure-activity regression (QSAR) models, are increasingly used in assisting the process of preclinical, small molecule drug discovery. The models are trained on data consisting of a finite dimensional representation of molecular structures and their corresponding target specific activities. These models can then be used to predict the acti…
▽ More
Supervised learning models, also known as quantitative structure-activity regression (QSAR) models, are increasingly used in assisting the process of preclinical, small molecule drug discovery. The models are trained on data consisting of a finite dimensional representation of molecular structures and their corresponding target specific activities. These models can then be used to predict the activity of previously unmeasured novel compounds. In this work we address two problems related to this approach. The first is to estimate the extent to which the quality of the model predictions degrades for compounds very different from the compounds in the training data. The second is to adjust for the screening dependent selection bias inherent in many training data sets. In the most extreme cases, only compounds which pass an activity-dependent screening are reported. By using a semi-supervised learning framework, we show that it is possible to make predictions which take into account the similarity of the testing compounds to those in the training data and adjust for the reporting selection bias. We illustrate this approach using publicly available structure-activity data on a large set of compounds reported by GlaxoSmithKline (the Tres Cantos AntiMalarial Set) to inhibit in vitro P. falciparum growth.
△ Less
Submitted 7 January, 2020;
originally announced January 2020.
-
P2L: Predicting Transfer Learning for Images and Semantic Relations
Authors:
Bishwaranjan Bhattacharjee,
John R. Kender,
Matthew Hill,
Parijat Dube,
Siyu Huo,
Michael R. Glass,
Brian Belgodere,
Sharath Pankanti,
Noel Codella,
Patrick Watson
Abstract:
Transfer learning enhances learning across tasks, by leveraging previously learned representations -- if they are properly chosen. We describe an efficient method to accurately estimate the appropriateness of a previously trained model for use in a new learning task. We use this measure, which we call "Predict To Learn" ("P2L"), in the two very different domains of images and semantic relations, w…
▽ More
Transfer learning enhances learning across tasks, by leveraging previously learned representations -- if they are properly chosen. We describe an efficient method to accurately estimate the appropriateness of a previously trained model for use in a new learning task. We use this measure, which we call "Predict To Learn" ("P2L"), in the two very different domains of images and semantic relations, where it predicts, from a set of "source" models, the one model most likely to produce effective transfer for training a given "target" model. We validate our approach thoroughly, by assembling a collection of candidate source models, then fine-tuning each candidate to perform each of a collection of target tasks, and finally measuring how well transfer has been enhanced. Across 95 tasks within multiple domains (images classification and semantic relations), the P2L approach was able to select the best transfer learning model on average, while the heuristic of choosing model trained with the largest data set selected the best model in only 55 cases. These results suggest that P2L captures important information in common between source and target tasks, and that this shared informational structure contributes to successful transfer learning more than simple data size.
△ Less
Submitted 15 October, 2020; v1 submitted 20 August, 2019;
originally announced August 2019.
-
The impact of stochastic physics on climate sensitivity in EC-Earth
Authors:
K. Strommen,
P. A. G. Watson,
T. N. Palmer
Abstract:
Stochastic schemes, designed to represent unresolved sub-grid scale variability, are frequently used in short and medium-range weather forecasts, where they are found to improve several aspects of the model. In recent years, the impact of stochastic physics has also been found to be beneficial for the model's long term climate. In this paper, we demonstrate for the first time that the inclusion of…
▽ More
Stochastic schemes, designed to represent unresolved sub-grid scale variability, are frequently used in short and medium-range weather forecasts, where they are found to improve several aspects of the model. In recent years, the impact of stochastic physics has also been found to be beneficial for the model's long term climate. In this paper, we demonstrate for the first time that the inclusion of a stochastic physics scheme can notably affect a model's projection of global warming, as well as its historical climatological global temperature. Specifically, we find that when including the 'stochastically perturbed parametrisation tendencies' scheme (SPPT) in the fully coupled climate model EC-Earth v3.1, the predicted level of global warming between 1850 and 2100 is reduced by 10% under an RCP8.5 forcing scenario. We link this reduction in climate sensitivity to a change in the cloud feedbacks with SPPT. In particular, the scheme appears to reduce the positive low cloud cover feedback, and increase the negative cloud optical feedback. A key role is played by a robust, rapid increase in cloud liquid water with SPPT, which we speculate is due to the scheme's non-linear interaction with condensation.
△ Less
Submitted 16 May, 2019;
originally announced May 2019.
-
Applying machine learning to improve simulations of a chaotic dynamical system using empirical error correction
Authors:
Peter A. G. Watson
Abstract:
Dynamical weather and climate prediction models underpin many studies of the Earth system and hold the promise of being able to make robust projections of future climate change based on physical laws. However, simulations from these models still show many differences compared with observations. Machine learning has been applied to solve certain prediction problems with great success, and recently…
▽ More
Dynamical weather and climate prediction models underpin many studies of the Earth system and hold the promise of being able to make robust projections of future climate change based on physical laws. However, simulations from these models still show many differences compared with observations. Machine learning has been applied to solve certain prediction problems with great success, and recently it's been proposed that this could replace the role of physically-derived dynamical weather and climate models to give better quality simulations. Here, instead, a framework using machine learning together with physically-derived models is tested, in which it is learnt how to correct the errors of the latter from timestep to timestep. This maintains the physical understanding built into the models, whilst allowing performance improvements, and also requires much simpler algorithms and less training data. This is tested in the context of simulating the chaotic Lorenz '96 system, and it is shown that the approach yields models that are stable and that give both improved skill in initialised predictions and better long-term climate statistics. Improvements in long-term statistics are smaller than for single time-step tendencies, however, indicating that it would be valuable to develop methods that target improvements on longer time scales. Future strategies for the development of this approach and possible applications to making progress on important scientific problems are discussed.
△ Less
Submitted 24 April, 2019;
originally announced April 2019.
-
Provenance and data differencing for workflow reproducibility analysis
Authors:
Paolo Missier,
Simon Woodman,
Hugo Hiden,
Paul Watson
Abstract:
One of the foundations of science is that researchers must publish the methodology used to achieve their results so that others can attempt to reproduce them. This has the added benefit of allowing methods to be adopted and adapted for other purposes. In the field of e-Science, services -- often choreographed through workflow, process data to generate results. The reproduction of results is often…
▽ More
One of the foundations of science is that researchers must publish the methodology used to achieve their results so that others can attempt to reproduce them. This has the added benefit of allowing methods to be adopted and adapted for other purposes. In the field of e-Science, services -- often choreographed through workflow, process data to generate results. The reproduction of results is often not straightforward as the computational objects may not be made available or may have been updated since the results were generated. For example, services are often updated to fix bugs or improve algorithms. This paper addresses these problems in three ways. Firstly, it introduces a new framework to clarify the range of meanings of "reproducibility". Secondly, it describes a new algorithm, \PDIFF, that uses a comparison of workflow provenance traces to determine whether an experiment has been reproduced; the main innovation is that if this is not the case then the specific point(s) of divergence are identified through graph analysis, assisting any researcher wishing to understand those differences. One key feature is support for user-defined, semantic data comparison operators. Finally, the paper describes an implementation of \PDIFF that leverages the power of the e-Science Central platform which enacts workflows in the cloud. As well as automatically generating a provenance trace for consumption by \PDIFF, the platform supports the storage and re-use of old versions of workflows, data and services; the paper shows how this can be powerfully exploited in order to achieve reproduction and re-use.
△ Less
Submitted 3 June, 2014;
originally announced June 2014.
-
A Flow Sensitive Security Model for Cloud Computing Systems
Authors:
Wen Zeng,
Chunyan Mu,
Maciej Koutny,
Paul Watson
Abstract:
The extent and importance of cloud computing is rapidly increasing due to the ever increasing demand for internet services and communications. Instead of building individual information technology infrastructure to host databases or software, a third party can host them in its large server clouds. Large organizations may wish to keep sensitive information on their more restricted servers rather th…
▽ More
The extent and importance of cloud computing is rapidly increasing due to the ever increasing demand for internet services and communications. Instead of building individual information technology infrastructure to host databases or software, a third party can host them in its large server clouds. Large organizations may wish to keep sensitive information on their more restricted servers rather than in the public cloud. This has led to the introduction of federated cloud computing (FCC) in which both public and private cloud computing resources are used.
△ Less
Submitted 30 April, 2014;
originally announced April 2014.
-
Probing the Repulsive Core of the Nucleon-Nucleon Interaction via the 4He(e,e'pN) Triple-Coincidence Reaction
Authors:
I. Korover N. Muangma,
O. Hen,
R. Shneor,
V. Sulkosky,
A. Kelleher,
S. Gilad,
D. W. Higinbotham,
E. Piasetzky J. Watson,
S. Wood,
Abdurahim Rakhman,
P. Aguilera,
Z. Ahmed,
H. Albataineh,
K. Allada,
B. Anderson,
D. Anez,
K. Aniol,
J. Annand,
W. Armstrong,
J. Arrington,
T. Averett,
T. Badman,
H. Baghdasaryan,
X. Bai,
A. Beck
, et al. (93 additional authors not shown)
Abstract:
We studied simultaneously the 4He(e,e'p), 4He(e,e'pp), and 4He(e,e'pn) reactions at Q^2=2 [GeV/c]2 and x_B>1, for a (e,e'p) missing-momentum range of 400 to 830 MeV/c. The knocked-out proton was detected in coincidence with a proton or neutron recoiling almost back to back to the missing momentum, leaving the residual A=2 system at low excitation energy. These data were used to identify two-nucleo…
▽ More
We studied simultaneously the 4He(e,e'p), 4He(e,e'pp), and 4He(e,e'pn) reactions at Q^2=2 [GeV/c]2 and x_B>1, for a (e,e'p) missing-momentum range of 400 to 830 MeV/c. The knocked-out proton was detected in coincidence with a proton or neutron recoiling almost back to back to the missing momentum, leaving the residual A=2 system at low excitation energy. These data were used to identify two-nucleon short-range correlated pairs and to deduce their isospin structure as a function of missing momentum in a region where the nucleon-nucleon force is expected to change from predominantly tensor to repulsive. Neutron-proton pairs dominate the high-momentum tail of the nucleon momentum distributions, but their abundance is reduced as the nucleon momentum increases beyond ~500 MeV/c. The extracted fraction of proton-proton pairs is small and almost independent of the missing momentum in the range we studied. Our data are compared with ab-initio calculations of two-nucleon momentum distributions in 4He.
△ Less
Submitted 11 July, 2014; v1 submitted 23 January, 2014;
originally announced January 2014.
-
Symmetries of quasiplatonic Riemann surfaces
Authors:
Gareth A. Jones,
David Singerman,
Paul D. Watson
Abstract:
We state and prove a corrected version of a theorem of Singerman, which relates the existence of symmetries (anticonformal involutions) of a quasiplatonic Riemann surface $\mathcal S$ (one uniformised by a normal subgroup $N$ of finite index in a cocompact triangle group $Δ$) to the properties of the group $G=Δ/N$. We give examples to illustrate the revised necessary and sufficient conditions for…
▽ More
We state and prove a corrected version of a theorem of Singerman, which relates the existence of symmetries (anticonformal involutions) of a quasiplatonic Riemann surface $\mathcal S$ (one uniformised by a normal subgroup $N$ of finite index in a cocompact triangle group $Δ$) to the properties of the group $G=Δ/N$. We give examples to illustrate the revised necessary and sufficient conditions for the existence of symmetries, and we relate them to properties of the associated dessins d'enfants, or hypermaps.
△ Less
Submitted 11 January, 2014;
originally announced January 2014.
-
Quarks and gluons in a magnetic field
Authors:
Peter Watson,
Hugo Reinhardt
Abstract:
The quark gap equation under the rainbow truncation, with two versions of a phenomenological one-gluon exchange interaction and in the presence of a uniform magnetic field is considered. It is argued that in order to describe the quark condensate in the limit of vanishing magnetic fields, one must sum over the Landau levels. The resulting chiral quark condensate rises quadratically for small magne…
▽ More
The quark gap equation under the rainbow truncation, with two versions of a phenomenological one-gluon exchange interaction and in the presence of a uniform magnetic field is considered. It is argued that in order to describe the quark condensate in the limit of vanishing magnetic fields, one must sum over the Landau levels. The resulting chiral quark condensate rises quadratically for small magnetic fields and linearly for large fields, in qualitative agreement with various recent lattice results. It is observed that when discussing quarks, the magnitude of the magnetic field must be considered relative to the scale of the strong interaction.
△ Less
Submitted 28 November, 2013;
originally announced November 2013.
-
Quark gap equation in an external magnetic field
Authors:
Peter Watson,
Hugo Reinhardt
Abstract:
The nonperturbative quark gap equation under the rainbow truncation and with two versions of a phenomenological one-gluon exchange interaction is studied in the presence of a uniform external magnetic field, with emphasis on the small field limit. The chiral quark condensate, magnetic moment and susceptibility are calculated and compared to recent lattice data.
The nonperturbative quark gap equation under the rainbow truncation and with two versions of a phenomenological one-gluon exchange interaction is studied in the presence of a uniform external magnetic field, with emphasis on the small field limit. The chiral quark condensate, magnetic moment and susceptibility are calculated and compared to recent lattice data.
△ Less
Submitted 22 October, 2013;
originally announced October 2013.
-
A Search for Astrophysical Burst Signals at the Sudbury Neutrino Observatory
Authors:
B. Aharmim,
S. N. Ahmed,
A. E. Anthony,
N. Barros,
E. W. Beier,
A. Bellerive,
B. Beltran,
M. Bergevin,
S. D. Biller,
K. Boudjemline,
M. G. Boulay,
B. Cai,
Y. D. Chan,
D. Chauhan,
M. Chen,
B. T. Cleveland,
G. A. Cox,
X. Dai,
H. Deng,
J. A. Detwiler,
M. DiMarco,
M. D. Diamond,
P. J. Doe,
G. Doucas,
P. -L. Drouin
, et al. (102 additional authors not shown)
Abstract:
The Sudbury Neutrino Observatory (SNO) has confirmed the standard solar model and neutrino oscillations through the observation of neutrinos from the solar core. In this paper we present a search for neutrinos associated with sources other than the solar core, such as gamma-ray bursters and solar flares. We present a new method for looking for temporal coincidences between neutrino events and astr…
▽ More
The Sudbury Neutrino Observatory (SNO) has confirmed the standard solar model and neutrino oscillations through the observation of neutrinos from the solar core. In this paper we present a search for neutrinos associated with sources other than the solar core, such as gamma-ray bursters and solar flares. We present a new method for looking for temporal coincidences between neutrino events and astrophysical bursts of widely varying intensity. No correlations were found between neutrinos detected in SNO and such astrophysical sources.
△ Less
Submitted 4 September, 2013;
originally announced September 2013.
-
Surface Morphology and Phase Stability of Titanium Foils Irradiated by 136 MeV 136Xe
Authors:
S. Sadi,
A. Paulenova,
W. Loveland,
P. R. Watson,
J. P. Greene,
S. Zhu,
G. Zinkann
Abstract:
A stack of titanium foils was irradiated with 136 MeV 136Xe to study microstructure damage and phase stability of titanium upon irradiation. X- ray diffraction, scanning electron microscopy/energy dispersive spectroscopy and atomic force microscopy were used to study the resulting microstructure damage and phase stability of titanium. We observed the phase transfor- mation of polycrystalline titan…
▽ More
A stack of titanium foils was irradiated with 136 MeV 136Xe to study microstructure damage and phase stability of titanium upon irradiation. X- ray diffraction, scanning electron microscopy/energy dispersive spectroscopy and atomic force microscopy were used to study the resulting microstructure damage and phase stability of titanium. We observed the phase transfor- mation of polycrystalline titanium from alpha-Ti (hexagonally closed packed (hcp)) to face centered cubic (fcc) after irradiation with 2.2 x 1015 ions/cm2. Irradiation of Ti with 1.8 x 1014-2.2 x 1015 ions/cm2 resulted in the forma- tion of voids, hillocks, dislocation loops, dislocation lines, as well as polygonal ridge networks.
△ Less
Submitted 7 January, 2013;
originally announced January 2013.
-
Gap and Bethe-Salpeter equations in Coulomb gauge
Authors:
Peter Watson,
Hugo Reinhardt
Abstract:
I will discuss the gap and Bethe-Salpeter equations for quantum chromodynamics in Coulomb gauge under a leading order truncation scheme. Within this scheme, the role of charge conservation and the cancellation of infrared divergences becomes particularly clear. The quark gap equation exhibits not only chiral symmetry breaking, but explicitly reproduces the heavy quark limit. The corresponding gluo…
▽ More
I will discuss the gap and Bethe-Salpeter equations for quantum chromodynamics in Coulomb gauge under a leading order truncation scheme. Within this scheme, the role of charge conservation and the cancellation of infrared divergences becomes particularly clear. The quark gap equation exhibits not only chiral symmetry breaking, but explicitly reproduces the heavy quark limit. The corresponding gluon equation has a massive solution with interesting nonperturbative renormalization properties. I will further discuss various aspects of the Bethe-Salpeter equation for pseudoscalar and vector mesons with arbitrary quark masses and present numerical results for the meson masses and leptonic decay constants.
△ Less
Submitted 19 December, 2012;
originally announced December 2012.
-
Bethe-Salpeter equation at leading order in Coulomb gauge
Authors:
Peter Watson,
Hugo Reinhardt
Abstract:
The Bethe-Salpeter equation and leptonic decay constants for pseudoscalar and vector quark-antiquark mesons with arbitrary quark masses are studied in Coulomb gauge, under a leading order truncation. As input, we use a pure linear rising potential, supplemented by a contact term arising from the conservation of total color charge. It is shown how the equations can be written in terms of manifestly…
▽ More
The Bethe-Salpeter equation and leptonic decay constants for pseudoscalar and vector quark-antiquark mesons with arbitrary quark masses are studied in Coulomb gauge, under a leading order truncation. As input, we use a pure linear rising potential, supplemented by a contact term arising from the conservation of total color charge. It is shown how the equations can be written in terms of manifestly finite functions, despite the infrared singular interaction. The resulting equations are solved numerically. Both the pattern of dynamical chiral symmetry breaking and the leading order heavy quark limit are visible.
△ Less
Submitted 19 November, 2012;
originally announced November 2012.
-
Secure & Rapid Composition of Infrastructure Services in the Cloud
Authors:
Pierre de Leusse,
Panos Periorellis,
Paul Watson,
Andreas Maierhofer
Abstract:
A fundamental ambition of grid and distributed systems is to be capable of sustaining evolution and allowing for adaptability ((F. Losavio et al., 2002), (S. Radhakrishnan, 2005)). Furthermore, as the complexity and sophistication of theses structures increases, so does the need for adaptability of each component. One of the primary benefits of service oriented architecture (SOA) is the ability to…
▽ More
A fundamental ambition of grid and distributed systems is to be capable of sustaining evolution and allowing for adaptability ((F. Losavio et al., 2002), (S. Radhakrishnan, 2005)). Furthermore, as the complexity and sophistication of theses structures increases, so does the need for adaptability of each component. One of the primary benefits of service oriented architecture (SOA) is the ability to compose applications, processes or more complex services from other services which increases the capacity for adaptation. This document proposes a novel infrastructure composition model that aims at increasing the adaptability of the capabilities exposed through it by dynamically managing their non functional requirements.
△ Less
Submitted 2 March, 2012;
originally announced March 2012.
-
Nonperturbative study of the four-point heavy quark Green's functions in Coulomb gauge
Authors:
Carina Popovici,
Peter Watson,
Hugo Reinhardt
Abstract:
The heavy quark sector of Coulomb gauge QCD is investigated, by making a heavy quark mass expansion of the QCD action and restricting to the leading order. With the truncation of the Yang-Mills sector to include only dressed two-point functions, we study the Dyson-Schwinger equations for the four-point quark Green's functions (proper and amputated). In this limit, we provide an exact solution for…
▽ More
The heavy quark sector of Coulomb gauge QCD is investigated, by making a heavy quark mass expansion of the QCD action and restricting to the leading order. With the truncation of the Yang-Mills sector to include only dressed two-point functions, we study the Dyson-Schwinger equations for the four-point quark Green's functions (proper and amputated). In this limit, we provide an exact solution for the four-point quark Green's functions and show that the corresponding poles relate to the bound state energy of the heavy quark system. Moreover, a natural separation between the physical and unphysical poles in the Green's functions emerges.
△ Less
Submitted 21 December, 2011;
originally announced December 2011.
-
Leading order QCD in Coulomb gauge
Authors:
Peter Watson,
Hugo Reinhardt
Abstract:
Coulomb gauge QCD in the first order formalism can be written in terms of a ghost-free, nonlocal action that ensures total color charge conservation via Gauss' law. Making an Ansatz whereby the nonlocal term (the Coulomb kernel) is replaced by its expectation value, the resulting Dyson-Schwinger equations can be derived. With a leading order truncation, these equations reduce to the gap equations…
▽ More
Coulomb gauge QCD in the first order formalism can be written in terms of a ghost-free, nonlocal action that ensures total color charge conservation via Gauss' law. Making an Ansatz whereby the nonlocal term (the Coulomb kernel) is replaced by its expectation value, the resulting Dyson-Schwinger equations can be derived. With a leading order truncation, these equations reduce to the gap equations for the static gluon and quark propagators obtained from a quasi-particle approximation to the canonical Hamiltonian approach. Moreover a connection to the heavy quark limit can be established, allowing an intuitive explanation for the charge constraint and infrared divergences.
△ Less
Submitted 30 November, 2011;
originally announced November 2011.
-
Leading order infrared quantum chromodynamics in Coulomb gauge
Authors:
Peter Watson,
Hugo Reinhardt
Abstract:
A truncation scheme for the Dyson-Schwinger equations of quantum chromodynamics in Coulomb gauge within the first order formalism is presented. The truncation is based on an Ansatz for the Coulomb kernel occurring in the action. Results at leading loop order and in the infrared are discussed for both the Yang-Mills and quark sectors. It is found that the resulting equations for the static gluon an…
▽ More
A truncation scheme for the Dyson-Schwinger equations of quantum chromodynamics in Coulomb gauge within the first order formalism is presented. The truncation is based on an Ansatz for the Coulomb kernel occurring in the action. Results at leading loop order and in the infrared are discussed for both the Yang-Mills and quark sectors. It is found that the resulting equations for the static gluon and quark propagators agree with those derived in a quasi-particle approximation to the canonical Hamiltonian approach. Moreover, a connection to the heavy quark limit is established. The equations are analyzed numerically and it is seen that in both the gluonic and quark sectors, a nontrivial dynamical infrared mass scale emerges.
△ Less
Submitted 25 November, 2011;
originally announced November 2011.
-
Linking Search Space Structure, Run-Time Dynamics, and Problem Difficulty: A Step Toward Demystifying Tabu Search
Authors:
A. E. Howe,
J. P. Watson,
L. D. Whitley
Abstract:
Tabu search is one of the most effective heuristics for locating high-quality solutions to a diverse array of NP-hard combinatorial optimization problems. Despite the widespread success of tabu search, researchers have a poor understanding of many key theoretical aspects of this algorithm, including models of the high-level run-time dynamics and identification of those search space features that i…
▽ More
Tabu search is one of the most effective heuristics for locating high-quality solutions to a diverse array of NP-hard combinatorial optimization problems. Despite the widespread success of tabu search, researchers have a poor understanding of many key theoretical aspects of this algorithm, including models of the high-level run-time dynamics and identification of those search space features that influence problem difficulty. We consider these questions in the context of the job-shop scheduling problem (JSP), a domain where tabu search algorithms have been shown to be remarkably effective. Previously, we demonstrated that the mean distance between random local optima and the nearest optimal solution is highly correlated with problem difficulty for a well-known tabu search algorithm for the JSP introduced by Taillard. In this paper, we discuss various shortcomings of this measure and develop a new model of problem difficulty that corrects these deficiencies. We show that Taillards algorithm can be modeled with high fidelity as a simple variant of a straightforward random walk. The random walk model accounts for nearly all of the variability in the cost required to locate both optimal and sub-optimal solutions to random JSPs, and provides an explanation for differences in the difficulty of random versus structured JSPs. Finally, we discuss and empirically substantiate two novel predictions regarding tabu search algorithm behavior. First, the method for constructing the initial solution is highly unlikely to impact the performance of tabu search. Second, tabu tenure should be selected to be as small as possible while simultaneously avoiding search stagnation; values larger than necessary lead to significant degradations in performance.
△ Less
Submitted 11 September, 2011;
originally announced September 2011.
-
Combined Analysis of all Three Phases of Solar Neutrino Data from the Sudbury Neutrino Observatory
Authors:
SNO Collaboration,
B. Aharmim,
S. N. Ahmed,
A. E. Anthony,
N. Barros,
E. W. Beier,
A. Bellerive,
B. Beltran,
M. Bergevin,
S. D. Biller,
K. Boudjemline,
M. G. Boulay,
B. Cai,
Y. D. Chan,
D. Chauhan,
M. Chen,
B. T. Cleveland,
G. A. Cox,
X. Dai,
H. Deng,
J. A. Detwiler,
M. DiMarco,
P. J. Doe,
G. Doucas,
P. -L. Drouin
, et al. (99 additional authors not shown)
Abstract:
We report results from a combined analysis of solar neutrino data from all phases of the Sudbury Neutrino Observatory. By exploiting particle identification information obtained from the proportional counters installed during the third phase, this analysis improved background rejection in that phase of the experiment. The combined analysis resulted in a total flux of active neutrino flavors from 8…
▽ More
We report results from a combined analysis of solar neutrino data from all phases of the Sudbury Neutrino Observatory. By exploiting particle identification information obtained from the proportional counters installed during the third phase, this analysis improved background rejection in that phase of the experiment. The combined analysis resulted in a total flux of active neutrino flavors from 8B decays in the Sun of (5.25 \pm 0.16(stat.)+0.11-0.13(syst.))\times10^6 cm^{-2}s^{-1}. A two-flavor neutrino oscillation analysis yielded \Deltam^2_{21} = (5.6^{+1.9}_{-1.4})\times10^{-5} eV^2 and tan^2θ_{12}= 0.427^{+0.033}_{-0.029}. A three-flavor neutrino oscillation analysis combining this result with results of all other solar neutrino experiments and the KamLAND experiment yielded \Deltam^2_{21} = (7.41^{+0.21}_{-0.19})\times10^{-5} eV^2, tan^2θ_{12} = 0.446^{+0.030}_{-0.029}, and sin^2θ_{13} = (2.5^{+1.8}_{-1.5})\times10^{-2}. This implied an upper bound of sin^2θ_{13} < 0.053 at the 95% confidence level (C.L.).
△ Less
Submitted 4 September, 2011;
originally announced September 2011.
-
Higher order heavy quark Green's functions in Coulomb gauge
Authors:
C. Popovici,
P. Watson,
H. Reinhardt
Abstract:
The Dyson-Schwinger equation for the 4-point quark Green's functions is studied. In the limit of the heavy quark mass and with the truncation to include only the dressed two point functions for the Yang-Mills sector, we provide an exact solution for the 4-point quark Green's functions, in both quark-antiquark and diquark channels, and show that the corresponding poles relate to the bound state ene…
▽ More
The Dyson-Schwinger equation for the 4-point quark Green's functions is studied. In the limit of the heavy quark mass and with the truncation to include only the dressed two point functions for the Yang-Mills sector, we provide an exact solution for the 4-point quark Green's functions, in both quark-antiquark and diquark channels, and show that the corresponding poles relate to the bound state energy of the heavy quark systems. Moreover, a natural separation between physical and unphysical poles in the Green's functions emerges.
△ Less
Submitted 24 March, 2011;
originally announced March 2011.
-
Heavy quarks, gluons and the confinement potential in Coulomb gauge
Authors:
Carina Popovici,
Peter Watson,
Hugo Reinhardt
Abstract:
We consider the heavy quark limit of Coulomb gauge QCD, with the truncation of the Yang-Mills sector to include only (dressed) two-point functions. We find that the rainbow-ladder approximation to the gap and Bethe-Salpeter equations is nonperturbatively exact and moreover, we provide a direct connection between the temporal gluon propagator and the quark confinement potential. Further, we show th…
▽ More
We consider the heavy quark limit of Coulomb gauge QCD, with the truncation of the Yang-Mills sector to include only (dressed) two-point functions. We find that the rainbow-ladder approximation to the gap and Bethe-Salpeter equations is nonperturbatively exact and moreover, we provide a direct connection between the temporal gluon propagator and the quark confinement potential. Further, we show that only bound states of color singlet quark-antiquark (meson) and quark-quark (SU(2) baryon) pairs are physically allowed.
△ Less
Submitted 9 November, 2010;
originally announced November 2010.
-
The ghost propagator in Coulomb gauge
Authors:
P. Watson,
H. Reinhardt
Abstract:
We present results for a numerical study of the ghost propagator in Coulomb gauge whereby lattice results for the spatial gluon propagator are used as input to solving the ghost Dyson-Schwinger equation. We show that in order to solve completely, the ghost equation must be supplemented by a boundary condition (the value of the inverse ghost propagator dressing function at zero momentum) which dete…
▽ More
We present results for a numerical study of the ghost propagator in Coulomb gauge whereby lattice results for the spatial gluon propagator are used as input to solving the ghost Dyson-Schwinger equation. We show that in order to solve completely, the ghost equation must be supplemented by a boundary condition (the value of the inverse ghost propagator dressing function at zero momentum) which determines if the solution is critical (zero value for the boundary condition) or subcritical (finite value). The various solutions exhibit a characteristic behavior where all curves follow the same (critical) solution when going from high to low momenta until `forced' to freeze out in the infrared to the value of the boundary condition. The boundary condition can be interpreted in terms of the Gribov gauge-fixing ambiguity; we also demonstrate that this is not connected to the renormalization. Further, the connection to the temporal gluon propagator and the infrared slavery picture of confinement is discussed.
△ Less
Submitted 9 November, 2010;
originally announced November 2010.
-
Three-quark confinement potential from the Faddeev equation
Authors:
C. Popovici,
P. Watson,
H. Reinhardt
Abstract:
In the heavy quark limit of Coulomb gauge QCD and by truncating the Yang-Mills sector to include only dressed two-point functions, an analytic nonperturbative solution to the Faddeev equation for three-quark bound states in the case of equal quark separations is presented. A direct connection between the temporal gluon propagator and the three-quark confinement potential is provided and it is show…
▽ More
In the heavy quark limit of Coulomb gauge QCD and by truncating the Yang-Mills sector to include only dressed two-point functions, an analytic nonperturbative solution to the Faddeev equation for three-quark bound states in the case of equal quark separations is presented. A direct connection between the temporal gluon propagator and the three-quark confinement potential is provided and it is shown that only color singlet qqq (baryon) states are physically allowed.
△ Less
Submitted 20 October, 2010;
originally announced October 2010.
-
Phase space spinor amplitudes for spin 1/2 systems
Authors:
P Watson,
A J Bracken
Abstract:
The concept of phase space amplitudes for systems with continuous degrees of freedom is generalized to finite-dimensional spin systems. Complex amplitudes are obtained on both a sphere and a finite lattice, in each case enabling a more fundamental description of pure spin states than that previously given by Wigner functions. In each case the Wigner function can be expressed as the star product of…
▽ More
The concept of phase space amplitudes for systems with continuous degrees of freedom is generalized to finite-dimensional spin systems. Complex amplitudes are obtained on both a sphere and a finite lattice, in each case enabling a more fundamental description of pure spin states than that previously given by Wigner functions. In each case the Wigner function can be expressed as the star product of the amplitude and its conjugate, so providing a generalized Born interpretation of amplitudes that emphasizes their more fundamental status. The ordinary product of the amplitude and its conjugate produces a (generalized) spin Husimi function. The case of spin-$\half$ is treated in detail, and it is shown that phase space amplitudes on the sphere transform correctly as spinors under under rotations, despite their expression in terms of spherical harmonics. Spin amplitudes on a lattice are also found to transform as spinors. Applications are given to the phase space description of state superposition, and to the evolution in phase space of the state of a spin-$\half$ magnetic dipole in a time-dependent magnetic field.
△ Less
Submitted 6 March, 2011; v1 submitted 18 October, 2010;
originally announced October 2010.
-
The Coulomb gauge ghost Dyson-Schwinger equation
Authors:
Peter Watson,
Hugo Reinhardt
Abstract:
A numerical study of the ghost Dyson-Schwinger equation in Coulomb gauge is performed and solutions for the ghost propagator found. As input, lattice results for the spatial gluon propagator are used. It is shown that in order to solve completely, the equation must be supplemented by a nonperturbative boundary condition (the value of the inverse ghost propagator dressing function at zero momentum)…
▽ More
A numerical study of the ghost Dyson-Schwinger equation in Coulomb gauge is performed and solutions for the ghost propagator found. As input, lattice results for the spatial gluon propagator are used. It is shown that in order to solve completely, the equation must be supplemented by a nonperturbative boundary condition (the value of the inverse ghost propagator dressing function at zero momentum) which determines if the solution is critical (zero value for the boundary condition) or subcritical (finite value). The various solutions exhibit a characteristic behavior where all curves follow the same (critical) solution when going from high to low momenta until `forced' to freeze out in the infrared to the value of the boundary condition. The renormalization is shown to be largely independent of the boundary condition. The boundary condition and the pattern of the solutions can be interpreted in terms of the Gribov gauge-fixing ambiguity. The connection to the temporal gluon propagator and the infrared slavery picture of confinement is explored.
△ Less
Submitted 21 November, 2010; v1 submitted 15 July, 2010;
originally announced July 2010.
-
The quantum state vector in phase space and Gabor's windowed Fourier transform
Authors:
A. J. Bracken,
P. Watson
Abstract:
Representations of quantum state vectors by complex phase space amplitudes, complementing the description of the density operator by the Wigner function, have been defined by applying the Weyl-Wigner transform to dyadic operators, linear in the state vector and anti-linear in a fixed `window state vector'. Here aspects of this construction are explored, with emphasis on the connection with Gabor's…
▽ More
Representations of quantum state vectors by complex phase space amplitudes, complementing the description of the density operator by the Wigner function, have been defined by applying the Weyl-Wigner transform to dyadic operators, linear in the state vector and anti-linear in a fixed `window state vector'. Here aspects of this construction are explored, with emphasis on the connection with Gabor's `windowed Fourier transform'. The amplitudes that arise for simple quantum states from various choices of window are presented as illustrations. Generalized Bargmann representations of the state vector appear as special cases, associated with Gaussian windows. For every choice of window, amplitudes lie in a corresponding linear subspace of square-integrable functions on phase space. A generalized Born interpretation of amplitudes is described, with both the Wigner function and a generalized Husimi function appearing as quantities linear in an amplitude and anti-linear in its complex conjugate. Schrödinger's time-dependent and time-independent equations are represented on phase space amplitudes, and their solutions described in simple cases.
△ Less
Submitted 13 October, 2010; v1 submitted 26 May, 2010;
originally announced June 2010.
-
Coulomb gauge confinement in the heavy quark limit
Authors:
C. Popovici,
P. Watson,
H. Reinhardt
Abstract:
The relationship between the nonperturbative Green's functions of Yang-Mills theory and the confinement potential is investigated. By rewriting the generating functional of quantum chromodynamics in terms of a heavy quark mass expansion in Coulomb gauge, restricting to leading order in this expansion and considering only the two-point functions of the Yang-Mills sector, the rainbow-ladder approxim…
▽ More
The relationship between the nonperturbative Green's functions of Yang-Mills theory and the confinement potential is investigated. By rewriting the generating functional of quantum chromodynamics in terms of a heavy quark mass expansion in Coulomb gauge, restricting to leading order in this expansion and considering only the two-point functions of the Yang-Mills sector, the rainbow-ladder approximation to the gap and Bethe-Salpeter equations is shown to be exact in this case and an analytic, nonperturbative solution is presented. It is found that there is a direct connection between the string tension and the temporal gluon propagator. Further, it is shown that for the 4-point quark correlation functions, only confined bound states of color-singlet quark-antiquark (meson) and quark-quark (baryon) pairs exist.
△ Less
Submitted 19 March, 2010;
originally announced March 2010.
-
Slavnov-Taylor Identities in Coulomb Gauge Yang-Mills Theory
Authors:
Peter Watson,
Hugo Reinhardt
Abstract:
Two aspects of the color charge in Coulomb gauge continuum Yang-Mills theory are discussed. The first aspect is the existence of a conserved and vanishing total charge exhibited within the first order functional formalism. The second aspect is the closure of the set of Slavnov-Taylor identities in the second order functional formalism, such that the exact solution for temporal Green's functions…
▽ More
Two aspects of the color charge in Coulomb gauge continuum Yang-Mills theory are discussed. The first aspect is the existence of a conserved and vanishing total charge exhibited within the first order functional formalism. The second aspect is the closure of the set of Slavnov-Taylor identities in the second order functional formalism, such that the exact solution for temporal Green's functions is in principle possible and thereby preserving the color charge.
△ Less
Submitted 11 December, 2008;
originally announced December 2008.
-
Slavnov-Taylor identities in Coulomb gauge Yang-Mills theory
Authors:
P. Watson,
H. Reinhardt
Abstract:
The Slavnov-Taylor identities of Coulomb gauge Yang-Mills theory are derived from the (standard, second order) functional formalism. It is shown how these identities form closed sets from which one can in principle fully determine the Green's functions involving the temporal component of the gauge field without approximation, given appropriate input.
The Slavnov-Taylor identities of Coulomb gauge Yang-Mills theory are derived from the (standard, second order) functional formalism. It is shown how these identities form closed sets from which one can in principle fully determine the Green's functions involving the temporal component of the gauge field without approximation, given appropriate input.
△ Less
Submitted 10 December, 2008;
originally announced December 2008.
-
Quarks in Coulomb gauge perturbation theory
Authors:
C. Popovici,
P. Watson,
H. Reinhardt
Abstract:
Coulomb gauge quantum chromodynamics within the first order functional formalism is considered. The quark contributions to the Dyson-Schwinger equations are derived and one-loop perturbative results for the two-point functions are presented.
Coulomb gauge quantum chromodynamics within the first order functional formalism is considered. The quark contributions to the Dyson-Schwinger equations are derived and one-loop perturbative results for the two-point functions are presented.
△ Less
Submitted 27 October, 2008;
originally announced October 2008.
-
Resolving temporal Gribov copies in Coulomb gauge Yang-Mills theory
Authors:
Hugo Reinhardt,
Peter Watson
Abstract:
The continuum Yang-Mills functional integral within the first order formalism and in Coulomb gauge is studied. In particular, the temporal zero-modes of the Faddeev-Popov operator are explicitly accounted for. It is shown that the treatment of these zero-modes results in the constraint that the total color charge of the system vanishes at all times. Further, it is argued that the functional inte…
▽ More
The continuum Yang-Mills functional integral within the first order formalism and in Coulomb gauge is studied. In particular, the temporal zero-modes of the Faddeev-Popov operator are explicitly accounted for. It is shown that the treatment of these zero-modes results in the constraint that the total color charge of the system vanishes at all times. Further, it is argued that the functional integral is effectively fully gauge-fixed once Gauss' law has been resolved in Coulomb gauge.
△ Less
Submitted 18 August, 2008;
originally announced August 2008.
-
Completing Continuum Coulomb Gauge in the Functional Formalism
Authors:
P. Watson,
H. Reinhardt
Abstract:
It is argued that within the continuum functional formalism, there is no need to supply a further (spatially independent) gauge constraint to complete the Coulomb gauge of Yang-Mills theory. It is shown explicitly that a natural completion of the gauge-fixing leads to a contradiction with the perturbative renormalizability of the theory.
It is argued that within the continuum functional formalism, there is no need to supply a further (spatially independent) gauge constraint to complete the Coulomb gauge of Yang-Mills theory. It is shown explicitly that a natural completion of the gauge-fixing leads to a contradiction with the perturbative renormalizability of the theory.
△ Less
Submitted 19 November, 2007;
originally announced November 2007.
-
Two-Point Functions of Coulomb Gauge Yang-Mills Theory
Authors:
Peter Watson,
Hugo Reinhardt
Abstract:
The functional approach to Coulomb gauge Yang-Mills theory is considered within the standard, second order, formalism. The Dyson-Schwinger equations and Slavnov-Taylor identities concerning the two-point functions are derived explicitly and one-loop perturbative results are presented.
The functional approach to Coulomb gauge Yang-Mills theory is considered within the standard, second order, formalism. The Dyson-Schwinger equations and Slavnov-Taylor identities concerning the two-point functions are derived explicitly and one-loop perturbative results are presented.
△ Less
Submitted 25 September, 2007;
originally announced September 2007.