Showing 1–50 of 65 results for author: Xia, F

Optical wavefront shaping in deep tissue using photoacoustic feedback

Authors: Fei Xia, Ivo Leite, Robert Prevedel, Thomas Chaigne

Abstract: Over the past decade, optical wavefront shaping has been developed to focus light through highly opaque scattering layers, opening new possibilities for biomedical applications. To probe light intensity deep inside soft scattering media such as biological tissues, internal guide-stars are required. Here, we give an overview of the main principles and describe in depth the use of a photoacoustic fe… ▽ More Over the past decade, optical wavefront shaping has been developed to focus light through highly opaque scattering layers, opening new possibilities for biomedical applications. To probe light intensity deep inside soft scattering media such as biological tissues, internal guide-stars are required. Here, we give an overview of the main principles and describe in depth the use of a photoacoustic feedback signal for this purpose. We further present first principles calculations and simulations to estimate important experimental parameters, and detailed instructions on designing and conducting these experiments. Finally, we provide guidance towards selecting suitable equipment for building a typical experimental setup, paving the way for further innovative biomedical imaging and therapy applications. △ Less

Submitted 5 July, 2024; originally announced July 2024.

Comments: 25 pages, 9 figures, 3 tables

PANDORA: Deep graph learning based COVID-19 infection risk level forecasting

Authors: Shuo Yu, Feng Xia, Yueru Wang, Shihao Li, Falih Febrinanto, Madhu Chetty

Abstract: COVID-19 as a global pandemic causes a massive disruption to social stability that threatens human life and the economy. Policymakers and all elements of society must deliver measurable actions based on the pandemic's severity to minimize the detrimental impact of COVID-19. A proper forecasting system is arguably important to provide an early signal of the risk of COVID-19 infection so that the au… ▽ More COVID-19 as a global pandemic causes a massive disruption to social stability that threatens human life and the economy. Policymakers and all elements of society must deliver measurable actions based on the pandemic's severity to minimize the detrimental impact of COVID-19. A proper forecasting system is arguably important to provide an early signal of the risk of COVID-19 infection so that the authorities are ready to protect the people from the worst. However, making a good forecasting model for infection risks in different cities or regions is not an easy task, because it has a lot of influential factors that are difficult to be identified manually. To address the current limitations, we propose a deep graph learning model, called PANDORA, to predict the infection risks of COVID-19, by considering all essential factors and integrating them into a geographical network. The framework uses geographical position relations and transportation frequency as higher-order structural properties formulated by higher-order network structures (i.e., network motifs). Moreover, four significant node attributes (i.e., multiple features of a particular area, including climate, medical condition, economy, and human mobility) are also considered. We propose three different aggregators to better aggregate node attributes and structural features, namely, Hadamard, Summation, and Connection. Experimental results over real data show that PANDORA outperforms the baseline method with higher accuracy and faster convergence speed, no matter which aggregator is chosen. We believe that PANDORA using deep graph learning provides a promising approach to get superior performance in infection risk level forecasting and help humans battle the COVID-19 crisis. △ Less

Submitted 7 June, 2024; originally announced June 2024.

Identifying L-H transition in HL-2A through deep learning

Authors: Meihuizi He, Songfen Liu, Fan Xia, Zongyu Yang, Wulyu Zhong

Abstract: During the operation of tokamak devices, addressing the thermal load issues caused by Edge Localized Modes (ELMs) eruption is crucial. Ideally, mitigation and suppression measures for ELMs should be promptly initiated as soon as the first low-to-high confinement (L-H) transition occurs, which necessitates the real-time monitoring and accurate identification of the L-H transition process. Motivated… ▽ More During the operation of tokamak devices, addressing the thermal load issues caused by Edge Localized Modes (ELMs) eruption is crucial. Ideally, mitigation and suppression measures for ELMs should be promptly initiated as soon as the first low-to-high confinement (L-H) transition occurs, which necessitates the real-time monitoring and accurate identification of the L-H transition process. Motivated by this, and by recent deep learning boom, we propose a deep learning-based L-H transition identification algorithm on HL-2A tokamak. In this work, we have constructed a neural network comprising layers of Residual Long Short-Term Memory (LSTM) and Temporal Convolutional Network (TCN). Unlike previous work based on recognition for ELMs by slice, this method implements recognition on L-H transition process before the first ELMs crash. Therefore the mitigation techniques can be triggered in time to suppress the initial ELMs bursts. In order to further explain the effectiveness of the algorithm, we developed a series of evaluation indicators by shots, and the results show that this algorithm can provide necessary reference for the mitigation and suppression system. △ Less

Submitted 17 May, 2024; originally announced May 2024.

Coherent XUV super continuum emission from atomic bound states

Authors: Jing Zhao, Xiaowei Wang, Li Wang, Jiacan Wang, Yalei Zhu, Fan Xiao, Wenkai Tao, Zhigang Zheng, Haizhong Wu, Xu Sun, Yue Lang, Congsen Meng, Dongwen Zhang, Zhihui Lv, Jinlei Liu, Zengxiu Zhao

Abstract: Coherent supercontinuum radiation in the extreme-ultraviolet (XUV) range is indispensable for synthesizing attosecond light pulses and for exploring transient atomic structures. Here, we report the striking observations of coherent XUV supercontinuum (XSC) extended from below to far above the ionization threshold, which exhibits completely different temporal and spatial properties comparing to the… ▽ More Coherent supercontinuum radiation in the extreme-ultraviolet (XUV) range is indispensable for synthesizing attosecond light pulses and for exploring transient atomic structures. Here, we report the striking observations of coherent XUV supercontinuum (XSC) extended from below to far above the ionization threshold, which exhibits completely different temporal and spatial properties comparing to the conventional rescattering induced high harmonic generation (HHG). We demonstrate that the strong-field created coherence among bound orbitals strongly distort the atomic transition energies during the pulse, leading to coherent emission spanning tens of electron-volts, in contrast to the line emission via free-induction decay occurring after the pulse. The supposed non-radiating bound dark states contribute as well by emitting dressed energy through dark-to-bright emission mechanism. All the processes modulated at sub-cycle time scale jointly form this new-type coherent XSC. This work achieves the strong-field attosecond control of the exotic atomic radiation dynamics and provides the means of simultaneous generation of separated attosecond sources, i.e., XSC and HHG, with potential advancing attosecond interferometry. △ Less

Submitted 3 May, 2024; originally announced May 2024.

Replica-assisted super-resolution fluorescence imaging in scattering media

Authors: Tengfei Wu, YoonSeok Baek, Fei Xia, Sylvain Gigan, Hilton B. de Aguiar

Abstract: Far-field super-resolution fluorescence microscopy has been rapidly developed for applications ranging from cell biology to nanomaterials. However, it remains a significant challenge to achieve super-resolution imaging at depth in opaque materials. In this study, we present a super-resolution microscopy technique for imaging hidden fluorescent objects through scattering media, started by exploitin… ▽ More Far-field super-resolution fluorescence microscopy has been rapidly developed for applications ranging from cell biology to nanomaterials. However, it remains a significant challenge to achieve super-resolution imaging at depth in opaque materials. In this study, we present a super-resolution microscopy technique for imaging hidden fluorescent objects through scattering media, started by exploiting the inherent object replica generation arising from the memory effect, i.e. the seemingly informationless emission speckle can be regarded as a random superposition of multiple object copies. Inspired by the concept of super-resolution optical fluctuation imaging, we use temporally-fluctuating speckles to excite fluorescent signals and perform high-order cumulant analysis on the fluctuation, which can not only improve the image resolution, but also increase the speckle contrast to isolate only the bright object replicas. A super-resolved image can be finally retrieved by simply unmixing the sparsely distributed replicas with their location map. This methodology allows to overcome scattering and achieve robust super-resolution fluorescence imaging, circumventing the need of heavy computational steps. △ Less

Submitted 30 April, 2024; originally announced April 2024.

Comments: 6 pages, 3 figures in the main text

Zincophilic armor: Phytate ammonium as a multifunctional additive for enhanced performance in aqueous zinc-ion batteries

Authors: Fangyuan Xiao, Xiaoke Wang, Kaitong Sun, Qian Zhao, Cuiping Han, Hai-Feng Li

Abstract: Corrosion and the formation of by-products resulting from parasitic side reactions, as well as random dendrite growth, pose significant challenges for aqueous zinc-ion batteries (AZIBs). In this study, phytate ammonium is introduced into the traditional dilute Zinc sulfate electrolyte as a multi-functional additive. Leveraging the inherent zincophilic nature of the phytic anion, a protective layer… ▽ More Corrosion and the formation of by-products resulting from parasitic side reactions, as well as random dendrite growth, pose significant challenges for aqueous zinc-ion batteries (AZIBs). In this study, phytate ammonium is introduced into the traditional dilute Zinc sulfate electrolyte as a multi-functional additive. Leveraging the inherent zincophilic nature of the phytic anion, a protective layer is formed on the surface of the zinc anode. This layer can effectively manipulate the deposition process, mitigate parasitic reactions, and reduce the accumulation of detrimental by-products. Additionally, the competitive deposition between dissociated ammonium ions and Zn2+ promotes uniform deposition, thereby alleviating dendrite growth. Consequently, the modified electrolyte with a lower volume addition exhibits superior performance. The zinc symmetric battery demonstrates much more reversible plating/stripping, sustaining over 2000 hours at 5 mA cm-2 and 1 mA h cm-2. A high average deposition/stripping efficiency of 99.83% is achieved, indicating the significant boosting effect and practical potential of our strategy for high-performance aqueous zinc-ion batteries. △ Less

Submitted 7 April, 2024; originally announced April 2024.

Neurophotonics beyond the Surface: Unmasking the Brain's Complexity Exploiting Optical Scattering

Authors: Fei Xia, Caio Vaz Rimoli, Walther Akemann, Cathie Ventalon, Laurent Bourdieu, Sylvain Gigan, Hilton B de Aguiar

Abstract: The intricate nature of the brain necessitates the application of advanced probing techniques to comprehensively study and understand its working mechanisms. Neurophotonics offers minimally invasive methods to probe the brain using optics at cellular and even molecular levels. However, multiple challenges persist, especially concerning imaging depth, field of view, speed, and biocompatibility. A m… ▽ More The intricate nature of the brain necessitates the application of advanced probing techniques to comprehensively study and understand its working mechanisms. Neurophotonics offers minimally invasive methods to probe the brain using optics at cellular and even molecular levels. However, multiple challenges persist, especially concerning imaging depth, field of view, speed, and biocompatibility. A major hindrance to solving these challenges in optics is the scattering nature of the brain. This perspective highlights the potential of complex media optics, a specialized area of study focused on light propagation in materials with intricate heterogeneous optical properties, in advancing and improving neuronal readouts for structural imaging and optical recordings of neuronal activity. Key strategies include wavefront shaping techniques and computational imaging and sensing techniques that exploit scattering properties for enhanced performance. We discuss the potential merger of the two fields as well as potential challenges and perspectives toward longer term in vivo applications. △ Less

