Showing 1–50 of 149 results for author: Schilling, K

Field-of-View Extension for Diffusion MRI via Deep Generative Models

Authors: Chenyu Gao, Shunxing Bao, Michael Kim, Nancy Newlin, Praitayini Kanakaraj, Tianyuan Yao, Gaurav Rudravaram, Yuankai Huo, Daniel Moyer, Kurt Schilling, Walter Kukull, Arthur Toga, Derek Archer, Timothy Hohman, Bennett Landman, Zhiyuan Li

Abstract: Purpose: In diffusion MRI (dMRI), the volumetric and bundle analyses of whole-brain tissue microstructure and connectivity can be severely impeded by an incomplete field-of-view (FOV). This work aims to develop a method for imputing the missing slices directly from existing dMRI scans with an incomplete FOV. We hypothesize that the imputed image with complete FOV can improve the whole-brain tracto… ▽ More Purpose: In diffusion MRI (dMRI), the volumetric and bundle analyses of whole-brain tissue microstructure and connectivity can be severely impeded by an incomplete field-of-view (FOV). This work aims to develop a method for imputing the missing slices directly from existing dMRI scans with an incomplete FOV. We hypothesize that the imputed image with complete FOV can improve the whole-brain tractography for corrupted data with incomplete FOV. Therefore, our approach provides a desirable alternative to discarding the valuable dMRI data, enabling subsequent tractography analyses that would otherwise be challenging or unattainable with corrupted data. Approach: We propose a framework based on a deep generative model that estimates the absent brain regions in dMRI scans with incomplete FOV. The model is capable of learning both the diffusion characteristics in diffusion-weighted images (DWI) and the anatomical features evident in the corresponding structural images for efficiently imputing missing slices of DWI outside of incomplete FOV. Results: For evaluating the imputed slices, on the WRAP dataset the proposed framework achieved PSNRb0=22.397, SSIMb0=0.905, PSNRb1300=22.479, SSIMb1300=0.893; on the NACC dataset it achieved PSNRb0=21.304, SSIMb0=0.892, PSNRb1300=21.599, SSIMb1300= 0.877. The proposed framework improved the tractography accuracy, as demonstrated by an increased average Dice score for 72 tracts (p < 0.001) on both the WRAP and NACC datasets. Conclusions: Results suggest that the proposed framework achieved sufficient imputation performance in dMRI data with incomplete FOV for improving whole-brain tractography, thereby repairing the corrupted data. Our approach achieved more accurate whole-brain tractography results with extended and complete FOV and reduced the uncertainty when analyzing bundles associated with Alzheimer's Disease. △ Less

Submitted 6 May, 2024; originally announced May 2024.

Comments: 20 pages, 11 figures

Tractography with T1-weighted MRI and associated anatomical constraints on clinical quality diffusion MRI

Authors: Tian Yu, Yunhe Li, Michael E. Kim, Chenyu Gao, Qi Yang, Leon Y. Cai, Susane M. Resnick, Lori L. Beason-Held, Daniel C. Moyer, Kurt G. Schilling, Bennett A. Landman

Abstract: Diffusion MRI (dMRI) streamline tractography, the gold standard for in vivo estimation of brain white matter (WM) pathways, has long been considered indicative of macroscopic relationships with WM microstructure. However, recent advances in tractography demonstrated that convolutional recurrent neural networks (CoRNN) trained with a teacher-student framework have the ability to learn and propagate… ▽ More Diffusion MRI (dMRI) streamline tractography, the gold standard for in vivo estimation of brain white matter (WM) pathways, has long been considered indicative of macroscopic relationships with WM microstructure. However, recent advances in tractography demonstrated that convolutional recurrent neural networks (CoRNN) trained with a teacher-student framework have the ability to learn and propagate streamlines directly from T1 and anatomical contexts. Training for this network has previously relied on high-resolution dMRI. In this paper, we generalize the training mechanism to traditional clinical resolution data, which allows generalizability across sensitive and susceptible study populations. We train CoRNN on a small subset of the Baltimore Longitudinal Study of Aging (BLSA), which better resembles clinical protocols. Then, we define a metric, termed the epsilon ball seeding method, to compare T1 tractography and traditional diffusion tractography at the streamline level. Under this metric, T1 tractography generated by CoRNN reproduces diffusion tractography with approximately two millimeters of error. △ Less

Submitted 27 March, 2024; originally announced March 2024.

AI-enabled Cyber-Physical In-Orbit Factory -- AI approaches based on digital twin technology for robotic small satellite production

Authors: Florian Leutert, David Bohlig, Florian Kempf, Klaus Schilling, Maximilian Mühlbauer, Bengisu Ayan, Thomas Hulin, Freek Stulp, Alin Albu-Schäffer, Vladimir Kutscher, Christian Plesker, Thomas Dasbach, Stephan Damm, Reiner Anderl, Benjamin Schleich

Abstract: With the ever increasing number of active satellites in space, the rising demand for larger formations of small satellites and the commercialization of the space industry (so-called New Space), the realization of manufacturing processes in orbit comes closer to reality. Reducing launch costs and risks, allowing for faster on-demand deployment of individually configured satellites as well as the pr… ▽ More With the ever increasing number of active satellites in space, the rising demand for larger formations of small satellites and the commercialization of the space industry (so-called New Space), the realization of manufacturing processes in orbit comes closer to reality. Reducing launch costs and risks, allowing for faster on-demand deployment of individually configured satellites as well as the prospect for possible on-orbit servicing for satellites makes the idea of realizing an in-orbit factory promising. In this paper, we present a novel approach to an in-orbit factory of small satellites covering a digital process twin, AI-based fault detection, and teleoperated robot-control, which are being researched as part of the "AI-enabled Cyber-Physical In-Orbit Factory" project. In addition to the integration of modern automation and Industry 4.0 production approaches, the question of how artificial intelligence (AI) and learning approaches can be used to make the production process more robust, fault-tolerant and autonomous is addressed. This lays the foundation for a later realisation of satellite production in space in the form of an in-orbit factory. Central aspect is the development of a robotic AIT (Assembly, Integration and Testing) system where a small satellite could be assembled by a manipulator robot from modular subsystems. Approaches developed to improving this production process with AI include employing neural networks for optical and electrical fault detection of components. Force sensitive measuring and motion training helps to deal with uncertainties and tolerances during assembly. An AI-guided teleoperated control of the robot arm allows for human intervention while a Digital Process Twin represents process data and provides supervision during the whole production process. Approaches and results towards automated satellite production are presented in detail. △ Less

Submitted 5 February, 2024; v1 submitted 31 January, 2024; originally announced January 2024.

Journal ref: Acta Astronautica (2024), vol. 217, page 1-17

Evaluation of Mean Shift, ComBat, and CycleGAN for Harmonizing Brain Connectivity Matrices Across Sites

Authors: Hanliang Xu, Nancy R. Newlin, Michael E. Kim, Chenyu Gao, Praitayini Kanakaraj, Aravind R. Krishnan, Lucas W. Remedios, Nazirah Mohd Khairi, Kimberly Pechman, Derek Archer, Timothy J. Hohman, Angela L. Jefferson, The BIOCARD Study Team, Ivana Isgum, Yuankai Huo, Daniel Moyer, Kurt G. Schilling, Bennett A. Landman

Abstract: Connectivity matrices derived from diffusion MRI (dMRI) provide an interpretable and generalizable way of understanding the human brain connectome. However, dMRI suffers from inter-site and between-scanner variation, which impedes analysis across datasets to improve robustness and reproducibility of results. To evaluate different harmonization approaches on connectivity matrices, we compared graph… ▽ More Connectivity matrices derived from diffusion MRI (dMRI) provide an interpretable and generalizable way of understanding the human brain connectome. However, dMRI suffers from inter-site and between-scanner variation, which impedes analysis across datasets to improve robustness and reproducibility of results. To evaluate different harmonization approaches on connectivity matrices, we compared graph measures derived from these matrices before and after applying three harmonization techniques: mean shift, ComBat, and CycleGAN. The sample comprises 168 age-matched, sex-matched normal subjects from two studies: the Vanderbilt Memory and Aging Project (VMAP) and the Biomarkers of Cognitive Decline Among Normal Individuals (BIOCARD). First, we plotted the graph measures and used coefficient of variation (CoV) and the Mann-Whitney U test to evaluate different methods' effectiveness in removing site effects on the matrices and the derived graph measures. ComBat effectively eliminated site effects for global efficiency and modularity and outperformed the other two methods. However, all methods exhibited poor performance when harmonizing average betweenness centrality. Second, we tested whether our harmonization methods preserved correlations between age and graph measures. All methods except for CycleGAN in one direction improved correlations between age and global efficiency and between age and modularity from insignificant to significant with p-values less than 0.05. △ Less

Submitted 24 January, 2024; v1 submitted 8 January, 2024; originally announced January 2024.

Comments: 11 pages, 5 figures, to be published in SPIE Medical Imaging 2024: Image Processing

Predicting Age from White Matter Diffusivity with Residual Learning

Authors: Chenyu Gao, Michael E. Kim, Ho Hin Lee, Qi Yang, Nazirah Mohd Khairi, Praitayini Kanakaraj, Nancy R. Newlin, Derek B. Archer, Angela L. Jefferson, Warren D. Taylor, Brian D. Boyd, Lori L. Beason-Held, Susan M. Resnick, The BIOCARD Study Team, Yuankai Huo, Katherine D. Van Schaik, Kurt G. Schilling, Daniel Moyer, Ivana Išgum, Bennett A. Landman

