-
Physics-informed motion registration of lung parenchyma across static CT images
Authors:
Sunder Neelakantan,
Tanmay Mukherjee,
Kyle J. Myers,
Rahim Rizi,
Reza Avazmohammadi
Abstract:
Lung injuries, such as ventilator-induced lung injury and radiation-induced lung injury, can lead to heterogeneous alterations in the biomechanical behavior of the lungs. While imaging methods, e.g., X-ray and static computed tomography (CT), can point to regional alterations in lung structure between healthy and diseased tissue, they fall short of delineating timewise kinematic variations between…
▽ More
Lung injuries, such as ventilator-induced lung injury and radiation-induced lung injury, can lead to heterogeneous alterations in the biomechanical behavior of the lungs. While imaging methods, e.g., X-ray and static computed tomography (CT), can point to regional alterations in lung structure between healthy and diseased tissue, they fall short of delineating timewise kinematic variations between the former and the latter. Image registration has gained recent interest as a tool to estimate the displacement experienced by the lungs during respiration via regional deformation metrics such as volumetric expansion and distortion. However, successful image registration commonly relies on a temporal series of image stacks with small displacements in the lungs across succeeding image stacks, which remains limited in static imaging. In this study, we have presented a finite element (FE) method to estimate strains from static images acquired at the end-expiration (EE) and end-inspiration (EI) timepoints, i.e., images with a large deformation between the two distant timepoints. Physiologically realistic loads were applied to the geometry obtained at EE to deform this geometry to match the geometry obtained at EI. The results indicated that the simulation could minimize the error between the two geometries. Using four-dimensional (4D) dynamic CT in a rat, the strain at an isolated transverse plane estimated by our method showed sufficient agreement with that estimated through non-rigid image registration that used all the timepoints. Through the proposed method, we can estimate the lung deformation at any timepoint between EE and EI. The proposed method offers a tool to estimate timewise regional deformation in the lungs using only static images acquired at EE and EI.
△ Less
Submitted 3 July, 2024;
originally announced July 2024.
-
Acceptance Tests of more than 10 000 Photomultiplier Tubes for the multi-PMT Digital Optical Modules of the IceCube Upgrade
Authors:
R. Abbasi,
M. Ackermann,
J. Adams,
S. K. Agarwalla,
J. A. Aguilar,
M. Ahlers,
J. M. Alameddine,
N. M. Amin,
K. Andeen,
C. Argüelles,
Y. Ashida,
S. Athanasiadou,
L. Ausborm,
S. N. Axani,
X. Bai,
A. Balagopal V.,
M. Baricevic,
S. W. Barwick,
S. Bash,
V. Basu,
R. Bay,
J. J. Beatty,
J. Becker Tjus,
J. Beise,
C. Bellenghi
, et al. (399 additional authors not shown)
Abstract:
More than 10,000 photomultiplier tubes (PMTs) with a diameter of 80 mm will be installed in multi-PMT Digital Optical Modules (mDOMs) of the IceCube Upgrade. These have been tested and pre-calibrated at two sites. A throughput of more than 1000 PMTs per week with both sites was achieved with a modular design of the testing facilities and highly automated testing procedures. The testing facilities…
▽ More
More than 10,000 photomultiplier tubes (PMTs) with a diameter of 80 mm will be installed in multi-PMT Digital Optical Modules (mDOMs) of the IceCube Upgrade. These have been tested and pre-calibrated at two sites. A throughput of more than 1000 PMTs per week with both sites was achieved with a modular design of the testing facilities and highly automated testing procedures. The testing facilities can easily be adapted to other PMTs, such that they can, e.g., be re-used for testing the PMTs for IceCube-Gen2. Single photoelectron response, high voltage dependence, time resolution, prepulse, late pulse, afterpulse probabilities, and dark rates were measured for each PMT. We describe the design of the testing facilities, the testing procedures, and the results of the acceptance tests.
△ Less
Submitted 20 June, 2024; v1 submitted 30 April, 2024;
originally announced April 2024.
-
Improved modeling of in-ice particle showers for IceCube event reconstruction
Authors:
R. Abbasi,
M. Ackermann,
J. Adams,
S. K. Agarwalla,
J. A. Aguilar,
M. Ahlers,
J. M. Alameddine,
N. M. Amin,
K. Andeen,
G. Anton,
C. Argüelles,
Y. Ashida,
S. Athanasiadou,
L. Ausborm,
S. N. Axani,
X. Bai,
A. Balagopal V.,
M. Baricevic,
S. W. Barwick,
S. Bash,
V. Basu,
R. Bay,
J. J. Beatty,
J. Becker Tjus,
J. Beise
, et al. (394 additional authors not shown)
Abstract:
The IceCube Neutrino Observatory relies on an array of photomultiplier tubes to detect Cherenkov light produced by charged particles in the South Pole ice. IceCube data analyses depend on an in-depth characterization of the glacial ice, and on novel approaches in event reconstruction that utilize fast approximations of photoelectron yields. Here, a more accurate model is derived for event reconstr…
▽ More
The IceCube Neutrino Observatory relies on an array of photomultiplier tubes to detect Cherenkov light produced by charged particles in the South Pole ice. IceCube data analyses depend on an in-depth characterization of the glacial ice, and on novel approaches in event reconstruction that utilize fast approximations of photoelectron yields. Here, a more accurate model is derived for event reconstruction that better captures our current knowledge of ice optical properties. When evaluated on a Monte Carlo simulation set, the median angular resolution for in-ice particle showers improves by over a factor of three compared to a reconstruction based on a simplified model of the ice. The most substantial improvement is obtained when including effects of birefringence due to the polycrystalline structure of the ice. When evaluated on data classified as particle showers in the high-energy starting events sample, a significantly improved description of the events is observed.
△ Less
Submitted 22 April, 2024; v1 submitted 4 March, 2024;
originally announced March 2024.
-
Measurement of Atmospheric Neutrino Mixing with Improved IceCube DeepCore Calibration and Data Processing
Authors:
IceCube Collaboration,
R. Abbasi,
M. Ackermann,
J. Adams,
S. K. Agarwalla,
J. A. Aguilar,
M. Ahlers,
J. M. Alameddine,
N. M. Amin,
K. Andeen,
G. Anton,
C. Argüelles,
Y. Ashida,
S. Athanasiadou,
S. N. Axani,
X. Bai,
A. Balagopal V.,
M. Baricevic,
S. W. Barwick,
V. Basu,
R. Bay,
J. J. Beatty,
K. -H. Becker,
J. Becker Tjus,
J. Beise
, et al. (383 additional authors not shown)
Abstract:
We describe a new data sample of IceCube DeepCore and report on the latest measurement of atmospheric neutrino oscillations obtained with data recorded between 2011-2019. The sample includes significant improvements in data calibration, detector simulation, and data processing, and the analysis benefits from a detailed treatment of systematic uncertainties, with significantly higher level of detai…
▽ More
We describe a new data sample of IceCube DeepCore and report on the latest measurement of atmospheric neutrino oscillations obtained with data recorded between 2011-2019. The sample includes significant improvements in data calibration, detector simulation, and data processing, and the analysis benefits from a detailed treatment of systematic uncertainties, with significantly higher level of detail since our last study. By measuring the relative fluxes of neutrino flavors as a function of their reconstructed energies and arrival directions we constrain the atmospheric neutrino mixing parameters to be $\sin^2θ_{23} = 0.51\pm 0.05$ and $Δm^2_{32} = 2.41\pm0.07\times 10^{-3}\mathrm{eV}^2$, assuming a normal mass ordering. The resulting 40\% reduction in the error of both parameters with respect to our previous result makes this the most precise measurement of oscillation parameters using atmospheric neutrinos. Our results are also compatible and complementary to those obtained using neutrino beams from accelerators, which are obtained at lower neutrino energies and are subject to different sources of uncertainties.
△ Less
Submitted 8 August, 2023; v1 submitted 24 April, 2023;
originally announced April 2023.
-
Graph Neural Networks for Low-Energy Event Classification & Reconstruction in IceCube
Authors:
R. Abbasi,
M. Ackermann,
J. Adams,
N. Aggarwal,
J. A. Aguilar,
M. Ahlers,
M. Ahrens,
J. M. Alameddine,
A. A. Alves Jr.,
N. M. Amin,
K. Andeen,
T. Anderson,
G. Anton,
C. Argüelles,
Y. Ashida,
S. Athanasiadou,
S. Axani,
X. Bai,
A. Balagopal V.,
M. Baricevic,
S. W. Barwick,
V. Basu,
R. Bay,
J. J. Beatty,
K. -H. Becker
, et al. (359 additional authors not shown)
Abstract:
IceCube, a cubic-kilometer array of optical sensors built to detect atmospheric and astrophysical neutrinos between 1 GeV and 1 PeV, is deployed 1.45 km to 2.45 km below the surface of the ice sheet at the South Pole. The classification and reconstruction of events from the in-ice detectors play a central role in the analysis of data from IceCube. Reconstructing and classifying events is a challen…
▽ More
IceCube, a cubic-kilometer array of optical sensors built to detect atmospheric and astrophysical neutrinos between 1 GeV and 1 PeV, is deployed 1.45 km to 2.45 km below the surface of the ice sheet at the South Pole. The classification and reconstruction of events from the in-ice detectors play a central role in the analysis of data from IceCube. Reconstructing and classifying events is a challenge due to the irregular detector geometry, inhomogeneous scattering and absorption of light in the ice and, below 100 GeV, the relatively low number of signal photons produced per event. To address this challenge, it is possible to represent IceCube events as point cloud graphs and use a Graph Neural Network (GNN) as the classification and reconstruction method. The GNN is capable of distinguishing neutrino events from cosmic-ray backgrounds, classifying different neutrino event types, and reconstructing the deposited energy, direction and interaction vertex. Based on simulation, we provide a comparison in the 1-100 GeV energy range to the current state-of-the-art maximum likelihood techniques used in current IceCube analyses, including the effects of known systematic uncertainties. For neutrino event classification, the GNN increases the signal efficiency by 18% at a fixed false positive rate (FPR), compared to current IceCube methods. Alternatively, the GNN offers a reduction of the FPR by over a factor 8 (to below half a percent) at a fixed signal efficiency. For the reconstruction of energy, direction, and interaction vertex, the resolution improves by an average of 13%-20% compared to current maximum likelihood techniques in the energy range of 1-30 GeV. The GNN, when run on a GPU, is capable of processing IceCube events at a rate nearly double of the median IceCube trigger rate of 2.7 kHz, which opens the possibility of using low energy neutrinos in online searches for transient events.
△ Less
Submitted 11 October, 2022; v1 submitted 7 September, 2022;
originally announced September 2022.
-
Low Energy Event Reconstruction in IceCube DeepCore
Authors:
R. Abbasi,
M. Ackermann,
J. Adams,
J. A. Aguilar,
M. Ahlers,
M. Ahrens,
J. M. Alameddine,
A. A. Alves Jr.,
N. M. Amin,
K. Andeen,
T. Anderson,
G. Anton,
C. Argüelles,
Y. Ashida,
S. Axani,
X. Bai,
A. Balagopal V.,
S. W. Barwick,
B. Bastian,
V. Basu,
S. Baur,
R. Bay,
J. J. Beatty,
K. -H. Becker,
J. Becker Tjus
, et al. (360 additional authors not shown)
Abstract:
The reconstruction of event-level information, such as the direction or energy of a neutrino interacting in IceCube DeepCore, is a crucial ingredient to many physics analyses. Algorithms to extract this high level information from the detector's raw data have been successfully developed and used for high energy events. In this work, we address unique challenges associated with the reconstruction o…
▽ More
The reconstruction of event-level information, such as the direction or energy of a neutrino interacting in IceCube DeepCore, is a crucial ingredient to many physics analyses. Algorithms to extract this high level information from the detector's raw data have been successfully developed and used for high energy events. In this work, we address unique challenges associated with the reconstruction of lower energy events in the range of a few to hundreds of GeV and present two separate, state-of-the-art algorithms. One algorithm focuses on the fast directional reconstruction of events based on unscattered light. The second algorithm is a likelihood-based multipurpose reconstruction offering superior resolutions, at the expense of larger computational cost.
△ Less
Submitted 4 March, 2022;
originally announced March 2022.
-
Stabilization method with Relativistic Configuration-interaction applied to two-electron resonances
Authors:
P. Amaro,
J. P. Santos,
S. Bhattacharyya,
T. K. Mukherjee,
J. K. Saha
Abstract:
We applied a relativistic configuration-interaction (CI) framework to the stabilization method as an approach for obtaining the autoionization resonance structure of heliumlike ions. In this method, the ion is confined within an impenetrable spherical cavity, the size of which determines the radial space available for electron wavefunctions and electron-electron interactions. By varying the size o…
▽ More
We applied a relativistic configuration-interaction (CI) framework to the stabilization method as an approach for obtaining the autoionization resonance structure of heliumlike ions. In this method, the ion is confined within an impenetrable spherical cavity, the size of which determines the radial space available for electron wavefunctions and electron-electron interactions. By varying the size of the cavity, one can obtain the autoionization resonance position and width. The applicability of this method is tested on the resonances of He atom while comparing with benchmark data available in the literature. The present method is further applied to the determination of the resonance structure of heliumlike uranium ion, where a relativistic framework is mandatory. In the strong-confinement region, the present method can be useful to simulate the properties of an atom or ion under extreme pressure. An exemplary application of the present method to determine the structure of ions embedded in a dense plasma environment is briefly discussed.
△ Less
Submitted 18 December, 2020;
originally announced December 2020.
-
Singly excited S-states of compressed two-electron ions
Authors:
J. K. Saha,
S. Bhattacharyya,
T. K. Mukherjee
Abstract:
A detailed analysis on the effect of spherical impenetrable confinement on the structural properties of two-electron ions in S-states have been done. The energy values of 1sns [n = 2-4] (3Se) states of helium-like ions (Z = 2-5) are estimated within the framework of Ritz variational method by using explicitly correlated Hylleraas-type basis sets. The correlated wave functions used here are consist…
▽ More
A detailed analysis on the effect of spherical impenetrable confinement on the structural properties of two-electron ions in S-states have been done. The energy values of 1sns [n = 2-4] (3Se) states of helium-like ions (Z = 2-5) are estimated within the framework of Ritz variational method by using explicitly correlated Hylleraas-type basis sets. The correlated wave functions used here are consistent with the finite boundary conditions due to spherical confinement. A comparative study between the singlet and triplet states originating from a particular electronic configuration shows incidental degeneracy and the subsequent level-crossing phenomena. The thermodynamic pressure felt by the ion inside the sphere pushes the energy levels towards continuum. Critical pressures for the transition to strong confinement regime (where the singly excited two-electron energy levels cross the corresponding one-electron threshold) as well as for the complete destabilization are also estimated.
△ Less
Submitted 25 November, 2015;
originally announced November 2015.
-
Electronic structure of helium atom in a quantum dot
Authors:
Jayanta K. Saha,
S. Bhattacharyya,
T. K. Mukherjee
Abstract:
Bound and resonance states of helium atom have been investigated inside a quantum dot by using explicitly correlated Hylleraas type basis set within the framework of stabilization method. To be specific, precise energy eigenvalues of bound 1sns (1Se) [n = 1-6] states and the resonance parameters i.e. positions and widths of 1Se states due to 2sns [n = 2-5] and 2pnp [n = 2-5] configuration of confi…
▽ More
Bound and resonance states of helium atom have been investigated inside a quantum dot by using explicitly correlated Hylleraas type basis set within the framework of stabilization method. To be specific, precise energy eigenvalues of bound 1sns (1Se) [n = 1-6] states and the resonance parameters i.e. positions and widths of 1Se states due to 2sns [n = 2-5] and 2pnp [n = 2-5] configuration of confined helium below N = 2 ionization threshold of He+ have been estimated. The two-parameter (Depth and Width) finite oscillator potential is used to represent the confining potential representing the quantum dot. It has been explicitly demonstrated that electronic structure properties become a sensitive function of the dot size. It is observed from the calculations of ionization potential that the stability of an impurity ion within quantum dot may be manipulated by varying the confinement parameters. A possibility of controlling the autoionization lifetime of doubly excited states of two-electron ions by tuning the width of the quantum cavity is also discussed here.
△ Less
Submitted 10 August, 2015;
originally announced August 2015.
-
Ab-initio calculations on two-electron ions in strongly coupled plasma environment
Authors:
S. Bhattacharyya,
J. K. Saha,
T. K. Mukherjee
Abstract:
In this work, the controversy between the interpretations of recent measurements on dense aluminum plasma created with Linac coherent light sources (LCLS) X-ray free electron laser (FEL) and Orion laser has been addressed. In both kind of experiments, helium-like and hydrogen-like spectral lines are used for plasma diagnostics . However, there exist no precise theoretical calculations for He-like…
▽ More
In this work, the controversy between the interpretations of recent measurements on dense aluminum plasma created with Linac coherent light sources (LCLS) X-ray free electron laser (FEL) and Orion laser has been addressed. In both kind of experiments, helium-like and hydrogen-like spectral lines are used for plasma diagnostics . However, there exist no precise theoretical calculations for He-like ions within dense plasma environment. The strong need for an accurate theoretical estimates for spectral properties of He-like ions in strongly coupled plasma environment leads us to perform ab initio calculations in the framework of Rayleigh-Ritz variation principle in Hylleraas coordinates where ion-sphere potential is used. An approach to resolve the long-drawn problem of numerical instability for evaluating two-electron integrals with extended basis inside a finite domain is presented here. The present values of electron densities corresponding to disappearance of different spectral lines obtained within the framework of ion-sphere potential show excellent agreement with Orion laser experiments in Al plasma and with recent theories. Moreover, this method is extended to predict the critical plasma densities at which the spectral lines of H-like and He-like carbon and argon ions disappear. Incidental degeneracy and level-crossing phenomena are being reported for the first time for two-electron ions embedded in strongly coupled plasma. Thermodynamic pressure experienced by the ions in their respective ground states inside the ion-spheres are also reported.
△ Less
Submitted 19 March, 2015;
originally announced March 2015.
-
Precise energy eigenvalues of hydrogen-like ion moving in quantum plasmas
Authors:
S. Dutta,
J. K. Saha,
T. K. Mukherjee
Abstract:
The analytic form of the electrostatic potential felt by a slowly moving test charge in quantum plasma is being derived. It has been shown that the potential composed of two parts: Debye-Huckel screening term and near-field wake potential which depends on the velocity of the test charge and the number density of the plasma electrons. Rayleigh-Ritz variational calculation has been done to estimate…
▽ More
The analytic form of the electrostatic potential felt by a slowly moving test charge in quantum plasma is being derived. It has been shown that the potential composed of two parts: Debye-Huckel screening term and near-field wake potential which depends on the velocity of the test charge and the number density of the plasma electrons. Rayleigh-Ritz variational calculation has been done to estimate precise energy eigenvalues of hydrogen-like ion under such plasma environment. A detailed analysis shows that the energy levels are gradually moves to the continuum with increasing plasma electron density while level crossing phenomenon have been observed with the variation of ion velocity.
△ Less
Submitted 1 December, 2014;
originally announced December 2014.
-
Precise Variational Calculation For The Doubly Excited State (2p^2)^3P^e of Helium
Authors:
Tapan K. Mukherjee,
Prasanta K. Mukherjee
Abstract:
Highly precise variational calculations of non-relativistic energies of the (2p^2)^3P^e state of Helium atom are presented.We get an upper bound energy E=-0.71050015565678 a.u.,the lowest yet obtained.
Highly precise variational calculations of non-relativistic energies of the (2p^2)^3P^e state of Helium atom are presented.We get an upper bound energy E=-0.71050015565678 a.u.,the lowest yet obtained.
△ Less
Submitted 15 June, 2004;
originally announced June 2004.