Submitted 21 March, 2024; originally announced March 2024.

Comments: 22 pages, 3 figures, 1 table

A Modelling Investigation for Solar Flare X-ray Stereoscopy with Solar Orbiter/STIX and Earth Orbiting Missions

Authors: Natasha L. S. Jeffrey, Säm Krucker, Morgan Stores, Eduard P. Kontar, Pascal Saint-Hilaire, Andrea F. Battaglia, Laura Hayes, Hannah Collier, Astrid Veronig, Yang Su, Srikar Paavan Tadepalli, Fanxiaoyu Xia

Abstract: The Spectrometer/Telescope for Imaging X-rays (STIX) on board Solar Orbiter (SolO) provides a unique opportunity to systematically perform stereoscopic X-ray observations of solar flares with current and upcoming X-ray missions at Earth. These observations will produce the first reliable measurements of hard X-ray (HXR) directivity in decades, providing a new diagnostic of the flare-accelerated el… ▽ More The Spectrometer/Telescope for Imaging X-rays (STIX) on board Solar Orbiter (SolO) provides a unique opportunity to systematically perform stereoscopic X-ray observations of solar flares with current and upcoming X-ray missions at Earth. These observations will produce the first reliable measurements of hard X-ray (HXR) directivity in decades, providing a new diagnostic of the flare-accelerated electron angular distribution and helping to constrain the processes that accelerate electrons in flares. However, such observations must be compared to modelling, taking into account electron and X-ray transport effects and realistic plasma conditions, all of which can change the properties of the measured HXR directivity. Here, we show how HXR directivity, defined as the ratio of X-ray spectra at different spacecraft viewing angles, varies with different electron and flare properties (e.g., electron angular distribution, highest energy electrons, and magnetic configuration), and how modelling can be used to extract these typically unknown properties from the data. Lastly, we present a preliminary HXR directivity analysis of two flares, observed by the Fermi Gamma-ray Burst Monitor (GBM) and SolO/STIX, demonstrating the feasibility and challenges associated with such observations, and how HXR directivity can be extracted by comparison with the modelling presented here. △ Less

Submitted 29 January, 2024; originally announced January 2024.

Comments: Accepted for publication in ApJ (January 2024)

Nonlinear Self-Calibrated Spectrometer with Single GeSe-InSe Heterojunction Device

Authors: Rana Darweesh, Rajesh Kumar Yadav, Elior Adler, Michal Poplinger, Adi Levi, Jea-Jung Lee, Amir Leshem, Ashwin Ramasubramaniam, Fengnian Xia, Doron Naveh

Abstract: Optical spectroscopy the measurement of electromagnetic spectra is fundamental to various scientific domains and serves as the building block of numerous technologies. Computational spectrometry is an emerging field that employs an array of photodetectors with different spectral responses or a single photodetector device with tunable spectral response, in conjunction with numerical algorithms, for… ▽ More Optical spectroscopy the measurement of electromagnetic spectra is fundamental to various scientific domains and serves as the building block of numerous technologies. Computational spectrometry is an emerging field that employs an array of photodetectors with different spectral responses or a single photodetector device with tunable spectral response, in conjunction with numerical algorithms, for spectroscopic measurements. Compact single photodetectors made from layered materials are particularly attractive, since they eliminate the need for bulky mechanical and optical components used in traditional spectrometers and can easily be engineered as heterostructures to optimize device performance. However, compact tunable photodetectors are typically nonlinear devices and this adds complexity to extracting optical spectra from the device response. Here, we report on the training of an artificial neural network (ANN) to recover the full nonlinear spectral photoresponse of a nonlinear problem of high dimensionality of a single GeSe-InSe p-n heterojunction device. We demonstrate the functionality of a calibrated spectrometer in the spectral range of 400-1100 nm, with a small device footprint of ~25X25 micrometers, and we achieve a mean reconstruction error of 0.0002 for the power-spectrum at a spectral resolution of 0.35 nm. Using our device, we demonstrate a solution to metamerism, an apparent matching of colors with different power spectral distributions, which is a fundamental problem in optical imaging. △ Less

Submitted 21 November, 2023; originally announced November 2023.

Comments: 15 pages, 5 figures

MSC Class: Nonlinear Self-Calibrated Spectrometer; Artificial neural network; Optical imaging

Large-scale photonic computing with nonlinear disordered media

Authors: Hao Wang, Jianqi Hu, Andrea Morandi, Alfonso Nardi, Fei Xia, Xuanchen Li, Romolo Savo, Qiang Liu, Rachel Grange, Sylvain Gigan

Abstract: Neural networks find widespread use in scientific and technological applications, yet their implementations in conventional computers have encountered bottlenecks due to ever-expanding computational needs. Photonic neuromorphic hardware, which manipulates information and represents data continuously in the optical domain, is one of the promising platforms with potential advantages of massive paral… ▽ More Neural networks find widespread use in scientific and technological applications, yet their implementations in conventional computers have encountered bottlenecks due to ever-expanding computational needs. Photonic neuromorphic hardware, which manipulates information and represents data continuously in the optical domain, is one of the promising platforms with potential advantages of massive parallelism, ultralow latency, and reduced energy consumption. While linear photonic neural networks are within reach, photonic computing with large-scale optical nonlinear nodes remains largely unexplored. Here, we demonstrate a large-scale, high-performance nonlinear photonic neural system based on a disordered polycrystalline slab composed of lithium niobate nanocrystals. Mediated by random quasi-phase-matching and multiple scattering, linear and nonlinear optical speckle features are generated as the interplay between the simultaneous linear random scattering and the second-harmonic generation, defining a complex neural network in which the second-order nonlinearity acts as internal nonlinear activation functions. Benchmarked against linear random projection, such nonlinear mapping embedded with rich physical computational operations shows improved performance across a large collection of machine learning tasks in image classification, regression, and graph classification with varying complexity. Demonstrating up to 27,648 input and 3,500 nonlinear output nodes, the combination of optical nonlinearity and random scattering serves as a scalable computing engine for diverse applications. △ Less

Submitted 11 October, 2023; originally announced October 2023.

Deep Learning with Passive Optical Nonlinear Mapping

Authors: Fei Xia, Kyungduk Kim, Yaniv Eliezer, Liam Shaughnessy, Sylvain Gigan, Hui Cao

Abstract: Deep learning has fundamentally transformed artificial intelligence, but the ever-increasing complexity in deep learning models calls for specialized hardware accelerators. Optical accelerators can potentially offer enhanced performance, scalability, and energy efficiency. However, achieving nonlinear mapping, a critical component of neural networks, remains challenging optically. Here, we introdu… ▽ More Deep learning has fundamentally transformed artificial intelligence, but the ever-increasing complexity in deep learning models calls for specialized hardware accelerators. Optical accelerators can potentially offer enhanced performance, scalability, and energy efficiency. However, achieving nonlinear mapping, a critical component of neural networks, remains challenging optically. Here, we introduce a design that leverages multiple scattering in a reverberating cavity to passively induce optical nonlinear random mapping, without the need for additional laser power. A key advantage emerging from our work is that we show we can perform optical data compression, facilitated by multiple scattering in the cavity, to efficiently compress and retain vital information while also decreasing data dimensionality. This allows rapid optical information processing and generation of low dimensional mixtures of highly nonlinear features. These are particularly useful for applications demanding high-speed analysis and responses such as in edge computing devices. Utilizing rapid optical information processing capabilities, our optical platforms could potentially offer more efficient and real-time processing solutions for a broad range of applications. We demonstrate the efficacy of our design in improving computational performance across tasks, including classification, image reconstruction, key-point detection, and object detection, all achieved through optical data compression combined with a digital decoder. Notably, we observed high performance, at an extreme compression ratio, for real-time pedestrian detection. Our findings pave the way for novel algorithms and architectural designs for optical computing. △ Less

Submitted 18 July, 2023; v1 submitted 17 July, 2023; originally announced July 2023.

Comments: 16 pages, 7 figures

Plasma heating and nanoflare caused by slow-mode wave in a coronal loop

Authors: Fanxiaoyu Xia, Tongjiang Wang, Yang Su, Jie Zhao, Qingmin Zhang, Astrid M. Veronig, Weiqun Gan

Abstract: We present a detailed analysis of a reflecting intensity perturbation in a large coronal loop that appeared as sloshing oscillation and lasted for at least one and a half periods. The perturbation is initiated by a microflare at one footpoint of the loop, propagates along the loop and is eventually reflected at the remote footpoint where significant brightenings are observed in all the AIA extreme… ▽ More We present a detailed analysis of a reflecting intensity perturbation in a large coronal loop that appeared as sloshing oscillation and lasted for at least one and a half periods. The perturbation is initiated by a microflare at one footpoint of the loop, propagates along the loop and is eventually reflected at the remote footpoint where significant brightenings are observed in all the AIA extreme-ultraviolet (EUV) channels. This unique observation provides us with the opportunity to better understand not only the thermal properties and damping mechanisms of the sloshing oscillation, but also the energy transfer at the remote footpoint. Based on differential emission measures (DEM) analysis and the technique of coronal seismology, we find that 1) the calculated local sound speed is consistent with the observed propagation speed of the perturbation during the oscillation, which is suggestive of a slow magnetoacoustic wave; 2) thermal conduction is the major damping mechanism of the wave but additional damping mechanism such as anomalous enhancement of compressive viscosity or wave leakage is also required to account for the rapid decay of the observed waves; 3) the wave produced a nanoflare at the remote footpoint, with a peak thermal energy of $\thicksim10^{24}-10^{25}$ erg. This work provides a consistent picture of the magnetoacoustic wave propagation and reflection in a coronal loop, and reports the first solid evidence of a wave-induced nanoflare. The results reveal new clues for further simulation studies and may help solving the coronal heating problem. △ Less

Submitted 21 August, 2022; originally announced August 2022.

Comments: 13pages, 5 figures; Accepted by ApJL

Sub-monolayer Biolasers: Lower Gain, Higher Sensitivity

Authors: C. Gong, X. Yang, S. J. Tang, Q. Q. Zhang, Y. Wang, Y. L. Liu, Y. C. Chen, G. D. Peng, X. Fan, Y. F. Xiao, Y. J. Rao, Y. Gong