Abstract: Imaging findings inconsistent with those expected at specific chronological age ranges may serve as early indicators of neurological disorders and increased mortality risk. Estimation of chronological age, and deviations from expected results, from structural MRI data has become an important task for developing biomarkers that are sensitive to such deviations. Complementary to structural analysis,… ▽ More Imaging findings inconsistent with those expected at specific chronological age ranges may serve as early indicators of neurological disorders and increased mortality risk. Estimation of chronological age, and deviations from expected results, from structural MRI data has become an important task for developing biomarkers that are sensitive to such deviations. Complementary to structural analysis, diffusion tensor imaging (DTI) has proven effective in identifying age-related microstructural changes within the brain white matter, thereby presenting itself as a promising additional modality for brain age prediction. Although early studies have sought to harness DTI's advantages for age estimation, there is no evidence that the success of this prediction is owed to the unique microstructural and diffusivity features that DTI provides, rather than the macrostructural features that are also available in DTI data. Therefore, we seek to develop white-matter-specific age estimation to capture deviations from normal white matter aging. Specifically, we deliberately disregard the macrostructural information when predicting age from DTI scalar images, using two distinct methods. The first method relies on extracting only microstructural features from regions of interest. The second applies 3D residual neural networks (ResNets) to learn features directly from the images, which are non-linearly registered and warped to a template to minimize macrostructural variations. When tested on unseen data, the first method yields mean absolute error (MAE) of 6.11 years for cognitively normal participants and MAE of 6.62 years for cognitively impaired participants, while the second method achieves MAE of 4.69 years for cognitively normal participants and MAE of 4.96 years for cognitively impaired participants. We find that the ResNet model captures subtler, non-macrostructural features for brain age prediction. △ Less

Submitted 21 January, 2024; v1 submitted 6 November, 2023; originally announced November 2023.

Comments: SPIE Medical Imaging: Image Processing. San Diego, CA. February 2024 (accepted as poster presentation)

BundleSeg: A versatile, reliable and reproducible approach to white matter bundle segmentation

Authors: Etienne St-Onge, Kurt G Schilling, Francois Rheault

Abstract: This work presents BundleSeg, a reliable, reproducible, and fast method for extracting white matter pathways. The proposed method combines an iterative registration procedure with a recently developed precise streamline search algorithm that enables efficient segmentation of streamlines without the need for tractogram clustering or simplifying assumptions. We show that BundleSeg achieves improved… ▽ More This work presents BundleSeg, a reliable, reproducible, and fast method for extracting white matter pathways. The proposed method combines an iterative registration procedure with a recently developed precise streamline search algorithm that enables efficient segmentation of streamlines without the need for tractogram clustering or simplifying assumptions. We show that BundleSeg achieves improved repeatability and reproducibility than state-of-the-art segmentation methods, with significant speed improvements. The enhanced precision and reduced variability in extracting white matter connections offer a valuable tool for neuroinformatic studies, increasing the sensitivity and specificity of tractography-based studies of white matter pathways. △ Less

Submitted 21 August, 2023; originally announced August 2023.

Robust Fiber ODF Estimation Using Deep Constrained Spherical Deconvolution for Diffusion MRI

Authors: Tianyuan Yao, Francois Rheault, Leon Y Cai, Vishwesh nath, Zuhayr Asad, Nancy Newlin, Can Cui, Ruining Deng, Karthik Ramadass, Andrea Shafer, Susan Resnick, Kurt Schilling, Bennett A. Landman, Yuankai Huo

Abstract: Diffusion-weighted magnetic resonance imaging (DW-MRI) is a critical imaging method for capturing and modeling tissue microarchitecture at a millimeter scale. A common practice to model the measured DW-MRI signal is via fiber orientation distribution function (fODF). This function is the essential first step for the downstream tractography and connectivity analyses. With recent advantages in data… ▽ More Diffusion-weighted magnetic resonance imaging (DW-MRI) is a critical imaging method for capturing and modeling tissue microarchitecture at a millimeter scale. A common practice to model the measured DW-MRI signal is via fiber orientation distribution function (fODF). This function is the essential first step for the downstream tractography and connectivity analyses. With recent advantages in data sharing, large-scale multi-site DW-MRI datasets are being made available for multi-site studies. However, measurement variabilities (e.g., inter- and intra-site variability, hardware performance, and sequence design) are inevitable during the acquisition of DW-MRI. Most existing model-based methods (e.g., constrained spherical deconvolution (CSD)) and learning based methods (e.g., deep learning (DL)) do not explicitly consider such variabilities in fODF modeling, which consequently leads to inferior performance on multi-site and/or longitudinal diffusion studies. In this paper, we propose a novel data-driven deep constrained spherical deconvolution method to explicitly constrain the scan-rescan variabilities for a more reproducible and robust estimation of brain microstructure from repeated DW-MRI scans. Specifically, the proposed method introduces a new 3D volumetric scanner-invariant regularization scheme during the fODF estimation. We study the Human Connectome Project (HCP) young adults test-retest group as well as the MASiVar dataset (with inter- and intra-site scan/rescan data). The Baltimore Longitudinal Study of Aging (BLSA) dataset is employed for external validation. From the experimental results, the proposed data-driven framework outperforms the existing benchmarks in repeated fODF estimation. The proposed method is assessing the downstream connectivity analysis and shows increased performance in distinguishing subjects with different biomarkers. △ Less

Submitted 5 June, 2023; originally announced June 2023.

Comments: 33 pages, 7 figures

A Unified Learning Model for Estimating Fiber Orientation Distribution Functions on Heterogeneous Multi-shell Diffusion-weighted MRI

Authors: Tianyuan Yao, Nancy Newlin, Praitayini Kanakaraj, Vishwesh nath, Leon Y Cai, Karthik Ramadass, Kurt Schilling, Bennett A. Landman, Yuankai Huo

Abstract: Diffusion-weighted (DW) MRI measures the direction and scale of the local diffusion process in every voxel through its spectrum in q-space, typically acquired in one or more shells. Recent developments in micro-structure imaging and multi-tissue decomposition have sparked renewed attention to the radial b-value dependence of the signal. Applications in tissue classification and micro-architecture… ▽ More Diffusion-weighted (DW) MRI measures the direction and scale of the local diffusion process in every voxel through its spectrum in q-space, typically acquired in one or more shells. Recent developments in micro-structure imaging and multi-tissue decomposition have sparked renewed attention to the radial b-value dependence of the signal. Applications in tissue classification and micro-architecture estimation, therefore, require a signal representation that extends over the radial as well as angular domain. Multiple approaches have been proposed that can model the non-linear relationship between the DW-MRI signal and biological microstructure. In the past few years, many deep learning-based methods have been developed towards faster inference speed and higher inter-scan consistency compared with traditional model-based methods (e.g., multi-shell multi-tissue constrained spherical deconvolution). However, a multi-stage learning strategy is typically required since the learning process relies on various middle representations, such as simple harmonic oscillator reconstruction (SHORE) representation. In this work, we present a unified dynamic network with a single-stage spherical convolutional neural network, which allows efficient fiber orientation distribution function (fODF) estimation through heterogeneous multi-shell diffusion MRI sequences. We study the Human Connectome Project (HCP) young adults with test-retest scans. From the experimental results, the proposed single-stage method outperforms prior multi-stage approaches in repeated fODF estimation with shell dropoff and single-shell DW-MRI sequences. △ Less

Submitted 29 January, 2024; v1 submitted 28 March, 2023; originally announced March 2023.

Biomedical image analysis competitions: The state of current participation practice

Authors: Matthias Eisenmann, Annika Reinke, Vivienn Weru, Minu Dietlinde Tizabi, Fabian Isensee, Tim J. Adler, Patrick Godau, Veronika Cheplygina, Michal Kozubek, Sharib Ali, Anubha Gupta, Jan Kybic, Alison Noble, Carlos Ortiz de Solórzano, Samiksha Pachade, Caroline Petitjean, Daniel Sage, Donglai Wei, Elizabeth Wilden, Deepak Alapatt, Vincent Andrearczyk, Ujjwal Baid, Spyridon Bakas, Niranjan Balu, Sophia Bano , et al. (331 additional authors not shown)

Abstract: The number of international benchmarking competitions is steadily increasing in various fields of machine learning (ML) research and practice. So far, however, little is known about the common practice as well as bottlenecks faced by the community in tackling the research questions posed. To shed light on the status quo of algorithm development in the specific field of biomedical imaging analysis,… ▽ More The number of international benchmarking competitions is steadily increasing in various fields of machine learning (ML) research and practice. So far, however, little is known about the common practice as well as bottlenecks faced by the community in tackling the research questions posed. To shed light on the status quo of algorithm development in the specific field of biomedical imaging analysis, we designed an international survey that was issued to all participants of challenges conducted in conjunction with the IEEE ISBI 2021 and MICCAI 2021 conferences (80 competitions in total). The survey covered participants' expertise and working environments, their chosen strategies, as well as algorithm characteristics. A median of 72% challenge participants took part in the survey. According to our results, knowledge exchange was the primary incentive (70%) for participation, while the reception of prize money played only a minor role (16%). While a median of 80 working hours was spent on method development, a large portion of participants stated that they did not have enough time for method development (32%). 25% perceived the infrastructure to be a bottleneck. Overall, 94% of all solutions were deep learning-based. Of these, 84% were based on standard architectures. 43% of the respondents reported that the data samples (e.g., images) were too large to be processed at once. This was most commonly addressed by patch-based training (69%), downsampling (37%), and solving 3D analysis tasks as a series of 2D tasks. K-fold cross-validation on the training set was performed by only 37% of the participants and only 50% of the participants performed ensembling based on multiple identical models (61%) or heterogeneous models (39%). 48% of the respondents applied postprocessing steps. △ Less

Submitted 12 September, 2023; v1 submitted 16 December, 2022; originally announced December 2022.

Towards Wireless Ranging and Synchronization using CubeSat Software-Defined Radio Subsystems

Authors: Markus Gardill, Dominik Pearson, Julian Scharnagl, Klaus Schilling

Abstract: An approach towards wireless ranging and synchronization using commercial of-the-shelf software-defined radio payloads for small satellites, esp. CubeSats is studied. The approach only relies on the programmable logic configuration and processing system software. No hardware modifications or additions to the payloads are necessary. Experimental evaluation of the initial implementation shows a stan… ▽ More An approach towards wireless ranging and synchronization using commercial of-the-shelf software-defined radio payloads for small satellites, esp. CubeSats is studied. The approach only relies on the programmable logic configuration and processing system software. No hardware modifications or additions to the payloads are necessary. Experimental evaluation of the initial implementation shows a standard deviation of time-of-flight-based ranging measurements in the order of 1 cm, which renders the concept of interest for distributed satellite system missions. △ Less

Submitted 1 December, 2022; originally announced December 2022.

Recommendations and guidelines from the ISMRM Diffusion Study Group for preclinical diffusion MRI: Part 2 -- Ex vivo imaging

