-
Fast Polynomial Approximation of Heat Kernel Convolution on Manifolds and Its Application to Brain Sulcal and Gyral Graph Pattern Analysis
Authors:
Shih-Gu Huang,
Ilwoo Lyu,
Anqi Qiu,
Moo K. Chung
Abstract:
Heat diffusion has been widely used in brain imaging for surface fairing, mesh regularization and cortical data smoothing. Motivated by diffusion wavelets and convolutional neural networks on graphs, we present a new fast and accurate numerical scheme to solve heat diffusion on surface meshes. This is achieved by approximating the heat kernel convolution using high degree orthogonal polynomials in…
▽ More
Heat diffusion has been widely used in brain imaging for surface fairing, mesh regularization and cortical data smoothing. Motivated by diffusion wavelets and convolutional neural networks on graphs, we present a new fast and accurate numerical scheme to solve heat diffusion on surface meshes. This is achieved by approximating the heat kernel convolution using high degree orthogonal polynomials in the spectral domain. We also derive the closed-form expression of the spectral decomposition of the Laplace-Beltrami operator and use it to solve heat diffusion on a manifold for the first time. The proposed fast polynomial approximation scheme avoids solving for the eigenfunctions of the Laplace-Beltrami operator, which is computationally costly for large mesh size, and the numerical instability associated with the finite element method based diffusion solvers. The proposed method is applied in localizing the male and female differences in cortical sulcal and gyral graph patterns obtained from MRI in an innovative way. The MATLAB code is available at http://www.stat.wisc.edu/~mchung/chebyshev.
△ Less
Submitted 17 January, 2020; v1 submitted 6 November, 2019;
originally announced November 2019.
-
Linear Beam Stability in Periodic Focusing Systems:Krein Signature and Band Structure
Authors:
Moses Chung,
Yoolim Cheon,
Hong Qin
Abstract:
The general question of how a beam becomes unstable has been one of the fundamental research topics among beam and accelerator physicists for several decades. In this study, we revisited the general problem of linear beam stability in periodic focusing systems by applying the concepts of Krein signature and band structure. We numerically calculated the eigenvalues and other associated characterist…
▽ More
The general question of how a beam becomes unstable has been one of the fundamental research topics among beam and accelerator physicists for several decades. In this study, we revisited the general problem of linear beam stability in periodic focusing systems by applying the concepts of Krein signature and band structure. We numerically calculated the eigenvalues and other associated characteristics of one-period maps, and discussed the stability properties of single-particle motions with skew quadrupoles and envelope perturbations in high-intensity beams on an equal footing. In particular, an application of the Krein theory to envelope instability analysis was newly attempted in this study. The appearance of instabilities is interpreted as the result of the collision between eigenmodes of opposite Krein signatures and the formation of a band gap.
△ Less
Submitted 19 September, 2019;
originally announced September 2019.
-
Deep learning of topological phase transitions from entanglement aspects
Authors:
Yuan-Hong Tsai,
Meng-Zhe Yu,
Yu-Hao Hsu,
Ming-Chiang Chung
Abstract:
The one-dimensional $p$-wave superconductor proposed by Kitaev has long been a classic example for understanding topological phase transitions through various methods, such as examining Berry phase, edge states of open chains and, in particular, aspects from quantum entanglement of ground states. In order to understand the amount of information carried in the entanglement-related quantities, here…
▽ More
The one-dimensional $p$-wave superconductor proposed by Kitaev has long been a classic example for understanding topological phase transitions through various methods, such as examining Berry phase, edge states of open chains and, in particular, aspects from quantum entanglement of ground states. In order to understand the amount of information carried in the entanglement-related quantities, here we study topological phase transitions of the model with emphasis of using the deep learning approach. We feed different quantities, including Majorana correlation matrices (MCMs), entanglement spectra (ES) or entanglement eigenvectors (EE) originated from Block correlation matrices (BCMs), into the deep neural networks for training, and investigate which one could be the most useful input format in this approach. We find that ES is indeed too compressed information compared to MCM or EE. MCM and EE can provide us abundant information to recognize not only the topological phase transitions in the model but also phases of matter with different $U$(1) gauges, which is not reachable by using ES only.
△ Less
Submitted 26 December, 2019; v1 submitted 10 September, 2019;
originally announced September 2019.
-
Resource Optimized Neural Architecture Search for 3D Medical Image Segmentation
Authors:
Woong Bae,
Seungho Lee,
Yeha Lee,
Beomhee Park,
Minki Chung,
Kyu-Hwan Jung
Abstract:
Neural Architecture Search (NAS), a framework which automates the task of designing neural networks, has recently been actively studied in the field of deep learning. However, there are only a few NAS methods suitable for 3D medical image segmentation. Medical 3D images are generally very large; thus it is difficult to apply previous NAS methods due to their GPU computational burden and long train…
▽ More
Neural Architecture Search (NAS), a framework which automates the task of designing neural networks, has recently been actively studied in the field of deep learning. However, there are only a few NAS methods suitable for 3D medical image segmentation. Medical 3D images are generally very large; thus it is difficult to apply previous NAS methods due to their GPU computational burden and long training time. We propose the resource-optimized neural architecture search method which can be applied to 3D medical segmentation tasks in a short training time (1.39 days for 1GB dataset) using a small amount of computation power (one RTX 2080Ti, 10.8GB GPU memory). Excellent performance can also be achieved without retraining(fine-tuning) which is essential in most NAS methods. These advantages can be achieved by using a reinforcement learning-based controller with parameter sharing and focusing on the optimal search space configuration of macro search rather than micro search. Our experiments demonstrate that the proposed NAS method outperforms manually designed networks with state-of-the-art performance in 3D medical image segmentation.
△ Less
Submitted 2 September, 2019;
originally announced September 2019.
-
Classical potential for general spinning bodies
Authors:
Ming-Zhi Chung,
Yu-tin Huang,
Jung-Wook Kim
Abstract:
In this paper we compute the spin-dependent terms of the gravitational potential for general spinning bodies at the leading Newton's constant $G$ and to all orders in spin. We utilize the on-shell approach, which extracts the classical potential directly from the scattering amplitude. For spinning particles, extra care is required due to the fact that the spin space of each particle is independent…
▽ More
In this paper we compute the spin-dependent terms of the gravitational potential for general spinning bodies at the leading Newton's constant $G$ and to all orders in spin. We utilize the on-shell approach, which extracts the classical potential directly from the scattering amplitude. For spinning particles, extra care is required due to the fact that the spin space of each particle is independent. Once the appropriate matching procedures are applied, taking the classical-spin limit we obtain the potential for general spinning bodies. When the Wilson coefficients are set to unity, we successfully reproduced the potential for the Kerr black hole. Interestingly, for finite spins, we find that the finite-spin deviations from Kerr Wilson coefficients cancel with that in the matching procedure, reproducing the Kerr potential without the need for taking the classical-spin limit. Finally, we find that when cast into the chiral basis, the spin-dependence of minimal coupling exhibits factorization, allowing us to take the classical-spin limit straight forwardly.
△ Less
Submitted 25 September, 2020; v1 submitted 22 August, 2019;
originally announced August 2019.
-
Automatic Registration between Cone-Beam CT and Scanned Surface via Deep-Pose Regression Neural Networks and Clustered Similarities
Authors:
Minyoung Chung,
Jingyu Lee,
Wisoo Song,
Youngchan Song,
Il-Hyung Yang,
Jeongjin Lee,
Yeong-Gil Shin
Abstract:
Computerized registration between maxillofacial cone-beam computed tomography (CT) images and a scanned dental model is an essential prerequisite in surgical planning for dental implants or orthognathic surgery. We propose a novel method that performs fully automatic registration between a cone-beam CT image and an optically scanned model. To build a robust and automatic initial registration metho…
▽ More
Computerized registration between maxillofacial cone-beam computed tomography (CT) images and a scanned dental model is an essential prerequisite in surgical planning for dental implants or orthognathic surgery. We propose a novel method that performs fully automatic registration between a cone-beam CT image and an optically scanned model. To build a robust and automatic initial registration method, our method applies deep-pose regression neural networks in a reduced domain (i.e., 2-dimensional image). Subsequently, fine registration is performed via optimal clusters. Majority voting system achieves globally optimal transformations while each cluster attempts to optimize local transformation parameters. The coherency of clusters determines their candidacy for the optimal cluster set. The outlying regions in the iso-surface are effectively removed based on the consensus among the optimal clusters. The accuracy of registration was evaluated by the Euclidean distance of 10 landmarks on a scanned model which were annotated by the experts in the field. The experiments show that the proposed method's registration accuracy, measured in landmark distance, outperforms other existing methods by 30.77% to 70%. In addition to achieving high accuracy, our proposed method requires neither human-interactions nor priors (e.g., iso-surface extraction). The main significance of our study is twofold: 1) the employment of light-weighted neural networks which indicates the applicability of neural network in extracting pose cues that can be easily obtained and 2) the introduction of an optimal cluster-based registration method that can avoid metal artifacts during the matching procedures.
△ Less
Submitted 29 July, 2019;
originally announced July 2019.
-
First demonstration of ionization cooling by the Muon Ionization Cooling Experiment
Authors:
M. Bogomilov,
R. Tsenov,
G. Vankova-Kirilova,
Y. P. Song,
J. Y. Tang,
Z. H. Li,
R. Bertoni,
M. Bonesini,
F. Chignoli,
R. Mazza,
V. Palladino,
A. de Bari,
D. Orestano,
L. Tortora,
Y. Kuno,
H. Sakamoto,
A. Sato,
S. Ishimoto,
M. Chung,
C. K. Sung,
F. Filthaut,
D. Jokovic,
D. Maletic,
M. Savic,
N. Jovancevic
, et al. (110 additional authors not shown)
Abstract:
High-brightness muon beams of energy comparable to those produced by state-of-the-art electron, proton and ion accelerators have yet to be realised. Such beams have the potential to carry the search for new phenomena in lepton-antilepton collisions to extremely high energy and also to provide uniquely well-characterised neutrino beams. A muon beam may be created through the decay of pions produced…
▽ More
High-brightness muon beams of energy comparable to those produced by state-of-the-art electron, proton and ion accelerators have yet to be realised. Such beams have the potential to carry the search for new phenomena in lepton-antilepton collisions to extremely high energy and also to provide uniquely well-characterised neutrino beams. A muon beam may be created through the decay of pions produced in the interaction of a proton beam with a target. To produce a high-brightness beam from such a source requires that the phase space volume occupied by the muons be reduced (cooled). Ionization cooling is the novel technique by which it is proposed to cool the beam. The Muon Ionization Cooling Experiment collaboration has constructed a section of an ionization cooling cell and used it to provide the first demonstration of ionization cooling. We present these ground-breaking measurements.
△ Less
Submitted 19 July, 2019;
originally announced July 2019.
-
Self-imitating Feedback Generation Using GAN for Computer-Assisted Pronunciation Training
Authors:
Seung Hee Yang,
Minhwa Chung
Abstract:
Self-imitating feedback is an effective and learner-friendly method for non-native learners in Computer-Assisted Pronunciation Training. Acoustic characteristics in native utterances are extracted and transplanted onto learner's own speech input, and given back to the learner as a corrective feedback. Previous works focused on speech conversion using prosodic transplantation techniques based on PS…
▽ More
Self-imitating feedback is an effective and learner-friendly method for non-native learners in Computer-Assisted Pronunciation Training. Acoustic characteristics in native utterances are extracted and transplanted onto learner's own speech input, and given back to the learner as a corrective feedback. Previous works focused on speech conversion using prosodic transplantation techniques based on PSOLA algorithm. Motivated by the visual differences found in spectrograms of native and non-native speeches, we investigated applying GAN to generate self-imitating feedback by utilizing generator's ability through adversarial training. Because this mapping is highly under-constrained, we also adopt cycle consistency loss to encourage the output to preserve the global structure, which is shared by native and non-native utterances. Trained on 97,200 spectrogram images of short utterances produced by native and non-native speakers of Korean, the generator is able to successfully transform the non-native spectrogram input to a spectrogram with properties of self-imitating feedback. Furthermore, the transformed spectrogram shows segmental corrections that cannot be obtained by prosodic transplantation. Perceptual test comparing the self-imitating and correcting abilities of our method with the baseline PSOLA method shows that the generative approach with cycle consistency loss is promising.
△ Less
Submitted 20 April, 2019;
originally announced April 2019.
-
A New Approach for Measuring the Muon Anomalous Magnetic Moment and Electric Dipole Moment
Authors:
M. Abe,
S. Bae,
G. Beer,
G. Bunce,
H. Choi,
S. Choi,
M. Chung,
W. da Silva,
S. Eidelman,
M. Finger,
Y. Fukao,
T. Fukuyama,
S. Haciomeroglu,
K. Hasegawa,
K. Hayasaka,
N. Hayashizaki,
H. Hisamatsu,
T. Iijima,
H. Iinuma,
K. Inami,
H. Ikeda,
M. Ikeno,
K. Ishida,
T. Itahashi,
M. Iwasaki
, et al. (71 additional authors not shown)
Abstract:
This paper introduces a new approach to measure the muon magnetic moment anomaly $a_μ = (g-2)/2$, and the muon electric dipole moment (EDM) $d_μ$ at the J-PARC muon facility. The goal of our experiment is to measure $a_μ$ and $d_μ$ using an independent method with a factor of 10 lower muon momentum, and a factor of 20 smaller diameter storage-ring solenoid compared with previous and ongoing muon…
▽ More
This paper introduces a new approach to measure the muon magnetic moment anomaly $a_μ = (g-2)/2$, and the muon electric dipole moment (EDM) $d_μ$ at the J-PARC muon facility. The goal of our experiment is to measure $a_μ$ and $d_μ$ using an independent method with a factor of 10 lower muon momentum, and a factor of 20 smaller diameter storage-ring solenoid compared with previous and ongoing muon $g-2$ experiments with unprecedented quality of the storage magnetic field. Additional significant differences from the present experimental method include a factor of 1,000 smaller transverse emittance of the muon beam (reaccelerated thermal muon beam), its efficient vertical injection into the solenoid, and tracking each decay positron from muon decay to obtain its momentum vector. The precision goal for $a_μ$ is statistical uncertainty of 450 part per billion (ppb), similar to the present experimental uncertainty, and a systematic uncertainty less than 70 ppb. The goal for EDM is a sensitivity of $1.5\times 10^{-21}~e\cdot\mbox{cm}$.
△ Less
Submitted 10 March, 2019; v1 submitted 10 January, 2019;
originally announced January 2019.
-
Statistical Model for Dynamically-Changing Correlation Matrices with Application to Brain Connectivity
Authors:
Shih-Gu Huang,
S. Balqis Samdin,
Chee-Ming Ting,
Hernando Ombao,
Moo K. Chung
Abstract:
Background: Recent studies have indicated that functional connectivity is dynamic even during rest. A common approach to modeling the dynamic functional connectivity in whole-brain resting-state fMRI is to compute the correlation between anatomical regions via sliding time windows. However, the direct use of the sample correlation matrices is not reliable due to the image acquisition and processin…
▽ More
Background: Recent studies have indicated that functional connectivity is dynamic even during rest. A common approach to modeling the dynamic functional connectivity in whole-brain resting-state fMRI is to compute the correlation between anatomical regions via sliding time windows. However, the direct use of the sample correlation matrices is not reliable due to the image acquisition and processing noises in resting-sate fMRI.
New method: To overcome these limitations, we propose a new statistical model that smooths out the noise by exploiting the geometric structure of correlation matrices. The dynamic correlation matrix is modeled as a linear combination of symmetric positive-definite matrices combined with cosine series representation. The resulting smoothed dynamic correlation matrices are clustered into disjoint brain connectivity states using the k-means clustering algorithm.
Results: The proposed model preserves the geometric structure of underlying physiological dynamic correlation, eliminates unwanted noise in connectivity and obtains more accurate state spaces. The difference in the estimated dynamic connectivity states between males and females is identified.
Comparison with existing methods: We demonstrate that the proposed statistical model has less rapid state changes caused by noise and improves the accuracy in identifying and discriminating different states.
Conclusions: We propose a new regression model on dynamically changing correlation matrices that provides better performance over existing windowed correlation and is more reliable for the modeling of dynamic connectivity.
△ Less
Submitted 3 November, 2019; v1 submitted 25 December, 2018;
originally announced December 2018.
-
The simplest massive S-matrix: from minimal coupling to Black Holes
Authors:
Ming-Zhi Chung,
Yu-tin Huang,
Jung-Wook Kim,
Sangmin Lee
Abstract:
In this paper, we explore the physics of electromagnetically and gravitationally coupled massive higher spin states from the on-shell point of view. Starting with the three-point amplitude, we focus on the simplest amplitude which is characterized by matching to minimal coupling in the UV. In the IR such amplitude leads to g = 2 for arbitrary charged spin states, and the best high energy behavior…
▽ More
In this paper, we explore the physics of electromagnetically and gravitationally coupled massive higher spin states from the on-shell point of view. Starting with the three-point amplitude, we focus on the simplest amplitude which is characterized by matching to minimal coupling in the UV. In the IR such amplitude leads to g = 2 for arbitrary charged spin states, and the best high energy behavior for a given spin. We proceed to construct the (gravitational) Compton amplitude for generic spins. We find that the leading deformation away from minimal coupling, in the gravitation sector, will lead to inconsistent factorizations and are thus forbidden. As the corresponding deformation in the gauge sector encodes the anomalous magnetic dipole moment, this leads to the prediction that for systems with gauge2 =gravity relations, such as perturbative string theory, all charged states must have g = 2. It is then natural to ask for generic spin, what is the theory that yields such minimal coupling. By matching to the one body effective action, remarkably we verify that for large spins, the answer is Kerr black holes. This identification is then an on-shell avatar of the no hair theorem. Finally using this identification as well as the newly constructed Compton amplitudes, we proceed to compute the spin dependent pieces for the classical potential at 2PM order up to degree four in spin operator of either black holes.
△ Less
Submitted 22 August, 2019; v1 submitted 20 December, 2018;
originally announced December 2018.
-
Rapid Acceleration of the Permutation Test via Slow Random Walks in the Permutation Group
Authors:
Moo K. Chung,
Yixian Wang,
Shih-Gu Huang,
Ilwoo Lyu
Abstract:
The permutation test is an often used test procedure in brain imaging. Unfortunately, generating every possible permutation for large-scale brain image datasets such as HCP and ADNI with hundreds images is not practical. Many previous attempts at speeding up the permutation test rely on various approximation strategies such as estimating the tail distribution with known parametric distributions. I…
▽ More
The permutation test is an often used test procedure in brain imaging. Unfortunately, generating every possible permutation for large-scale brain image datasets such as HCP and ADNI with hundreds images is not practical. Many previous attempts at speeding up the permutation test rely on various approximation strategies such as estimating the tail distribution with known parametric distributions. In this study, we show how to rapidly accelerate the permutation test without any type of approximate strategies by exploiting the underlying algebraic structure of the permutation group. The method is applied to large number of MRIs in two applications: (1) localizing the male and female differences and (2) localizing the regions of high genetic heritability in the sulcal and gyral pattern of the human cortical brain.
△ Less
Submitted 9 April, 2019; v1 submitted 17 December, 2018;
originally announced December 2018.
-
Sampled Tikhonov Regularization for Large Linear Inverse Problems
Authors:
J. Tanner Slagel,
Julianne Chung,
Matthias Chung,
David Kozak,
Luis Tenorio
Abstract:
In this paper, we investigate iterative methods that are based on sampling of the data for computing Tikhonov-regularized solutions. We focus on very large inverse problems where access to the entire data set is not possible all at once (e.g., for problems with streaming or massive datasets). Row-access methods provide an ideal framework for solving such problems, since they only require access to…
▽ More
In this paper, we investigate iterative methods that are based on sampling of the data for computing Tikhonov-regularized solutions. We focus on very large inverse problems where access to the entire data set is not possible all at once (e.g., for problems with streaming or massive datasets). Row-access methods provide an ideal framework for solving such problems, since they only require access to "blocks" of the data at any given time. However, when using these iterative sampling methods to solve inverse problems, the main challenges include a proper choice of the regularization parameter, appropriate sampling strategies, and a convergence analysis. To address these challenges, we first describe a family of sampled iterative methods that can incorporate data as they become available (e.g., randomly sampled). We consider two sampled iterative methods, where the iterates can be characterized as solutions to a sequence of approximate Tikhonov problems. The first method requires the regularization parameter to be fixed a priori and converges asymptotically to an unregularized solution for randomly sampled data. This is undesirable for inverse problems. Thus, we focus on the second method where the main benefits are that the regularization parameter can be updated during the iterative process and the iterates converge asymptotically to a Tikhonov-regularized solution. We describe adaptive approaches to update the regularization parameter that are based on sampled residuals, and we describe a limited-memory variant for larger problems. Numerical examples, including a large-scale super-resolution imaging example, demonstrate the potential for these methods.
△ Less
Submitted 14 December, 2018;
originally announced December 2018.
-
Harmonic holes as the submodules of brain network and network dissimilarity
Authors:
Hyekyoung Lee,
Moo K. Chung,
Hongyoon Choi,
Hyejin Kang,
Seunggyun Ha,
Yu Kyeong Kim,
Dong Soo Lee
Abstract:
Persistent homology has been applied to brain network analysis for finding the shape of brain networks across multiple thresholds. In the persistent homology, the shape of networks is often quantified by the sequence of $k$-dimensional holes and Betti numbers.The Betti numbers are more widely used than holes themselves in topological brain network analysis. However, the holes show the local connec…
▽ More
Persistent homology has been applied to brain network analysis for finding the shape of brain networks across multiple thresholds. In the persistent homology, the shape of networks is often quantified by the sequence of $k$-dimensional holes and Betti numbers.The Betti numbers are more widely used than holes themselves in topological brain network analysis. However, the holes show the local connectivity of networks, and they can be very informative features in analysis. In this study, we propose a new method of measuring network differences based on the dissimilarity measure of harmonic holes (HHs). The HHs, which represent the substructure of brain networks, are extracted by the Hodge Laplacian of brain networks. We also find the most contributed HHs to the network difference based on the HH dissimilarity. We applied our proposed method to clustering the networks of 4 groups, normal control (NC), stable and progressive mild cognitive impairment (sMCI and pMCI), and Alzheimer's disease (AD). The results showed that the clustering performance of the proposed method was better than that of network distances based on only the global change of topology.
△ Less
Submitted 11 November, 2018;
originally announced November 2018.
-
First particle-by-particle measurement of emittance in the Muon Ionization Cooling Experiment
Authors:
The MICE Collaboration,
D. Adams,
D. Adey,
R. Asfandiyarov,
G. Barber,
A. de Bari,
R. Bayes,
V. Bayliss,
R. Bertoni,
V. Blackmore,
A. Blondel,
J. Boehm,
M. Bogomilov,
M. Bonesini,
C. N. Booth,
D. Bowring,
S. Boyd,
T. W. Bradshaw,
A. D. Bross,
C. Brown,
L. Coney,
G. Charnley,
G. T. Chatzitheodoridis,
F. Chignoli,
M. Chung
, et al. (111 additional authors not shown)
Abstract:
The Muon Ionization Cooling Experiment (MICE) collaboration seeks to demonstrate the feasibility of ionization cooling, the technique by which it is proposed to cool the muon beam at a future neutrino factory or muon collider. The emittance is measured from an ensemble of muons assembled from those that pass through the experiment. A pure muon ensemble is selected using a particle-identification s…
▽ More
The Muon Ionization Cooling Experiment (MICE) collaboration seeks to demonstrate the feasibility of ionization cooling, the technique by which it is proposed to cool the muon beam at a future neutrino factory or muon collider. The emittance is measured from an ensemble of muons assembled from those that pass through the experiment. A pure muon ensemble is selected using a particle-identification system that can reject efficiently both pions and electrons. The position and momentum of each muon are measured using a high-precision scintillating-fibre tracker in a 4\,T solenoidal magnetic field. This paper presents the techniques used to reconstruct the phase-space distributions and reports the first particle-by-particle measurement of the emittance of the MICE Muon Beam as a function of muon-beam momentum.
△ Less
Submitted 26 March, 2019; v1 submitted 31 October, 2018;
originally announced October 2018.
-
Hyper $b$-ary expansions and Stern polynomials
Authors:
Tanay Wakhare,
Caleb Kendrick,
Matthew Chung,
Catherine Cassell,
Stefano Santini,
William Colin Mosley,
Anand Raghu,
Robert Morrison,
Iman Schurman,
Timothy Kevin Beal,
Matthew Patrick
Abstract:
We study a recently introduced base $b$ polynomial analog of Stern's diatomic sequence, which generalizes Stern polynomials of Klavar, Dilcher, Ericksen, Mansour, Stolarsky, and others. We lift some basic properties of base $2$ Stern polynomials to arbitrary base, and introduce a matrix characterization of Stern polynomials. By specializing, we recover some new number theoretic results about hyper…
▽ More
We study a recently introduced base $b$ polynomial analog of Stern's diatomic sequence, which generalizes Stern polynomials of Klavar, Dilcher, Ericksen, Mansour, Stolarsky, and others. We lift some basic properties of base $2$ Stern polynomials to arbitrary base, and introduce a matrix characterization of Stern polynomials. By specializing, we recover some new number theoretic results about hyper $b$-ary partitions, which count partitions of $n$ into powers of $b$.
△ Less
Submitted 25 October, 2018;
originally announced October 2018.
-
Experimental observation of proton bunch modulation in a plasma, at varying plasma densities
Authors:
E. Adli,
A. Ahuja,
O. Apsimon,
R. Apsimon,
A. -M. Bachmann,
D. Barrientos,
M. M. Barros,
J. Batkiewicz,
F. Batsch,
J. Bauche,
V. K. Berglyd Olsen,
M. Bernardini,
B. Biskup,
A. Boccardi,
T. Bogey,
T. Bohl,
C. Bracco,
F. Braunmüller,
S. Burger,
G. Burt,
S. Bustamante,
B. Buttenschön,
A. Caldwell,
M. Cascella,
J. Chappell
, et al. (87 additional authors not shown)
Abstract:
We give direct experimental evidence for the observation of the full transverse self-modulation of a relativistic proton bunch propagating through a dense plasma. The bunch exits the plasma with a density modulation resulting from radial wakefield effects with a period reciprocal to the plasma frequency. We show that the modulation is seeded by using an intense laser pulse co-propagating with the…
▽ More
We give direct experimental evidence for the observation of the full transverse self-modulation of a relativistic proton bunch propagating through a dense plasma. The bunch exits the plasma with a density modulation resulting from radial wakefield effects with a period reciprocal to the plasma frequency. We show that the modulation is seeded by using an intense laser pulse co-propagating with the proton bunch which creates a relativistic ionization front within the bunch. We show by varying the plasma density over one order of magnitude that the modulation period scales with the expected dependence on the plasma density.
△ Less
Submitted 1 April, 2019; v1 submitted 12 September, 2018;
originally announced September 2018.
-
Topological Brain Network Distances
Authors:
Moo K. Chung,
Hyekyoung Lee,
Andrey Gritsenko,
Alex DiChristofano,
Dustin Pluta,
Hernando Ombao,
Victor Solo
Abstract:
Existing brain network distances are often based on matrix norms. The element-wise differences in the existing matrix norms may fail to capture underlying topological differences. Further, matrix norms are sensitive to outliers. A major disadvantage to element-wise distance calculations is that it could be severely affected even by a small number of extreme edge weights. Thus it is necessary to de…
▽ More
Existing brain network distances are often based on matrix norms. The element-wise differences in the existing matrix norms may fail to capture underlying topological differences. Further, matrix norms are sensitive to outliers. A major disadvantage to element-wise distance calculations is that it could be severely affected even by a small number of extreme edge weights. Thus it is necessary to develop network distances that recognize topology. In this paper, we provide a survey of bottleneck, Gromov-Hausdorff (GH) and Kolmogorov-Smirnov (KS) distances that are adapted for brain networks, and compare them against matrix-norm based network distances. Bottleneck and GH-distances are often used in persistent homology. However, they were rarely utilized to measure similarity between brain networks. KS-distance is recently introduced to measure the similarity between networks across different filtration values. The performance analysis was conducted using the random network simulations with the ground truths. Using a twin imaging study, which provides biological ground truth, we demonstrate that the KS distance has the ability to determine heritability.
△ Less
Submitted 9 September, 2018;
originally announced September 2018.
-
Acceleration of electrons in the plasma wakefield of a proton bunch
Authors:
The AWAKE Collaboration,
E. Adli,
A. Ahuja,
O. Apsimon,
R. Apsimon,
A. -M. Bachmann,
D. Barrientos,
F. Batsch,
J. Bauche,
V. K. Berglyd Olsen,
M. Bernardini,
T. Bohl,
C. Bracco,
F. Braunmueller,
G. Burt,
B. Buttenschoen,
A. Caldwell,
M. Cascella,
J. Chappell,
E. Chevallay,
M. Chung,
D. Cooke,
H. Damerau,
L. Deacon,
L. H. Deubner
, et al. (69 additional authors not shown)
Abstract:
High energy particle accelerators have been crucial in providing a deeper understanding of fundamental particles and the forces that govern their interactions. In order to increase the energy or reduce the size of the accelerator, new acceleration schemes need to be developed. Plasma wakefield acceleration, in which the electrons in a plasma are excited, leading to strong electric fields, is one s…
▽ More
High energy particle accelerators have been crucial in providing a deeper understanding of fundamental particles and the forces that govern their interactions. In order to increase the energy or reduce the size of the accelerator, new acceleration schemes need to be developed. Plasma wakefield acceleration, in which the electrons in a plasma are excited, leading to strong electric fields, is one such promising novel acceleration technique. Pioneering experiments have shown that an intense laser pulse or electron bunch traversing a plasma, drives electric fields of 10s GV/m and above. These values are well beyond those achieved in conventional RF accelerators which are limited to ~0.1 GV/m. A limitation of laser pulses and electron bunches is their low stored energy, which motivates the use of multiple stages to reach very high energies. The use of proton bunches is compelling, as they have the potential to drive wakefields and accelerate electrons to high energy in a single accelerating stage. The long proton bunches currently available can be used, as they undergo self-modulation, a particle-plasma interaction which longitudinally splits the bunch into a series of high density microbunches, which then act resonantly to create large wakefields. The AWAKE experiment at CERN uses intense bunches of protons, each of energy 400 GeV, with a total bunch energy of 19 kJ, to drive a wakefield in a 10 m long plasma. Bunches of electrons are injected into the wakefield formed by the proton microbunches. This paper presents measurements of electrons accelerated up to 2 GeV at AWAKE. This constitutes the first demonstration of proton-driven plasma wakefield acceleration. The potential for this scheme to produce very high energy electron bunches in a single accelerating stage means that the results shown here are a significant step towards the development of future high energy particle accelerators.
△ Less
Submitted 11 October, 2018; v1 submitted 29 August, 2018;
originally announced August 2018.
-
Relation between particle beam eigen-emittances and canonical vorticity flux
Authors:
Lars Groening,
Moses Chung
Abstract:
Particle beams provided by accelerators occupy a finite volume of the four dimensional transverse phase space. The latter is spanned by the four degrees of freedom, i.e., horizontal/vertical position and momentum. This volume is referred to as emittance. Horizontal and vertical emittances are obtained through projections onto the two transverse sub-phase spaces. Eigen-emittances are obtained from…
▽ More
Particle beams provided by accelerators occupy a finite volume of the four dimensional transverse phase space. The latter is spanned by the four degrees of freedom, i.e., horizontal/vertical position and momentum. This volume is referred to as emittance. Horizontal and vertical emittances are obtained through projections onto the two transverse sub-phase spaces. Eigen-emittances are obtained from the latter by removing all horizontal-vertical correlations through an appropriate beam optics section. Canonical vorticity flux is used for instance for modelling dynamics of tubes formed by magnetic field lines and particle currents embedded into plasmas. This report is on the relation of eigen-emittances and canonical vorticity flux. Change of beam eigen-emittances is equivalent to change of beam canonical vorticity flux.
△ Less
Submitted 24 December, 2020; v1 submitted 2 August, 2018;
originally announced August 2018.
-
Deeply Self-Supervised Contour Embedded Neural Network Applied to Liver Segmentation
Authors:
Minyoung Chung,
Jingyu Lee,
Minkyung Lee,
Jeongjin Lee,
Yeong-Gil Shin
Abstract:
Objective: Herein, a neural network-based liver segmentation algorithm is proposed, and its performance was evaluated using abdominal computed tomography (CT) images. Methods: A fully convolutional network was developed to overcome the volumetric image segmentation problem. To guide a neural network to accurately delineate a target liver object, the network was deeply supervised by applying the ad…
▽ More
Objective: Herein, a neural network-based liver segmentation algorithm is proposed, and its performance was evaluated using abdominal computed tomography (CT) images. Methods: A fully convolutional network was developed to overcome the volumetric image segmentation problem. To guide a neural network to accurately delineate a target liver object, the network was deeply supervised by applying the adaptive self-supervision scheme to derive the essential contour, which acted as a complement with the global shape. The discriminative contour, shape, and deep features were internally merged for the segmentation results. Results and Conclusion: 160 abdominal CT images were used for training and validation. The quantitative evaluation of the proposed network was performed through an eight-fold cross-validation. The result showed that the method, which uses the contour feature, segmented the liver more accurately than the state-of-the-art with a 2.13% improvement in the dice score. Significance: In this study, a new framework was introduced to guide a neural network and learn complementary contour features. The proposed neural network demonstrates that the guided contour features can significantly improve the performance of the segmentation task.
△ Less
Submitted 17 October, 2019; v1 submitted 2 August, 2018;
originally announced August 2018.
-
Development of a globally optimised model of the cerebral arteries
Authors:
Jonathan Keelan,
Emma M. L. Chung,
James P. Hague
Abstract:
The cerebral arteries are difficult to reproduce from first principles, featuring interwoven territories, and intricate layers of grey and white matter with differing metabolic demand. The aim of this study was to identify the ideal configuration of arteries required to sustain an entire brain hemisphere based on minimisation of the energy required to supply the tissue. The 3D distribution of grey…
▽ More
The cerebral arteries are difficult to reproduce from first principles, featuring interwoven territories, and intricate layers of grey and white matter with differing metabolic demand. The aim of this study was to identify the ideal configuration of arteries required to sustain an entire brain hemisphere based on minimisation of the energy required to supply the tissue. The 3D distribution of grey and white matter within a healthy human brain was first segmented from Magnetic Resonance Images. A novel simulated annealing algorithm was then applied to determine the optimal configuration of arteries required to supply brain tissue. The model is validated through comparison of this ideal, entirely optimised, brain vasculature with the known structure of real arteries. This establishes that the human cerebral vasculature is highly optimised; closely resembling the most energy efficient arrangement of vessels. In addition to local adherence to fluid dynamics optimisation principles, the optimised vasculature reproduces global brain perfusion territories with well defined boundaries between anterior, middle and posterior regions. This validated brain vascular model and algorithm can be used for patient-specific modelling of stroke and cerebral haemodynamics, identification of sub-optimal conditions associated with vascular disease, and optimising vascular structures for tissue engineering and artificial organ design.
△ Less
Submitted 11 December, 2018; v1 submitted 20 July, 2018;
originally announced July 2018.
-
Operation of normal-conducting RF cavities in multi-tesla magnetic fields for muon ionization cooling: a feasibility demonstration
Authors:
Daniel Bowring,
Alexey Kochemirovskiy,
Yagmur Torun,
Chris Adolphsen,
Alan Bross,
Moses Chung,
Ben Freemire,
Lixin Ge,
Andrew Haase,
Peter Lane,
Maria Leonova,
Derun Li,
Zenghai Li,
Ao Liu,
Tianhuan Luo,
David Martin,
Alfred Moretti,
David Neuffer,
Ralph Pasquinelli,
Mark Palmer,
David Peterson,
Milorad Popovic,
Diktys Stratakis,
Katsuya Yonehara
Abstract:
Ionization cooling is the preferred method for producing bright muon beams. This cooling technique requires the operation of normal conducting, radio-frequency (RF) accelerating cavities within the multi-tesla fields of DC solenoid magnets. Under these conditions, cavities exhibit increased susceptibility to RF breakdown, which can damage channel components and imposes limits on channel length and…
▽ More
Ionization cooling is the preferred method for producing bright muon beams. This cooling technique requires the operation of normal conducting, radio-frequency (RF) accelerating cavities within the multi-tesla fields of DC solenoid magnets. Under these conditions, cavities exhibit increased susceptibility to RF breakdown, which can damage channel components and imposes limits on channel length and transmission efficiency. We present a solution to the problem of breakdown in strong magnetic fields. We report, for the first time, stable high-vacuum, copper cavity operation at gradients above 50 MV/m and in an external magnetic field of three tesla. This eliminates a significant technical risk that has previously been inherent in ionization cooling channel designs.
△ Less
Submitted 10 July, 2018;
originally announced July 2018.
-
Exact Combinatorial Inference for Brain Images
Authors:
Moo K. Chung,
Zhan Luo,
Alex D. Leow,
Andrew L. Alexander,
Richard J. Davidson,
H. Hill Goldsmith
Abstract:
The permutation test is known as the exact test procedure in statistics. However, often it is not exact in practice and only an approximate method since only a small fraction of every possible permutation is generated. Even for a small sample size, it often requires to generate tens of thousands permutations, which can be a serious computational bottleneck. In this paper, we propose a novel combin…
▽ More
The permutation test is known as the exact test procedure in statistics. However, often it is not exact in practice and only an approximate method since only a small fraction of every possible permutation is generated. Even for a small sample size, it often requires to generate tens of thousands permutations, which can be a serious computational bottleneck. In this paper, we propose a novel combinatorial inference procedure that enumerates all possible permutations combinatorially without any resampling. The proposed method is validated against the standard permutation test in simulation studies with the ground truth. The method is further applied in twin DTI study in determining the genetic contribution of the minimum spanning tree of the structural brain connectivity.
△ Less
Submitted 8 July, 2018;
originally announced July 2018.
-
Determination of lattice functions and 2nd-order transfer matrix for High Power Cyclotron
Authors:
B-H. Oh,
C. Choi,
G-R. Hahn,
M. Chung,
S. Shin
Abstract:
This paper describes the development of beam dynamic simulation code for cyclotron. Starting from a description of beam dynamics in the cyclotron, lattice functions were determined and the solutions for the 2nd-order nonlinear Hamiltonian were revised. Based on the description of beam dynamics in the cyclotron, simulation code was also developed for cyclotron design.
This paper describes the development of beam dynamic simulation code for cyclotron. Starting from a description of beam dynamics in the cyclotron, lattice functions were determined and the solutions for the 2nd-order nonlinear Hamiltonian were revised. Based on the description of beam dynamics in the cyclotron, simulation code was also developed for cyclotron design.
△ Less
Submitted 3 July, 2018;
originally announced July 2018.
-
Automatic Identification of Twin Zygosity in Resting-State Functional MRI
Authors:
Andrey Gritsenko,
Martin A. Lindquist,
Gregory R. Kirk,
Moo K. Chung
Abstract:
A key strength of twin studies arises from the fact that there are two types of twins, monozygotic and dizygotic, that share differing amounts of genetic information. Accurate differentiation of twin types allows efficient inference on genetic influences in a population. However, identification of zygosity is often prone to errors without genotying. In this study, we propose a novel pairwise featu…
▽ More
A key strength of twin studies arises from the fact that there are two types of twins, monozygotic and dizygotic, that share differing amounts of genetic information. Accurate differentiation of twin types allows efficient inference on genetic influences in a population. However, identification of zygosity is often prone to errors without genotying. In this study, we propose a novel pairwise feature representation to classify the zygosity of twin pairs of resting state functional magnetic resonance images (rs-fMRI). For this, we project an fMRI signal to a set of basis functions and use the projection coefficients as the compact and discriminative feature representation of noisy fMRI. We encode the relationship between twins as the correlation between the new feature representations across brain regions. We employ hill climbing variable selection to identify brain regions that are the most genetically affected. The proposed framework was applied to 208 twin pairs and achieved 94.19% classification accuracy in automatically identifying the zygosity of paired images.
△ Less
Submitted 26 October, 2018; v1 submitted 30 June, 2018;
originally announced July 2018.
-
CRAM: Clued Recurrent Attention Model
Authors:
Minki Chung,
Sungzoon Cho
Abstract:
To overcome the poor scalability of convolutional neural network, recurrent attention model(RAM) selectively choose what and where to look on the image. By directing recurrent attention model how to look the image, RAM can be even more successful in that the given clue narrow down the scope of the possible focus zone. In this perspective, this work proposes clued recurrent attention model (CRAM) w…
▽ More
To overcome the poor scalability of convolutional neural network, recurrent attention model(RAM) selectively choose what and where to look on the image. By directing recurrent attention model how to look the image, RAM can be even more successful in that the given clue narrow down the scope of the possible focus zone. In this perspective, this work proposes clued recurrent attention model (CRAM) which add clue or constraint on the RAM better problem solving. CRAM follows encoder-decoder framework, encoder utilizes recurrent attention model with spatial transformer network and decoder which varies depending on the task. To ensure the performance, CRAM tackles two computer vision task. One is the image classification task, with clue given as the binary image saliency which indicates the approximate location of object. The other is the inpainting task, with clue given as binary mask which indicates the occluded part. In both tasks, CRAM shows better performance than existing methods showing the successful extension of RAM.
△ Less
Submitted 28 April, 2018;
originally announced April 2018.
-
An MCMC Algorithm for Estimating the Q-matrix in a Bayesian Framework
Authors:
Mengta Chung,
Matthew S. Johnson
Abstract:
The purpose of this research is to develop an MCMC algorithm for estimating the Q-matrix. Based on the DINA model, the algorithm starts with estimating correlated attributes. Using a saturated model and a binary decimal conversion, the algorithm transforms possible attribute patterns to a Multinomial distribution. Along with the likelihood of an attribute pattern, a Dirichlet distribution, constru…
▽ More
The purpose of this research is to develop an MCMC algorithm for estimating the Q-matrix. Based on the DINA model, the algorithm starts with estimating correlated attributes. Using a saturated model and a binary decimal conversion, the algorithm transforms possible attribute patterns to a Multinomial distribution. Along with the likelihood of an attribute pattern, a Dirichlet distribution, constructed using Gamma distributions, is used as the prior to sample from the posterior. Correlated attributes of examinees are generated using inverse transform sampling. Closed form posteriors for sampling guess and slip parameters are found. A distribution for sampling the Q-matrix is derived. A relabeling algorithm that accounts for potential label switching is presented. A method for simulating data with correlated attributes for the DINA model is offered. Three simulation studies are conducted to evaluate the performance of the algorithm. An empirical study using the ECPE data is performed. The algorithm is implemented using customized R codes.
△ Less
Submitted 6 February, 2018;
originally announced February 2018.
-
Parameter and Uncertainty Estimation for Dynamical Systems Using Surrogate Stochastic Processes
Authors:
M. Chung,
M. Binois,
R. B. Gramacy,
D. J. Moquin,
A. P. Smith,
A. M. Smith
Abstract:
Inference on unknown quantities in dynamical systems via observational data is essential for providing meaningful insight, furnishing accurate predictions, enabling robust control, and establishing appropriate designs for future experiments. Merging mathematical theory with empirical measurements in a statistically coherent way is critical and challenges abound, e.g.,: ill-posedness of the paramet…
▽ More
Inference on unknown quantities in dynamical systems via observational data is essential for providing meaningful insight, furnishing accurate predictions, enabling robust control, and establishing appropriate designs for future experiments. Merging mathematical theory with empirical measurements in a statistically coherent way is critical and challenges abound, e.g.,: ill-posedness of the parameter estimation problem, proper regularization and incorporation of prior knowledge, and computational limitations on full uncertainty qualification. To address these issues, we propose a new method for learning parameterized dynamical systems from data. In many ways, our proposal turns the canonical framework on its head. We first fit a surrogate stochastic process to observational data, enforcing prior knowledge (e.g., smoothness), and coping with challenging data features like heteroskedasticity, heavy tails and censoring. Then, samples of the stochastic process are used as "surrogate data" and point estimates are computed via ordinary point estimation methods in a modular fashion. An attractive feature of this approach is that it is fully Bayesian and simultaneously parallelizable. We demonstrate the advantages of our new approach on a predator prey simulation study and on a real world application involving within-host influenza virus infection data paired with a viral kinetic model.
△ Less
Submitted 2 February, 2018;
originally announced February 2018.
-
Pressurized rf cavities in ionizing beams
Authors:
B. Freemire,
A. V. Tollestrup,
K. Yonehara,
M. Chung,
Y. Torun,
R. P. Johnson,
G. Flanagan,
P. M. Hanlet,
M. G. Collura,
M. R. Jana,
M. Leonova,
A. Moretti,
T. Schwarz
Abstract:
A muon collider or Higgs factory requires significant reduction of the six dimensional emittance of the beam prior to acceleration. One method to accomplish this involves building a cooling channel using high pressure gas filled radio frequency cavities. The performance of such a cavity when subjected to an intense particle beam must be investigated before this technology can be validated. To this…
▽ More
A muon collider or Higgs factory requires significant reduction of the six dimensional emittance of the beam prior to acceleration. One method to accomplish this involves building a cooling channel using high pressure gas filled radio frequency cavities. The performance of such a cavity when subjected to an intense particle beam must be investigated before this technology can be validated. To this end, a high pressure gas filled radio frequency (rf) test cell was built and placed in a 400 MeV beam line from the Fermilab linac to study the plasma evolution and its effect on the cavity. Hydrogen, deuterium, helium and nitrogen gases were studied. Additionally, sulfur hexafluoride and dry air were used as dopants to aid in the removal of plasma electrons. Measurements were made using a variety of beam intensities, gas pressures, dopant concentrations, and cavity rf electric fields, both with and without a 3 T external solenoidal magnetic field. Energy dissipation per electron-ion pair, electron-ion recombination rates, ion-ion recombination rates, and electron attachment times to $SF_6$ and $O_2$ were measured.
△ Less
Submitted 5 January, 2018;
originally announced January 2018.
-
Squeezed Convolutional Variational AutoEncoder for Unsupervised Anomaly Detection in Edge Device Industrial Internet of Things
Authors:
Dohyung Kim,
Hyochang Yang,
Minki Chung,
Sungzoon Cho
Abstract:
In this paper, we propose Squeezed Convolutional Variational AutoEncoder (SCVAE) for anomaly detection in time series data for Edge Computing in Industrial Internet of Things (IIoT). The proposed model is applied to labeled time series data from UCI datasets for exact performance evaluation, and applied to real world data for indirect model performance comparison. In addition, by comparing the mod…
▽ More
In this paper, we propose Squeezed Convolutional Variational AutoEncoder (SCVAE) for anomaly detection in time series data for Edge Computing in Industrial Internet of Things (IIoT). The proposed model is applied to labeled time series data from UCI datasets for exact performance evaluation, and applied to real world data for indirect model performance comparison. In addition, by comparing the models before and after applying Fire Modules from SqueezeNet, we show that model size and inference times are reduced while similar levels of performance is maintained.
△ Less
Submitted 18 December, 2017;
originally announced December 2017.
-
The Experimental Program for High Pressure Gas Filled Radio Frequency Cavities for Muon Cooling Channels
Authors:
Ben Freemire,
Moses Chung,
Pierrick M. Hanlet,
Rolland P. Johnson,
Alfred Moretti,
Yagmur Torun,
Katsuya Yonehara
Abstract:
An intense beam of muons is needed to provide a luminosity on the order of 10$^{34}$ cm$^{-2}$s$^{-1}$ for a multi-TeV collider. Because muons produced by colliding a multi-MW proton beam with a target made of carbon or mercury have a large phase space, significant six dimensional cooling is required. Through ionization cooling - the only cooling method that works within the lifetime of the muon -…
▽ More
An intense beam of muons is needed to provide a luminosity on the order of 10$^{34}$ cm$^{-2}$s$^{-1}$ for a multi-TeV collider. Because muons produced by colliding a multi-MW proton beam with a target made of carbon or mercury have a large phase space, significant six dimensional cooling is required. Through ionization cooling - the only cooling method that works within the lifetime of the muon - and emittance exchange, the desired emittances for a Higgs Factory or higher energy collider are attainable. A cooling channel utilizing gas filled radio frequency cavities has been designed to deliver the requisite cool muon beam. Technology development of these RF cavities has progressed from breakdown studies, through beam tests, to dielectric loaded and reentrant cavity designs. The results of these experiments are summarized.
△ Less
Submitted 22 December, 2017; v1 submitted 26 October, 2017;
originally announced October 2017.
-
An MCMC Algorithm for Estimating the Reduced RUM
Authors:
Meng-Ta Chung,
Matthew S. Johnson
Abstract:
The RRUM is a model that is frequently seen in language assessment studies. The objective of this research is to advance an MCMC algorithm for the Bayesian RRUM. The algorithm starts with estimating correlated attributes. Using a saturated model and a binary decimal conversion, the algorithm transforms possible attribute patterns to a Multinomial distribution. Along with the likelihood of an attri…
▽ More
The RRUM is a model that is frequently seen in language assessment studies. The objective of this research is to advance an MCMC algorithm for the Bayesian RRUM. The algorithm starts with estimating correlated attributes. Using a saturated model and a binary decimal conversion, the algorithm transforms possible attribute patterns to a Multinomial distribution. Along with the likelihood of an attribute pattern, a Dirichlet distribution is used as the prior to sample from the posterior. The Dirichlet distribution is constructed using Gamma distributions. Correlated attributes of examinees are generated using the inverse transform sampling. Model parameters are estimated using the Metropolis within Gibbs sampler sequentially. Two simulation studies are conducted to evaluate the performance of the algorithm. The first simulation uses a complete and balanced Q-matrix that measures 5 attributes. Comprised of 28 items and 9 attributes, the Q-matrix for the second simulation is incomplete and imbalanced. The empirical study uses the ECPE data obtained from the CDM R package. Parameter estimates from the MCMC algorithm and from the CDM R package are presented and compared. The algorithm developed in this research is implemented in R.
△ Less
Submitted 23 October, 2017;
originally announced October 2017.
-
Heat Kernel Smoothing in Irregular Image Domains
Authors:
Moo K. Chung,
Yanli Wang,
Gurong Wu
Abstract:
We present the discrete version of heat kernel smoothing on graph data structure. The method is used to smooth data in an irregularly shaped domains in 3D images.
New statistical properties are derived. As an application, we show how to filter out data in the lung blood vessel trees obtained from computed tomography. The method can be further used in representing the complex vessel trees paramet…
▽ More
We present the discrete version of heat kernel smoothing on graph data structure. The method is used to smooth data in an irregularly shaped domains in 3D images.
New statistical properties are derived. As an application, we show how to filter out data in the lung blood vessel trees obtained from computed tomography. The method can be further used in representing the complex vessel trees parametrically and extracting the skeleton representation of the trees.
△ Less
Submitted 21 October, 2017;
originally announced October 2017.
-
Statistical Challenges of Big Brain Network Data
Authors:
Moo K. Chung
Abstract:
We explore the main characteristics of big brain network data that offer unique statistical challenges. The brain networks are biologically expected to be both sparse and hierarchical. Such unique characterizations put specific topological constraints onto statistical approaches and models we can use effectively. We explore the limitations of the current models used in the field and offer alternat…
▽ More
We explore the main characteristics of big brain network data that offer unique statistical challenges. The brain networks are biologically expected to be both sparse and hierarchical. Such unique characterizations put specific topological constraints onto statistical approaches and models we can use effectively. We explore the limitations of the current models used in the field and offer alternative approaches and explain new challenges.
△ Less
Submitted 23 December, 2017; v1 submitted 27 September, 2017;
originally announced September 2017.
-
Extension of Busch's Theorem to Particle Beams
Authors:
L. Groening,
M. Chung,
C. Xiao
Abstract:
In 1926, H. Busch formulated a theorem for one single charged particle moving along a region with a longitudinal magnetic field [H. Busch, Berechnung der Bahn von Kathodenstrahlen in axial symmetrischen electromagnetischen Felde, Z. Phys. 81 (5) p. 974, (1926)]. The theorem relates particle angular momentum to the amount of field lines being enclosed by the particle cyclotron motion. This paper ex…
▽ More
In 1926, H. Busch formulated a theorem for one single charged particle moving along a region with a longitudinal magnetic field [H. Busch, Berechnung der Bahn von Kathodenstrahlen in axial symmetrischen electromagnetischen Felde, Z. Phys. 81 (5) p. 974, (1926)]. The theorem relates particle angular momentum to the amount of field lines being enclosed by the particle cyclotron motion. This paper extends the theorem to many particles forming a beam without cylindrical symmetry. A quantity being preserved is derived, which represents the sum of difference of eigen-emittances, magnetic flux through the beam area, and beam rms-vorticity multiplied by the magnetic flux. Tracking simulations and analytical calculations using the generalized Courant-Snyder formalism confirm the validity of the extended theorem. The new theorem has been applied for fast modelling of experiments with electron and ion beams on transverse emittance re-partitioning conducted at FERMILAB and at GSI.
△ Less
Submitted 27 September, 2017;
originally announced September 2017.
-
Optimal Experimental Design for Constrained Inverse Problems
Authors:
Lars Ruthotto,
Julianne Chung,
Matthias Chung
Abstract:
In this paper, we address the challenging problem of optimal experimental design (OED) of constrained inverse problems. We consider two OED formulations that allow reducing the experimental costs by minimizing the number of measurements. The first formulation assumes a fine discretization of the design parameter space and uses sparsity promoting regularization to obtain an efficient design. The se…
▽ More
In this paper, we address the challenging problem of optimal experimental design (OED) of constrained inverse problems. We consider two OED formulations that allow reducing the experimental costs by minimizing the number of measurements. The first formulation assumes a fine discretization of the design parameter space and uses sparsity promoting regularization to obtain an efficient design. The second formulation parameterizes the design and seeks optimal placement for these measurements by solving a small-dimensional optimization problem. We consider both problems in a Bayes risk as well as an empirical Bayes risk minimization framework. For the unconstrained inverse state problem, we exploit the closed form solution for the inner problem to efficiently compute derivatives for the outer OED problem. The empirical formulation does not require an explicit solution of the inverse problem and therefore allows to integrate constraints efficiently. A key contribution is an efficient optimization method for solving the resulting, typically high-dimensional, bilevel optimization problem using derivative-based methods. To overcome the lack of non-differentiability in active set methods for inequality constraints problems, we use a relaxed interior point method. To address the growing computational complexity of empirical Bayes OED, we parallelize the computation over the training models. Numerical examples and illustrations from tomographic reconstruction, for various data sets and under different constraints, demonstrate the impact of constraints on the optimal design and highlight the importance of OED for constrained problems.
△ Less
Submitted 15 August, 2017;
originally announced August 2017.
-
Large Deviation Principle for $S$-unimodal maps with flat critical point
Authors:
Yong Moo Chung,
Hiroki Takahasi
Abstract:
We study a topologically exact, negative Schwarzian unimodal map whose critical point is non-recurrent and flat. Assuming the critical order is either logarithmic or polynomial, we establish the Large Deviation Principle and give a partial description of the zeros of the corresponding rate functions. We apply our main results to a certain parametrized family of unimodal maps in the same topologica…
▽ More
We study a topologically exact, negative Schwarzian unimodal map whose critical point is non-recurrent and flat. Assuming the critical order is either logarithmic or polynomial, we establish the Large Deviation Principle and give a partial description of the zeros of the corresponding rate functions. We apply our main results to a certain parametrized family of unimodal maps in the same topological conjugacy class, and give a complete description of the zeros of the rate functions. We observe a qualitative change at a transition parameter, and show that the sets of zeros depend continuously on the parameter even at the transition.
△ Less
Submitted 17 December, 2017; v1 submitted 11 August, 2017;
originally announced August 2017.
-
AWAKE readiness for the study of the seeded self-modulation of a 400\,GeV proton bunch
Authors:
P. Muggli,
E. Adli,
R. Apsimon,
F. Asmus,
R. Baartman,
A. -M. Bachmann,
M. Barros Marin,
F. Batsch,
J. Bauche,
V. K. Berglyd Olsen,
M. Bernardini,
B. Biskup,
A. Boccardi,
T. Bogey,
T. Bohl,
C. Bracco,
F. Braunmuller,
S. Burger,
G. Burt,
S. Bustamante,
B. Buttenschon,
A. Butterworth,
A. Caldwell,
M. Cascella,
E. Chevallay
, et al. (82 additional authors not shown)
Abstract:
AWAKE is a proton-driven plasma wakefield acceleration experiment. % We show that the experimental setup briefly described here is ready for systematic study of the seeded self-modulation of the 400\,GeV proton bunch in the 10\,m-long rubidium plasma with density adjustable from 1 to 10$\times10^{14}$\,cm$^{-3}$. % We show that the short laser pulse used for ionization of the rubidium vapor propag…
▽ More
AWAKE is a proton-driven plasma wakefield acceleration experiment. % We show that the experimental setup briefly described here is ready for systematic study of the seeded self-modulation of the 400\,GeV proton bunch in the 10\,m-long rubidium plasma with density adjustable from 1 to 10$\times10^{14}$\,cm$^{-3}$. % We show that the short laser pulse used for ionization of the rubidium vapor propagates all the way along the column, suggesting full ionization of the vapor. % We show that ionization occurs along the proton bunch, at the laser time and that the plasma that follows affects the proton bunch. %
△ Less
Submitted 3 August, 2017;
originally announced August 2017.
-
Comparison of the Deep-Learning-Based Automated Segmentation Methods for the Head Sectioned Images of the Virtual Korean Human Project
Authors:
Mohammad Eshghi,
Holger R. Roth,
Masahiro Oda,
Min Suk Chung,
Kensaku Mori
Abstract:
This paper presents an end-to-end pixelwise fully automated segmentation of the head sectioned images of the Visible Korean Human (VKH) project based on Deep Convolutional Neural Networks (DCNNs). By converting classification networks into Fully Convolutional Networks (FCNs), a coarse prediction map, with smaller size than the original input image, can be created for segmentation purposes. To refi…
▽ More
This paper presents an end-to-end pixelwise fully automated segmentation of the head sectioned images of the Visible Korean Human (VKH) project based on Deep Convolutional Neural Networks (DCNNs). By converting classification networks into Fully Convolutional Networks (FCNs), a coarse prediction map, with smaller size than the original input image, can be created for segmentation purposes. To refine this map and to obtain a dense pixel-wise output, standard FCNs use deconvolution layers to upsample the coarse map. However, upsampling based on deconvolution increases the number of network parameters and causes loss of detail because of interpolation. On the other hand, dilated convolution is a new technique introduced recently that attempts to capture multi-scale contextual information without increasing the network parameters while keeping the resolution of the prediction maps high. We used both a standard FCN and a dilated convolution based FCN for semantic segmentation of the head sectioned images of the VKH dataset. Quantitative results showed approximately 20% improvement in the segmentation accuracy when using FCNs with dilated convolutions.
△ Less
Submitted 15 March, 2017;
originally announced March 2017.
-
Stochastic Newton and Quasi-Newton Methods for Large Linear Least-squares Problems
Authors:
Julianne Chung,
Matthias Chung,
J. Tanner Slagel,
Luis Tenorio
Abstract:
We describe stochastic Newton and stochastic quasi-Newton approaches to efficiently solve large linear least-squares problems where the very large data sets present a significant computational burden (e.g., the size may exceed computer memory or data are collected in real-time). In our proposed framework, stochasticity is introduced in two different frameworks as a means to overcome these computat…
▽ More
We describe stochastic Newton and stochastic quasi-Newton approaches to efficiently solve large linear least-squares problems where the very large data sets present a significant computational burden (e.g., the size may exceed computer memory or data are collected in real-time). In our proposed framework, stochasticity is introduced in two different frameworks as a means to overcome these computational limitations, and probability distributions that can exploit structure and/or sparsity are considered. Theoretical results on consistency of the approximations for both the stochastic Newton and the stochastic quasi-Newton methods are provided. The results show, in particular, that stochastic Newton iterates, in contrast to stochastic quasi-Newton iterates, may not converge to the desired least-squares solution. Numerical examples, including an example from extreme learning machines, demonstrate the potential applications of these methods.
△ Less
Submitted 23 February, 2017;
originally announced February 2017.
-
Topological Distances between Networks and Its Application to Brain Imaging
Authors:
Hyekyoung Lee,
Zhiwei Ma,
Yuan Wang,
Moo K. Chung
Abstract:
This paper surveys various distance measures for networks and graphs that were introduced in persistent homology. The scope of the paper is limited to network distances that were actually used in brain networks but the methods can be easily adapted to any weighted graph in other fields. The network version of Gromov-Hausdorff, bottleneck, kernel distances are introduced. We also introduce a recent…
▽ More
This paper surveys various distance measures for networks and graphs that were introduced in persistent homology. The scope of the paper is limited to network distances that were actually used in brain networks but the methods can be easily adapted to any weighted graph in other fields. The network version of Gromov-Hausdorff, bottleneck, kernel distances are introduced. We also introduce a recently developed KS-test like distance based on monotonic topology features such as the zeroth Betti number. Numerous toy examples and the result of applying many different distances to the brain networks of different clinical status and populations are given.
△ Less
Submitted 11 July, 2017; v1 submitted 16 January, 2017;
originally announced January 2017.
-
Compact Full Duplex MIMO Radios in D2D Underlaid Cellular Networks: From System Design to Prototype Results
Authors:
MinKeun Chung,
Min Soo Sim,
Dong Ku Kim,
Chan-Byoung Chae
Abstract:
This paper considers the implementation and application possibilities of a compact full duplex multiple-input multiple-output (MIMO) architecture where direct communication exists between users, e.g., device-to-device (D2D) and cellular link coexisting on the same spectrum. For the architecture of the compact full duplex radio, we combine an analog self-interference canceler based dual-polarizatio…
▽ More
This paper considers the implementation and application possibilities of a compact full duplex multiple-input multiple-output (MIMO) architecture where direct communication exists between users, e.g., device-to-device (D2D) and cellular link coexisting on the same spectrum. For the architecture of the compact full duplex radio, we combine an analog self-interference canceler based dual-polarization with high cross-polarization discrimination (XPD) and Long Term Evolution (LTE)-based per-subcarrier digital self-interference canceler. While we consider the compactness and power efficiency of an analog solution, we focus on the digital canceler design with robustness to a frequency-selective channel and high compatibility with a conventional LTE system. For an over-the-air wireless experiment of full duplex testbed with a two-user-pair, we implement a full duplex MIMO physical layer (PHY), supporting 20 MHz bandwidth, on an FPGA-based software-defined radio platform. Further, we propose a novel timing synchronization method to construct a more viable full duplex MIMO link. By having the full duplex link prototype fully operating in real-time, we present the first characterization of the proposed compact full duplex MIMO performance depending on the transmit power of the full duplex node. We also show the link quality between nodes. One of the crucial insights of this work is that the full duplex operation of a user is capable of acquiring the throughput gain if the user has self-interference capability with guaranteed performance.
△ Less
Submitted 17 March, 2017; v1 submitted 19 December, 2016;
originally announced December 2016.
-
Half-Chain Entanglement Entropy in the One-Dimensional Spinless Fermion Model
Authors:
Myung-Hoon Chung
Abstract:
We calculate the half-chain entanglement entropy of the ground state in the one-dimensional spinless fermion model. Considering a tiny corner of the Hilbert space represented by matrix product states, we efficiently find the ground state by the infinite time-evolving block decimation. The Schmidt coefficients are used to determine the half-chain entanglement entropy. Using the bond dimension scali…
▽ More
We calculate the half-chain entanglement entropy of the ground state in the one-dimensional spinless fermion model. Considering a tiny corner of the Hilbert space represented by matrix product states, we efficiently find the ground state by the infinite time-evolving block decimation. The Schmidt coefficients are used to determine the half-chain entanglement entropy. Using the bond dimension scaling of the half-chain entanglement entropy, we find the critical region, which is consistent with the previous results.
△ Less
Submitted 17 October, 2016;
originally announced October 2016.
-
The Broadband and Spectrally-Resolved H-band Eclipse of KELT-1b and the Role of Surface Gravity in Stratospheric Inversions in Hot Jupiters
Authors:
Thomas G. Beatty,
Nikku Madhusudhan,
Richard Pogge,
Sun Mi Chung,
Allyson Bierlya,
B. Scott Gaudi,
David W. Latham
Abstract:
We present a high precision H-band emission spectrum of the transiting brown dwarf KELT-1b, which we spectrophotometrically observed during a single secondary eclipse using the LUCI1 multi-object spectrograph on the Large Binocular Telescope. Using a Gaussian-process regression model, we are able to clearly measure the broadband eclipse depth as Delta-H=1418+/-94ppm. We are also able to spectrally…
▽ More
We present a high precision H-band emission spectrum of the transiting brown dwarf KELT-1b, which we spectrophotometrically observed during a single secondary eclipse using the LUCI1 multi-object spectrograph on the Large Binocular Telescope. Using a Gaussian-process regression model, we are able to clearly measure the broadband eclipse depth as Delta-H=1418+/-94ppm. We are also able to spectrally-resolve the H-band into five separate wavechannels and measure the eclipse spectrum of KELT-1b at R~50 with an average precision of +/-135ppm. We find that the day side has an average brightness temperature of 3250+/-50K, with significant variation as a function of wavelength. Based on our observations, and previous measurements of KELT-1b's eclipse at other wavelengths, we find that KELT-1b's day side appears identical to an isolated 3200K brown dwarf, and our modeling of the atmospheric emission shows a monotonically decreasing temperature-pressure profile. This is in contrast to hot Jupiters with similar day side brightness temperatures near 3000K, all of which appear to be either isothermal or posses a stratospheric temperature inversion. We hypothesize that the lack of an inversion in KELT-1b is due to its high surface gravity, which we argue could be caused by the increased efficiency of cold-trap processes within its atmosphere.
△ Less
Submitted 23 August, 2017; v1 submitted 11 October, 2016;
originally announced October 2016.
-
Relaxation of the entanglement spectrum in quench dynamics of topological systems
Authors:
Yi-Hao Jhu,
Pochung Chen,
Ming-Chiang Chung
Abstract:
We study how the entanglement spectrum relaxes to its steady state in one-dimensional quadratic systems after a quantum quench. In particular we apply the saddle point expansion to the dimerized chains and 1-D p-wave superconductors. We find that the entanglement spectrum always exhibits a power-law relaxation superimposed with oscillations at certain characteristic angular frequencies. For the di…
▽ More
We study how the entanglement spectrum relaxes to its steady state in one-dimensional quadratic systems after a quantum quench. In particular we apply the saddle point expansion to the dimerized chains and 1-D p-wave superconductors. We find that the entanglement spectrum always exhibits a power-law relaxation superimposed with oscillations at certain characteristic angular frequencies. For the dimerized chains, we find that the exponent $ν$ of the power-law decay is always $3/2$. For 1-D p-wave superconductors, however, we find that depending on the initial and final Hamiltonian, the exponent $ν$ can take value from a limited list of values. The smallest possible value is $ν=1/2$, which leads to a very slow convergence to its steady state value.
△ Less
Submitted 8 March, 2017; v1 submitted 4 October, 2016;
originally announced October 2016.
-
Large deviation principle in one-dimensional dynamics
Authors:
Yong Moo Chung,
Juan Rivera-Letelier,
Hiroki Takahasi
Abstract:
We study the dynamics of smooth interval maps with non-flat critical points. For every such a map that is topologically exact, we establish the full (level-2) Large Deviation Principle for empirical means. In particular, the Large Deviation Principle holds for every non\nobreakdash-renormalizable quadratic map. This includes the maps without physical measure found by Hofbauer and Keller, and chall…
▽ More
We study the dynamics of smooth interval maps with non-flat critical points. For every such a map that is topologically exact, we establish the full (level-2) Large Deviation Principle for empirical means. In particular, the Large Deviation Principle holds for every non\nobreakdash-renormalizable quadratic map. This includes the maps without physical measure found by Hofbauer and Keller, and challenges the widely-shared view of the Large Deviation Principle as a refinement of laws of large numbers.
△ Less
Submitted 17 July, 2019; v1 submitted 3 October, 2016;
originally announced October 2016.
-
Generalized Kapchinskij-Vladimirskij Distribution and Beam Matrix for Phase-Space Manipulations of High-Intensity Beams
Authors:
Moses Chung,
Hong Qin,
Ronald C. Davidson,
Lars Groening,
Chen Xiao
Abstract:
In an uncoupled linear lattice system, the Kapchinskij-Vladimirskij (KV) distribution, formulated on the basis of the single-particle Courant-Snyder (CS) invariants, has served as a fundamental theoretical basis for the analyses of the equilibrium, stability, and transport properties of high-intensity beams for the past several decades. Recent applications of high-intensity beams, however, require…
▽ More
In an uncoupled linear lattice system, the Kapchinskij-Vladimirskij (KV) distribution, formulated on the basis of the single-particle Courant-Snyder (CS) invariants, has served as a fundamental theoretical basis for the analyses of the equilibrium, stability, and transport properties of high-intensity beams for the past several decades. Recent applications of high-intensity beams, however, require beam phase-space manipulations by intentionally introducing strong coupling. In this Letter, we report the full generalization of the KV model by including all of the linear (both external and space-charge) coupling forces, beam energy variations, and arbitrary emittance partition, which all form essential elements for phase-space manipulations. The new generalized KV model yields spatially uniform density profiles and corresponding linear self-field forces as desired. The corresponding matrix envelope equations and beam matrix for the generalized KV model provide important new theoretical tools for the detailed design and analysis of high-intensity beam manipulations, for which previous theoretical models are not easily applicable.
△ Less
Submitted 1 October, 2016;
originally announced October 2016.
-
The Importance of Being Negative: A serious treatment of non-trivial edges in brain functional connectome
Authors:
Liang Zhan,
Lisanne M. Jenkins,
Ouri E. Wolfson,
Johnson J. GadElkarim,
Kevin Nocito,
Paul M. Thompson,
Olusola A. Ajilore,
Moo K. Chung,
Alex D. Leow
Abstract:
Understanding the modularity of fMRI-derived brain networks or connectomes can inform the study of brain function organization. However, fMRI connectomes additionally involve negative edges, which are not rigorously accounted for by existing approaches to modularity that either ignores or arbitrarily weight these connections. Furthermore, most Q maximization-based modularity algorithms yield varia…
▽ More
Understanding the modularity of fMRI-derived brain networks or connectomes can inform the study of brain function organization. However, fMRI connectomes additionally involve negative edges, which are not rigorously accounted for by existing approaches to modularity that either ignores or arbitrarily weight these connections. Furthermore, most Q maximization-based modularity algorithms yield variable results with suboptimal reproducibility. Here we present an alternative, reproducible approach that exploits how frequent the BOLD-signal correlation between two nodes is negative. We validated this novel probability-based modularity approach on two independent publicly-available resting-state connectome dataset (the Human Connectome Project and the 1000 Functional Connectomes) and demonstrated that negative correlations alone are sufficient in understanding resting-state modularity. In fact, this approach a) permits a dual formulation, leading to equivalent solutions regardless of whether one considers positive or negative edges; b) is theoretically linked to the Ising model defined on the connectome, thus yielding modularity result that maximizes data likelihood. We additionally were able to detect sex differences in modularity that the most widely utilized methods did not. Results confirmed the superiority of our approach in that: a) correlations with the highest probability of being negative are consistently placed between modules, b) due to the equivalent dual forms, no arbitrary weighting factor is required to balance the influence between negative and positive correlations, unlike existing Q maximization-based modularity approaches. As datasets like HCP become widely available for analysis by the neuroscience community at large, appropriate computational tools to understand the neurobiological information of negative edges in fMRI connectomes are increasingly important.
△ Less
Submitted 5 June, 2017; v1 submitted 6 September, 2016;
originally announced September 2016.
-
Nonlinear Self-Interference Cancellation for Full-Duplex Radios: From Link- and System-Level Performance Perspectives
Authors:
Min Soo Sim,
MinKeun Chung,
Dongkyu Kim,
Jaehoon Chung,
Dong Ku Kim,
Chan-Byoung Chae
Abstract:
One of the promising technologies for LTE Evolution is full-duplex radio, an innovation is expected to double the spectral efficiency. To realize full-duplex in practice, the main challenge is overcoming self-interference, and to do so, researchers have developed self-interference cancellation techniques. Since most wireless transceivers use power amplifiers, especially in cellular systems, resear…
▽ More
One of the promising technologies for LTE Evolution is full-duplex radio, an innovation is expected to double the spectral efficiency. To realize full-duplex in practice, the main challenge is overcoming self-interference, and to do so, researchers have developed self-interference cancellation techniques. Since most wireless transceivers use power amplifiers, especially in cellular systems, researchers have revealed the importance of nonlinear self-interference cancellation. In this article, we first explore several nonlinear digital self-interference cancellation techniques. We then propose a low complexity pre-calibration-based nonlinear digital self-interference cancellation technique. Next we discuss issues about reference signal allocation and the overhead of each technique. For performance evaluations, we carry out extensive measurements through a real-time prototype and link-/system-level simulations. For link-level analysis, we measure the amount of cancelled self-interference for each technique. We also evaluate system-level performances through 3D ray-tracing-based simulations. Numerical results confirm the significant performance improvement over a half-duplex system even in interference-limited indoor environments.
△ Less
Submitted 14 February, 2017; v1 submitted 7 July, 2016;
originally announced July 2016.