Abstract: Biomarker detection is the key to identifying health risks. However, designing sensitive biosensors in a single-use mode for disease diagnosis remains a major challenge. Here, we report sub-monolayer biolasers with remarkable repeatability for ultrasensitive and disposable biomarker detection. The biolaser sensors are designed by employing the telecom optical fibers as distributed optical microcav… ▽ More Biomarker detection is the key to identifying health risks. However, designing sensitive biosensors in a single-use mode for disease diagnosis remains a major challenge. Here, we report sub-monolayer biolasers with remarkable repeatability for ultrasensitive and disposable biomarker detection. The biolaser sensors are designed by employing the telecom optical fibers as distributed optical microcavities and pushing the gain molecules down to the sub-monolayer level. We observe a status transition from the monolayer biolaser to the sub-monolayer biolaser by tuning the specific conjugation. By reducing the fluorophores down to the threshold density (~ 3.2 x 10-13 mol/cm2), we demonstrate an ultimate sensitivity of sub-monolayer biolaser with six orders of magnitude enhancement compared with the monolayer biolasers. We further achieved ultrasensitive immunoassay for Parkinson's disease biomarker, alpha-synuclein, with a lower limit of detection of 0.32 pM in serum. This biosensor with massive fabrication capability at ultralow cost provides a general method for the ultrasensitive disposable biodetection of disease biomarkers. △ Less

Submitted 5 July, 2022; originally announced July 2022.

Comments: 27 pages, 15 figures

MSC Class: 78A70

A clinically relevant online patient QA solution with daily CT scans and EPID-based in vivo dosimetry: A feasible study on rectal cancer

Authors: Liyuan Chen, Zhiyuan Zhang, Lei Yu, Jiyou Peng, Bin Feng, Jun Zhao, Yanfang Liu, Fan Xia, Zhen Zhang, Weigang Hu, Jiazhou Wang

Abstract: Adaptive radiation therapy (ART) could protect organs at risk (OARs) while maintain high dose coverage to targets. However, there still lack efficient online patient QA methods. We aim to develop a clinically relevant online patient quality assurance (QA) solution for ART using daily CT scans and electronic portal imaging device (EPID)-based in vivo dosimetry. Ten patients with rectal cancer at ou… ▽ More Adaptive radiation therapy (ART) could protect organs at risk (OARs) while maintain high dose coverage to targets. However, there still lack efficient online patient QA methods. We aim to develop a clinically relevant online patient quality assurance (QA) solution for ART using daily CT scans and electronic portal imaging device (EPID)-based in vivo dosimetry. Ten patients with rectal cancer at our center were included. Patients' daily CT scans and portal images were collected to generate reconstructed 3D dose distributions. Contours of targets and OARs were recontoured on these daily CT scans by a clinician or an auto-segmentation algorithm, then dose-volume indices were calculated, and the percent deviation of these indices to their original plans were determined. This deviation was regarded as the metric for clinically relevant patient QA. The tolerance level was obtained using a 95% interval of the QA metric distribution. These deviations could be further divided into anatomically relevant or delivery relevant indicators for error source analysis. Finally, our QA solution was validated on an additional six clinical patients. In rectal cancer, the lower and upper tolerance of the QA metric for PTV ΔD95 (%) were [-3.11%, 2.35%], and for PTV ΔD2 (%) were [-0.78%, 3.23%]. In validation, the 68% for PTV ΔD95 (%) and the 79% for PTV ΔD2 ({%)of the 28 fractions are within tolerances of the QA metrics. By using four or more out-of-tolerance QA metrics as an action level, there were 5 fractions (18%) have four or more out-of-tolerance QA metrics in validation patient dataset. The online patient QA solution using daily CT scans and EPID-based in vivo dosimetry is clinically feasible. Source of error analysis has the potential for distinguishing sources of error and guiding ART for future treatments. △ Less

Submitted 27 June, 2022; originally announced June 2022.

Physics-based neural network for non-invasive control of coherent light in scattering media

Authors: Alexandra d'Arco, Fei Xia, Antoine Boniface, Jonathan Dong, Sylvain Gigan

Abstract: Optical imaging through complex media, such as biological tissues or fog, is challenging due to light scattering. In the multiple scattering regime, wavefront shaping provides an effective method to retrieve information; it relies on measuring how the propagation of different optical wavefronts are impacted by scattering. Based on this principle, several wavefront shaping techniques were successfu… ▽ More Optical imaging through complex media, such as biological tissues or fog, is challenging due to light scattering. In the multiple scattering regime, wavefront shaping provides an effective method to retrieve information; it relies on measuring how the propagation of different optical wavefronts are impacted by scattering. Based on this principle, several wavefront shaping techniques were successfully developed, but most of them are highly invasive and limited to proof-of-principle experiments. Here, we propose to use a neural network approach to non-invasively characterize and control light scattering inside the medium and also to retrieve information of hidden objects buried within it. Unlike most of the recently-proposed approaches, the architecture of our neural network with its layers, connected nodes and activation functions has a true physical meaning as it mimics the propagation of light in our optical system. It is trained with an experimentally-measured input/output dataset built from a series of incident light patterns and corresponding camera snapshots. We apply our physics-based neural network to a fluorescence microscope in epi-configuration and demonstrate its performance through numerical simulations and experiments. This flexible method can include physical priors and we show that it can be applied to other systems as, for example, non-linear or coherent contrast mechanisms. △ Less

Submitted 1 June, 2022; originally announced June 2022.

Comments: 15 pages, 11 figures

Denoising of Three-Dimensional Fast Spin Echo Magnetic Resonance Images of Knee Joints using Spatial-Variant Noise-Relevant Residual Learning of Convolution Neural Network

Authors: Shutian Zhao, Donal G. Cahill, Siyue Li, Fan Xiao, Thierry Blu, James F Griffith, Weitian Chen

Abstract: Two-dimensional (2D) fast spin echo (FSE) techniques play a central role in the clinical magnetic resonance imaging (MRI) of knee joints. Moreover, three-dimensional (3D) FSE provides high-isotropic-resolution magnetic resonance (MR) images of knee joints, but it has a reduced signal-to-noise ratio compared to 2D FSE. Deep-learning denoising methods are a promising approach for denoising MR images… ▽ More Two-dimensional (2D) fast spin echo (FSE) techniques play a central role in the clinical magnetic resonance imaging (MRI) of knee joints. Moreover, three-dimensional (3D) FSE provides high-isotropic-resolution magnetic resonance (MR) images of knee joints, but it has a reduced signal-to-noise ratio compared to 2D FSE. Deep-learning denoising methods are a promising approach for denoising MR images, but they are often trained using synthetic noise due to challenges in obtaining true noise distributions for MR images. In this study, inherent true noise information from 2-NEX acquisition was used to develop a deep-learning model based on residual learning of convolutional neural network (CNN), and this model was used to suppress the noise in 3D FSE MR images of knee joints. The proposed CNN used two-step residual learning over parallel transporting and residual blocks and was designed to comprehensively learn real noise features from 2-NEX training data. The results of an ablation study validated the network design. The new method achieved improved denoising performance of 3D FSE knee MR images compared with current state-of-the-art methods, based on the peak signal-to-noise ratio and structural similarity index measure. The improved image quality after denoising using the new method was verified by radiological evaluation. A deep CNN using the inherent spatial-varying noise information in 2-NEX acquisitions was developed. This method showed promise for clinical MRI assessments of the knee, and has potential applications for the assessment of other anatomical structures. △ Less

Submitted 20 April, 2022; originally announced April 2022.

Comments: 6 figures, abstract accepted by Joint Annual Meeting ISMRM-ESMRMB & ISMRT 31st Annual Meeting

Journal ref: Computers in Biology and Medicine, Volume 151, Part A, 2022, 106295, ISSN 0010-4825

A nonhydrostatic atmospheric dynamical core on cubed sphere using hybrid multi-moment finite-volume/finite difference methods: formulations and preliminary tests

Authors: Chungang Chen, Xingliang Li, Feng Xiao, Xueshun Shen

Abstract: A nonhydrostatic dynamical core has been developed by using the multi-moment finite volume method that ensures the rigorous numerical conservation. To represent the spherical geometry free of polar problems, the cubed-sphere grid is adopted. A fourth-order multi-moment discretization formulation is applied to solve the governing equations cast in the local curvilinear coordinates on each patch of… ▽ More A nonhydrostatic dynamical core has been developed by using the multi-moment finite volume method that ensures the rigorous numerical conservation. To represent the spherical geometry free of polar problems, the cubed-sphere grid is adopted. A fourth-order multi-moment discretization formulation is applied to solve the governing equations cast in the local curvilinear coordinates on each patch of cubed sphere through a gnomonic projection. In vertical direction, the height-based terrain-following grid is used to deal with the topography and a conservative finite difference scheme is adopted for the spatial discretization. The dynamical core adopts the nonhydrostatic governing equations. To get around the CFL stability restriction imposed by sound wave propagation and relatively small grid spacing in the vertical direction, the dimensional-splitting time integration algorithm using the HEVI (horizontally-explicit and vertically-implicit) strategy is implemented by applying the IMEX (implicit-explicit) Runge-Kutta method. The proposed model was checked by the widely-used benchmark tests in this study. The numerical results show that the multi-moment model has superior solution quality and great practical potential as a numerical platform for development of the atmospheric general circulation models. △ Less

Submitted 27 February, 2022; originally announced February 2022.

Comments: 37 pages, 5 figures. arXiv admin note: text overlap with arXiv:2004.06290

First implementation of full-workflow automation in radiotherapy: the All-in-One solution on rectal cancer

Authors: Lei Yu, Jun Zhao, Fan Xia, Zhiyuan Zhang, Yanfang Liu, Wei Zhang, Jingjie Zhou, Jiazhou Wang, Weigang Hu, Zhen Zhang

Abstract: The aim of this work is to describe the technical characteristics of an AI-powered radiotherapy workflow that enables full-process automation (All-in-One), evaluate its performance implemented for on-couch initial treatment of rectal cancer, and provide insight into the behavior of full-workflow automation in the specialty of radiotherapy. The All-in-One workflow was developed based on a CT-integr… ▽ More The aim of this work is to describe the technical characteristics of an AI-powered radiotherapy workflow that enables full-process automation (All-in-One), evaluate its performance implemented for on-couch initial treatment of rectal cancer, and provide insight into the behavior of full-workflow automation in the specialty of radiotherapy. The All-in-One workflow was developed based on a CT-integrated linear accelerator. It incorporates routine radiotherapy procedures from simulation, autosegmentation, autoplanning, image guidance, beam delivery, and in vivo quality assurance (QA) into one scheme, with critical decision points involved, while the patient is on the treatment couch during the whole process. For the enrolled ten patients with rectal cancer, minor modifications of the autosegmented target volumes were required, and the Dice similarity coefficient and 95% Hausdorff distance before and after modifications were 0.892{\pm}0.061 and 18.2{\pm}13.0 mm, respectively. The autosegmented normal tissues and automatic plans were clinically acceptable without any modifications or reoptimization. The pretreatment IGRT corrections were within 2 mm in all directions, and the EPID-based in vivo QA showed a γ passing rate better than 97{\%} (3{\%}/3 mm/10{\%} threshold). The duration of the whole process was 23.2{\pm}3.5 minutes, depending mostly on the time required for manual modification and plan evaluation. The All-in-One workflow enables full automation of the entire radiotherapy process by seamlessly integrating multiple routine procedures. The one-stop solution shortens the time scale it takes to ready the first treatment from days to minutes, significantly improving the patient experience and the efficiency of the workflow, and shows potential to facilitate the clinical application of online adaptive replanning. △ Less