Authors: Kurt G Schilling, Francesco Grussu, Andrada Ianus, Brian Hansen, Rachel L C Barrett, Manisha Aggarwal, Stijn Michielse, Fatima Nasrallah, Warda Syeda, Nian Wang, Jelle Veraart, Alard Roebroeck, Andrew F Bagdasarian, Cornelius Eichner, Farshid Sepehrband, Jan Zimmermann, Lucas Soustelle, Christien Bowman, Benjamin C Tendler, Andreea Hertanu, Ben Jeurissen, Lucio Frydman, Yohan van de Looij, David Hike, Jeff F Dunn , et al. (31 additional authors not shown)

Abstract: The value of preclinical diffusion MRI (dMRI) is substantial. While dMRI enables in vivo non-invasive characterization of tissue, ex vivo dMRI is increasingly being used to probe tissue microstructure and brain connectivity. Ex vivo dMRI has several experimental advantages including higher signal-to-noise ratio and spatial resolution compared to in vivo studies, and more advanced diffusion contras… ▽ More The value of preclinical diffusion MRI (dMRI) is substantial. While dMRI enables in vivo non-invasive characterization of tissue, ex vivo dMRI is increasingly being used to probe tissue microstructure and brain connectivity. Ex vivo dMRI has several experimental advantages including higher signal-to-noise ratio and spatial resolution compared to in vivo studies, and more advanced diffusion contrasts for improved microstructure and connectivity characterization. Another major advantage is direct comparison with histological data as a crucial methodological validation. However, there are a number of considerations that must be made when performing ex vivo experiments. The steps from tissue preparation, image acquisition and processing, and interpretation of results are complex, with many decisions that not only differ dramatically from in vivo imaging, but ultimately affect what questions can be answered using the data. This work represents 'Part 2' of a series of recommendations and considerations for preclinical dMRI, where we focus on best practices for dMRI of ex vivo tissue. We first describe the value that ex vivo imaging adds to the field of dMRI, followed by general considerations and foundational knowledge that must be considered when designing experiments. We then give guidelines for ex vivo protocols, including tissue preparation, imaging sequences and data processing including pre-processing, model-fitting, and tractography. Finally, we provide an online resource which lists publicly available ex vivo dMRI datasets and dedicated software packages. In each section, we attempt to provide guidelines and recommendations, but also highlight areas for which no guidelines exist, and where future work should lie. An overarching goal herein is to enhance the rigor and reproducibility of ex vivo dMRI acquisitions and analyses, and thereby advance biomedical knowledge. △ Less

Submitted 7 February, 2023; v1 submitted 27 September, 2022; originally announced September 2022.

Comments: 59 pages, 12 figures, part of ongoing efforts on ISMRM Diffusion Study Group initiative 'Best Practices (Consensus) for diffusion MRI'. arXiv admin note: text overlap with arXiv:2209.12994

Recommendations and guidelines from the ISMRM Diffusion Study Group for preclinical diffusion MRI: Part 1 -- In vivo small-animal imaging

Authors: Ileana O Jelescu, Francesco Grussu, Andrada Ianus, Brian Hansen, Rachel L C Barrett, Manisha Aggarwal, Stijn Michielse, Fatima Nasrallah, Warda Syeda, Nian Wang, Jelle Veraart, Alard Roebroeck, Andrew F Bagdasarian, Cornelius Eichner, Farshid Sepehrband, Jan Zimmermann, Lucas Soustelle, Christien Bowman, Benjamin C Tendler, Andreea Hertanu, Ben Jeurissen, Marleen Verhoye, Lucio Frydman, Yohan van de Looij, David Hike , et al. (32 additional authors not shown)

Abstract: The value of in vivo preclinical diffusion MRI (dMRI) is substantial. Small-animal dMRI has been used for methodological development and validation, characterizing the biological basis of diffusion phenomena, and comparative anatomy. Many of the influential works in this field were first performed in small animals or ex vivo samples. The steps from animal setup and monitoring, to acquisition, anal… ▽ More The value of in vivo preclinical diffusion MRI (dMRI) is substantial. Small-animal dMRI has been used for methodological development and validation, characterizing the biological basis of diffusion phenomena, and comparative anatomy. Many of the influential works in this field were first performed in small animals or ex vivo samples. The steps from animal setup and monitoring, to acquisition, analysis, and interpretation are complex, with many decisions that may ultimately affect what questions can be answered using the data. This work aims to serve as a reference, presenting selected recommendations and guidelines from the diffusion community, on best practices for preclinical dMRI of in vivo animals. In each section, we also highlight areas for which no guidelines exist (and why), and where future work should focus. We first describe the value that small animal imaging adds to the field of dMRI, followed by general considerations and foundational knowledge that must be considered when designing experiments. We briefly describe differences in animal species and disease models and discuss how they are appropriate for different studies. We then give guidelines for in vivo acquisition protocols, including decisions on hardware, animal preparation, imaging sequences and data processing, including pre-processing, model-fitting, and tractography. Finally, we provide an online resource which lists publicly available preclinical dMRI datasets and software packages, to promote responsible and reproducible research. An overarching goal herein is to enhance the rigor and reproducibility of small animal dMRI acquisitions and analyses, and thereby advance biomedical knowledge. △ Less

Submitted 21 April, 2023; v1 submitted 26 September, 2022; originally announced September 2022.

Comments: 69 pages, 6 figures, 1 table

Joint analysis of structural connectivity and cortical surface features: correlates with mild traumatic brain injury

Authors: Cailey I. Kerley, Leon Y. Cai, Chang Yu, Logan M. Crawford, Jason M. Elenberger, Eden S. Singh, Kurt G. Schilling, Katherine S. Aboud, Bennett A. Landman, Tonia S. Rex

Abstract: Mild traumatic brain injury (mTBI) is a complex syndrome that affects up to 600 per 100,000 individuals, with a particular concentration among military personnel. About half of all mTBI patients experience a diverse array of chronic symptoms which persist long after the acute injury. Hence, there is an urgent need for better understanding of the white matter and gray matter pathologies associated… ▽ More Mild traumatic brain injury (mTBI) is a complex syndrome that affects up to 600 per 100,000 individuals, with a particular concentration among military personnel. About half of all mTBI patients experience a diverse array of chronic symptoms which persist long after the acute injury. Hence, there is an urgent need for better understanding of the white matter and gray matter pathologies associated with mTBI to map which specific brain systems are impacted and identify courses of intervention. Previous works have linked mTBI to disruptions in white matter pathways and cortical surface abnormalities. Herein, we examine these hypothesized links in an exploratory study of joint structural connectivity and cortical surface changes associated with mTBI and its chronic symptoms. Briefly, we consider a cohort of 12 mTBI and 26 control subjects. A set of 588 cortical surface metrics and 4,753 structural connectivity metrics were extracted from cortical surface regions and diffusion weighted magnetic resonance imaging in each subject. Principal component analysis (PCA) was used to reduce the dimensionality of each metric set. We then applied independent component analysis (ICA) both to each PCA space individually and together in a joint ICA approach. We identified a stable independent component across the connectivity-only and joint ICAs which presented significant group differences in subject loadings (p<0.05, corrected). Additionally, we found that two mTBI symptoms, slowed thinking and forgetfulness, were significantly correlated (p<0.05, corrected) with mTBI subject loadings in a surface-only ICA. These surface-only loadings captured an increase in bilateral cortical thickness. △ Less

Submitted 15 December, 2020; v1 submitted 18 November, 2020; originally announced December 2020.

Comments: To be published in Proc SPIE Int Soc Opt Eng. 2021 Feb

Deep Learning Estimation of Multi-Tissue Constrained Spherical Deconvolution with Limited Single Shell DW-MRI

Authors: Vishwesh Nath, Sudhir K. Pathak, Kurt G. Schilling, Walt Schneider, Bennett A. Landman

Abstract: Diffusion-weighted magnetic resonance imaging (DW-MRI) is the only non-invasive approach for estimation of intra-voxel tissue microarchitecture and reconstruction of in vivo neural pathways for the human brain. With improvement in accelerated MRI acquisition technologies, DW-MRI protocols that make use of multiple levels of diffusion sensitization have gained popularity. A well-known advanced meth… ▽ More Diffusion-weighted magnetic resonance imaging (DW-MRI) is the only non-invasive approach for estimation of intra-voxel tissue microarchitecture and reconstruction of in vivo neural pathways for the human brain. With improvement in accelerated MRI acquisition technologies, DW-MRI protocols that make use of multiple levels of diffusion sensitization have gained popularity. A well-known advanced method for reconstruction of white matter microstructure that uses multi-shell data is multi-tissue constrained spherical deconvolution (MT-CSD). MT-CSD substantially improves the resolution of intra-voxel structure over the traditional single shell version, constrained spherical deconvolution (CSD). Herein, we explore the possibility of using deep learning on single shell data (using the b=1000 s/mm2 from the Human Connectome Project (HCP)) to estimate the information content captured by 8th order MT-CSD using the full three shell data (b=1000, 2000, and 3000 s/mm2 from HCP). Briefly, we examine two network architectures: 1.) Sequential network of fully connected dense layers with a residual block in the middle (ResDNN), 2.) Patch based convolutional neural network with a residual block (ResCNN). For both networks an additional output block for estimation of voxel fraction was used with a modified loss function. Each approach was compared against the baseline of using MT-CSD on all data on 15 subjects from the HCP divided into 5 training, 2 validation, and 8 testing subjects with a total of 6.7 million voxels. The fiber orientation distribution function (fODF) can be recovered with high correlation (0.77 vs 0.74 and 0.65) as compared to the ground truth of MT-CST, which was derived from the multi-shell DW-MRI acquisitions. Source code and models have been made publicly available. △ Less

Submitted 20 February, 2020; originally announced February 2020.

Comments: 10 pages, 7 figures

MRI correlates of chronic symptoms in mild traumatic brain injury

Authors: Cailey I. Kerley, Kurt G. Schilling, Justin Blaber, Beth Miller, Allen Newton, Adam W. Anderson, Bennett A. Landman, Tonia S. Rex