Submitted 24 February, 2022; originally announced February 2022.

Exciting Magnetic Dipole Mode of Split Ring Plasmonic Nano Resonator by Photonic Crystal Nanocavity

Authors: Yingke Ji, Binbin Wang, Liang Fang, Qiang Zhao, Fajun Xiao, Xuetao Gan

Abstract: On chip exciting electric modes in individual plasmonic nanostructures are realized widely; nevertheless, the excitation of their magnetic counterparts is seldom reported. Here, we propose a highly efficient on chip excitation approach of the magnetic dipole mode of an individual split ring resonator (SRR) by integrating it onto a photonic crystal nanocavity (PCNC). A high excitation efficiency of… ▽ More On chip exciting electric modes in individual plasmonic nanostructures are realized widely; nevertheless, the excitation of their magnetic counterparts is seldom reported. Here, we propose a highly efficient on chip excitation approach of the magnetic dipole mode of an individual split ring resonator (SRR) by integrating it onto a photonic crystal nanocavity (PCNC). A high excitation efficiency of up to 58% is realized through the resonant coupling between the modes of the SRR and PCNC. A further fine adjustment of the excited magnetic dipole mode is demonstrated by tuning the relative position and twist angle between the SRR and PCNC. Finally, a structure with a photonic crystal waveguide side coupled with the hybrid SRR PCNC is illustrated, which could excite the magnetic dipole mode with an in plane coupling geometry and potentially facilitate the future device application. Our result may open a way for developing chip integrated photonic devices employing a magnetic field component in the optical field. △ Less

Submitted 2 December, 2021; originally announced December 2021.

Hourly Warning for Strong Earthquakes

Authors: T. Chen, L. Li, X. -X. Zhang, C. Wang, X. -B. Jin, Q. -M. Ma, J. -Y. Xu, Z. -H. He, H. Li, S. -G. Xiao, X. -Z. Wang, X. -H. Shen, X. -M. Zhang, H. -B. Li, Z. -M. Zeren, J. -P. Huang, F. -Q. Huang, S. Che, Z. -M. Zou, P. Xiong, J. Liu, L. -Q. Zhang, Q. Guo, I. Roth, V. S. Makhmutov , et al. (32 additional authors not shown)

Abstract: A promising perspective is presented that humans can provide hourly warning for strong land earthquakes (EQs, Ms6). Two important atmospheric electrostatic signal features are described. A table that lists 9 strong land EQs with shock time, epicenter, magnitude, weather in the region near the epicenter, precursor beginning time, and precursor duration demonstrates that at approximately several hou… ▽ More A promising perspective is presented that humans can provide hourly warning for strong land earthquakes (EQs, Ms6). Two important atmospheric electrostatic signal features are described. A table that lists 9 strong land EQs with shock time, epicenter, magnitude, weather in the region near the epicenter, precursor beginning time, and precursor duration demonstrates that at approximately several hours to one day before a strong land EQ, the weather conditions are fair near the epicenter, and an abnormal negative atmospheric electrostatic signal is very obvious. Moreover, the mechanism is explained. A method by which someone could determine the epicenter and the magnitude of a forthcoming strong EQ is suggested. Finally, the possibility of realizing hourly warning for strong land EQs in the near future is pointed out. △ Less

Submitted 23 June, 2021; originally announced June 2021.

Symmetry-preserving enforcement of low-dissipation method based on boundary variation diminishing principle

Authors: Hiro Wakimura, Shinichi Takagi, Feng Xiao

Abstract: A class of high-order shock-capturing schemes, P$_n$T$_m$-BVD (Deng et al., J. Comp. Phys., 386:323-349, 2019; Comput. & Fluids, 200:104433, 2020.) schemes, have been devised to solve the Euler equations with substantially reduced numerical dissipation, which enable high-resolution simulations to resolve flow structures of wider range scales. In such simulations with low dissipation, errors of rou… ▽ More A class of high-order shock-capturing schemes, P$_n$T$_m$-BVD (Deng et al., J. Comp. Phys., 386:323-349, 2019; Comput. & Fluids, 200:104433, 2020.) schemes, have been devised to solve the Euler equations with substantially reduced numerical dissipation, which enable high-resolution simulations to resolve flow structures of wider range scales. In such simulations with low dissipation, errors of round-off level might grow and contaminate the numerical solutions. A typical example of such problems is the loss of symmetry in the numerical solutions for physical problems of symmetric configurations even if the schemes are mathematically in line with the symmetry rules. In this study, the mechanisms of symmetry-breaking in a finite volume framework with the P$_4$T$_2$-BVD reconstruction scheme are thoroughly examined. Particular attention has been paid to remove the possible causes due to the lack of associativity in floating-point arithmetic which is associated with round-off errors. Modifications and new techniques are proposed to completely remove the possible causes for symmetry breaking in different components of the P$_4$T$_2$-BVD finite volume solver. Benchmark tests that have symmetric solution structures are used to verify the proposed methods. The numerical results demonstrate the perfect symmetric solution structures. △ Less

Submitted 3 June, 2021; originally announced June 2021.

Comments: 32 pages, 8 figures

Turing Number: How Far Are You to A. M. Turing Award?

Authors: Feng Xia, Jiaying Liu, Jing Ren, Wei Wang, Xiangjie Kong

Abstract: The ACM A.M. Turing Award is commonly acknowledged as the highest distinction in the realm of computer science. Since 1960s, it has been awarded to computer scientists who made outstanding contributions. The significance of this award is far-reaching to the laureates as well as their research teams. However, unlike the Nobel Prize that has been extensively investigated, little research has been do… ▽ More The ACM A.M. Turing Award is commonly acknowledged as the highest distinction in the realm of computer science. Since 1960s, it has been awarded to computer scientists who made outstanding contributions. The significance of this award is far-reaching to the laureates as well as their research teams. However, unlike the Nobel Prize that has been extensively investigated, little research has been done to explore this most important award. To this end, we propose the Turing Number (TN) index to measure how far a specific scholar is to this award. Inspired by previous works on Erdos Number and Bacon Number, this index is defined as the shortest path between a given scholar to any Turing Award Laureate. Experimental results suggest that TN can reflect the closeness of collaboration between scholars and Turing Award Laureates. With the correlation analysis between TN and metrics from the bibliometric-level and network-level, we demonstrate that TN has the potential of reflecting a scholar's academic influence and reputation. △ Less

Submitted 3 May, 2021; originally announced May 2021.

Comments: 8 pages, 3 figures

MSC Class: 68P20 ACM Class: H.3.3; H.3.7

Journal ref: ACM SIGWEB Newsletter (2020)

Pandemic Drugs at Pandemic Speed: Infrastructure for Accelerating COVID-19 Drug Discovery with Hybrid Machine Learning- and Physics-based Simulations on High Performance Computers

Authors: Agastya P. Bhati, Shunzhou Wan, Dario Alfè, Austin R. Clyde, Mathis Bode, Li Tan, Mikhail Titov, Andre Merzky, Matteo Turilli, Shantenu Jha, Roger R. Highfield, Walter Rocchia, Nicola Scafuri, Sauro Succi, Dieter Kranzlmüller, Gerald Mathias, David Wifling, Yann Donon, Alberto Di Meglio, Sofia Vallecorsa, Heng Ma, Anda Trifan, Arvind Ramanathan, Tom Brettin, Alexander Partin , et al. (4 additional authors not shown)

Abstract: The race to meet the challenges of the global pandemic has served as a reminder that the existing drug discovery process is expensive, inefficient and slow. There is a major bottleneck screening the vast number of potential small molecules to shortlist lead compounds for antiviral drug development. New opportunities to accelerate drug discovery lie at the interface between machine learning methods… ▽ More The race to meet the challenges of the global pandemic has served as a reminder that the existing drug discovery process is expensive, inefficient and slow. There is a major bottleneck screening the vast number of potential small molecules to shortlist lead compounds for antiviral drug development. New opportunities to accelerate drug discovery lie at the interface between machine learning methods, in this case developed for linear accelerators, and physics-based methods. The two in silico methods, each have their own advantages and limitations which, interestingly, complement each other. Here, we present an innovative infrastructural development that combines both approaches to accelerate drug discovery. The scale of the potential resulting workflow is such that it is dependent on supercomputing to achieve extremely high throughput. We have demonstrated the viability of this workflow for the study of inhibitors for four COVID-19 target proteins and our ability to perform the required large-scale calculations to identify lead antiviral compounds through repurposing on a variety of supercomputers. △ Less

Submitted 4 September, 2021; v1 submitted 4 March, 2021; originally announced March 2021.

Journal ref: Interface Focus. 2021. 11 (6): 20210018

The upgrade of EAST Safety and Interlock system

Authors: Z. C. Zhang, B. J. Xiao, Z. S. Ji, Y. Wang, F. Xia, Z. H. Xu

Abstract: The Experimental Advanced Superconducting Tokamak (EAST), a nation-level large-scale scientific project of China, plays a key role for the research of peaceful utilizations of fusion energy. The safety and interlock system (SIS) is in charge of the supervision and control of all the EAST components involved in the protection of human and tokamak from potential accidents. With the development of ph… ▽ More The Experimental Advanced Superconducting Tokamak (EAST), a nation-level large-scale scientific project of China, plays a key role for the research of peaceful utilizations of fusion energy. The safety and interlock system (SIS) is in charge of the supervision and control of all the EAST components involved in the protection of human and tokamak from potential accidents. With the development of physical experiment, the SIS had come close to reaching its limits for expandability. Therefore, a prototype for upgrading EAST SIS has been designed, and a fast architecture based on COTS FPGA is absorbed into the new SIS. This paper presents EAST machine and human protection mechanism and the architecture of the upgrading safety and interlock system. △ Less

Submitted 30 October, 2020; originally announced November 2020.

Quantifying Success in Science: An Overview

Authors: Xiaomei Bai, Hanxiao Pan, Jie Hou, Teng Guo, Ivan Lee, Feng Xia

Abstract: Quantifying success in science plays a key role in guiding funding allocations, recruitment decisions, and rewards. Recently, a significant amount of progresses have been made towards quantifying success in science. This lack of detailed analysis and summary continues a practical issue. The literature reports the factors influencing scholarly impact and evaluation methods and indices aimed at over… ▽ More Quantifying success in science plays a key role in guiding funding allocations, recruitment decisions, and rewards. Recently, a significant amount of progresses have been made towards quantifying success in science. This lack of detailed analysis and summary continues a practical issue. The literature reports the factors influencing scholarly impact and evaluation methods and indices aimed at overcoming this crucial weakness. We focus on categorizing and reviewing the current development on evaluation indices of scholarly impact, including paper impact, scholar impact, and journal impact. Besides, we summarize the issues of existing evaluation methods and indices, investigate the open issues and challenges, and provide possible solutions, including the pattern of collaboration impact, unified evaluation standards, implicit success factor mining, dynamic academic network embedding, and scholarly impact inflation. This paper should help the researchers obtaining a broader understanding of quantifying success in science, and identifying some potential research directions. △ Less

Submitted 10 August, 2020; originally announced August 2020.

Random Walks: A Review of Algorithms and Applications

Authors: Feng Xia, Jiaying Liu, Hansong Nie, Yonghao Fu, Liangtian Wan, Xiangjie Kong

Abstract: A random walk is known as a random process which describes a path including a succession of random steps in the mathematical space. It has increasingly been popular in various disciplines such as mathematics and computer science. Furthermore, in quantum mechanics, quantum walks can be regarded as quantum analogues of classical random walks. Classical random walks and quantum walks can be used to c… ▽ More A random walk is known as a random process which describes a path including a succession of random steps in the mathematical space. It has increasingly been popular in various disciplines such as mathematics and computer science. Furthermore, in quantum mechanics, quantum walks can be regarded as quantum analogues of classical random walks. Classical random walks and quantum walks can be used to calculate the proximity between nodes and extract the topology in the network. Various random walk related models can be applied in different fields, which is of great significance to downstream tasks such as link prediction, recommendation, computer vision, semi-supervised learning, and network embedding. In this paper, we aim to provide a comprehensive review of classical random walks and quantum walks. We first review the knowledge of classical random walks and quantum walks, including basic concepts and some typical algorithms. We also compare the algorithms based on quantum walks and classical random walks from the perspective of time complexity. Then we introduce their applications in the field of computer science. Finally we discuss the open issues from the perspectives of efficiency, main-memory volume, and computing time of existing algorithms. This study aims to contribute to this growing area of research by exploring random walks and quantum walks together. △ Less

Submitted 8 August, 2020; originally announced August 2020.

Comments: 13 pages, 4 figures

Journal ref: IEEE Transactions on Emerging Topics in Computational Intelligence, 4(2): 95-107, April 2020

Big Networks: A Survey

Authors: Hayat Dino Bedru, Shuo Yu, Xinru Xiao, Da Zhang, Liangtian Wan, He Guo, Feng Xia

Abstract: A network is a typical expressive form of representing complex systems in terms of vertices and links, in which the pattern of interactions amongst components of the network is intricate. The network can be static that does not change over time or dynamic that evolves through time. The complication of network analysis is different under the new circumstance of network size explosive increasing. In… ▽ More A network is a typical expressive form of representing complex systems in terms of vertices and links, in which the pattern of interactions amongst components of the network is intricate. The network can be static that does not change over time or dynamic that evolves through time. The complication of network analysis is different under the new circumstance of network size explosive increasing. In this paper, we introduce a new network science concept called big network. Big networks are generally in large-scale with a complicated and higher-order inner structure. This paper proposes a guideline framework that gives an insight into the major topics in the area of network science from the viewpoint of a big network. We first introduce the structural characteristics of big networks from three levels, which are micro-level, meso-level, and macro-level. We then discuss some state-of-the-art advanced topics of big network analysis. Big network models and related approaches, including ranking methods, partition approaches, as well as network embedding algorithms are systematically introduced. Some typical applications in big networks are then reviewed, such as community detection, link prediction, recommendation, etc. Moreover, we also pinpoint some critical open issues that need to be investigated further. △ Less

Submitted 8 August, 2020; originally announced August 2020.

Comments: 69 pages, 4 figures

Journal ref: Computer Science Review, Volume 37, August 2020, 100247

A new paradigm of dissipation-controllable, multi-scale resolving schemes for compressible flows

Authors: Xi Deng, Zhen-hua Jiang, Peter Vincent, Feng Xiao, Chao Yan

Abstract: The scale-resolving simulation of high speed compressible flow through direct numerical simulation (DNS) or large eddy simulation (LES) requires shock-capturing schemes to be more accurate for resolving broadband turbulence and robust for capturing strong shock waves. In this work, we develop a new paradigm of dissipation-controllable, shock capturing scheme to resolve multi-scale flow structures… ▽ More The scale-resolving simulation of high speed compressible flow through direct numerical simulation (DNS) or large eddy simulation (LES) requires shock-capturing schemes to be more accurate for resolving broadband turbulence and robust for capturing strong shock waves. In this work, we develop a new paradigm of dissipation-controllable, shock capturing scheme to resolve multi-scale flow structures in high speed compressible flow. This novel paradigm of shock-capturing scheme is named as PnTm-BVD-CD. The proposed PnTm-BVD-CD scheme has following desirable properties. First, it can capture large-scale discontinuous structures such as strong shock waves without obvious non-physical oscillations while resolving sharp contact, material interface and shear layer. Secondly, the numerical dissipation property of PnTm-BVD-CD can be effectively controlled between n+1 order upwind-biased scheme and non-dissipative n+2 order central scheme through a simple tunable parameter $λ$. Thirdly, with $λ=0.5$ the scheme can recover to n+2 order non-dissipative central interpolation for smooth solution over all wavenumber, which is preferable for solving small-scale structures in DNS as well as resolvable-scale in explicit LES. Finally, the under-resolved small-scale can be solved with dissipation controllable algorithm through so-called implicit LES (ILES) approach. △ Less

Submitted 14 July, 2020; originally announced July 2020.

Positivity-preserving hybrid DG/FV method with subcell resolution for compressible Euler equations with stiff source terms

Authors: Zhen-Hua Jiang, Xi Deng, Chao Yan, Feng Xiao, Jian Yu

Abstract: In numerical investigations of stringent problems in science and engineering research, non-physical negative density or pressure may emerge and cause blow-ups of the computation. We build on the realization that the positivity of density and pressure can be preserved in an efficient way of exploiting the advantage of discontinuous Galerkin and finite volume (DG/FV) hybrid computation framework. We… ▽ More In numerical investigations of stringent problems in science and engineering research, non-physical negative density or pressure may emerge and cause blow-ups of the computation. We build on the realization that the positivity of density and pressure can be preserved in an efficient way of exploiting the advantage of discontinuous Galerkin and finite volume (DG/FV) hybrid computation framework. We thus present an effective and simple method with great practical significance to maintain the positivity of density and pressure in solving the reactive Euler equations. The approach is able to maximize the multiscale capability of the DG method by achieving both positivity-preserving and oscillation-free solution in the subgrid level. In the designed scheme, a priori detection and computation with hyperbolic tangent function prevent the occurrence of negativity in the flux evaluation, while a posteriori detection and computation with the first-order Godunov scheme guarantee the positivity of subcell solution. The a priori computation achieves bounded reconstruction and less oscillatory solution so that the a posteriori computation can be active in few cells where the extremely complex computation condition appears. Furthermore, the valuable information from the reconstruction process is utilized to design the indication strategy that identifies the DG and FV cells, and the technique of adaptively choosing reconstruction candidates is adopted to overcome the excessive numerical dissipation in the shock-capturing scheme. Numerical tests, including demanding examples in stiff detonation simulation, demonstrate the positivity-preserving, non-oscillatory and subcell resolution property of the present method. △ Less

Submitted 11 July, 2020; originally announced July 2020.

A new correction method for quasi-Keplerian orbits

Authors: Yue Chen, Da-Zhu Ma, Fang Xia

Abstract: A pure two-body problem has seven integrals including the Kepler energy, the Laplace vector, and the angular momentum vector. However, only five of them are independent. When the five independent integrals are preserved, the two other dependent integrals are naturally preserved from a theoretical viewpoint; but they may not be either from a numerical computational viewpoint. Because of this, we us… ▽ More A pure two-body problem has seven integrals including the Kepler energy, the Laplace vector, and the angular momentum vector. However, only five of them are independent. When the five independent integrals are preserved, the two other dependent integrals are naturally preserved from a theoretical viewpoint; but they may not be either from a numerical computational viewpoint. Because of this, we use seven scale factors to adjust the integrated positions and velocities so that the adjusted solutions strictly satisfy the seven constraints. Noticing the existence of the two dependent integrals, we adopt the Newton iterative method combined with the singular value decomposition to calculate these factors. This correction scheme can be applied to perturbed two-body and N-body problems in the solar system. In this case, the seven quantities of each planet slowly vary with time. More accurate values can be given to the seven slowly-varying quantities by integrating the integral invariant relations of these quantities and the equations of motion. They should be satisfied with the adjusted solutions. Numerical tests show that the new method can significantly reduce the rapid growth of numerical errors of all orbital elements. △ Less

Submitted 24 May, 2020; originally announced May 2020.

Comments: 16 pages, 48 figures

Proofreading through spatial gradients

Authors: Vahe Galstyan, Kabir Husain, Fangzhou Xiao, Arvind Murugan, Rob Phillips

Abstract: Key enzymatic processes in biology use the nonequilibrium error correction mechanism called kinetic proofreading to enhance their specificity. Kinetic proofreading typically requires several dedicated structural features in the enzyme, such as a nucleotide hydrolysis site and multiple enzyme-substrate conformations that delay product formation. Such requirements limit the applicability and the ada… ▽ More Key enzymatic processes in biology use the nonequilibrium error correction mechanism called kinetic proofreading to enhance their specificity. Kinetic proofreading typically requires several dedicated structural features in the enzyme, such as a nucleotide hydrolysis site and multiple enzyme-substrate conformations that delay product formation. Such requirements limit the applicability and the adaptability of traditional proofreading schemes. Here, we explore an alternative conceptual mechanism of error correction that achieves delays between substrate binding and subsequent product formation by having these events occur at distinct physical locations. The time taken by the enzyme-substrate complex to diffuse from one location to another is leveraged to discard wrong substrates. This mechanism does not require dedicated structural elements on the enzyme, making it easier to overlook in experiments but also making proofreading tunable on the fly. We discuss how tuning the length scales of enzyme or substrate concentration gradients changes the fidelity, speed and energy dissipation, and quantify the performance limitations imposed by realistic diffusion and reaction rates in the cell. Our work broadens the applicability of kinetic proofreading, and sets the stage for the study of spatial gradients as a possible route to specificity. △ Less