Abstract: Veterans with mild traumatic brain injury (mTBI) have reported auditory and visual dysfunction that persists beyond the acute incident. The etiology behind these symptoms is difficult to characterize with current clinical imaging. These functional deficits may be caused by shear injury or micro-bleeds, which can be detected with special imaging modalities. We explore these hypotheses in a pilot st… ▽ More Veterans with mild traumatic brain injury (mTBI) have reported auditory and visual dysfunction that persists beyond the acute incident. The etiology behind these symptoms is difficult to characterize with current clinical imaging. These functional deficits may be caused by shear injury or micro-bleeds, which can be detected with special imaging modalities. We explore these hypotheses in a pilot study of multi-parametric MRI. We extract over 1,000 imaging and clinical metrics and project them to a low-dimensional space, where we can discriminate between healthy controls and patients with mTBI. We also show correlations between the metric representations and patient symptoms. △ Less

Submitted 22 June, 2020; v1 submitted 6 December, 2019; originally announced December 2019.

Comments: SPIE Medical Imaging 2020

Deep Learning Captures More Accurate Diffusion Fiber Orientations Distributions than Constrained Spherical Deconvolution

Authors: Vishwesh Nath, Kurt G. Schilling, Colin B. Hansen, Prasanna Parvathaneni, Allison E. Hainline, Camilo Bermudez, Andrew J. Plassard, Vaibhav Janve, Yurui Gao, Justin A. Blaber, Iwona Stępniewska, Adam W. Anderson, Bennett A. Landman

Abstract: Confocal histology provides an opportunity to establish intra-voxel fiber orientation distributions that can be used to quantitatively assess the biological relevance of diffusion weighted MRI models, e.g., constrained spherical deconvolution (CSD). Here, we apply deep learning to investigate the potential of single shell diffusion weighted MRI to explain histologically observed fiber orientation… ▽ More Confocal histology provides an opportunity to establish intra-voxel fiber orientation distributions that can be used to quantitatively assess the biological relevance of diffusion weighted MRI models, e.g., constrained spherical deconvolution (CSD). Here, we apply deep learning to investigate the potential of single shell diffusion weighted MRI to explain histologically observed fiber orientation distributions (FOD) and compare the derived deep learning model with a leading CSD approach. This study (1) demonstrates that there exists additional information in the diffusion signal that is not currently exploited by CSD, and (2) provides an illustrative data-driven model that makes use of this information. △ Less

Submitted 13 November, 2019; originally announced November 2019.

Comments: 2 pages, 4 figures. This work was accepted and published as an abstract at ISMRM 2018 held in Paris, France

Enabling Multi-Shell b-Value Generalizability of Data-Driven Diffusion Models with Deep SHORE

Authors: Vishwesh Nath, Ilwoo Lyu, Kurt G. Schilling, Prasanna Parvathaneni, Colin B. Hansen, Yucheng Tang, Yuankai Huo, Vaibhav A. Janve, Yurui Gao, Iwona Stepniewska, Adam W. Anderson, Bennett A. Landman

Abstract: Intra-voxel models of the diffusion signal are essential for interpreting organization of the tissue environment at micrometer level with data at millimeter resolution. Recent advances in data driven methods have enabled direct compari-son and optimization of methods for in-vivo data with externally validated histological sections with both 2-D and 3-D histology. Yet, all existing methods make lim… ▽ More Intra-voxel models of the diffusion signal are essential for interpreting organization of the tissue environment at micrometer level with data at millimeter resolution. Recent advances in data driven methods have enabled direct compari-son and optimization of methods for in-vivo data with externally validated histological sections with both 2-D and 3-D histology. Yet, all existing methods make limiting assumptions of either (1) model-based linkages between b-values or (2) limited associations with single shell data. We generalize prior deep learning models that used single shell spherical harmonic transforms to integrate the re-cently developed simple harmonic oscillator reconstruction (SHORE) basis. To enable learning on the SHORE manifold, we present an alternative formulation of the fiber orientation distribution (FOD) object using the SHORE basis while rep-resenting the observed diffusion weighted data in the SHORE basis. To ensure consistency of hyper-parameter optimization for SHORE, we present our Deep SHORE approach to learn on a data-optimized manifold. Deep SHORE is evalu-ated with eight-fold cross-validation of a preclinical MRI-histology data with four b-values. Generalizability of in-vivo human data is evaluated on two separate 3T MRI scanners. Specificity in terms of angular correlation (ACC) with the preclinical data improved on single shell: 0.78 relative to 0.73 and 0.73, multi-shell: 0.80 relative to 0.74 (p < 0.001). In the in-vivo human data, Deep SHORE was more consistent across scanners with 0.63 relative to other multi-shell methods 0.39, 0.52 and 0.57 in terms of ACC. In conclusion, Deep SHORE is a promising method to enable data driven learning with DW-MRI under conditions with varying b-values, number of diffusion shells, and gradient directions per shell. △ Less

Submitted 22 February, 2020; v1 submitted 14 July, 2019; originally announced July 2019.

Novel data and a new parametrization of the electric dipole strength in nuclei with 88 < A < 116

Authors: Eckart Grosse, Arnd R. Junghans, Krasimir Kosev, Gencho Rusev, Klaus D. Schilling, Ronald Schwengner, Andreas Wagner

Abstract: A hitherto unexplored method for the experimental determination of the photon strength function up to the neutron separation energy was developed at the Radiation Source ELBE in Dresden. It was applied to various heavy nuclei, preferentially to nuclides with increasing distance to the N=50 neutron shell, and it covers the high level density excitation energy range above 4 MeV. The observed quasi-c… ▽ More A hitherto unexplored method for the experimental determination of the photon strength function up to the neutron separation energy was developed at the Radiation Source ELBE in Dresden. It was applied to various heavy nuclei, preferentially to nuclides with increasing distance to the N=50 neutron shell, and it covers the high level density excitation energy range above 4 MeV. The observed quasi-continuous spectra of scattered photons can be -- after a proper correction for multi-step processes -- directly combined to nuclear photo effect data from literature. A remarkably good match of the photon strengths as measured below and above the neutron emission threshold is observed. The wide energy coverage of the combined data forms an excellent basis to derive a parameterization for the dipole strength function fully covering the range across the nucleon separation energies. In addition to the parameters defining the deformation of the nuclear ground states only one additional constant is needed to describe the dipole strength in the nuclei with 88<A<116. The new parameterization differs significantly from the prescriptions generally used in network calculations, e.g. those of interest for the cosmic nucleo-synthesis. △ Less

Submitted 10 January, 2019; originally announced January 2019.

Comments: 12 pages, 6 figures,

Journal ref: Proceedings of Science(PSF)004 (2007)

Inter-Scanner Harmonization of High Angular Resolution DW-MRI using Null Space Deep Learning

Authors: Vishwesh Nath, Prasanna Parvathaneni, Colin B. Hansen, Allison E. Hainline, Camilo Bermudez, Samuel Remedios, Justin A. Blaber, Kurt G. Schilling, Ilwoo Lyu, Vaibhav Janve, Yurui Gao, Iwona Stepniewska, Baxter P. Rogers, Allen T. Newton, L. Taylor Davis, Jeff Luci, Adam W. Anderson, Bennett A. Landman

Abstract: Diffusion-weighted magnetic resonance imaging (DW-MRI) allows for non-invasive imaging of the local fiber architecture of the human brain at a millimetric scale. Multiple classical approaches have been proposed to detect both single (e.g., tensors) and multiple (e.g., constrained spherical deconvolution, CSD) fiber population orientations per voxel. However, existing techniques generally exhibit l… ▽ More Diffusion-weighted magnetic resonance imaging (DW-MRI) allows for non-invasive imaging of the local fiber architecture of the human brain at a millimetric scale. Multiple classical approaches have been proposed to detect both single (e.g., tensors) and multiple (e.g., constrained spherical deconvolution, CSD) fiber population orientations per voxel. However, existing techniques generally exhibit low reproducibility across MRI scanners. Herein, we propose a data-driven tech-nique using a neural network design which exploits two categories of data. First, training data were acquired on three squirrel monkey brains using ex-vivo DW-MRI and histology of the brain. Second, repeated scans of human subjects were acquired on two different scanners to augment the learning of the network pro-posed. To use these data, we propose a new network architecture, the null space deep network (NSDN), to simultaneously learn on traditional observed/truth pairs (e.g., MRI-histology voxels) along with repeated observations without a known truth (e.g., scan-rescan MRI). The NSDN was tested on twenty percent of the histology voxels that were kept completely blind to the network. NSDN significantly improved absolute performance relative to histology by 3.87% over CSD and 1.42% over a recently proposed deep neural network approach. More-over, it improved reproducibility on the paired data by 21.19% over CSD and 10.09% over a recently proposed deep approach. Finally, NSDN improved gen-eralizability of the model to a third in vivo human scanner (which was not used in training) by 16.08% over CSD and 10.41% over a recently proposed deep learn-ing approach. This work suggests that data-driven approaches for local fiber re-construction are more reproducible, informative and precise and offers a novel, practical method for determining these models. △ Less

Submitted 9 October, 2018; originally announced October 2018.

Comments: 10 pages, 5 figures

Nuclear deformation and neutron excess as competing effects for pygmy dipole strength

Authors: R. Massarczyk, R. Schwengner, F. Dönau, S. Frauendorf, M. Anders, D. Bemmerer, R. Beyer, C. Bhatia, E. Birgersson, M. Butterling, Z. Elekes, A. Ferrari, M. E. Gooden, R. Hannaske, A. R. Junghans, M. Kempe, J. H. Kelley, T. Kögler, A. Matic, M. L. Menzel, S. Müller, T. P. Reinhardt, M. Röder, G. Rusev, K. D. Schilling , et al. (5 additional authors not shown)