Submitted 23 May, 2020; originally announced May 2020.

Electrical molecular switch addressed by chemical stimuli

Authors: H. Audi, Y. Viero, N. Alwhaibi, Z. Chen, M. Iazykov, A. Heynderickx, F. Xiao, D. Guerin, C Krzeminski, I. M. Grace, C. J. Lambert, O. Siri, D. Vuillaume, S Lenfant, H. Klein

Abstract: We demonstrate that the conductance switching of benzo-bis(imidazole) molecules upon protonation depends on the lateral functional groups. The protonated H-substituted molecule shows a higher conductance than the neutral one (Gpro>Gneu), while the opposite (Gneu>Gpro) is observed for a molecule laterally functionalized by amino-phenyl groups. These results are demonstrated at various scale lengths… ▽ More We demonstrate that the conductance switching of benzo-bis(imidazole) molecules upon protonation depends on the lateral functional groups. The protonated H-substituted molecule shows a higher conductance than the neutral one (Gpro>Gneu), while the opposite (Gneu>Gpro) is observed for a molecule laterally functionalized by amino-phenyl groups. These results are demonstrated at various scale lengths : self-assembled monolayer, tiny nanodot-molecule junction and single molecules. From ab-initio theoretical calculations, we conclude that for the H-substituted molecule, the result Gpro>Gneu is correctly explained by a reduction of the LUMO-HOMO gap, while for the amino-phenyl functionnalized molecule, the result Gneu>Gpro is consistent with a shift of HOMO, which reduces the density of states at the Fermi energy. △ Less

Submitted 30 April, 2020; originally announced May 2020.

Comments: Full paper and supporting information

Semimetals for high performance photodetection

Authors: Jing Liu, Fengnian Xia, Di Xiao, Javier García de Abajo, Dong Sun

Abstract: Semimetals are being explored for their unique advantages in low-energy high-speed photodetection, although they suffer from serious drawbacks such as an intrinsically high dark current. In this perspective, we envision the exploitation of topological effects in the photoresponse of these materials as a promising route to circumvent these problems. We overview recent studies on photodetection base… ▽ More Semimetals are being explored for their unique advantages in low-energy high-speed photodetection, although they suffer from serious drawbacks such as an intrinsically high dark current. In this perspective, we envision the exploitation of topological effects in the photoresponse of these materials as a promising route to circumvent these problems. We overview recent studies on photodetection based on graphene and other semimetals, and further discuss the exciting opportunities created by the topological effects, along with the additional requirements that they impose on photodetector designs. △ Less

Submitted 24 March, 2020; originally announced April 2020.

Comments: Perspective on photodetection based on topological semimetals

A global nonhydrostatic dynamical core on cubed sphere using multi-moment finite volume method: formulation and preliminary test

Authors: Chungang Chen, Xingliang Li, Feng Xiao, Xueshun Shen

Abstract: A nonhydrostatic dynamical core has been developed by using the multi-moment finite volume method that ensures the rigorous numerical conservation. To represent the spherical geometry free of polar problems, the cubed-sphere grid is adopted. A fourth-order multi-moment discretization formulation is applied to the nonhydrostatic governing equations cast in local curvilinear coordinates on each patc… ▽ More A nonhydrostatic dynamical core has been developed by using the multi-moment finite volume method that ensures the rigorous numerical conservation. To represent the spherical geometry free of polar problems, the cubed-sphere grid is adopted. A fourth-order multi-moment discretization formulation is applied to the nonhydrostatic governing equations cast in local curvilinear coordinates on each patch of cubed sphere through a gnomonic projection. In vertical direction, the height-based terrain-following grid is used to represent the topography. To get around the CFL stability restriction imposed by relatively small grid spacing in the vertical direction, the dimensional-splitting time integration using the HEVI (Horizontal Explicit and Vertical Implicit) strategy is implemented by applying the IMEX Runge-Kutta scheme. The proposed dynamical core preserves the fourth-order accuracy in spherical geometry and has been verified by the widely-used benchmark tests. The results of our numerical experiments show that the present numerical core has superior solution quality and great practical potential as a platform for atmospheric models. A new unified model for numerical weather prediction and global atmospheric circulation simulation based on this dynamical core is under development. △ Less

Submitted 13 April, 2020; originally announced April 2020.

Review on development of a scalable high-order nonhydrostatic multi-moment constrained finite volume dynamical core

Authors: Xingliang Li, Chungang Chen, Xueshun Shen, Feng Xiao

Abstract: This report summarizes the major progresses to develop the dynamic core for next-generation atmospherical model for both numerical weather prediction and climate simulation. The numerical framework is based on a general formulation, so-called multi-moment constrained finite volume (MCV) method, which is well-balanced among solution quality (accuracy and robustness), algorithmic simplicity, computa… ▽ More This report summarizes the major progresses to develop the dynamic core for next-generation atmospherical model for both numerical weather prediction and climate simulation. The numerical framework is based on a general formulation, so-called multi-moment constrained finite volume (MCV) method, which is well-balanced among solution quality (accuracy and robustness), algorithmic simplicity, computational efficiency and flexibility for model configuration. A local high-order limiting projection is also devised to remove spurious oscillations and noises in numerical solutions, which allows the numerical model working well alone without artificial diffusion or filter. The resulted numerical schemes are very simple, efficient and easy to implement for both structured and unstructured grids, which provide a promising plateform of great practical significance. We have implemented the MCV method to shallow water equations on various spherical grids, including Yin-Yang overset grid, cubed sphere grid and geodesic icosahedral grid, non-hydrostatic compressible atmosherical model under complex topographic boundary conditions. In addition, the moist dynamics simulation like moist thermal bubble has been validated by using direct microphysical feedback. We have also constructed a prototype of 3D global non-hydrostatic compressible atmosherical model on cubed sphere grid with an explicit/implicit-hybrid time integration scheme, which can be used as the base to develop global atmospheric GCM. All the MCV models have been verified with widely used benchmark tests. The numerical results show that the present MCV models have solution quality competitive to other exiting high order models. Parallelization of the MCV shallow water model on cubed sphere grid reveals its suitability for large scale parallel processing with desirable scalability. △ Less

Submitted 13 April, 2020; originally announced April 2020.

Comments: 49 pages, 33 figures

A positivity-preserving conservative Semi-Lagrangian Multi-moment Global Transport Model on the Cubed Sphere

Authors: Jie Tang, Chungang Chen, Xueshun Shen, Feng Xiao, Xingliang Li

Abstract: A positivity-preserving conservative semi-Lagrangian transport model by multi-moment finite volume method has been developed on the cubed-sphere grid. In this paper, two kinds of moments, i.e. point values (PV moment) at cell boundaries and volume integrated average (VIA) value, are defined within a single cell. The PV moment is updated by a conventional semi-Lagrangian method, while the VIA momen… ▽ More A positivity-preserving conservative semi-Lagrangian transport model by multi-moment finite volume method has been developed on the cubed-sphere grid. In this paper, two kinds of moments, i.e. point values (PV moment) at cell boundaries and volume integrated average (VIA) value, are defined within a single cell. The PV moment is updated by a conventional semi-Lagrangian method, while the VIA moment is cast by the flux form formulation that assures the exact numerical conservation. Different from the spatial approximation used in CSL2 (conservative semi-Lagrangian scheme with second order polynomial function) scheme, a monotonic rational function which can effectively remove non-physical oscillations and preserve the shape, is reconstructed in a single cell by the PV moment and VIA moment. The resulting scheme is inherently conservative and can allow a CFL number larger than one. Moreover, the scheme uses only one cell for spatial reconstruction, which is very easy for practical implementation. The proposed model is evaluated by several widely used benchmark tests on cubed-sphere geometry. Numerical results show that the proposed transport model can effectively remove unphysical oscillations compared with the CSL2 scheme and preserve the numerical non-negativity, and it has the potential to transport the tracers accurately in real atmospheric model. △ Less

Submitted 13 April, 2020; originally announced April 2020.

Comments: 16 pages, 11 figures

Development of a unified high-order nonhydrostatic multi-moment constrained finite volume dynamical core: derivation of flux-form governing equations in the general curvilinear coordinate system

Authors: Xingliang Li, Chungang Chen, Xueshun Shen, Feng Xiao

Abstract: In the manuscript we have derived the flux-form atmospheric governing equations in the general curvilinear coordinate system which is used by a high-order nonhydrostatic multi-moment constrained finite volume (MCV) dynamical core, and given the explicit formulations in the shallow-atmosphere approximation. In general curvilinear coordinate x^i(i = 1,2,3), unlike the Cartesian coordinate, the base… ▽ More In the manuscript we have derived the flux-form atmospheric governing equations in the general curvilinear coordinate system which is used by a high-order nonhydrostatic multi-moment constrained finite volume (MCV) dynamical core, and given the explicit formulations in the shallow-atmosphere approximation. In general curvilinear coordinate x^i(i = 1,2,3), unlike the Cartesian coordinate, the base vectors are not constants either in magnitude or direction. Following the representations such as base vectors, vector and tensor and so on in general curvilinear coordinate, we can obtain the differential relations of base vectors, the gradient and divergence operator etc. which are the component parts of the atmospheric governing equation. Then we apply them in the two specific curvilinear coordinate system: the spherical polar and cubed-sphere coordinates that are adopted in high-order nonhydrostatic MCV dynamical core. By switching the geometrics such as the metric tensors (covariant and contravariant), Jacobian of the transformation, the Christoffel symbol of the second kind between the spherical polar and cubed-sphere coordinates, the resulting flux-form governing equations in the specific coordinate system can be easily achieved. Of course, the Cartesian coordinate can be recovered. Noted that the projection metric tensors like spherical polar system and Cartesian coordinate become simple due to orthogonal properties of coordinate. △ Less

Submitted 13 April, 2020; originally announced April 2020.

Discontinuity-resolving shock-capturing schemes on unstructured grids

Authors: Lidong Cheng, Xi Deng, Bin Xie, Yi Jiang, Feng Xiao