Abstract: The electromagnetic dipole strength below the neutron-separation energy has been studied for the xenon isotopes with mass numbers A = 124, 128, 132, and 134 in nuclear resonance fluorescence experiments using the ELBE bremsstrahlung facility at Helmholtz-Zentrum Dresden-Rossendorf and the HIgS facility at Triangle Universities Nuclear Laboratory Durham. The systematic study gained new information… ▽ More The electromagnetic dipole strength below the neutron-separation energy has been studied for the xenon isotopes with mass numbers A = 124, 128, 132, and 134 in nuclear resonance fluorescence experiments using the ELBE bremsstrahlung facility at Helmholtz-Zentrum Dresden-Rossendorf and the HIgS facility at Triangle Universities Nuclear Laboratory Durham. The systematic study gained new information about the influence of the neutron excess as well as of nuclear deformation on the strength in the region of the pygmy dipole resonance. The results are compared with those obtained for the chain of molybdenum isotopes and with predictions of a random-phase approximation in a deformed basis. It turned out that the effect of nuclear deformation plays a minor role compared with the one caused by neutron excess. A global parametrization of the strength in terms of neutron and proton numbers allowed us to derive a formula capable of predicting the summed E1 strengths in the pygmy region for a wide mass range of nuclides. △ Less

Submitted 8 November, 2013; originally announced November 2013.

Comments: 5 pages, subimtted to Physical Review Letters

Dipole strength in 144Sm studied via (gamma,n), (gamma,p) and (gamma,alpha) reactions

Authors: Chithra Nair, Arnd R. Junghans, Martin Erhard, Daniel Bemmerer, Roland Beyer, Eckart Grosse, Krasimir Kosev, Michele Marta, Gencho Rusev, Klaus-Dieter Schilling, Ronald Schwengner, Andreas Wagner

Abstract: Photoactivation measurements on 144Sm have been performed with bremsstrahlung endpoint energies from 10.0 to 15.5 MeV at the bremsstrahlung facility of the superconducting electron accelerator ELBE of Forschungszentrum Dresden-Rossendorf. The measured activation yield for the 144Sm(gamma,n) reaction is compared with the calculated yield using cross sections from previous photoneutron experiments.… ▽ More Photoactivation measurements on 144Sm have been performed with bremsstrahlung endpoint energies from 10.0 to 15.5 MeV at the bremsstrahlung facility of the superconducting electron accelerator ELBE of Forschungszentrum Dresden-Rossendorf. The measured activation yield for the 144Sm(gamma,n) reaction is compared with the calculated yield using cross sections from previous photoneutron experiments. The activation yields measured for all disintegration channels 144Sm(gamma,n), (gamma,p) and (gamma,alpha) are compared to the yield calculated by using Hauser-Feshbach statistical models. A new parametrization of the photon strength function is presented and the yield simulated by using the modified photon strength parameters are compared to the experimental data. △ Less

Submitted 12 May, 2010; originally announced May 2010.

Comments: Phys Rev C, Accepted for publication

Journal ref: Phys.Rev.C81:055806,2010

Photoactivation experiment on 197Au and its implications for the dipole strength in heavy nuclei

Authors: C. Nair, M. Erhard, A. R. Junghans, D. Bemmerer, R. Beyer, E. Grosse, J. Klug, K. Kosev, G. Rusev, K. D. Schilling, R. Schwengner, A. Wagner

Abstract: The 197Au(gamma,n) reaction is used as an activation standard for photodisintegration studies on astrophysically relevant nuclei. At the bremsstrahlung facility of the superconducting electron accelerator ELBE (Electron Linear accelerator of high Brilliance and low Emittance) of Forschungszentrum Dresden-Rossendorf, photoactivation measurements on 197Au have been performed with bremsstrahlung en… ▽ More The 197Au(gamma,n) reaction is used as an activation standard for photodisintegration studies on astrophysically relevant nuclei. At the bremsstrahlung facility of the superconducting electron accelerator ELBE (Electron Linear accelerator of high Brilliance and low Emittance) of Forschungszentrum Dresden-Rossendorf, photoactivation measurements on 197Au have been performed with bremsstrahlung endpoint energies from 8.0 to 15.5 MeV. The measured activation yield is compared with previous experiments as well as with calculations using Hauser-Feshbach statistical models. It is shown that the experimental data are best described by a two-Lorentzian parametrization with taking the axial deformation of 197Au into account. The experimental 197Au(gamma,n) reaction yield measured at ELBE via the photoactivation method is found to be consistent with previous experimental data using photon scattering or neutron detection methods. △ Less

Submitted 28 November, 2008; originally announced November 2008.

Comments: 9 pages

Journal ref: Phys.Rev.C78:055802,2008

Pygmy dipole strength in 90Zr

Authors: R. Schwengner, G. Rusev, N. Tsoneva, N. Benouaret, R. Beyer, M. Erhard, E. Grosse, A. R. Junghans, J. Klug, K. Kosev, H. Lenske, C. Nair, K. D. Schilling, A. Wagner

Abstract: The dipole response of the N=50 nucleus 90Zr was studied in photon-scattering experiments at the electron linear accelerator ELBE with bremsstrahlung produced at kinetic electron energies of 7.9, 9.0, and 13.2 MeV. We identified 189 levels up to an excitation energy of 12.9 MeV. Statistical methods were applied to estimate intensities of inelastic transitions and to correct the intensities of th… ▽ More The dipole response of the N=50 nucleus 90Zr was studied in photon-scattering experiments at the electron linear accelerator ELBE with bremsstrahlung produced at kinetic electron energies of 7.9, 9.0, and 13.2 MeV. We identified 189 levels up to an excitation energy of 12.9 MeV. Statistical methods were applied to estimate intensities of inelastic transitions and to correct the intensities of the ground-state transitions for their branching ratios. In this way we derived the photoabsorption cross section up to the neutron-separation energy. This cross section matches well the photoabsorption cross section obtained from (gamma,n) data and thus provides information about the extension of the dipole-strength distribution toward energies below the neutron-separation energy. An enhancement of E1 strength has been found in the range of 6 MeV to 11 MeV. Calculations within the framework of the quasiparticle-phonon model ascribe this strength to a vibration of the excessive neutrons against the N = Z neutron-proton core, giving rise to a pygmy dipole resonance. △ Less

Submitted 27 November, 2008; originally announced November 2008.

Comments: 27 pages, 15 figures, to be published in PRC

Journal ref: Phys.Rev.C78:064314,2008

Photodisintegration studies on p-nuclei: The case of Mo and Sm isotopes

Authors: C. Nair, A. R. Junghans, M. Erhard, D. Bemmerer, R. Beyer, P. Crespo, E. Grosse, M. Fauth, K. Kosev, G. Rusev, K. D. Schilling, R. Schwengner, A. Wagner

Abstract: In explosive stellar environments like supernovae, the temperatures are high enough for the production of heavy neutron-deficient nuclei, the socalled p-nuclei. Up to now, the knowledge of the reaction rates of p-nuclei is based on theoretical parameterizations using statistical model calculations. At the bremsstrahlung facility of the superconducting electron accelerator ELBE of FZ Dresden-Ross… ▽ More In explosive stellar environments like supernovae, the temperatures are high enough for the production of heavy neutron-deficient nuclei, the socalled p-nuclei. Up to now, the knowledge of the reaction rates of p-nuclei is based on theoretical parameterizations using statistical model calculations. At the bremsstrahlung facility of the superconducting electron accelerator ELBE of FZ Dresden-Rossendorf, we aim to measure the photodisintegration rates of heavy nuclei experimentally. Photoactivation measurements on the astrophysically relevant p-nuclei 92Mo and 144Sm have been performed with bremsstrahlung end-point energies from 10.0 to 16.5 MeV. First experiments on the short-lived decays following the reaction 144Sm(gamma,n) are carried out using a pneumatic delivery system for rapid transport of activated samples. The activation yields are compared with calculations using cross sections from recent Hauser-Feshbach models. △ Less

Submitted 26 October, 2007; originally announced October 2007.

Comments: 6 pages, 3 figures

Journal ref: J.Phys.G35:014036,2008

Photon strength distributions in stable even-even molybdenum isotopes

Authors: A. Wagner, R. Beyer, M. Erhard, E. Grosse, A. R. Junghans, J. Klug, K. Kosev, C. Nair, N. Nankov, G. Rusev, K. D. Schilling, R. Schwengner

Abstract: Electromagnetic dipole-strength distributions up to the particle separation energies are studied for the stable even-even nuclides $^{92,94,96,98,100}$Mo in photon scattering experiments at the superconducting electron accelerator ELBE of the Forschungszentrum Dresden-Rossendorf. The influence of inelastic transitions to low-lying excited states has been corrected by a simulation of $γ$ cascades… ▽ More Electromagnetic dipole-strength distributions up to the particle separation energies are studied for the stable even-even nuclides $^{92,94,96,98,100}$Mo in photon scattering experiments at the superconducting electron accelerator ELBE of the Forschungszentrum Dresden-Rossendorf. The influence of inelastic transitions to low-lying excited states has been corrected by a simulation of $γ$ cascades using a statistical model. After corrections for branching ratios of ground-state transitions, the photon-scattering cross-sections smoothly connect to data obtained from $(γ,n)$-reactions. With the newly determined electromagnetic dipole response of nuclei well below the particle separation energies the parametrisation of the isovector giant-dipole resonance is done with improved precision. △ Less

Submitted 26 July, 2007; originally announced July 2007.

Comments: Proceedings Nuclear Physics in Astrophysics 3, March 2007, Dresden Journal of Physics G, IOP Publishing

Journal ref: J.Phys.G35:014035,2008

Pygmy dipole strength close to particle-separation energies - the case of the Mo isotopes

Authors: Gencho Rusev, Eckart Grosse, Martin Erhard, Arnd Junghans, Krasimir Kosev, Klaus-Dieter Schilling, Ronald Schwengner, Andreas Wagner

Abstract: The distribution of electromagnetic dipole strength in 92, 98, 100 Mo has been investigated by photon scattering using bremsstrahlung from the new ELBE facility. The experimental data for well separated nuclear resonances indicate a transition from a regular to a chaotic behaviour above 4 MeV of excitation energy. As the strength distributions follow a Porter-Thomas distribution much of the dipo… ▽ More The distribution of electromagnetic dipole strength in 92, 98, 100 Mo has been investigated by photon scattering using bremsstrahlung from the new ELBE facility. The experimental data for well separated nuclear resonances indicate a transition from a regular to a chaotic behaviour above 4 MeV of excitation energy. As the strength distributions follow a Porter-Thomas distribution much of the dipole strength is found in weak and in unresolved resonances appearing as fluctuating cross section. An analysis of this quasi-continuum - here applied to nuclear resonance fluorescence in a novel way - delivers dipole strength functions, which are combining smoothly to those obtained from (g,n)-data. Enhancements at 6.5 MeV and at ~9 MeV are linked to the pygmy dipole resonances postulated to occur in heavy nuclei. △ Less

Submitted 20 December, 2005; originally announced December 2005.

Comments: 6 pages, 5 figures, proceedings Nuclear Physics in Astrophysics II, May 16-20, Debrecen, Hungary. The original publication is available at www.eurphysj.org

Photodissociation of p-process nuclei studied by bremsstrahlung induced activation

Authors: M. Erhard, A. R. Junghans, R. Beyer, E. Grosse, J. Klug, K. Kosev, C. Nair, N. Nankov, G. Rusev, K. D. Schilling, R. Schwengner, A. Wagner

Abstract: A research program has been started to study experimentally the near-threshold photodissociation of nuclides in the chain of cosmic heavy element production with bremsstrahlung from the ELBE accelerator. An important prerequisite for such studies is good knowledge of the bremsstrahlung distribution which was determined by measuring the photodissociation of the deuteron and by comparison with mod… ▽ More A research program has been started to study experimentally the near-threshold photodissociation of nuclides in the chain of cosmic heavy element production with bremsstrahlung from the ELBE accelerator. An important prerequisite for such studies is good knowledge of the bremsstrahlung distribution which was determined by measuring the photodissociation of the deuteron and by comparison with model calculations. First data were obtained for the astrophysically important target nucleus 92-Mo by observing the radioactive decay of the nuclides produced by bremsstrahlung irradiation at end-point energies between 11.8 MeV and 14.0 MeV. The results are compared to recent statistical model calculations. △ Less

Submitted 20 December, 2005; originally announced December 2005.

Comments: 6 pages, 8 figures, Proceedings Nuclear Physics in Astrophysics II, May 16-20, 2005, Debrecen, Hungary. The original publication is available at www.eurphysj.org

Journal ref: Eur.Phys.J. A27 (2006) 135-140

String breaking

Authors: Gunnar S Bali, Thomas Dussel, Thomas Lippert, Hartmut Neff, Zdravko Prkacin, Klaus Schilling

Abstract: We numerically investigate the transition of the static quark-antiquark string into a static-light meson-antimeson system. Improving noise reduction techniques, we are able to resolve the signature of string breaking dynamics for Nf=2 lattice QCD at zero temperature. We discuss the lattice techniques used and present results on energy levels and mixing angle of the static two-state system. We vi… ▽ More We numerically investigate the transition of the static quark-antiquark string into a static-light meson-antimeson system. Improving noise reduction techniques, we are able to resolve the signature of string breaking dynamics for Nf=2 lattice QCD at zero temperature. We discuss the lattice techniques used and present results on energy levels and mixing angle of the static two-state system. We visualize the action density distribution in the region of string breaking as a function of the static colour source-antisource separation. The results can be related to properties of quarkonium systems. △ Less

Submitted 14 December, 2005; originally announced December 2005.

Comments: 8 pages, Talk given at the Workshop on Computational Hadron Physics, Nicosia, Cyprus, 14--17 September 2005

Journal ref: Nucl.Phys.Proc.Suppl. 153 (2006) 9-16

Anatomy of String Breaking in QCD

Authors: Zdravko Prkacin, Gunnar S Bali, Thomas Dussel, Thomas Lippert, Hartmut Neff, Klaus Schilling

Abstract: We investigate the string breaking mechanism in n_f=2 QCD. We discuss the lattice techniques used and present results on energy levels and mixing angle of the static BBbar|QbarQ two-state system. The string breaking is visualized, by means of an animation of the action density distribution as a function of the static colour source-antisource separation. We investigate the string breaking mechanism in n_f=2 QCD. We discuss the lattice techniques used and present results on energy levels and mixing angle of the static BBbar|QbarQ two-state system. The string breaking is visualized, by means of an animation of the action density distribution as a function of the static colour source-antisource separation. △ Less

Submitted 7 October, 2005; originally announced October 2005.

Comments: Talk presented at Lattice 2005 (Topology and Confinement), Dublin, July 25-30, 2005; 6 pages, 5 figures, 2 mpeg animations not included because of file size problems but href linked, uses PoS.cls; to appear in Proceedings of Science

Journal ref: PoS LAT2005 (2005) 308

Observation of String Breaking in QCD

Authors: SESAM Collaboration, Gunnar S. Bali, Hartmut Neff, Thomas Duessel, Thomas Lippert, Klaus Schilling

Abstract: We numerically investigate the transition of the static quark-antiquark string into a static-light meson-antimeson system. Improving noise reduction techniques, we are able to resolve the signature of string breaking dynamics for n_f=2 lattice QCD at zero temperature. This result can be related to properties of quarkonium systems. We also study short-distance interactions between two static-ligh… ▽ More We numerically investigate the transition of the static quark-antiquark string into a static-light meson-antimeson system. Improving noise reduction techniques, we are able to resolve the signature of string breaking dynamics for n_f=2 lattice QCD at zero temperature. This result can be related to properties of quarkonium systems. We also study short-distance interactions between two static-light mesons. △ Less

Submitted 1 July, 2005; v1 submitted 13 May, 2005; originally announced May 2005.

Comments: 27 pages, 22 figures, changed decimal place of errors in 3 entries of Table, corrected references

Journal ref: Phys.Rev.D71:114513,2005

Finite-Size Effects in Lattice QCD with Dynamical Wilson Fermions

Authors: B. Orth, Th. Lippert, K. Schilling

Abstract: As computing resources are limited, choosing the parameters for a full Lattice QCD simulation always amounts to a compromise between the competing objectives of a lattice spacing as small, quarks as light, and a volume as large as possible. Aiming to push unquenched simulations with the Wilson action towards the computationally expensive regime of small quark masses we address the question wheth… ▽ More As computing resources are limited, choosing the parameters for a full Lattice QCD simulation always amounts to a compromise between the competing objectives of a lattice spacing as small, quarks as light, and a volume as large as possible. Aiming to push unquenched simulations with the Wilson action towards the computationally expensive regime of small quark masses we address the question whether one can possibly save computing time by extrapolating results from small lattices to the infinite volume, prior to the usual chiral and continuum extrapolations. In the present work the systematic volume dependence of simulated pion and nucleon masses is investigated and compared with a long-standing analytic formula by Luescher and with results from Chiral Perturbation Theory. We analyze data from Hybrid Monte Carlo simulations with the standard (unimproved) two-flavor Wilson action at two different lattice spacings of a=0.08fm and 0.13fm. The quark masses considered correspond to approximately 85 and 50% (at the smaller a) and 36% (at the larger a) of the strange quark mass. At each quark mass we study at least three different lattices with L/a=10 to 24 sites in the spatial directions (L=0.85-2.08fm). △ Less

Submitted 11 April, 2005; v1 submitted 15 March, 2005; originally announced March 2005.

Comments: 21 pages, 20 figures, REVTeX 4; v2: caption of Fig.7 corrected, one reference added

Report number: WUB 05-03, FZJ-ZAM-IB-2005-01

Journal ref: Phys.Rev. D72 (2005) 014503

Numerical Methods for the QCD Overlap Operator IV: Hybrid Monte Carlo

Authors: N. Cundy, S. Krieg, G. Arnold, A. Frommer, Th. Lippert, K. Schilling

Abstract: The extreme computational costs of calculating the sign of the Wilson matrix within the overlap operator have so far prevented four dimensional dynamical overlap simulations on realistic lattice sizes, because the computational power required to invert the overlap operator, the time consuming part of the Hybrid Monte Carlo algorithm, is too high. In this series of papers we introduced the optima… ▽ More The extreme computational costs of calculating the sign of the Wilson matrix within the overlap operator have so far prevented four dimensional dynamical overlap simulations on realistic lattice sizes, because the computational power required to invert the overlap operator, the time consuming part of the Hybrid Monte Carlo algorithm, is too high. In this series of papers we introduced the optimal approximation of the sign function and have been developing preconditioning and relaxation techniques which reduce the time needed for the inversion of the overlap operator by over a factor of four, bringing the simulation of dynamical overlap fermions on medium-size lattices within the range of Teraflop-computers. In this paper we adapt the HMC algorithm to overlap fermions. We approximate the matrix sign function using the Zolotarev rational approximation, treating the smallest eigenvalues of the Wilson operator exactly within the fermionic force. We then derive the fermionic force for the overlap operator, elaborating on the problem of Dirac delta-function terms from zero crossings of eigenvalues of the Wilson operator. The crossing scheme proposed shows energy violations which are better than O($Δτ^2$) and thus are comparable with the violations of the standard leapfrog algorithm over the course of a trajectory. We explicitly prove that our algorithm satisfies reversibility and area conservation. Finally, we test our algorithm on small $4^4$, $6^4$, and $8^4$ lattices at large masses. △ Less

Submitted 18 December, 2008; v1 submitted 15 February, 2005; originally announced February 2005.

Comments: v2 60 pages; substantial changes to all parts of the article; v3 minor revsions

Journal ref: Comput.Phys.Commun.180:26-54,2009

Helicity Dependent and Independent Generalized Parton Distributions of the Nucleon in Lattice QCD

Authors: Ph. Hagler, J. W. Negele, D. B. Renner, W. Schroers, T. Lippert, K. Schilling

Abstract: A complete description of the nucleon structure in terms of generalized parton distributions (GPDs) at twist 2 level requires the measurement/computation of the eight functions H, E, \tilde H, \tilde E, H_T, E_T, \tilde H_T and \tilde E_T, all depending on the three variables x, ξand t. In this talk, we present and discuss our first steps in the framework of lattice QCD towards this enormous tas… ▽ More A complete description of the nucleon structure in terms of generalized parton distributions (GPDs) at twist 2 level requires the measurement/computation of the eight functions H, E, \tilde H, \tilde E, H_T, E_T, \tilde H_T and \tilde E_T, all depending on the three variables x, ξand t. In this talk, we present and discuss our first steps in the framework of lattice QCD towards this enormous task. Dynamical lattice QCD results for the lowest three Mellin moments of the helicity dependent and independent GPDs are shown in terms of their corresponding generalized form factors. Implications for the transverse coordinate space structure of the nucleon as well as the orbital angular momentum (OAM) contribution of quarks to the nucleon spin are discussed in some detail. △ Less