Abstract: Solving compressible flows containing discontinuities remains a major challenge for numerical methods especially on unstructured grids. Thus in this work, we make contributions to shock capturing schemes on unstructured grids with aim of resolving discontinuities with low numerical dissipation. Different from conventional shock capturing schemes which only use polynomials as interpolation function… ▽ More Solving compressible flows containing discontinuities remains a major challenge for numerical methods especially on unstructured grids. Thus in this work, we make contributions to shock capturing schemes on unstructured grids with aim of resolving discontinuities with low numerical dissipation. Different from conventional shock capturing schemes which only use polynomials as interpolation functions on unstructured grids, the proposed scheme employs the linear polynomial as well as non-polynomial as reconstruction candidates. For linear polynomial, the second order MUSCL scheme with the MLP (Multi-dimensional Limiting Process) slope limiter is adopted. The multi-dimensional THINC (Tangent of Hyperbola for INterface Capturing) function with quadratic surface representation and Gaussian quadrature, so-called THINC/QQ, is used as the non-polynomial reconstruction candidate. With these reconstruction candidates, a multi-stage boundary variation diminishing (BVD) algorithm which aims to minimize numerical dissipation is designed on unstructured grids to select the final reconstruction function. The resulted shock capturing scheme is named as MUSCL-THINC/QQ-BVD. The performance of the proposed scheme is demonstrated through solving compressible single-phase and multi-phase problems where the discontinuity is the typical flow structure. The numerical results show that the proposed scheme is capable of capturing sharp discontinuous profiles without numerical oscillations as well as resolving vortices associated with Kelvin-Helmholtz instabilities along shear layers and material interfaces. In comparison with schemes only replying on high order polynomials, the proposed scheme shows significant improvement of resolution across discontinuities. Thus, this work provides an accurate and robust shock-capturing scheme to resolve discontinuities in compressible flows. △ Less

Submitted 20 March, 2020; originally announced March 2020.

Comments: 18 pages, 18 figures

MSC Class: 76M12 (Primary) 76N15; 76L05; 76T99 (Secondary) ACM Class: J.2

The dominance of big teams in China's scientific output

Authors: Linlin Liu, Jianfei Yu, Junming Huang, Feng Xia, Tao Jia

Abstract: Modern science is dominated by scientific productions from teams. A recent finding shows that teams with both large and small sizes are essential in research, prompting us to analyze the extent to which a country's scientific work is carried out by big/small teams. Here, using over 26 million publications from Web of Science, we find that China's research output is more dominated by big teams than… ▽ More Modern science is dominated by scientific productions from teams. A recent finding shows that teams with both large and small sizes are essential in research, prompting us to analyze the extent to which a country's scientific work is carried out by big/small teams. Here, using over 26 million publications from Web of Science, we find that China's research output is more dominated by big teams than the rest of the world, which is particularly the case in fields of natural science. Despite the global trend that more papers are done by big teams, China's drop in small team output is much steeper. As teams in China shift from small to large size, the team diversity that is essential for innovative works does not increase as much as that in other countries. Using the national average as the baseline, we find that the National Natural Science Foundation of China (NSFC) supports fewer small team works than the National Science Foundation of U.S. (NSF) does, implying that big teams are more preferred by grant agencies in China. Our finding provides new insights into the concern of originality and innovation in China, which urges a need to balance small and big teams. △ Less

Submitted 16 June, 2020; v1 submitted 29 February, 2020; originally announced March 2020.

The near surface vertical atmospheric electric field abnormality could be as a promising imminent precursor of major earthquakes

Authors: T. Chen, H. Wu, X. -X. Zhang, C. Wang, X. -B. Jin, Q. -M. Ma, J. -Y. Xu, S. -P. Duan, Z. -H. He, H. Li, S. -G. Xiao, X. -Z. Wang, X. -H Shen, Q. Guo, I. Roth, V. S. Makhmutov, Y. Liu, J. Luo, X. -J. Jiang, L. Dai, X. -D. Peng, X. Hu, L. Li, C. Zeng, J. -J. Song , et al. (6 additional authors not shown)

Abstract: A promising short term precursor of major earthquakes (EQ) is very crucial in saving people and preventing huge losses. Ez, atmospheric electrostatic field vertical component, under fair air conditions, is generally oriented downwards (positive). Anomalous negative Ez signals could be used as an indicator of a great number of radioactive gases which are released from great number of rock clefts ju… ▽ More A promising short term precursor of major earthquakes (EQ) is very crucial in saving people and preventing huge losses. Ez, atmospheric electrostatic field vertical component, under fair air conditions, is generally oriented downwards (positive). Anomalous negative Ez signals could be used as an indicator of a great number of radioactive gases which are released from great number of rock clefts just before major earthquakes. Enhanced emission of radon radioactive decay will produce an anomalously large number of ion pairs. The positive particles will be transported downward by the fair weather electrostatic field and pile up near the surface. Finally, obviously and abnormally, an oriented upward atmospheric electric field Ez near the ground could be formed. Therefore, monitoring this Ez may be applied effectively in earthquake warning. △ Less

Submitted 20 February, 2020; originally announced February 2020.

THINC-scaling scheme that unifies VOF and level set methods

Authors: Ronit Kumar, Lidong Cheng, Bin Xie, Feng Xiao

Abstract: We present a novel interface-capturing scheme, THINC-scaling, to unify the VOF (volume of fluid) and the level set methods, which have been developed as two completely different approaches widely used in various applications. The THINC-scaling scheme preserves at the samectime the advantages of both VOF and level set methods, i.e. the mass/volume conservation of the VOF method and the geometrical… ▽ More We present a novel interface-capturing scheme, THINC-scaling, to unify the VOF (volume of fluid) and the level set methods, which have been developed as two completely different approaches widely used in various applications. The THINC-scaling scheme preserves at the samectime the advantages of both VOF and level set methods, i.e. the mass/volume conservation of the VOF method and the geometrical faithfulness of the level set method. THINC-scaling scheme allows to represent interface with high-order polynomials, and has algorithmic simplicity which eases its implementation in unstructured grids. △ Less

Submitted 13 August, 2019; originally announced August 2019.

MSC Class: 65-XX; 76Txx ACM Class: G.1; J.2; J.6

Hour-scale persistent negative anomaly of atmospheric electrostatic field near the epicenter before earthquake

Authors: Tao Chen, Han Wu, Chi Wang, Xiaoxin Zhang, Xiaobin Jin, Qiming Ma, Jiyao Xu, Suping Duan, Zhaohai He, Hui Li, Saiguan Xiao, Xizhen Wang, Xuhui Shen, Quan Guo, Ilan Roth, Vladimir Makhmutov, Yong Liu, Jing Luo, Xiujie Jiang, Lei Dai, Xiaodong Peng, Xiong Hu, Lei Li, Chen Zeng, Jiajun Song , et al. (6 additional authors not shown)

Abstract: Although earthquake prediction is a big challenge in the world, some simple observational tools can capture many physical signals and demonstrate that an earthquake (EQ) may be forthcoming in short period. Many researchers have studied the significant variation of atmospheric electrostatic field related to the forthcoming earthquake. However, until now, there is not a compelling physical mechanism… ▽ More Although earthquake prediction is a big challenge in the world, some simple observational tools can capture many physical signals and demonstrate that an earthquake (EQ) may be forthcoming in short period. Many researchers have studied the significant variation of atmospheric electrostatic field related to the forthcoming earthquake. However, until now, there is not a compelling physical mechanism which can explain why atmospheric electrostatic abnormal signal could appear just before an earthquake. Here we present a precursor signal and propose a brief physical interpretation. Under fair air conditions, if the near-surface atmospheric electrostatic field Ez (oriented down when it is positive) presents a very stable negative anomaly (from -100 V/m to -5000 V/m), it will forebode that an earthquake (seismicity from 3-8) would take place in the next several to tens of hours within a distance less than 100 km. We name this prediction technique as "DC Ez Determination"(DED). In addition, the mechanism of such abnormal quasi-static electric field before a forthcoming earthquake has been proposed here: (1) Radon gas releases from the rock clefts near the epicenter during seismogenic process. (2) α particles are produced due to the radioactive decay of Radon gas in the air. (3) α particle ionizing radiation creates more positive particles in the air, which is much more effective than that β and γ particles produced during Radon radioactive decay. (4) The new positive particles change formal positive atmospheric electric field (Ez) into stably negative near the earth-surface. (5) The closer the instrument is to the epicenter, the more likely it is to observe the negative Ez signal related to the earthquake. It is recommended to establish an instrument network to capture the reliable precursor and make some warnings before the earthquake disaster. △ Less

Submitted 18 July, 2019; originally announced July 2019.

Comments: 18 pages, 3 figures, 2 supplementary figures

A New Injection and Acceleration Scheme of Positrons in the Laser-Plasma Bubble Regime

Authors: Z. Y. Xu, C. F. Xiao, H. Y. Lu, R. H. Hu, J. Q. Yu, Z. Gong, Y. R. Shou, J. X. Liu, C. Z. Xie, S. Y. Chen, H. G. Lu, T. Q. Xu, R. X. Li, N. Hafz, Z. Najmudin, P. P. Rajeev, D. Neely, X. Q. Yan

Abstract: A novel approach for positron injection and acceleration in laser driven plasma wakefield is proposed. A theoretical model is developed and confirmed through PIC simulation. One ring-shaped beam and one co-axially propagating Gaussian beam drive wakefields in a preformed plasma volume filled with both electrons and positrons. The laser's ponderomotive force as well as the charge separation force i… ▽ More A novel approach for positron injection and acceleration in laser driven plasma wakefield is proposed. A theoretical model is developed and confirmed through PIC simulation. One ring-shaped beam and one co-axially propagating Gaussian beam drive wakefields in a preformed plasma volume filled with both electrons and positrons. The laser's ponderomotive force as well as the charge separation force in the front bucket of the first bubble are utilized to provide the transverse momenta of injected positrons and those positrons can be trapped by the focusing field and then accelerated by the wakefield. The simulation shows that a relatively high-charge, quasi-monoenergetic positron beams can be obtained. △ Less

Submitted 8 December, 2019; v1 submitted 13 May, 2019; originally announced May 2019.

Comments: 6 pages, 4 figures, submit to Physical Review Letters

Journal ref: Phys. Rev. Accel. Beams 23, 091301 (2020)

Constructing high-order discontinuity-capturing schemes with linear-weight polynomials and boundary variation diminishing algorithm

Authors: Xi Deng, Yuya Shimizu, Feng Xiao