Submitted 14 October, 2004; v1 submitted 1 October, 2004; originally announced October 2004.

Comments: 5 pages, 5 figures, Talk presented by Ph.H. at Electron-Nucleus Scattering VIII, Elba, Italy, June 21-25, 2004; typos corrected, minor change in wording on p.4&5

Journal ref: Eur.Phys.J.A24s1:29-33,2005

String breaking with dynamical Wilson fermions

Authors: G. S. Bali, Th. Duessel, Th. Lippert, H. Neff, Z. Prkacin, K. Schilling

Abstract: We present results of our ongoing determination of string breaking in full QCD with N_f=2 Wilson fermions. Our investigation of the fission of the static quark-antiquark string into a static-light meson-antimeson system is based on dynamical configurations of size 24^3 x 40 produced by the TxL collaboration. Combining various optimization methods we determine the matrix elements of the two-by-tw… ▽ More We present results of our ongoing determination of string breaking in full QCD with N_f=2 Wilson fermions. Our investigation of the fission of the static quark-antiquark string into a static-light meson-antimeson system is based on dynamical configurations of size 24^3 x 40 produced by the TxL collaboration. Combining various optimization methods we determine the matrix elements of the two-by-two system with so far unprecedented accuracy. The all-to-all light quark propagators occurring in the transition element are computed from eigenmodes of the Hermitian Wilson-Dirac matrix complemented by stochastic estimates in the orthogonal subspace. We observe a clear signature for level-splitting between ground state and excited potential. Thus, for the first time, string breaking induced by sea quarks is observed in a simulation of 4-dimensional lattice-QCD. △ Less

Submitted 22 September, 2004; originally announced September 2004.

Comments: 3 pages, 3 figures, contribution to Lattice 2004

Journal ref: Nucl.Phys.Proc.Suppl.140:609-611,2005

Dynamical overlap simulations using HMC

Authors: N. Cundy, S. Krieg, A. Frommer, Th. Lippert, K. Schilling

Abstract: We apply the Hybrid Monte Carlo method to the simulation of overlap fermions. We give the fermionic force for the molecular dynamics update. We present early results on a small dynamical chiral ensemble. We apply the Hybrid Monte Carlo method to the simulation of overlap fermions. We give the fermionic force for the molecular dynamics update. We present early results on a small dynamical chiral ensemble. △ Less

Submitted 28 September, 2004; v1 submitted 8 September, 2004; originally announced September 2004.

Comments: Lattice2004(machines), 3 pages; references updated, minor changes to text

Insight into nucleon structure from lattice calculations of moments of parton and generalized parton distributions

Authors: J. W. Negele, R. C. Brower, P. Dreher, R. Edwards, G. Fleming, Ph. Hagler, U. M. Heller, Th. Lippert, A. V. Pochinsky, D. B. Renner, D. Richards, K. Schilling, W. Schroers

Abstract: This talk presents recent calculations in full QCD of the lowest three moments of generalized parton distributions and the insight they provide into the behavior of nucleon electromagnetic form factors, the origin of the nucleon spin, and the transverse structure of the nucleon. In addition, new exploratory calculations in the chiral regime of full QCD are discussed. This talk presents recent calculations in full QCD of the lowest three moments of generalized parton distributions and the insight they provide into the behavior of nucleon electromagnetic form factors, the origin of the nucleon spin, and the transverse structure of the nucleon. In addition, new exploratory calculations in the chiral regime of full QCD are discussed. △ Less

Submitted 5 April, 2004; originally announced April 2004.

Comments: 9 pages, 9 figures, Cairns Workshop on Lattice Hadron Physics

Journal ref: Nucl.Phys.Proc.Suppl.128:170-178,2004

Computing the eta and eta' Mesons in Lattice QCD

Authors: K. Schilling, H. Neff, Th. Lippert

Abstract: It has been known for a long time that the large experimental singlet-octet mass gap in the pseudoscalar meson mass spectrum originates from the anomaly of the axial vector current, i.e. from nonperturbative effects and the nontrivial topological structure of the QCD vacuum. In the N_colour -> infinity limit of the theory, this connection elucidates in the famous Witten-Veneziano relation betwee… ▽ More It has been known for a long time that the large experimental singlet-octet mass gap in the pseudoscalar meson mass spectrum originates from the anomaly of the axial vector current, i.e. from nonperturbative effects and the nontrivial topological structure of the QCD vacuum. In the N_colour -> infinity limit of the theory, this connection elucidates in the famous Witten-Veneziano relation between the eta'-mass and the topological susceptibility of the quenched QCD vacuum.While lattice QCD has by now produced impressive high precision results on the flavour nonsinglet hadron spectrum, the determination of the pseudoscalar singlet mesons from direct correlator studies is markedly lagging behind, due to the computational complexity in handling observables that include OZI-rule violating diagrams, like the eta' propagator. In this article, we report on some recent progress in dealing with the numerical bottleneck problem. △ Less

Submitted 2 January, 2004; originally announced January 2004.

Comments: 29 pages, 20 figures, Introductory review based on a talk presented by K. Schilling at the 1st Workshop on Lattice Hadron Physics, Cairns, Australia

A partially quenched analysis of the eta-eta' system in N_f=2 QCD

Authors: H. Neff, Th. Lippert, J. Negele, K. Schilling

Abstract: We report on a first, comprehensive partially quenched study of the eta-eta' problem, based on SESAM configurations on a 16^3x32 lattice at beta=5.6 QCD with two (mass degenerate) active sea quark flavours. By means of the spectral approximation of the two-loop (hairpin) diagrams, we find clear plateau formation in the effective masses which enables us both to determine the eta-eta' mass matrix… ▽ More We report on a first, comprehensive partially quenched study of the eta-eta' problem, based on SESAM configurations on a 16^3x32 lattice at beta=5.6 QCD with two (mass degenerate) active sea quark flavours. By means of the spectral approximation of the two-loop (hairpin) diagrams, we find clear plateau formation in the effective masses which enables us both to determine the eta-eta' mass matrix and the alpha-parameter in the effective chiral Lagrangian for the flavour singlet sector, alpha=0.028 +- 0.013. △ Less

Submitted 2 January, 2004; originally announced January 2004.

Comments: 3 pages, 5 figures, contribution to Lattice2003

Journal ref: Nucl.Phys.Proc.Suppl.129:218-220,2004

Transverse Structure of Nucleon Parton Distributions from Lattice QCD

Authors: LHPC, SESAM Collaborations, :, Ph. Hagler, J. W. Negele, D. B. Renner, W. Schroers, Th. Lippert, K. Schilling

Abstract: This work presents the first calculation in lattice QCD of three moments of spin-averaged and spin-polarized generalized parton distributions in the proton. It is shown that the slope of the associated generalized form factors decreases significantly as the moment increases, indicating that the transverse size of the light-cone quark distribution decreases as the momentum fraction of the struck… ▽ More This work presents the first calculation in lattice QCD of three moments of spin-averaged and spin-polarized generalized parton distributions in the proton. It is shown that the slope of the associated generalized form factors decreases significantly as the moment increases, indicating that the transverse size of the light-cone quark distribution decreases as the momentum fraction of the struck parton increases. △ Less

Submitted 8 December, 2003; originally announced December 2003.

Comments: 4 pages, 1 figure

Report number: CTP-3452

Journal ref: Phys.Rev.Lett.93:112001,2004

Volume dependence of light hadron masses in full lattice QCD

Authors: B. Orth, T. Lippert, K. Schilling

Abstract: The aim of the GRAL project is to simulate full QCD with standard Wilson fermions at light quark masses on small to medium-sized lattices and to obtain infinite-volume results by extrapolation. In order to establish the functional form of the volume dependence we study systematically the finite-size effects in the light hadron spectrum. We give an update on the status of the GRAL project and sho… ▽ More The aim of the GRAL project is to simulate full QCD with standard Wilson fermions at light quark masses on small to medium-sized lattices and to obtain infinite-volume results by extrapolation. In order to establish the functional form of the volume dependence we study systematically the finite-size effects in the light hadron spectrum. We give an update on the status of the GRAL project and show that our simulation data for the light hadron masses depend exponentially on the lattice size. △ Less

Submitted 15 September, 2003; originally announced September 2003.

Comments: 3 pages, 1 figure, Lattice2003(spectrum)

Moments of nucleon spin-dependent generalized parton distributions

Authors: W. Schroers, R. C. Brower, P. Dreher, R. Edwards, G. Fleming, Ph. Hagler, U. M. Heller, Th. Lippert, J. W. Negele, A. V. Pochinsky, D. B. Renner, D. Richards, K. Schilling

Abstract: We present a lattice measurement of the first two moments of the spin-dependent GPD H-tilde(x,xi,t). From these we obtain the axial coupling constant and the second moment of the spin-dependent forward parton distribution. The measurements are done in full QCD using Wilson fermions. In addition, we also present results from a first exploratory study of full QCD using Asqtad sea and domain-wall v… ▽ More We present a lattice measurement of the first two moments of the spin-dependent GPD H-tilde(x,xi,t). From these we obtain the axial coupling constant and the second moment of the spin-dependent forward parton distribution. The measurements are done in full QCD using Wilson fermions. In addition, we also present results from a first exploratory study of full QCD using Asqtad sea and domain-wall valence fermions. △ Less

Submitted 14 September, 2003; originally announced September 2003.

Comments: Lattice2003(Theory), 3 pages, 3 figures, to appear in the Proceedings of Lattice 2003

Journal ref: Nucl.Phys.Proc.Suppl.129:907-909,2004

Insight into nucleon structure from generalized parton distributions

Authors: LHPC Collaboration, SESAM Collaboration, J. W. Negele, R. C. Brower, P. Dreher, R. Edwards, G. Fleming, Ph. Hagler, Th. Lippert, A. V. Pochinsky, D. B. Renner, D. Richards, K. Schilling, W. Schroers

Abstract: The lowest three moments of generalized parton distributions are calculated in full QCD and provide new insight into the behavior of nucleon electromagnetic form factors, the origin of the nucleon spin, and the transverse structure of the nucleon. The lowest three moments of generalized parton distributions are calculated in full QCD and provide new insight into the behavior of nucleon electromagnetic form factors, the origin of the nucleon spin, and the transverse structure of the nucleon. △ Less