Abstract: In this study, a new framework of constructing very high order discontinuity-capturing schemes is proposed for finite volume method. These schemes, so-called $\mathrm{P}_{n}\mathrm{T}_{m}-\mathrm{BVD}$ (polynomial of $n$-degree and THINC function of $m$-level reconstruction based on BVD algorithm), are designed by employing high-order linear-weight polynomials and THINC (Tangent of Hyperbola for I… ▽ More In this study, a new framework of constructing very high order discontinuity-capturing schemes is proposed for finite volume method. These schemes, so-called $\mathrm{P}_{n}\mathrm{T}_{m}-\mathrm{BVD}$ (polynomial of $n$-degree and THINC function of $m$-level reconstruction based on BVD algorithm), are designed by employing high-order linear-weight polynomials and THINC (Tangent of Hyperbola for INterface Capturing) functions with adaptive steepness as the reconstruction candidates. The final reconstruction function in each cell is determined with a multi-stage BVD (Boundary Variation Diminishing) algorithm so as to effectively control numerical oscillation and dissipation. We devise the new schemes up to eleventh order in an efficient way by directly increasing the order of the underlying upwind scheme using linear-weight polynomials. The analysis of the spectral property and accuracy tests show that the new reconstruction strategy well preserves the low-dissipation property of the underlying upwind schemes with high-order linear-weight polynomials for smooth solution over all wave numbers and realizes $n+1$ order convergence rate. The performance of new schemes is examined through widely used benchmark tests, which demonstrate that the proposed schemes are capable of simultaneously resolving small-scale flow features with high resolution and capturing discontinuities with low dissipation. With outperforming results and simplicity in algorithm, the new reconstruction strategy shows great potential as an alternative numerical framework for computing nonlinear hyperbolic conservation laws that have discontinuous and smooth solutions of different scales. △ Less

Submitted 18 November, 2018; originally announced November 2018.

Comments: arXiv admin note: text overlap with arXiv:1811.01379

A fifth-order shock capturing scheme with BVD algorithm

Authors: Xi Deng, Yuya Shimizu, Feng Xiao

Abstract: A novel 5th-order shock capturing scheme is presented in this paper. The scheme, so-called P4-THINC-BVD (4th degree polynomial and THINC reconstruction based on BVD algorithm), is formulated as a two-stage cascade BVD (Boundary Variation Diminishing) algorithm following the BVD principle that minimizes the jumps of reconstructed values at cell boundaries. In the P4-THINC-BVD scheme, polynomial of… ▽ More A novel 5th-order shock capturing scheme is presented in this paper. The scheme, so-called P4-THINC-BVD (4th degree polynomial and THINC reconstruction based on BVD algorithm), is formulated as a two-stage cascade BVD (Boundary Variation Diminishing) algorithm following the BVD principle that minimizes the jumps of reconstructed values at cell boundaries. In the P4-THINC-BVD scheme, polynomial of degree four and THINC (Tangent of Hyperbola for INterface Capturing) functions with adaptive steepness are used as the candidate reconstruction functions. The final reconstruction function is selected from the candidate functions by a two-stage cascade BVD algorithm so as to effectively control numerical oscillation and dissipation. Spectral analysis and numerical verifications show that the P4-THINC-BVD scheme possesses the following desirable properties: 1) it effectively suppresses spurious numerical oscillation in the presence of strong shock or discontinuity; 2) it substantially reduces numerical dissipation errors; 3) it automatically retrieves the underlying linear 5th-order upwind scheme for smooth solution over all wave numbers; 4) it is able to resolve both smooth and discontinuous flow structures of all scales with substantially improved solution quality in comparison to other existing methods; and 5) it faithfully maintains the free-mode solutions in long term computation. P4-THINC-BVD, as well as the underlying idea presented in this paper, provides an innovative and practical approach to design high-fidelity numerical schemes for compressible flows involving strong discontinuities and flow structures of wide range scales. △ Less

Submitted 6 November, 2018; v1 submitted 4 November, 2018; originally announced November 2018.

Coupled THINC and level set method: a conservative interface capturing scheme with arbitrary-order surface representations

Authors: Longgen Qian, Yanhong Wei, Feng Xiao

Abstract: In this paper, we propose a simple and accurate numerical method for capturing moving interfaces on fixed Eulerian grids by coupling the Tangent of Hyperbola Interface Capturing (THINC) method and Level Set (LS) method. The innovative and practically-significant aspects of the proposed method, so-called THINC/LS method, lie in (1) representing the interface with polynomial of arbitrary-order rathe… ▽ More In this paper, we propose a simple and accurate numerical method for capturing moving interfaces on fixed Eulerian grids by coupling the Tangent of Hyperbola Interface Capturing (THINC) method and Level Set (LS) method. The innovative and practically-significant aspects of the proposed method, so-called THINC/LS method, lie in (1) representing the interface with polynomial of arbitrary-order rather than the plane representation commonly used in the Piecewise Linear Interface Calculation (PLIC) Volume of fluid (VOF) methods, (2) conserving rigorously the mass of the transported VOF field, and (3) providing a straightforward and easy-to-code algorithm for implementation. We verified the proposed scheme with the widely used benchmark tests. Numerical results show that this new method can obtain arbitrary-order accuracy for interface reconstruction and can produce more accurate results than the classical VOF methods. △ Less

Submitted 7 January, 2018; originally announced January 2018.

Some practical versions of boundary variation diminishing (BVD) algorithm

Authors: Xi Deng, Bin Xie, Feng Xiao

Abstract: This short note presents some variant schemes of boundary variation diminishing (BVD) algorithm in one dimension with the results of numerical tests for linear advection equation to facilitate practical use. In spite of being presented in 1D fashion, all the schemes are simple and easy to implement in multi-dimensions on structured and unstructured grids for nonlinear and system equations. This short note presents some variant schemes of boundary variation diminishing (BVD) algorithm in one dimension with the results of numerical tests for linear advection equation to facilitate practical use. In spite of being presented in 1D fashion, all the schemes are simple and easy to implement in multi-dimensions on structured and unstructured grids for nonlinear and system equations. △ Less

Submitted 1 August, 2017; originally announced August 2017.

A new shock-capturing scheme for stiff detonation waves problems

Authors: Xi Deng, Honghui Teng, Bin Xie, Feng Xiao

Abstract: A new approach to prevent spurious behavior caused by conventional shock-capturing schemes when solving stiff detonation waves problems is introduced in the present work. Due to smearing of discontinuous solution by the excessive numerical dissipation, conventional shock-capturing schemes have difficulties to obtain the correct location of detonation front without enough grids resolution. To overc… ▽ More A new approach to prevent spurious behavior caused by conventional shock-capturing schemes when solving stiff detonation waves problems is introduced in the present work. Due to smearing of discontinuous solution by the excessive numerical dissipation, conventional shock-capturing schemes have difficulties to obtain the correct location of detonation front without enough grids resolution. To overcome the excessive numerical errors around discontinuities by traditional discretized schemes used in non-reacting high speed compressible flow, we introduce a new shock-capturing scheme in which besides linear function constructed by MUSCL (Monotone Upstream-centered Schemes for Conservation Law) scheme, a step like THINC (Tangent of Hyperbola for INterface Capturing) function is also employed as another candidate in the reconstruction process. The final reconstruction function is determined by boundary variation diminishing (BVD) algorithm by which numerical dissipation around discontinuities can be reduced significantly. The new resulted shock-capturing scheme is named MUSCL-THINC-BVD. One- and two-dimensional comparative numerical tests about stiff detonation waves problems are conducted with the 5th order WENO (Weighted Essentially Non-Oscillatory) and MUSCL-THINC-BVD scheme respectively, which show MUSCL-THINC-BVD scheme can capture the correct position of detonation waves with improved resolution while WENO scheme, in spite of higher order, produces spurious waves. Compared with other existing methods which involves extra treatments by accepting the smeared out discontinuities profiles, the current method obtain the correct but also sharp detonation front by fundamentally reducing numerical dissipation errors from shock-capturing schemes. Thus the proposed approach is an effective but simple method to solve stiff detonation problems. △ Less

Submitted 3 August, 2017; originally announced August 2017.

Implementation of BVD (boundary variation diminishing) algorithm in simulations of compressible multiphase flows

Authors: Xi Deng, Satoshi Inaba, Bin Xie, Keh-Ming Shyue, Feng Xiao

Abstract: We present in this work a new reconstruction scheme, so-called MUSCL-THINC-BVD scheme, to solve the five-equation model for interfacial two phase flows. This scheme employs the traditional shock capturing MUSCL (Monotone Upstream-centered Schemes for Conservation Law) scheme as well as the interface sharpening THINC (Tangent of Hyperbola for INterface Capturing) scheme as two building-blocks of sp… ▽ More We present in this work a new reconstruction scheme, so-called MUSCL-THINC-BVD scheme, to solve the five-equation model for interfacial two phase flows. This scheme employs the traditional shock capturing MUSCL (Monotone Upstream-centered Schemes for Conservation Law) scheme as well as the interface sharpening THINC (Tangent of Hyperbola for INterface Capturing) scheme as two building-blocks of spatial reconstruction using the BVD (boundary variation diminishing) principle that minimizes the variations (jumps) of the reconstructed variables at cell boundaries, and thus effectively reduces the numerical dissipations in numerical solutions. The MUSCL-THINC-BVD scheme is implemented to all state variables and volume fraction, which realizes the consistency among volume fraction and other physical variables. Benchmark tests are carried out to verify the capability of the present method in capturing the material interface as a well-defined sharp jump in volume fraction, as well as significant improvement in solution quality. The proposed scheme is a simple and effective method of practical significance for simulating compressible interfacial multiphase flows. △ Less

Submitted 1 May, 2017; v1 submitted 26 April, 2017; originally announced April 2017.

A Simple, Fast and Fully Automated Approach for Midline Shift Measurement on Brain Computed Tomography

Authors: Huan-Chih Wang, Shih-Hao Ho, Furen Xiao, Jen-Hai Chou

Abstract: Brain CT has become a standard imaging tool for emergent evaluation of brain condition, and measurement of midline shift (MLS) is one of the most important features to address for brain CT assessment. We present a simple method to estimate MLS and propose a new alternative parameter to MLS: the ratio of MLS over the maximal width of intracranial region (MLS/ICWMAX). Three neurosurgeons and our aut… ▽ More Brain CT has become a standard imaging tool for emergent evaluation of brain condition, and measurement of midline shift (MLS) is one of the most important features to address for brain CT assessment. We present a simple method to estimate MLS and propose a new alternative parameter to MLS: the ratio of MLS over the maximal width of intracranial region (MLS/ICWMAX). Three neurosurgeons and our automated system were asked to measure MLS and MLS/ICWMAX in the same sets of axial CT images obtained from 41 patients admitted to ICU under neurosurgical service. A weighted midline (WML) was plotted based on individual pixel intensities, with higher weighted given to the darker portions. The MLS could then be measured as the distance between the WML and ideal midline (IML) near the foramen of Monro. The average processing time to output an automatic MLS measurement was around 10 seconds. Our automated system achieved an overall accuracy of 90.24% when the CT images were calibrated automatically, and performed better when the calibrations of head rotation were done manually (accuracy: 92.68%). MLS/ICWMAX and MLS both gave results in same confusion matrices and produced similar ROC curve results. We demonstrated a simple, fast and accurate automated system of MLS measurement and introduced a new parameter (MLS/ICWMAX) as a good alternative to MLS in terms of estimating the degree of brain deformation, especially when non-DICOM images (e.g. JPEG) are more easily accessed. △ Less

Submitted 2 March, 2017; originally announced March 2017.