Submitted 15 September, 2003; originally announced September 2003.

Comments: 3 pages, Lattice2003(Theoretical developments)

Journal ref: Nucl.Phys.Proc.Suppl.129:910-912,2004

$γN \to Δ$ transition form factors in Quenched and $N_F=2$ QCD

Authors: C. Alexandrou, Ph. de Forcrand, Th. Lippert, H. Neff, J. W. Negele, K. Schilling, W. Schroers, A. Tsapalis

Abstract: Calculations of the magnetic dipole, electric quadrupole and Coulomb quadrupole amplitudes for the transition $γN\to Δ$ are presented both in quenched QCD and with two flavours of degenerate dynamical quarks. Calculations of the magnetic dipole, electric quadrupole and Coulomb quadrupole amplitudes for the transition $γN\to Δ$ are presented both in quenched QCD and with two flavours of degenerate dynamical quarks. △ Less

Submitted 10 September, 2003; originally announced September 2003.

Comments: Lattice2003(Matrix), 3 pages

Journal ref: Nucl.Phys.Proc.Suppl.129:302-304,2004

N to Delta electromagnetic transition form factors from Lattice QCD

Authors: C. Alexandrou, Ph. de Forcrand, Th. Lippert, H. Neff, J. W. Negele, K. Schilling, W. Schroers, A. Tsapalis

Abstract: The magnetic dipole, the electric quadrupole and the Coulomb quadrupole amplitudes for the transition $γN\to Δ$ are evaluated both in quenched lattice QCD at $β=6.0$ and using two dynamical Wilson fermions simulated at $β=5.6$. The dipole transition form factor is accurately determined at several values of momentum transfer. On the lattices studied in this work, the electric quadrupole amplitude… ▽ More The magnetic dipole, the electric quadrupole and the Coulomb quadrupole amplitudes for the transition $γN\to Δ$ are evaluated both in quenched lattice QCD at $β=6.0$ and using two dynamical Wilson fermions simulated at $β=5.6$. The dipole transition form factor is accurately determined at several values of momentum transfer. On the lattices studied in this work, the electric quadrupole amplitude is found to be non-zero yielding a negative value for the ratio, $ R_{EM}$, of electric quadrupole to magnetic dipole amplitudes at three values of momentum transfer. △ Less

Submitted 3 April, 2004; v1 submitted 10 July, 2003; originally announced July 2003.

Comments: Discussion of systematic errors and chiral extrapolation extended; Typos corrected; Version to appear in Phys. Rev. D.; 19 pages including 21 figures

Journal ref: Phys.Rev.D69:114506,2004

On the scaling of computational particle physics codes on cluster computers

Authors: Z. Sroczynski, N. Eicker, Th. Lippert, B. Orth, K. Schilling

Abstract: Many appplications in computational science are sufficiently compute-intensive that they depend on the power of parallel computing for viability. For all but the "embarrassingly parallel" problems, the performance depends upon the level of granularity that can be achieved on the computer platform. Our computational particle physics applications require machines that can support a wide range of… ▽ More Many appplications in computational science are sufficiently compute-intensive that they depend on the power of parallel computing for viability. For all but the "embarrassingly parallel" problems, the performance depends upon the level of granularity that can be achieved on the computer platform. Our computational particle physics applications require machines that can support a wide range of granularities, but in general, compute-intensive state-of-the-art projects will require finely grained distributions. Of the different types of machines available for the task, we consider cluster computers. The use of clusters of commodity computers in high performance computing has many advantages including the raw price/performance ratio and the flexibility of machine configuration and upgrade. Here we focus on what is usually considered the weak point of cluster technology; the scaling behaviour when faced with a numerically intensive parallel computation. To this end we examine the scaling of our own applications from numerical quantum field theory on a cluster and infer conclusions about the more general case. △ Less

Submitted 10 July, 2003; v1 submitted 9 July, 2003; originally announced July 2003.

Comments: 26pp. LaTeX2e using package graphicx. 16 PostScript figures

Report number: LTH 583

Moments of Nucleon Generalized Parton Distributions in Lattice QCD

Authors: Philipp Hagler, John Negele, Dru Renner, Wolfram Schroers, Thomas Lippert, Klaus Schilling

Abstract: Calculation of moments of generalized parton distributions in lattice QCD requires more powerful techniques than those previously used to calculate moments of structure functions. Hence, we present a novel approach that exploits the full information content from a given lattice configuration by measuring an overdetermined set of lattice observables to provide maximal statistical constraints on t… ▽ More Calculation of moments of generalized parton distributions in lattice QCD requires more powerful techniques than those previously used to calculate moments of structure functions. Hence, we present a novel approach that exploits the full information content from a given lattice configuration by measuring an overdetermined set of lattice observables to provide maximal statistical constraints on the generalized form factors at a given virtuality, t. In an exploratory investigation using unquenched QCD configurations at intermediate sea quark masses, we demonstrate that our new technique is superior to conventional methods and leads to reliable numerical signals for the n=2 flavor singlet generalized form factors up to 3 GeV^2. The contribution from connected diagrams in the flavor singlet sector to the total quark angular momentum is measured to an accuracy of the order of one percent. △ Less

Submitted 4 July, 2003; v1 submitted 26 April, 2003; originally announced April 2003.

Comments: 16 pages, 8 figures, LaTeX, minor elaboration of formalism and singular value decomposition for non-specialists and addition of several references

Report number: MIT-3353

Journal ref: Phys.Rev.D68:034505,2003

Fast Parallel I/O on Cluster Computers

Authors: Thomas Duessel, Norbert Eicker, Florin Isaila, Thomas Lippert, Thomas Moschny, Hartmut Neff, Klaus Schilling, Walter Tichy

Abstract: Today's cluster computers suffer from slow I/O, which slows down I/O-intensive applications. We show that fast disk I/O can be achieved by operating a parallel file system over fast networks such as Myrinet or Gigabit Ethernet. In this paper, we demonstrate how the ParaStation3 communication system helps speed-up the performance of parallel I/O on clusters using the open source parallel virtua… ▽ More Today's cluster computers suffer from slow I/O, which slows down I/O-intensive applications. We show that fast disk I/O can be achieved by operating a parallel file system over fast networks such as Myrinet or Gigabit Ethernet. In this paper, we demonstrate how the ParaStation3 communication system helps speed-up the performance of parallel I/O on clusters using the open source parallel virtual file system (PVFS) as testbed and production system. We will describe the set-up of PVFS on the Alpha-Linux-Cluster-Engine (ALiCE) located at Wuppertal University, Germany. Benchmarks on ALiCE achieve write-performances of up to 1 GB/s from a 32-processor compute-partition to a 32-processor PVFS I/O-partition, outperforming known benchmark results for PVFS on the same network by more than a factor of 2. Read-performance from buffer-cache reaches up to 2.2 GB/s. Our benchmarks are giant, I/O-intensive eigenmode problems from lattice quantum chromodynamics, demonstrating stability and performance of PVFS over Parastation in large-scale production runs. △ Less

Submitted 19 March, 2003; originally announced March 2003.

Comments: 22 pages, 10 figures

ACM Class: B.4.3; C.1.2; C.2.2; D.4.3

Continuum Extrapolation of Moments of Nucleon Quark Distributions in Full QCD

Authors: LHPC, SESAM Collaborations, :, P. Dreher, R. Brower, S. Capitani, D. Dolgov, R. Edwards, N. Eicker, U. M. Heller, Th. Lipert, J. W. Negele, A. Pochinsky, D. B. Renner, K. Schilling

Abstract: Moments of light cone quark density, helicity, and transversity distributions are calculated in unquenched lattice QCD at $β= 5.5$ and $β= 5.3$ using Wilson fermions on $ 16^3 \times 32 $ lattices. These results are combined with earlier calculations at $β= 5.6$ using SESAM configurations to study the continuum limit. Moments of light cone quark density, helicity, and transversity distributions are calculated in unquenched lattice QCD at $β= 5.5$ and $β= 5.3$ using Wilson fermions on $ 16^3 \times 32 $ lattices. These results are combined with earlier calculations at $β= 5.6$ using SESAM configurations to study the continuum limit. △ Less

Submitted 11 November, 2002; originally announced November 2002.

Comments: Lattice 2002 (hadronic matrix elements) 3 pages, 1 figure

Journal ref: Nucl.Phys.Proc.Suppl.119:392-394,2003

Calculation of Nucleon Electromagnetic Form Factors

Authors: LHPC, SESAM Collaborations, :, D. B. Renner, R. Brower, D. Dolgov, N. Eicker, Th. Lippert, J. W. Negele, A. Pochinsky, K. Schilling

Abstract: The fomalism is developed to express nucleon matrix elements of the electromagnetic current in terms of form factors consistent with the translational, rotational, and parity symmetries of a cubic lattice. We calculate the number of these form factors and show how appropriate linear combinations approach the continuum limit. The fomalism is developed to express nucleon matrix elements of the electromagnetic current in terms of form factors consistent with the translational, rotational, and parity symmetries of a cubic lattice. We calculate the number of these form factors and show how appropriate linear combinations approach the continuum limit. △ Less

Submitted 26 November, 2002; v1 submitted 11 November, 2002; originally announced November 2002.

Comments: Lattice 2002 (hadronic matrix elements) 3 pages

Journal ref: Nucl.Phys.Proc.Suppl.119:395-397,2003

Compact QED under scrutiny: it's first order

Authors: G. Arnold, B. Bunk, Th. Lippert, K. Schilling

Abstract: We report new results from our finite size scaling analysis of 4d compact pure U(1) gauge theory with Wilson action. Investigating several cumulants of the plaquette energy within the Borgs-Kotecky finite size scaling scheme we find strong evidence for a first-order phase transition and present a high precision value for the critical coupling in the thermodynamic limit. We report new results from our finite size scaling analysis of 4d compact pure U(1) gauge theory with Wilson action. Investigating several cumulants of the plaquette energy within the Borgs-Kotecky finite size scaling scheme we find strong evidence for a first-order phase transition and present a high precision value for the critical coupling in the thermodynamic limit. △ Less

Submitted 7 October, 2002; originally announced October 2002.

Comments: Lattice2002(Spin)