Showing 1–50 of 154 results for author: Wu, K

Large-scale quantum reservoir learning with an analog quantum computer

Authors: Milan Kornjača, Hong-Ye Hu, Chen Zhao, Jonathan Wurtz, Phillip Weinberg, Majd Hamdan, Andrii Zhdanov, Sergio H. Cantu, Hengyun Zhou, Rodrigo Araiza Bravo, Kevin Bagnall, James I. Basham, Joseph Campo, Adam Choukri, Robert DeAngelo, Paige Frederick, David Haines, Julian Hammett, Ning Hsu, Ming-Guang Hu, Florian Huber, Paul Niklas Jepsen, Ningyuan Jia, Thomas Karolyshyn, Minho Kwon , et al. (28 additional authors not shown)

Abstract: Quantum machine learning has gained considerable attention as quantum technology advances, presenting a promising approach for efficiently learning complex data patterns. Despite this promise, most contemporary quantum methods require significant resources for variational parameter optimization and face issues with vanishing gradients, leading to experiments that are either limited in scale or lac… ▽ More Quantum machine learning has gained considerable attention as quantum technology advances, presenting a promising approach for efficiently learning complex data patterns. Despite this promise, most contemporary quantum methods require significant resources for variational parameter optimization and face issues with vanishing gradients, leading to experiments that are either limited in scale or lack potential for quantum advantage. To address this, we develop a general-purpose, gradient-free, and scalable quantum reservoir learning algorithm that harnesses the quantum dynamics of neutral-atom analog quantum computers to process data. We experimentally implement the algorithm, achieving competitive performance across various categories of machine learning tasks, including binary and multi-class classification, as well as timeseries prediction. Effective and improving learning is observed with increasing system sizes of up to 108 qubits, demonstrating the largest quantum machine learning experiment to date. We further observe comparative quantum kernel advantage in learning tasks by constructing synthetic datasets based on the geometric differences between generated quantum and classical data kernels. Our findings demonstrate the potential of utilizing classically intractable quantum correlations for effective machine learning. We expect these results to stimulate further extensions to different quantum hardware and machine learning paradigms, including early fault-tolerant hardware and generative machine learning tasks. △ Less

Submitted 2 July, 2024; originally announced July 2024.

Comments: 10 + 14 pages, 4 + 7 figures

An efficient singlet-triplet spin qubit to fiber interface assisted by a photonic crystal cavity

Authors: Kui Wu, Sebastian Kindel, Thomas Descamps, Tobias Hangleiter, Jan Christoph Müller, Rebecca Rodrigo, Florian Merget, Hendrik Bluhm, Jeremy Witzens

Abstract: We introduce a novel optical interface between a singlet-triplet spin qubit and a photonic qubit which would offer new prospects for future quantum communication applications. The interface is based on a 220 nm thick GaAs/Al-GaAs heterostructure membrane and features a gate-defined singlet-triplet qubit, a gate-defined optically active quantum dot, a photonic crystal cavity and a bot-tom gold refl… ▽ More We introduce a novel optical interface between a singlet-triplet spin qubit and a photonic qubit which would offer new prospects for future quantum communication applications. The interface is based on a 220 nm thick GaAs/Al-GaAs heterostructure membrane and features a gate-defined singlet-triplet qubit, a gate-defined optically active quantum dot, a photonic crystal cavity and a bot-tom gold reflector. All essential components can be lithographically defined and deterministically fabricated, which greatly increases the scalability of on-chip in-tegration. According to our FDTD simulations, the interface provides an overall coupling efficiency of 28.7% into a free space Gaussian beam, assuming an SiO2 interlayer filling the space between the reflector and the membrane. The performance can be further increased to 48.5% by undercutting this SiO2 interlayer below the photonic crystal. △ Less

Submitted 20 June, 2024; originally announced June 2024.

Journal ref: The 25th European Conference on Integrated Optics, Springer Proceedings in Physics 402, pp. 365-372, 2024

Structure-Preserving Oscillation-Eliminating Discontinuous Galerkin Schemes for Ideal MHD Equations: Locally Divergence-Free and Positivity-Preserving

Authors: Mengqing Liu, Kailiang Wu

Abstract: Numerically simulating magnetohydrodynamics (MHD) poses notable challenges, including the suppression of spurious oscillations near discontinuities (e.g., shocks) and preservation of essential physical structures (e.g., the divergence-free constraint of magnetic field and the positivity of density and pressure). This paper develops structure-preserving oscillation-eliminating discontinuous Galerki… ▽ More Numerically simulating magnetohydrodynamics (MHD) poses notable challenges, including the suppression of spurious oscillations near discontinuities (e.g., shocks) and preservation of essential physical structures (e.g., the divergence-free constraint of magnetic field and the positivity of density and pressure). This paper develops structure-preserving oscillation-eliminating discontinuous Galerkin (OEDG) schemes for ideal MHD. The schemes leverage a locally divergence-free (LDF) oscillation-eliminating (OE) procedure to suppress spurious oscillations while retaining the LDF property of magnetic field and many desirable attributes of original DG schemes, such as conservation, local compactness, and optimal convergence rates. The OE procedure is based on the solution operator of a novel damping equation, a linear system of ordinary differential equations that are exactly solvable without any discretization. The OE procedure is performed after each Runge-Kutta stage and does not impact DG spatial discretization, facilitating its easy integration into existing DG codes as an independent module. Moreover, this paper presents a rigorous positivity-preserving (PP) analysis of the LDF OEDG schemes on Cartesian meshes, utilizing the optimal convex decomposition technique and the geometric quasi-linearization (GQL) approach. Efficient PP LDF OEDG schemes are derived by incorporating appropriate discretization of Godunov-Powell source terms into only the discrete equations of cell averages, under a condition achievable through a simple PP limiter. Several one- and two-dimensional MHD tests verify the accuracy, effectiveness, and robustness of the proposed structure-preserving OEDG schemes. △ Less

Submitted 2 May, 2024; v1 submitted 25 April, 2024; originally announced April 2024.

Comments: 55 pages

Diamond surfaces with lateral gradients for systematic optimization of surface chemistry for relaxometry -- A low pressure plasma-based approach

Authors: Yuchen Tian, Ari R. Ortiz Moreno, Mayeul Chipaux, Kaiqi Wu, Felipe P. Perona Martinez, Hoda Shirzad, Thamir Hamoh, Aldona Mzyk, Patrick van Rijn, Romana Schirhagl

Abstract: Diamond is increasingly popular because of its unique material properties. Diamond defects called nitrogen vacancy (NV) centers allow measurements with unprecedented sensitivity. However, to achieve ideal sensing performance NV centers need to be within nanometers from the surface and are thus strongly dependent on the local surface chemistry. Several attempts have been made to compare diamond sur… ▽ More Diamond is increasingly popular because of its unique material properties. Diamond defects called nitrogen vacancy (NV) centers allow measurements with unprecedented sensitivity. However, to achieve ideal sensing performance NV centers need to be within nanometers from the surface and are thus strongly dependent on the local surface chemistry. Several attempts have been made to compare diamond surfaces. However, due to the high price of diamond crystals with shallow NV centers, a limited number of chemical modifications have been studied. Here, we developed a systematic method to investigate a continuity of different local environments with a varying density and nature of surface groups in a single experiment on a single diamond plate. To achieve this goal, we used diamonds with a shallow ensemble of NV centers and introduced a chemical gradient across the surface. More specifically we used air and hydrogen plasma. The gradients were formed by low pressure plasma treatment after masking with a right-angled triangular prism shield. As a result, the surface contained gradually more oxygen/hydrogen towards the open end of the shield. We then performed widefield relaxometry to determine the effect of surface chemistry on the sensing performance. As expected, relaxation times and thus sensing performance indeed varies along the gradient. △ Less

Submitted 18 April, 2024; originally announced April 2024.

Comments: Main text: 18 pages and 9 figures Supplementary materials: 6 pages and 6 figures

Polarized Adding Method of Discrete Ordinate Approximation for Ultraviolet-Visible and Near-Infrared Radiative Transfer

Authors: Kun Wu, Feng Zhang, Wenwen Li, Fengzi Bao, Yi-ning Shi

Abstract: The polarization characteristics of atmospheric scattering are important and should not be ignored in radiative transfer simulations. In this study, a new vector radiative transfer model called the polarized adding method of discrete ordinate approximation (POLDDA) is proposed for use in remote sensing applications for ultraviolet-visible and near-infrared spectra. The single-layer radiative trans… ▽ More The polarization characteristics of atmospheric scattering are important and should not be ignored in radiative transfer simulations. In this study, a new vector radiative transfer model called the polarized adding method of discrete ordinate approximation (POLDDA) is proposed for use in remote sensing applications for ultraviolet-visible and near-infrared spectra. The single-layer radiative transfer process and inhomogeneous multi-layer connection are solved using the discrete ordinate method (DOM) and adding methods, respectively. By combining the advantages of DOM and the adding method, the Stokes vector (including the I-, Q-, U-, and V-components) calculated using the new method conforms to the results of PolRadtran/RT3, whether in a Rayleigh scattering atmosphere or the water cloud case. Moreover, the relative root-mean-square error (RMSE) values of the Stokes vector for the test cases between MYSTIC and the new method or RT3 prove the accuracy of the proposed method. Meanwhile, the new method has a higher computational efficiency than RT3, particularly for an atmosphere with a large scattering optical depth. Unlike RT3, the computation time of the proposed method does not increase with the optical depth of each layer. △ Less

Submitted 16 April, 2024; originally announced April 2024.

Provably Convergent and Robust Newton-Raphson Method: A New Dawn in Primitive Variable Recovery for Relativistic MHD

Authors: Chaoyi Cai, Jianxian Qiu, Kailiang Wu

Abstract: A long-standing and formidable challenge faced by all conservative schemes for relativistic magnetohydrodynamics (RMHD) is the recovery of primitive variables from conservative ones. This process involves solving highly nonlinear equations subject to physical constraints. An ideal solver should be "robust, accurate, and fast -- it is at the heart of all conservative RMHD schemes," as emphasized in… ▽ More A long-standing and formidable challenge faced by all conservative schemes for relativistic magnetohydrodynamics (RMHD) is the recovery of primitive variables from conservative ones. This process involves solving highly nonlinear equations subject to physical constraints. An ideal solver should be "robust, accurate, and fast -- it is at the heart of all conservative RMHD schemes," as emphasized in [S.C. Noble et al., ApJ, 641:626-637, 2006]. Despite over three decades of research, seeking efficient solvers that can provably guarantee stability and convergence remains an open problem. This paper presents the first theoretical analysis for designing a robust, physical-constraint-preserving (PCP), and provably (quadratically) convergent Newton-Raphson (NR) method for primitive variable recovery in RMHD. Our key innovation is a unified approach for the initial guess, devised based on sophisticated analysis. It ensures that the NR iteration consistently converges and adheres to physical constraints. Given the extreme nonlinearity and complexity of the iterative function, the theoretical analysis is highly nontrivial and technical. We discover a pivotal inequality for delineating the convexity and concavity of the iterative function and establish theories to guarantee the PCP property and convergence. We also develop theories to determine a computable initial guess within a theoretical "safe" interval. Intriguingly, we find that the unique positive root of a cubic polynomial always falls within this interval. Our PCP NR method is versatile and can be seamlessly integrated into any RMHD scheme that requires the recovery of primitive variables, potentially leading to a broad impact in this field. As an application, we incorporate it into a discontinuous Galerkin method, resulting in fully PCP schemes. Several numerical experiments demonstrate the efficiency and robustness of the PCP NR method. △ Less

Submitted 8 April, 2024; originally announced April 2024.

Comments: 26 pages, 7 figures

Silicon Photonic Microresonator-Based High-Resolution Line-by-Line Pulse Shaping

Authors: Lucas M. Cohen, Kaiyi Wu, Karthik V. Myilswamy, Saleha Fatema, Navin B. Lingaraju, Andrew M. Weiner

Abstract: Optical pulse shaping stands as a formidable technique in ultrafast optics, radio-frequency photonics, and quantum communications. While existing systems rely on bulk optics or integrated platforms with planar waveguide sections for spatial dispersion, they face limitations in achieving finer (few- or sub-GHz) spectrum control. These methods either demand considerable space or suffer from pronounc… ▽ More Optical pulse shaping stands as a formidable technique in ultrafast optics, radio-frequency photonics, and quantum communications. While existing systems rely on bulk optics or integrated platforms with planar waveguide sections for spatial dispersion, they face limitations in achieving finer (few- or sub-GHz) spectrum control. These methods either demand considerable space or suffer from pronounced phase errors and optical losses when assembled to achieve fine resolution. Addressing these challenges, we present a foundry-fabricated six-channel silicon photonic shaper using microresonator filter banks with inline phase control and high spectral resolution. Leveraging existing comb-based spectroscopic techniques, we devise a novel system to mitigate thermal crosstalk and enable the versatile use of our on-chip shaper. Our results demonstrate the shaper's ability to phase-compensate six comb lines at tunable channel spacings of 3, 4, and 5 GHz. Specifically, at a 3 GHz channel spacing, we showcase the generation of arbitrary waveforms in the time domain. This scalable design and control scheme holds promise in meeting future demands for high-precision spectral shaping capabilities. △ Less

Submitted 1 March, 2024; originally announced March 2024.

Comments: 21 pages, 12 figures

High-order accurate positivity-preserving and well-balanced discontinuous Galerkin schemes for ten-moment Gaussian closure equations with source terms

Authors: Jiangfu Wang, Huazhong Tang, Kailiang Wu

Abstract: This paper proposes novel high-order accurate discontinuous Galerkin (DG) schemes for the one- and two-dimensional ten-moment Gaussian closure equations with source terms defined by a known potential function. Our DG schemes exhibit the desirable capability of being well-balanced (WB) for a known hydrostatic equilibrium state while simultaneously preserving positive density and positive-definite a… ▽ More This paper proposes novel high-order accurate discontinuous Galerkin (DG) schemes for the one- and two-dimensional ten-moment Gaussian closure equations with source terms defined by a known potential function. Our DG schemes exhibit the desirable capability of being well-balanced (WB) for a known hydrostatic equilibrium state while simultaneously preserving positive density and positive-definite anisotropic pressure tensor. The well-balancedness is built on carefully modifying the solution states in the Harten-Lax-van Leer-contact (HLLC) flux, and appropriate reformulation and discretization of the source terms. Our novel modification technique overcomes the difficulties posed by the anisotropic effects, maintains the high-order accuracy, and ensures that the modified solution state remains within the physically admissible state set. Positivity-preserving analyses of our WB DG schemes are conducted by using several key properties of the admissible state set, the HLLC flux and the HLLC solver, as well as the geometric quasilinearization (GQL) approach in [Wu & Shu, SIAM Review, 65: 1031-1073, 2023], which was originally applied to analyze the admissible state set and physical-constraints-preserving schemes for the relativistic magnetohydrodynamics in [Wu & Tang, M3AS, 27: 1871-1928, 2017], to address the difficulties arising from the nonlinear constraints on pressure tensor. Moreover, the proposed WB DG schemes satisfy the weak positivity for the cell averages, implying the use of a scaling limiter to enforce the physical admissibility of the DG solution polynomials at certain points of interest. Extensive numerical experiments are conducted to validate the preservation of equilibrium states, accuracy in capturing small perturbations to such states, robustness in solving problems involving low density or low pressure, and high resolution for both smooth and discontinuous solutions. △ Less

Submitted 23 February, 2024; originally announced February 2024.

Comments: 45 pages, 11 figures

GQL-Based Bound-Preserving and Locally Divergence-Free Central Discontinuous Galerkin Schemes for Relativistic Magnetohydrodynamics

Authors: Shengrong Ding, Kailiang Wu

Abstract: This paper develops novel and robust central discontinuous Galerkin (CDG) schemes of arbitrarily high-order accuracy for special relativistic magnetohydrodynamics (RMHD) with a general equation of state (EOS). These schemes are provably bound-preserving (BP), i.e., consistently preserve the upper bound for subluminal fluid velocity and the positivity of density and pressure, while also (locally) m… ▽ More This paper develops novel and robust central discontinuous Galerkin (CDG) schemes of arbitrarily high-order accuracy for special relativistic magnetohydrodynamics (RMHD) with a general equation of state (EOS). These schemes are provably bound-preserving (BP), i.e., consistently preserve the upper bound for subluminal fluid velocity and the positivity of density and pressure, while also (locally) maintaining the divergence-free (DF) constraint for the magnetic field. For 1D RMHD, the standard CDG method is exactly DF, and its BP property is proven under a condition achievable by BP limiter. For 2D RMHD, we design provably BP and locally DF CDG schemes based on the suitable discretization of a modified RMHD system. A key novelty in our schemes is the discretization of additional source terms in the modified RMHD equations, so as to precisely counteract the influence of divergence errors on the BP property across overlapping meshes. We provide rigorous proofs of the BP property for our CDG schemes and first establish the theoretical connection between BP and discrete DF properties on overlapping meshes for RMHD. Owing to the absence of explicit expressions for primitive variables in terms of conserved variables, the constraints of physical bounds are strongly nonlinear, making the BP proofs highly nontrivial. We overcome these challenges through technical estimates within the geometric quasilinearization (GQL) framework, which converts the nonlinear constraints into linear ones. Furthermore, we introduce a new 2D cell average decomposition on overlapping meshes, which relaxes the theoretical BP CFL constraint and reduces the number of internal nodes, thereby enhancing the efficiency of the 2D BP CDG method. We implement the proposed CDG schemes for extensive RMHD problems with various EOSs, demonstrating their robustness and effectiveness in challenging scenarios. △ Less

Submitted 23 February, 2024; originally announced February 2024.

Comments: 47 pages, 14 figures

A robust near-field body area network based on coaxially-shielded textile metamaterial

Authors: Xia Zhu, Ke Wu, Xiaohang Xie, Stephan W. Anderson, Xin Zhang

Abstract: A body area network (BAN) involving wearable sensors populated around the human body can continuously monitor physiological signals, finding applications in personal healthcare and athletic evaluation. Existing near-field communication (NFC)-enabled BAN solutions, while facilitating reliable and secure interconnection among battery-free sensors, face challenges such as limited spectral stability a… ▽ More A body area network (BAN) involving wearable sensors populated around the human body can continuously monitor physiological signals, finding applications in personal healthcare and athletic evaluation. Existing near-field communication (NFC)-enabled BAN solutions, while facilitating reliable and secure interconnection among battery-free sensors, face challenges such as limited spectral stability against external interference. Here we demonstrate a textile metamaterial featuring a coaxially-shielded internal structure designed to mitigate interference from extraneous loadings. The metamaterial can be patterned onto clothing to form a scalable, customizable network, enabling communication between NFC-enabled devices and developed battery-free textile NFC sensing nodes placed within the network. Proof of concept demonstration shows the metamaterial's robustness against mechanical deformation and exposure to lossy, conductive saline solutions, underscoring its potential applications in wet environments, particularly in athletic activities involving water or significant perspiration, offering insights for the future development of radio frequency components for a robust BAN at the system level. △ Less

Submitted 21 February, 2024; originally announced February 2024.

Ion Transport through Short Nanopores Modulated by Charged Exterior Surfaces

Authors: Long Ma, Zhe Liu, Bowen Ai, Jia Man, Jianyong Li, Kechen Wu, Yinghua Qiu

Abstract: Short nanopores find extensive applications capitalizing on their high throughput and detection resolution. Ionic behaviors through long nanopores are mainly determined by charged inner-pore walls. When pore lengths decrease to sub-200 nm, charged exterior surfaces provide considerable modulation to ion current. We find that the charge status of inner-pore walls affects the modulation of ion curre… ▽ More Short nanopores find extensive applications capitalizing on their high throughput and detection resolution. Ionic behaviors through long nanopores are mainly determined by charged inner-pore walls. When pore lengths decrease to sub-200 nm, charged exterior surfaces provide considerable modulation to ion current. We find that the charge status of inner-pore walls affects the modulation of ion current from charged exterior surfaces. For 50-nm-long nanopores with neutral inner-pore walls, charged exterior surfaces on the voltage (surfaceV) and ground (surfaceG) sides enhance and inhibit ion transport by forming ion enrichment and depletion zones inside nanopores, respectively. For nanopores with both charged inner-pore and exterior surfaces, continuous electric double layers enhance ion transport through nanopores significantly. The charged surfaceV results in higher ion current by simultaneously weakening ion depletion at pore entrances and enhancing the intra-pore ion enrichment. The charged surfaceG expedites the exit of ions from nanopores, resulting in a decrease in ion enrichment at pore exits. Through adjustment in the width of charged-ring regions near pore boundaries, the effective charged width of the charged exterior is explored at ~20nm. Our results may provide a theoretical guide for further optimizing the performance of nanopore-based applications, like seawater desalination, biosensing, and osmotic energy conversion. △ Less

Submitted 15 February, 2024; originally announced February 2024.

Comments: 18 pages, 5 figures

Journal ref: J. Chem. Phys. 160, 074701 (2024)

Chiral switching of many-body steady states in a dissipative Rydberg gas

Authors: Chongwu Xie, Konghao Sun, Kang-Da Wu, Chuan-Feng Li, Guang-Can Guo, Wei Yi, Guo-Yong Xiang

Abstract: Dissipative Rydberg gases are an outstanding platform for the investigation of many-body quantum open systems. Despite the wealth of existing studies, the non-equilibrium dynamics of dissipative Rydberg gases are rarely examined or harnessed from the perspective of non-Hermitian physics, which is but intrinsic to open systems. Here we report the experimental observation of a chiral switching betwe… ▽ More Dissipative Rydberg gases are an outstanding platform for the investigation of many-body quantum open systems. Despite the wealth of existing studies, the non-equilibrium dynamics of dissipative Rydberg gases are rarely examined or harnessed from the perspective of non-Hermitian physics, which is but intrinsic to open systems. Here we report the experimental observation of a chiral switching between many-body steady states in a dissipative thermal Rydberg vapor, where the interplay of many-body effects and non-Hermiticity plays a key role. Specifically, as the parameters are adiabatically varied around a closed contour, depending on the chirality of the parameter modulation, the Rydberg vapor can change between two collective steady states with distinct Rydberg excitations and optical transmissions. Adopting a mean-field description, we reveal that both the existence of the bistable steady states and chiral dynamics derive from an exceptional structure in the parameter space, where multiple steady states of the many-body Liouvillian superoperator coalesce. We demonstrate that both the exceptional structure and the resulting state-switching dynamics are tunable through microwave dressing and temperature variations, confirming their reliance on the many-body dissipative nature of the Rydberg vapor. △ Less

Submitted 5 February, 2024; originally announced February 2024.

Nonlinearity-Enhanced Continuous Microwave Detection Based on Stochastic Resonance

Authors: Kang-Da Wu, Chongwu Xie, Chuan-Feng Li, Guang-Can Guo, Chang-Ling Zou, Guo-Yong Xiang

Abstract: In practical sensing tasks, noise is usually regarded as an obstruction to better performance and will degrade the sensitivity. Fortunately, \textit{stochastic resonance} (SR), a counterintuitive concept, can utilize noise to greatly enhance the detected signal amplitude. Although fundamentally important as a mechanism of weak signal amplification, and has been continually explored in geological,… ▽ More In practical sensing tasks, noise is usually regarded as an obstruction to better performance and will degrade the sensitivity. Fortunately, \textit{stochastic resonance} (SR), a counterintuitive concept, can utilize noise to greatly enhance the detected signal amplitude. Although fundamentally important as a mechanism of weak signal amplification, and has been continually explored in geological, biological, and physical science, both theoretically and experimentally, SR has yet to be demonstrated in realistic sensing tasks. Here we develop a novel SR-based nonlinear sensor using a thermal ensemble of interacting Rydberg atoms. With the assistance of stochastic noise (either inherently in the system or added externally) and strong nonlinearity in the Rydberg ensembles, the signal encoded in a weak MW field is greatly enhanced (over 25 dB). Moreover, we show that the SR-based atomic sensor can achieve a better sensitivity in our system, which is over 6.6 dB compared to a heterodye atomic sensor. Our results show potential advantage of SR-based MW sensors in commercial or defense-related applications. △ Less

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

Liquid-State Semiconductor Lasers Based on Type-(I+II) Colloidal Quantum Dots

Authors: Donghyo Hahm, Valerio Pinchetti, Clément Livache, Namyoung Ahn, Jungchul Noh, Xueyang Li, Jun Du, Kaifeng Wu, Victor I. Klimov

Abstract: Present-day liquid-state lasers are based on organic dyes. Here we demonstrate an alternative class of liquid lasers that employ solutions of colloidal quantum dots (QDs). Previous efforts to realize such devices have been hampered by fast nonradiative Auger recombination of multi-carrier states needed for optical gain. We overcome this challenge using type-(I+II) QDs that feature a trion-like opt… ▽ More Present-day liquid-state lasers are based on organic dyes. Here we demonstrate an alternative class of liquid lasers that employ solutions of colloidal quantum dots (QDs). Previous efforts to realize such devices have been hampered by fast nonradiative Auger recombination of multi-carrier states needed for optical gain. We overcome this challenge using type-(I+II) QDs that feature a trion-like optical-gain state with strongly suppressed Auger recombination. When combined with a Littrow optical cavity, static (non-circulated) solutions of these QDs exhibit stable lasing tunable from 634 nm to 594 nm. These results point towards the feasibility of technologically viable dye-like QD lasers that feature wide spectral tunability and, importantly, allow for stable operation without the need for a bulky circulation system, a standard attribute of traditional dye lasers. △ Less

Submitted 16 January, 2024; originally announced January 2024.

Comments: Main text (29 pages, 4 Figs) + Supplementary Info (15 pages, 10 Figures)

Report number: LA-UR-23-33109

Wireless, Customizable Coaxially-shielded Coils for Magnetic Resonance Imaging

Authors: Ke Wu, Xia Zhu, Stephan W. Anderson, Xin Zhang

Abstract: Anatomy-specific RF receive coil arrays routinely adopted in magnetic resonance imaging (MRI) for signal acquisition, are commonly burdened by their bulky, fixed, and rigid configurations, which may impose patient discomfort, bothersome positioning, and suboptimal sensitivity in certain situations. Herein, leveraging coaxial cables' inherent flexibility and electric field confining property, for t… ▽ More Anatomy-specific RF receive coil arrays routinely adopted in magnetic resonance imaging (MRI) for signal acquisition, are commonly burdened by their bulky, fixed, and rigid configurations, which may impose patient discomfort, bothersome positioning, and suboptimal sensitivity in certain situations. Herein, leveraging coaxial cables' inherent flexibility and electric field confining property, for the first time, we present wireless, ultra-lightweight, coaxially-shielded MRI coils achieving a signal-to-noise ratio (SNR) comparable to or surpassing that of commercially available cutting-edge receive coil arrays with the potential for improved patient comfort, ease of implementation, and significantly reduced costs. The proposed coils demonstrate versatility by functioning both independently in form-fitting configurations, closely adapting to relatively small anatomical sites, and collectively by inductively coupling together as metamaterials, allowing for extension of the field-of-view of their coverage to encompass larger anatomical regions without compromising coil sensitivity. The wireless, coaxially-shielded MRI coils reported herein pave the way toward next generation MRI coils. △ Less

Submitted 19 December, 2023; originally announced December 2023.

Node-downloadable frequency transfer system based on a mode-locked laser with over 100 km of fiber

Authors: Ziyi Jin, Ziyang Chen, Kai Wu, Dongrui Yu, Guohua Wu, Song Yu, Bin Luo, Hong Guo

Abstract: To meet the requirements of time-frequency networks and enable frequency downloadability for nodes along the link, we demonstrated the extraction of stable frequency signals at nodes using a mode-locked laser under the condition of 100 km laboratory fiber. The node consists of a simple structure that utilizes widely used optoelectronic devices and enables plug-and-play applications. In addition, t… ▽ More To meet the requirements of time-frequency networks and enable frequency downloadability for nodes along the link, we demonstrated the extraction of stable frequency signals at nodes using a mode-locked laser under the condition of 100 km laboratory fiber. The node consists of a simple structure that utilizes widely used optoelectronic devices and enables plug-and-play applications. In addition, the node can recover frequency signals with multiple frequencies, which are useful for scenarios that require different frequencies. Here, we experimentally demonstrated a short-term frequency instability of $2.83\times {{10}^{-13}}$@1 s and a long-term frequency instability of $1.18\times {{10}^{-15}}$@10,000 s at the node, which is similar to that at the remote site of the frequency transfer system. At the same time, frequency signals with different frequencies also achieved stable extraction with the same performance at the node. Our results can support the distributed application under large-scale time-frequency networks. △ Less

Submitted 16 December, 2023; originally announced December 2023.

Comments: 11 pages, 9 figures

Journal ref: Opt. Express 31, 39681-39694 (2023)

Wearable Coaxially-shielded Metamaterial for Magnetic Resonance Imaging

Authors: Xia Zhu, Ke Wu, Stephan W. Anderson, Xin Zhang

Abstract: Recent advancements in metamaterials have yielded the possibility of a wireless solution to improve signal-to-noise ratio (SNR) in magnetic resonance imaging (MRI). Unlike traditional closely packed local coil arrays with rigid designs and numerous components, these lightweight, cost-effective metamaterials eliminate the need for radio frequency (RF) cabling, baluns, adapters, and interfaces. Howe… ▽ More Recent advancements in metamaterials have yielded the possibility of a wireless solution to improve signal-to-noise ratio (SNR) in magnetic resonance imaging (MRI). Unlike traditional closely packed local coil arrays with rigid designs and numerous components, these lightweight, cost-effective metamaterials eliminate the need for radio frequency (RF) cabling, baluns, adapters, and interfaces. However, their clinical adoption has been limited by their low sensitivity, bulky physical footprint, and limited, specific use cases. Herein, we introduce a wearable metamaterial developed using commercially available coaxial cable, designed for a 3.0 T MRI system. This metamaterial inherits the coaxially-shielded structure of its constituent coaxial cable, effectively containing the electric field within the cable, thereby mitigating the electric coupling to its loading while ensuring safer clinical adoption, lower signal loss, and resistance to frequency shifts. Weighing only 50g, the metamaterial maximizes its sensitivity by conforming to the anatomical region of interest. MRI images acquired using this metamaterial with various pulse sequences demonstrate an up to 2-fold SNR enhancement when compared to a state-of-the-art 16-channel knee coil. This work introduces a novel paradigm for constructing metamaterials in the MRI environment, paving the way for the development of next-generation wireless MRI technology. △ Less

Submitted 15 December, 2023; originally announced December 2023.

Computational-design Enabled Wearable and Tunable Metamaterials via Freeform Auxetics for Magnetic Resonance Imaging

Authors: Ke Wu, Xia Zhu, Thomas G. Bifano, Stephan W. Anderson, Xin Zhang

Abstract: Metamaterials hold significant promise for enhancing the imaging capabilities of MRI machines as an additive technology, due to their unique ability to enhance local magnetic fields. However, despite their potential, the metamaterials reported in the context of MRI applications have often been impractical. This impracticality arises from their predominantly flat configurations and their susceptibi… ▽ More Metamaterials hold significant promise for enhancing the imaging capabilities of MRI machines as an additive technology, due to their unique ability to enhance local magnetic fields. However, despite their potential, the metamaterials reported in the context of MRI applications have often been impractical. This impracticality arises from their predominantly flat configurations and their susceptibility to shifts in resonance frequencies, preventing them from realizing their optimal performance. Here, we introduce a computational method for designing wearable and tunable metamaterials via freeform auxetics. The proposed computational-design tools yield an approach to solving the complex circle packing problems in an interactive and efficient manner, thus facilitating the development of deployable metamaterials configured in freeform shapes. With such tools, the developed metamaterials may readily conform to a patient's kneecap, ankle, head, or any part of the body in need of imaging, and while ensuring an optimal resonance frequency, thereby paving the way for the widespread adoption of metamaterials in clinical MRI applications. △ Less

Submitted 21 November, 2023; originally announced November 2023.

Unraveling Diffusion in Fusion Plasma: A Case Study of In Situ Processing and Particle Sorting

Authors: Junmin Gu, Paul Lin, Kesheng Wu, Seung-Hoe Ku, C. S. Chang, R. Michael Churchill, Jong Choi, Norbert Podhorszki, Scott Klasky

Abstract: This work starts an in situ processing capability to study a certain diffusion process in magnetic confinement fusion. This diffusion process involves plasma particles that are likely to escape confinement. Such particles carry a significant amount of energy from the burning plasma inside the tokamak to the diverter and damaging the diverter plate. This study requires in situ processing because of… ▽ More This work starts an in situ processing capability to study a certain diffusion process in magnetic confinement fusion. This diffusion process involves plasma particles that are likely to escape confinement. Such particles carry a significant amount of energy from the burning plasma inside the tokamak to the diverter and damaging the diverter plate. This study requires in situ processing because of the fast changing nature of the particle diffusion process. However, the in situ processing approach is challenging because the amount of data to be retained for the diffusion calculations increases over time, unlike in other in situ processing cases where the amount of data to be processed is constant over time. Here we report our preliminary efforts to control the memory usage while ensuring the necessary analysis tasks are completed in a timely manner. Compared with an earlier naive attempt to directly computing the same diffusion displacements in the simulation code, this in situ version reduces the memory usage from particle information by nearly 60% and computation time by about 20%. △ Less

Submitted 2 November, 2023; originally announced November 2023.

Modeling of an efficient singlet-triplet spin qubit to photon interface assisted by a photonic crystal cavity

Authors: Kui Wu, Sebastian Kindel, Thomas Descamps, Tobias Hangleiter, Jan Christoph Müller, Rebecca Rodrigo, Florian Merget, Hendrik Bluhm, Jeremy Witzens

Abstract: Efficient interconnection between distant semiconductor spin qubits with the help of photonic qubits would offer exciting new prospects for future quantum communication applications. In this paper, we optimize the extraction efficiency of a novel interface between a singlet-triplet spin qubit and a photonic qubit. The interface is based on a 220 nm thick GaAs/AlGaAs heterostructure membrane and co… ▽ More Efficient interconnection between distant semiconductor spin qubits with the help of photonic qubits would offer exciting new prospects for future quantum communication applications. In this paper, we optimize the extraction efficiency of a novel interface between a singlet-triplet spin qubit and a photonic qubit. The interface is based on a 220 nm thick GaAs/AlGaAs heterostructure membrane and consists of a gate-defined double quantum dot (GDQD) supporting a singlet-triplet qubit, an optically active quantum dot (OAQD) consisting of a gate-defined exciton trap, a photonic crystal cavity providing in-plane optical confinement and efficient out-coupling to an ideal free space Gaussian beam while accommodating the gate wiring of the GDQD and OAQD, and a bottom gold reflector to recycle photons and increase the optical extraction efficiency. All essential components can be lithographically defined and deterministically fabricated on the GaAs/AlGaAs heterostructure membrane, which greatly increases the scalability of on-chip integration. According to our simulations, the interface provides an overall coupling efficiency of 28.7% into a free space Gaussian beam, assuming an SiO2 interlayer filling the space between the reflector and the membrane. The performance can be further increased by undercutting this SiO2 interlayer below the photonic crystal. In this case, the overall efficiency is calculated to be 48.5%. △ Less

Submitted 28 October, 2023; originally announced October 2023.

Single-pixel imaging based on deep learning

Authors: Kai Song, Yaoxing Bian, Ku Wu, Hongrui Liu, Shuangping Han, Jiaming Li, Jiazhao Tian, Chengbin Qin, Jianyong Hu, Liantuan Xiao

Abstract: Single-pixel imaging can collect images at the wavelengths outside the reach of conventional focal plane array detectors. However, the limited image quality and lengthy computational times for iterative reconstruction still impede the practical application of single-pixel imaging. Recently, deep learning has been introduced into single-pixel imaging, which has attracted a lot of attention due to i… ▽ More Single-pixel imaging can collect images at the wavelengths outside the reach of conventional focal plane array detectors. However, the limited image quality and lengthy computational times for iterative reconstruction still impede the practical application of single-pixel imaging. Recently, deep learning has been introduced into single-pixel imaging, which has attracted a lot of attention due to its exceptional reconstruction quality, fast reconstruction speed, and the potential to complete advanced sensing tasks without reconstructing images. Here, this advance is discussed and some opinions are offered. Firstly, based on the fundamental principles of single-pixel imaging and deep learning, the principles and algorithms of single-pixel imaging based on deep learning are described and analyzed. Subsequently, the implementation technologies of single-pixel imaging based on deep learning are reviewed. They are divided into super-resolution single-pixel imaging, single-pixel imaging through scattering media, photon-level single-pixel imaging, optical encryption based on single-pixel imaging, color single-pixel imaging, and image-free sensing according to diverse application fields. Finally, major challenges and corresponding feasible approaches are discussed, as well as more possible applications in the future. △ Less

Submitted 16 November, 2023; v1 submitted 25 October, 2023; originally announced October 2023.

Results and Limits of Time Division Multiplexing for the BICEP Array High Frequency Receivers

Authors: S. Fatigoni, P. A. R. Ade, Z. Ahmed, M. Amiri, D. Barkats, R. Basu Thakur, C. A. Bischoff, D. Beck, J. J. Bock, V. Buza, J. Cheshire, J. Connors, J. Cornelison, M. Crumrine, A. J. Cukierman, E. V. Denison, M. I. Dierickx, L. Duband, M. Eiben, J. P. Filippini, A. Fortes, M. Gao, C. Giannakopoulos, N. Goeckner-Wald, D. C. Goldfinger , et al. (62 additional authors not shown)

Abstract: Time-Division Multiplexing is the readout architecture of choice for many ground and space experiments, as it is a very mature technology with proven outstanding low-frequency noise stability, which represents a central challenge in multiplexing. Once fully populated, each of the two BICEP Array high frequency receivers, observing at 150GHz and 220/270GHz, will have 7776 TES detectors tiled on the… ▽ More Time-Division Multiplexing is the readout architecture of choice for many ground and space experiments, as it is a very mature technology with proven outstanding low-frequency noise stability, which represents a central challenge in multiplexing. Once fully populated, each of the two BICEP Array high frequency receivers, observing at 150GHz and 220/270GHz, will have 7776 TES detectors tiled on the focal plane. The constraints set by these two receivers required a redesign of the warm readout electronics. The new version of the standard Multi Channel Electronics, developed and built at the University of British Columbia, is presented here for the first time. BICEP Array operates Time Division Multiplexing readout technology to the limits of its capabilities in terms of multiplexing rate, noise and crosstalk, and applies them in rigorously demanding scientific application requiring extreme noise performance and systematic error control. Future experiments like CMB-S4 plan to use TES bolometers with Time Division/SQUID-based readout for an even larger number of detectors. △ Less

Submitted 24 October, 2023; v1 submitted 16 October, 2023; originally announced October 2023.

Comments: 10 pages, 7 figures, Submitted to Journal of Low Temperature Physics

Journal ref: Journal of Low Temperature Physics (2024)

FEcMD: A multi-physics and multi-scale computational program for electron emission characteristics dynamically coupled with atomic structure in metal nano-emitters

Authors: Nan Li, Xinyu Gao, Xianghui Feng, Kai Wu, Yonghong Cheng, Bing Xiao

Abstract: Field emission coupled with molecular dynamics simulation (FEcMD) software package is a computational tool for studying the electron emission characteristics and the atomic structure evolution of micro- and nano-protrusions made of pure metals or multi-component alloys by means of multi-physics and multi-scale methodology. The implementations of molecular dynamics, the electrodynamics, and the hea… ▽ More Field emission coupled with molecular dynamics simulation (FEcMD) software package is a computational tool for studying the electron emission characteristics and the atomic structure evolution of micro- and nano-protrusions made of pure metals or multi-component alloys by means of multi-physics and multi-scale methodology. The implementations of molecular dynamics, the electrodynamics, and the heat conduction in FEcMD program are addressed. For molecular dynamics simulation, the Lennard-Jones potentials, embedded atomic method (EAM), and moment tensor potentials (MTP) are fully supported for both alloys and pure metals. In the electrodynamics, the FEcMD program incorporates the space charge fields (space charge potential and exchange-correlation effects) in the Wentzel-Kramers-Brillouin-Jeffreys (WKBJ) approximation to evaluate the field emission current density more reliably for nano-gaps between two metal electrodes. Additionally, the advanced two-temperature heat conduction model is implemented in FEcMD program, and which provides more reliable descriptions for the temperature evolutions of electron and phonon subsystems under the radiofrequency (RF) or pulse electric fields. Comprehensive benchmark tests are performed for each module in FEcMD software to validate the numerical results, and also to access the accuracy and efficiency of the implemented algorithms. Finally, some typical applications of FEcMD program are also demonstrated for understanding the evolution of temperature and electric field coupled with the dynamic changing of atomic structures for metal micro- and nano-protrusions. △ Less

Submitted 7 October, 2023; originally announced October 2023.

Conformal Metamaterials with Active Tunability and Self-adaptivity for Magnetic Resonance Imaging

Authors: Ke Wu, Xia Zhu, Xiaoguang Zhao, Stephan W. Anderson, Xin Zhang

Abstract: Ongoing effort has been devoted to applying metamaterials to boost the imaging performance of magnetic resonance imaging owing to their unique capacity for electromagnetic field confinement and enhancement. However, there are still major obstacles to widespread clinical adoption of conventional metamaterials due to several notable restrictions, namely: their typically bulky and rigid structures, d… ▽ More Ongoing effort has been devoted to applying metamaterials to boost the imaging performance of magnetic resonance imaging owing to their unique capacity for electromagnetic field confinement and enhancement. However, there are still major obstacles to widespread clinical adoption of conventional metamaterials due to several notable restrictions, namely: their typically bulky and rigid structures, deviations in their optimal resonance frequency, and their inevitable interference with the transmission RF field in MRI. Herein, we address these restrictions and report a conformal, smart metamaterial, which may not only be readily tuned to achieve the desired, precise frequency match with MRI by a controlling circuit, but is also capable of selectively amplifying the magnetic field during the RF reception phase by sensing the excitation signal strength passively, thereby remaining off during the RF transmission phase and thereby ensuring its optimal performance when applied to MRI as an additive technology. By addressing a host of current technological challenges, the metamaterial presented herein paves the way toward the wide-ranging utilization of metamaterials in clinical MRI, thereby translating this promising technology to the MRI bedside. △ Less

Submitted 29 September, 2023; originally announced October 2023.

Comments: 21 pages, 7 figures

Vernier Microcombs for Integrated Optical Atomic Clocks

Authors: Kaiyi Wu, Nathan P. O'Malley, Saleha Fatema, Cong Wang, Marcello Girardi, Mohammed S. Alshaykh, Zhichao Ye, Daniel E. Leaird, Minghao Qi, Victor Torres-Company, Andrew M. Weiner

Abstract: CMOS-compatible Kerr microcombs have drawn substantial interest as mass-manufacturable, compact alternatives to bulk frequency combs. This could enable deployment of many comb-reliant applications previously confined to laboratories. Particularly enticing is the prospect of microcombs performing optical frequency division in compact optical atomic clocks. Unfortunately, it is difficult to meet the… ▽ More CMOS-compatible Kerr microcombs have drawn substantial interest as mass-manufacturable, compact alternatives to bulk frequency combs. This could enable deployment of many comb-reliant applications previously confined to laboratories. Particularly enticing is the prospect of microcombs performing optical frequency division in compact optical atomic clocks. Unfortunately, it is difficult to meet the self-referencing requirement of microcombs in these systems due to the $\sim$THz repetition rates typically required for octave-spanning comb generation. Additionally, it is challenging to spectrally engineer a microcomb system to align a comb mode with an atomic clock transition with sufficient signal-to-noise ratio. Here, we adopt a Vernier dual-microcomb scheme for optical frequency division of a stabilized ultranarrow-linewidth continuous-wave laser at 871 nm to a $\sim$235 MHz output frequency. In addition to enabling measurement of the comb repetition rates, this scheme brings the freedom to pick comb lines from either or both of the combs. We exploit this flexibility to shift an ultra-high-frequency ($\sim$100 GHz) carrier-envelope offset beat down to frequencies where detection is possible and to place a comb line close to the 871 nm laser - tuned so that if frequency-doubled it would fall close to the clock transition in $^{171}$Yb$^+$. Moreover, we introduce a novel scheme which suppresses frequency noise arising from interferometric phase fluctuations in our dual-comb system and reduces the frequency instability down to our measurement limit. Our dual-comb system can potentially combine with an integrated ion trap toward future chip-scale optical atomic clocks. △ Less

Submitted 21 September, 2023; v1 submitted 17 August, 2023; originally announced August 2023.

Realization of a two-dimensional checkerboard lattice in monolayer Cu$_2$N

Authors: Xuegao Hu, Run-Wu Zhang, Da-Shuai Ma, Zhihao Cai, Daiyu Geng, Zhenyu Sun, Qiaoxiao Zhao, Jisong Gao, Peng Cheng, Lan Chen, Kehui Wu, Yugui Yao, Baojie Feng

Abstract: Two-dimensional checkerboard lattice, the simplest line-graph lattice, has been intensively studied as a toy model, while material design and synthesis remain elusive. Here, we report theoretical prediction and experimental realization of the checkerboard lattice in monolayer Cu$_2$N. Experimentally, monolayer Cu$_2$N can be realized in the well-known N/Cu(100) and N/Cu(111) systems that were prev… ▽ More Two-dimensional checkerboard lattice, the simplest line-graph lattice, has been intensively studied as a toy model, while material design and synthesis remain elusive. Here, we report theoretical prediction and experimental realization of the checkerboard lattice in monolayer Cu$_2$N. Experimentally, monolayer Cu$_2$N can be realized in the well-known N/Cu(100) and N/Cu(111) systems that were previously mistakenly believed to be insulators. Combined angle-resolved photoemission spectroscopy measurements, first-principles calculations, and tight-binding analysis show that both systems host checkerboard-derived hole pockets near the Fermi level. In addition, monolayer Cu$_2$N has outstanding stability in air and organic solvents, which is crucial for further device applications. △ Less

Submitted 8 June, 2023; originally announced June 2023.

Comments: Nano Letters, in press

Bayesian model calibration for diblock copolymer thin film self-assembly using power spectrum of microscopy data and machine learning surrogate

Authors: Lianghao Cao, Keyi Wu, J. Tinsley Oden, Peng Chen, Omar Ghattas

Abstract: Identifying parameters of computational models from experimental data, or model calibration, is fundamental for assessing and improving the predictability and reliability of computer simulations. In this work, we propose a method for Bayesian calibration of models that predict morphological patterns of diblock copolymer (Di-BCP) thin film self-assembly while accounting for various sources of uncer… ▽ More Identifying parameters of computational models from experimental data, or model calibration, is fundamental for assessing and improving the predictability and reliability of computer simulations. In this work, we propose a method for Bayesian calibration of models that predict morphological patterns of diblock copolymer (Di-BCP) thin film self-assembly while accounting for various sources of uncertainties in pattern formation and data acquisition. This method extracts the azimuthally-averaged power spectrum (AAPS) of the top-down microscopy characterization of Di-BCP thin film patterns as summary statistics for Bayesian inference of model parameters via the pseudo-marginal method. We derive the analytical and approximate form of a conditional likelihood for the AAPS of image data. We demonstrate that AAPS-based image data reduction retains the mutual information, particularly on important length scales, between image data and model parameters while being relatively agnostic to the aleatoric uncertainties associated with the random long-range disorder of Di-BCP patterns. Additionally, we propose a phase-informed prior distribution for Bayesian model calibration. Furthermore, reducing image data to AAPS enables us to efficiently build surrogate models to accelerate the proposed Bayesian model calibration procedure. We present the formulation and training of two multi-layer perceptrons for approximating the parameter-to-spectrum map, which enables fast integrated likelihood evaluations. We validate the proposed Bayesian model calibration method through numerical examples, for which the neural network surrogate delivers a fivefold reduction of the number of model simulations performed for a single calibration task. △ Less

Submitted 3 August, 2023; v1 submitted 8 June, 2023; originally announced June 2023.

Comments: Minor changes from the original submission, including a change in the title

A new discretely divergence-free positivity-preserving high-order finite volume method for ideal MHD equations

Authors: Shengrong Ding, Kailiang Wu

Abstract: This paper proposes and analyzes a novel efficient high-order finite volume method for the ideal magnetohydrodynamics (MHD). As a distinctive feature, the method simultaneously preserves a discretely divergence-free (DDF) constraint on the magnetic field and the positivity-preserving (PP) property, which ensures the positivity of density, pressure, and internal energy. To enforce the DDF condition… ▽ More This paper proposes and analyzes a novel efficient high-order finite volume method for the ideal magnetohydrodynamics (MHD). As a distinctive feature, the method simultaneously preserves a discretely divergence-free (DDF) constraint on the magnetic field and the positivity-preserving (PP) property, which ensures the positivity of density, pressure, and internal energy. To enforce the DDF condition, we design a new discrete projection approach that projects the reconstructed point values at the cell interface into a DDF space, without using any approximation polynomials. This projection method is highly efficient, easy to implement, and particularly suitable for standard high-order finite volume WENO methods, which typically return only the point values in the reconstruction. Moreover, we also develop a new finite volume framework for constructing provably PP schemes for the ideal MHD system. The framework comprises the discrete projection technique, a suitable approximation to the Godunov--Powell source terms, and a simple PP limiter. We provide rigorous analysis of the PP property of the proposed finite volume method, demonstrating that the DDF condition and the proper approximation to the source terms eliminate the impact of magnetic divergence terms on the PP property. The analysis is challenging due to the internal energy function's nonlinearity and the intricate relationship between the DDF and PP properties. To address these challenges, the recently developed geometric quasilinearization approach is adopted, which transforms a nonlinear constraint into a family of linear constraints. Finally, we validate the effectiveness of the proposed method through several benchmark and demanding numerical examples. The results demonstrate that the proposed method is robust, accurate, and highly effective, confirming the significance of the proposed DDF projection and PP techniques. △ Less

Submitted 24 May, 2023; originally announced May 2023.

Comments: 26 pages

Provably convergent Newton-Raphson methods for recovering primitive variables with applications to physical-constraint-preserving Hermite WENO schemes for relativistic hydrodynamics

Authors: Chaoyi Cai, Jianxian Qiu, Kailiang Wu

Abstract: The relativistic hydrodynamics (RHD) equations have three crucial intrinsic physical constraints on the primitive variables: positivity of pressure and density, and subluminal fluid velocity. However, numerical simulations can violate these constraints, leading to nonphysical results or even simulation failure. Designing genuinely physical-constraint-preserving (PCP) schemes is very difficult, as… ▽ More The relativistic hydrodynamics (RHD) equations have three crucial intrinsic physical constraints on the primitive variables: positivity of pressure and density, and subluminal fluid velocity. However, numerical simulations can violate these constraints, leading to nonphysical results or even simulation failure. Designing genuinely physical-constraint-preserving (PCP) schemes is very difficult, as the primitive variables cannot be explicitly reformulated using conservative variables due to relativistic effects. In this paper, we propose three efficient Newton--Raphson (NR) methods for robustly recovering primitive variables from conservative variables. Importantly, we rigorously prove that these NR methods are always convergent and PCP, meaning they preserve the physical constraints throughout the NR iterations. The discovery of these robust NR methods and their PCP convergence analyses are highly nontrivial and technical. As an application, we apply the proposed NR methods to design PCP finite volume Hermite weighted essentially non-oscillatory (HWENO) schemes for solving the RHD equations. Our PCP HWENO schemes incorporate high-order HWENO reconstruction, a PCP limiter, and strong-stability-preserving time discretization. We rigorously prove the PCP property of the fully discrete schemes using convex decomposition techniques. Moreover, we suggest the characteristic decomposition with rescaled eigenvectors and scale-invariant nonlinear weights to enhance the performance of the HWENO schemes in simulating large-scale RHD problems. Several demanding numerical tests are conducted to demonstrate the robustness, accuracy, and high resolution of the proposed PCP HWENO schemes and to validate the efficiency of our NR methods. △ Less

Submitted 24 May, 2023; originally announced May 2023.

Comments: 49 pages

Experimental observation of Kerr-Raman solitons in a normal-dispersion FP resonator

Authors: Tieying Li, Kan Wu, Xujia Zhang, Minglu Cai, Jianping Chen

Abstract: Different from the Kerr effect,stimulated Raman scattering (SRS) is a delayed response to molecular vibrations in materials. In microcavities, when driven in an anomalous group velocity dispersion (GVD) regime, SRS typically leads to self-frequency shift of solitons and generation of breather solitons which have been verified both theoretically and experimentally. However, when driven in a normal… ▽ More Different from the Kerr effect,stimulated Raman scattering (SRS) is a delayed response to molecular vibrations in materials. In microcavities, when driven in an anomalous group velocity dispersion (GVD) regime, SRS typically leads to self-frequency shift of solitons and generation of breather solitons which have been verified both theoretically and experimentally. However, when driven in a normal GVD regime, recent theoretical work predicts that SRS can cause the locking of switching waves (SWs) and thus support bright moving localized structures (LSs), which we term as Kerr-Raman solitons (KRSs). Limited by the design of suitable experimental parameters, experimental observation of the KRSs is not achieved yet. Here, we provide numerical investigation, and to our knowledge, the first experimental observation of these SRS enabled KRSs in a fiber Fabry-Perot (FP) resonator with ultra-low normal GVD. Such Kerr-Raman solitons exhibit localized temporal features with strong oscillations at ~13 THz local frequency on the top of a flat-top pulse. The corresponding spectrum is a low-noise and broadband Kerr comb with typical platicon-like spectrum in the center and two Raman Stokes and anti-Stokes peaks located near 13 THz away from the center. With such SRS enabled broadband Kerr comb, we have achieved a KRS spectrum with a repetition rate of ~3.68 GHz and a -40 dB spectral width of 260 nm. The corresponding comb tooth count is >9000, covering the S+C+L telecommunication bands. Moreover, the formation process of such KRSs is also revealed, and it is found that the GVD plays a key role in its generation. Our work will help to advance the study of the dynamics of optical frequency combs under the influence of SRS, as well as providing a broadband coherent mode-locked optical source for wide applications. △ Less

Submitted 23 May, 2023; originally announced May 2023.

Self-passivated freestanding superconducting oxide film for flexible electronics

Authors: Zhuoyue Jia, Chi Sin Tang, Jing Wu, Changjian Li, Wanting Xu, Kairong Wu, Difan Zhou, Ping Yang, Shengwei Zeng, Zhigang Zeng, Dengsong Zhang, Ariando Ariando, Mark B. H. Breese, Chuanbing Cai, Xinmao Yin

Abstract: The integration of high-temperature superconducting YBa2Cu3O6+x (YBCO) into flexible electronic devices has the potential to revolutionize the technology industry. The effective preparation of high-quality flexible YBCO films therefore plays a key role in this development. We present a novel approach for transferring water-sensitive YBCO films onto flexible substrates without any buffer layer. Fre… ▽ More The integration of high-temperature superconducting YBa2Cu3O6+x (YBCO) into flexible electronic devices has the potential to revolutionize the technology industry. The effective preparation of high-quality flexible YBCO films therefore plays a key role in this development. We present a novel approach for transferring water-sensitive YBCO films onto flexible substrates without any buffer layer. Freestanding YBCO film on a polydimethylsiloxane substrate is extracted by etching the Sr3Al2O6 sacrificial layer from the LaAlO3 substrate. In addition to the obtained freestanding YBCO thin film having a Tc of 89.1 K, the freestanding YBCO thin films under inward and outward bending conditions have Tc of 89.6 K and 88.9 K, respectively. A comprehensive characterization involving multiple experimental techniques including high-resolution transmission electron microscopy, scanning electron microscopy, Raman and X-ray Absorption Spectroscopy is conducted to investigate the morphology, structural and electronic properties of the YBCO film before and after the extraction process where it shows the preservation of the structural and superconductive properties of the freestanding YBCO virtually in its pristine state. Further investigation reveals the formation of a YBCO passivated layer serves as a protective layer which effectively preserves the inner section of the freestanding YBCO during the etching process. This work plays a key role in actualizing the fabrication of flexible oxide thin films and opens up new possibilities for a diverse range of device applications involving thin-films and low-dimensional materials. △ Less

Submitted 6 July, 2023; v1 submitted 8 May, 2023; originally announced May 2023.

Comments: 22 pages,4 figures,references added

Journal ref: Applied Physics Reviews 10, 031401 (2023)

A neural network-based scale-adaptive cloud-fraction scheme for GCMs

Authors: Guoxing Chen, Wei-Chyung Wang, Shixi Yang, Yixin Wang, Feng Zhang, Kun Wu

Abstract: Cloud fraction significantly affects the short- and long-wave radiation. Its realistic representation in general circulation models (GCMs) still poses great challenges in modeling the atmosphere. Here, we present a neural network-based diagnostic scheme that uses the grid-mean temperature, pressure, liquid and ice water mixing ratios, and relative humidity to simulate the sub-grid cloud fraction.… ▽ More Cloud fraction significantly affects the short- and long-wave radiation. Its realistic representation in general circulation models (GCMs) still poses great challenges in modeling the atmosphere. Here, we present a neural network-based diagnostic scheme that uses the grid-mean temperature, pressure, liquid and ice water mixing ratios, and relative humidity to simulate the sub-grid cloud fraction. The scheme, trained using CloudSat data with explicit consideration of grid sizes, realistically simulates the observed cloud fraction with a correlation coefficient (r) > 0.9 for liquid-, mixed-, and ice-phase clouds. The scheme also captures the observed non-monotonic relationship between cloud fraction and relative humidity and is computationally efficient, and robust for GCMs with a variety of horizontal and vertical resolutions. For illustrative purposes, we conducted comparative analyses of the 2006-2019 climatological-mean cloud fractions among CloudSat, and simulations from the new scheme and the Xu-Randall scheme (optimized the same way as the new scheme). The network-based scheme improves not only the spatial distribution of the total cloud fraction but also the cloud vertical structure (r > 0.99). For example, the biases of too-many high-level clouds over the tropics and too-many low-level clouds over regions around 60°S and 60°N in the Xu-Randall scheme are significantly reduced. These improvements are also found to be insensitive to the spatio-temporal variability of large-scale meteorology conditions, implying that the scheme can be used in different climate regimes. △ Less

Submitted 18 June, 2023; v1 submitted 4 April, 2023; originally announced April 2023.

Journal ref: Journal of Advances in Modeling Earth Systems, 2023

Deep-OSG: Deep Learning of Operators in Semigroup

Authors: Junfeng Chen, Kailiang Wu

Abstract: This paper proposes a novel deep learning approach for learning operators in semigroup, with applications to modeling unknown autonomous dynamical systems using time series data collected at varied time lags. It is a sequel to the previous flow map learning (FML) works [T. Qin, K. Wu, and D. Xiu, J. Comput. Phys., 395:620--635, 2019], [K. Wu and D. Xiu, J. Comput. Phys., 408:109307, 2020], and [Z.… ▽ More This paper proposes a novel deep learning approach for learning operators in semigroup, with applications to modeling unknown autonomous dynamical systems using time series data collected at varied time lags. It is a sequel to the previous flow map learning (FML) works [T. Qin, K. Wu, and D. Xiu, J. Comput. Phys., 395:620--635, 2019], [K. Wu and D. Xiu, J. Comput. Phys., 408:109307, 2020], and [Z. Chen, V. Churchill, K. Wu, and D. Xiu, J. Comput. Phys., 449:110782, 2022], which focused on learning single evolution operator with a fixed time step. This paper aims to learn a family of evolution operators with variable time steps, which constitute a semigroup for an autonomous system. The semigroup property is very crucial and links the system's evolutionary behaviors across varying time scales, but it was not considered in the previous works. We propose for the first time a framework of embedding the semigroup property into the data-driven learning process, through a novel neural network architecture and new loss functions. The framework is very feasible, can be combined with any suitable neural networks, and is applicable to learning general autonomous ODEs and PDEs. We present the rigorous error estimates and variance analysis to understand the prediction accuracy and robustness of our approach, showing the remarkable advantages of semigroup awareness in our model. Moreover, our approach allows one to arbitrarily choose the time steps for prediction and ensures that the predicted results are well self-matched and consistent. Extensive numerical experiments demonstrate that embedding the semigroup property notably reduces the data dependency of deep learning models and greatly improves the accuracy, robustness, and stability for long-time prediction. △ Less

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

Discover governing differential equations from evolving systems

Authors: Yuanyuan Li, Kai Wu, Jing Liu

Abstract: Discovering the governing equations of evolving systems from available observations is essential and challenging. In this paper, we consider a new scenario: discovering governing equations from streaming data. Current methods struggle to discover governing differential equations with considering measurements as a whole, leading to failure to handle this task. We propose an online modeling method c… ▽ More Discovering the governing equations of evolving systems from available observations is essential and challenging. In this paper, we consider a new scenario: discovering governing equations from streaming data. Current methods struggle to discover governing differential equations with considering measurements as a whole, leading to failure to handle this task. We propose an online modeling method capable of handling samples one by one sequentially by modeling streaming data instead of processing the entire dataset. The proposed method performs well in discovering ordinary differential equations (ODEs) and partial differential equations (PDEs) from streaming data. Evolving systems are changing over time, which invariably changes with system status. Thus, finding the exact change points is critical. The measurement generated from a changed system is distributed dissimilarly to before; hence, the difference can be identified by the proposed method. Our proposal is competitive in identifying the change points and discovering governing differential equations in three hybrid systems and two switching linear systems. △ Less

Submitted 16 July, 2023; v1 submitted 18 January, 2023; originally announced January 2023.

Comments: 13 pages, 5 figures. Accepted by Physical Review Research

Frequency and Amplitude Optimizations for Magnetic Particle Spectroscopy Applications

Authors: Vinit Kumar Chugh, Arturo di Girolamo, Venkatramana D. Krishna, Kai Wu, Maxim C-J Cheeran, Jian-Ping Wang

Abstract: Nowadays, there is a growing interest in the field of magnetic particle spectroscopy (MPS)-based bioassays. MPS monitors the dynamic magnetic response of surface-functionalized magnetic nanoparticles (MNPs) upon excitation by an alternating magnetic field (AMF) to detect various target analytes. This technology has flourished in the past decade due to its low cost, low background magnetic noise in… ▽ More Nowadays, there is a growing interest in the field of magnetic particle spectroscopy (MPS)-based bioassays. MPS monitors the dynamic magnetic response of surface-functionalized magnetic nanoparticles (MNPs) upon excitation by an alternating magnetic field (AMF) to detect various target analytes. This technology has flourished in the past decade due to its low cost, low background magnetic noise interference from biomatrix, and fast response time. A large number of MPS variants have been reported by different groups around the world, with applications ranging from disease diagnosis to foodborne pathogen detection, and virus detection. However, there is an urgent need for guidance on how to optimize the sensitivity of MPS detection by choosing different types of MNPs, AMF modalities, and MPS assay strategies (i.e., volume- and surface-based assays). In this work, we systematically study the effect of AMF frequencies and amplitudes on the responses of single- and multi-core MNPs under two extreme conditions, namely, the bound and unbound states. Our results show that some modalities such as dual-frequency MPS utilizing multicore MNPs are more suitable for surface-based bioassay applications, whereas, single-frequency MPS systems using single- or multi-core MNPs are better suited for volumetric bioassay applications. Furthermore, the bioassay sensitivities for these modalities can be further improved by careful selection of AMF frequencies and amplitudes. △ Less

Submitted 10 November, 2022; originally announced November 2022.

Hamiltonian Transformation for Band Structure Calculations

Authors: Kai Wu, Yingzhou Li, Wentiao Wu, Lin Lin, Wei Hu, Jinlong Yang

Abstract: First-principles electronic band structure calculations are essential for understanding periodic systems in condensed matter physics and materials science. We propose an accurate and parameter-free method, called Hamiltonian transformation (HT), to calculate band structures in both density functional theory (DFT) and post-DFT calculations with plane-wave basis sets. The cost of HT is independent o… ▽ More First-principles electronic band structure calculations are essential for understanding periodic systems in condensed matter physics and materials science. We propose an accurate and parameter-free method, called Hamiltonian transformation (HT), to calculate band structures in both density functional theory (DFT) and post-DFT calculations with plane-wave basis sets. The cost of HT is independent of the choice of the density functional and scales as $\mathcal{O}(N_e^3N_k\log N_k)$, where $N_e$ and $N_k$ are the number of electrons and the number of $\mathbf{k}$-points. Compared to the widely used Wannier interpolation (WI), HT adopts an eigenvalue transformation to construct a spatial localized representation of the spectrally truncated Hamiltonian. HT also uses a non-iterative algorithm to change the basis sets to circumvent the construction of the maximally localized Wannier functions. As a result, HT can significantly outperform WI in terms of the accuracy of the band structure calculation. We also find that the eigenvalue transformation can be of independent interest, and can be used to improve the accuracy of the WI for systems with entangled bands. △ Less

Submitted 5 November, 2022; originally announced November 2022.

Comments: 5 pages, 4 figures

DIICAN: Dual Time-scale State-Coupled Co-estimation of SOC, SOH and RUL for Lithium-Ion Batteries

Authors: Ningbo Cai, Yuwen Qin, Xin Chen, Kai Wu

Abstract: Accurate co-estimations of battery states, such as state-of-charge (SOC), state-of-health (SOH,) and remaining useful life (RUL), are crucial to the battery management systems to assure safe and reliable management. Although the external properties of the battery charge with the aging degree, batteries' degradation mechanism shares similar evolving patterns. Since batteries are complicated chemica… ▽ More Accurate co-estimations of battery states, such as state-of-charge (SOC), state-of-health (SOH,) and remaining useful life (RUL), are crucial to the battery management systems to assure safe and reliable management. Although the external properties of the battery charge with the aging degree, batteries' degradation mechanism shares similar evolving patterns. Since batteries are complicated chemical systems, these states are highly coupled with intricate electrochemical processes. A state-coupled co-estimation method named Deep Inter and Intra-Cycle Attention Network (DIICAN) is proposed in this paper to estimate SOC, SOH, and RUL, which organizes battery measurement data into the intra-cycle and inter-cycle time scales. And to extract degradation-related features automatically and adapt to practical working conditions, the convolutional neural network is applied. The state degradation attention unit is utilized to extract the battery state evolution pattern and evaluate the battery degradation degree. To account for the influence of battery aging on the SOC estimation, the battery degradation-related state is incorporated in the SOC estimation for capacity calibration. The DIICAN method is validated on the Oxford battery dataset. The experimental results show that the proposed method can achieve SOH and RUL co-estimation with high accuracy and effectively improve SOC estimation accuracy for the whole lifespan. △ Less

Submitted 20 October, 2022; originally announced October 2022.

5-Minute Magnetic Particle Spectroscopy (MPS)-based Bioassay: Towards Ultra-fast, Convenient, and Large-scale Disease Diagnosis

Authors: Kai Wu, Vinit Kumar Chugh, Venkatramana D. Krishna, Yongqiang Andrew Wang, Timothy D. Gordon, Maxim C-J Cheeran, Jian-Ping Wang

Abstract: In this work, we report a 5-minute magnetic particle spectroscopy (MPS)-based bioassay strategy. In our approach, surface functionalized magnetic nanoparticles (MNPs) are incubated with target analytes at 37 °C with agitation for 3 minutes and the MPS reading is then taken at the 5th minute. We have proved the feasibility of 5-minute ultra-fast detection of SARS-CoV-2 spike protein with a detectio… ▽ More In this work, we report a 5-minute magnetic particle spectroscopy (MPS)-based bioassay strategy. In our approach, surface functionalized magnetic nanoparticles (MNPs) are incubated with target analytes at 37 °C with agitation for 3 minutes and the MPS reading is then taken at the 5th minute. We have proved the feasibility of 5-minute ultra-fast detection of SARS-CoV-2 spike protein with a detection limit below 5 nM (0.2 pmole). Our proposed 5-minute bioassay strategy may be applied to reduce assay time for other liquid phase, volumetric biosensors such as nuclear magnetic resonance (NMR), quantum dots (QD), fluorescent biosensors, etc. △ Less

Submitted 25 September, 2022; originally announced October 2022.

Comments: 10 pages, 5 figures

Observation of Gapped Dirac Cones in a Two-Dimensional Su-Schrieffer-Heeger Lattice

Authors: Daiyu Geng, Hui Zhou, Shaosheng Yue, Zhenyu Sun, Peng Cheng, Lan Chen, Sheng Meng, Kehui Wu, Baojie Feng

Abstract: The Su-Schrieffer-Heeger (SSH) model in a two-dimensional rectangular lattice features gapless or gapped Dirac cones with topological edge states along specific peripheries. While such a simple model has been recently realized in photonic/acoustic lattices and electric circuits, its material realization in condensed matter systems is still lacking. Here, we study the atomic and electronic structur… ▽ More The Su-Schrieffer-Heeger (SSH) model in a two-dimensional rectangular lattice features gapless or gapped Dirac cones with topological edge states along specific peripheries. While such a simple model has been recently realized in photonic/acoustic lattices and electric circuits, its material realization in condensed matter systems is still lacking. Here, we study the atomic and electronic structure of a rectangular Si lattice on Ag(001) by angle-resolved photoemission spectroscopy and theoretical calculations. We demonstrate that the Si lattice hosts gapped Dirac cones at the Brillouin zone corners. Our tight-binding analysis reveals that the Dirac bands can be described by a 2D SSH model with anisotropic polarizations. The gap of the Dirac cone is driven by alternative hopping amplitudes in one direction and staggered potential energies in the other one and hosts topological edge states. Our results establish an ideal platform to explore the rich physical properties of the 2D SSH model. △ Less

Submitted 12 October, 2022; originally announced October 2022.

Comments: Accepted for publication in Nature Communications

Journal ref: Nat. Commun. 13, 7000 (2022)

CBLab: Supporting the Training of Large-scale Traffic Control Policies with Scalable Traffic Simulation

Authors: Chumeng Liang, Zherui Huang, Yicheng Liu, Zhanyu Liu, Guanjie Zheng, Hanyuan Shi, Kan Wu, Yuhao Du, Fuliang Li, Zhenhui Li

Abstract: Traffic simulation provides interactive data for the optimization of traffic control policies. However, existing traffic simulators are limited by their lack of scalability and shortage in input data, which prevents them from generating interactive data from traffic simulation in the scenarios of real large-scale city road networks. In this paper, we present \textbf{C}ity \textbf{B}rain \textbf{… ▽ More Traffic simulation provides interactive data for the optimization of traffic control policies. However, existing traffic simulators are limited by their lack of scalability and shortage in input data, which prevents them from generating interactive data from traffic simulation in the scenarios of real large-scale city road networks. In this paper, we present \textbf{C}ity \textbf{B}rain \textbf{Lab}, a toolkit for scalable traffic simulation. CBLab consists of three components: CBEngine, CBData, and CBScenario. CBEngine is a highly efficient simulator supporting large-scale traffic simulation. CBData includes a traffic dataset with road network data of 100 cities all around the world. We also develop a pipeline to conduct a one-click transformation from raw road networks to input data of our traffic simulation. Combining CBEngine and CBData allows researchers to run scalable traffic simulations in the road network of real large-scale cities. Based on that, CBScenario implements an interactive environment and a benchmark for two scenarios of traffic control policies respectively, with which traffic control policies adaptable for large-scale urban traffic can be trained and tuned. To the best of our knowledge, CBLab is the first infrastructure supporting traffic control policy optimization in large-scale urban scenarios. CBLab has supported the City Brain Challenge @ KDD CUP 2021. The project is available on GitHub:~\url{https://github.com/CityBrainLab/CityBrainLab.git}. △ Less

Submitted 4 June, 2023; v1 submitted 3 October, 2022; originally announced October 2022.

Comments: Accepted by KDD2023 (Applied Data Science Track)

Implementing quantum dimensionality reduction for non-Markovian stochastic simulation

Authors: Kang-Da Wu, Chengran Yang, Ren-Dong He, Mile Gu, Guo-Yong Xiang, Chuan-Feng Li, Guang-Can Guo, Thomas J. Elliott

Abstract: Complex systems are embedded in our everyday experience. Stochastic modelling enables us to understand and predict the behaviour of such systems, cementing its utility across the quantitative sciences. Accurate models of highly non-Markovian processes -- where the future behaviour depends on events that happened far in the past -- must track copious amounts of information about past observations,… ▽ More Complex systems are embedded in our everyday experience. Stochastic modelling enables us to understand and predict the behaviour of such systems, cementing its utility across the quantitative sciences. Accurate models of highly non-Markovian processes -- where the future behaviour depends on events that happened far in the past -- must track copious amounts of information about past observations, requiring high-dimensional memories. Quantum technologies can ameliorate this cost, allowing models of the same processes with lower memory dimension than corresponding classical models. Here we implement such memory-efficient quantum models for a family of non-Markovian processes using a photonic setup. We show that with a single qubit of memory our implemented quantum models can attain higher precision than possible with any classical model of the same memory dimension. This heralds a key step towards applying quantum technologies in complex systems modelling. △ Less

Submitted 18 October, 2023; v1 submitted 26 August, 2022; originally announced August 2022.

Comments: 11+2 pages, 5 figures

Journal ref: Nature Communications 14, 2624 (2023)

Positivity-Preserving Well-Balanced Central Discontinuous Galerkin Schemes for the Euler Equations under Gravitational Fields

Authors: Haili Jiang, Huazhong Tang, Kailiang Wu

Abstract: This paper designs and analyzes positivity-preserving well-balanced (WB) central discontinuous Galerkin (CDG) schemes for the Euler equations with gravity. A distinctive feature of these schemes is that they not only are WB for a general known stationary hydrostatic solution, but also can preserve the positivity of the fluid density and pressure. The standard CDG method does not possess this featu… ▽ More This paper designs and analyzes positivity-preserving well-balanced (WB) central discontinuous Galerkin (CDG) schemes for the Euler equations with gravity. A distinctive feature of these schemes is that they not only are WB for a general known stationary hydrostatic solution, but also can preserve the positivity of the fluid density and pressure. The standard CDG method does not possess this feature, while directly applying some existing WB techniques to the CDG framework may not accommodate the positivity and keep other important properties at the same time. In order to obtain the WB and positivity-preserving properties simultaneously while also maintaining the conservativeness and stability of the schemes, a novel spatial discretization is devised in the CDG framework based on suitable modifications to the numerical dissipation term and the source term approximation. The modifications are based on a crucial projection operator for the stationary hydrostatic solution, which is proposed for the first time in this work. This novel projection has the same order of accuracy as the standard $L^2$-projection, can be explicitly calculated, and is easy to implement without solving any optimization problems. More importantly, it ensures that the projected stationary solution has the same cell averages on both the primal and dual meshes, which is a key to achieve the desired properties of our schemes. Based on some convex decomposition techniques, rigorous positivity-preserving analyses for the resulting WB CDG schemes are carried out. Several one- and two-dimensional numerical examples are performed to illustrate the desired properties of these schemes, including the high-order accuracy, the WB property, the robustness for simulations involving the low pressure or density, high resolution for the discontinuous solutions and the small perturbations around the equilibrium state. △ Less

Submitted 19 July, 2022; originally announced July 2022.

Comments: 57 pages

A Physical-Constraint-Preserving Finite Volume WENO Method for Special Relativistic Hydrodynamics on Unstructured Meshes

Authors: Yaping Chen, Kailiang Wu

Abstract: This paper presents a highly robust third-order accurate finite volume weighted essentially non-oscillatory (WENO) method for special relativistic hydrodynamics on unstructured triangular meshes. We rigorously prove that the proposed method is physical-constraint-preserving (PCP), namely, always preserves the positivity of the pressure and the rest-mass density as well as the subluminal constraint… ▽ More This paper presents a highly robust third-order accurate finite volume weighted essentially non-oscillatory (WENO) method for special relativistic hydrodynamics on unstructured triangular meshes. We rigorously prove that the proposed method is physical-constraint-preserving (PCP), namely, always preserves the positivity of the pressure and the rest-mass density as well as the subluminal constraint on the fluid velocity. The method is built on a highly efficient compact WENO reconstruction on unstructured meshes, a simple PCP limiter, the provably PCP property of the Harten--Lax--van Leer flux, and third-order strong-stability-preserving time discretization. Due to the relativistic effects, the primitive variables (namely, the rest-mass density, velocity, and pressure) are highly nonlinear implicit functions in terms of the conservative variables, making the design and analysis of our method nontrivial. To address the difficulties arising from the strong nonlinearity, we adopt a novel quasilinear technique for the theoretical proof of the PCP property. Three provable convergence-guaranteed iterative algorithms are also introduced for the robust recovery of primitive quantities from admissible conservative variables. We also propose a slight modification to an existing WENO reconstruction to ensure the scaling invariance of the nonlinear weights and thus to accommodate the homogeneity of the evolution operator, leading to the advantages of the modified WENO reconstruction in resolving multi-scale wave structures. Extensive numerical examples are presented to demonstrate the robustness, expected accuracy, and high resolution of the proposed method. △ Less

Submitted 19 July, 2022; originally announced July 2022.

Comments: 56 pages, 18 figures

Journal ref: Journal of Computational Physics, 466: 111398, 2022

Free-space point-to-multiplepoint optical frequency transfer with lens assisted integrated beam steering

Authors: Liang Hu, Ruimin Xue, Xianyi Cao, Jiao Liu, Kan Wu, Guiling Wu, Jianping Chen

Abstract: We report on the realization of high-performance silica integrated two-dimensional lens assisted beam-steering (LABS) arrays along with the first-of-their-kind point-to-multiplepoint optical frequency transfer. {The LABS equips with $N$ antennas} and has the capability to produce arbitrary number of output beams with different output angles with the simple control complexity. We demonstrate that t… ▽ More We report on the realization of high-performance silica integrated two-dimensional lens assisted beam-steering (LABS) arrays along with the first-of-their-kind point-to-multiplepoint optical frequency transfer. {The LABS equips with $N$ antennas} and has the capability to produce arbitrary number of output beams with different output angles with the simple control complexity. We demonstrate that the LABS has 16 scanning angles, which can support {the access capability for the maximum of simultaneous 16 user nodes.} The coaxial configuration for transmitting and receiving the light as a monolithic transceiver allows us to reduce the out-of-loop phase noise significantly. Finally, the LABS-based non-blocking point-to-multiplepoint in-door free-space optical frequency transfer links with 24 m and 50 m free-space links are shown. After being compensated for the free-space link up to 50 m, the fractional frequency instability of $4.5\times10^{-17}$ and $7.7\times10^{-20}$ at the averaging time of 1 s and 20,000 s, respectively, can be achieved. The present work proves the potential application of the 2D LABS in free-space optical time-frequency transfer and provides a guidance for developing a chip-scale optical time-frequency transfer system. △ Less

Submitted 28 June, 2022; originally announced July 2022.

Comments: 10 pages, 6 figures, accepted published in IEEE Transactions on Instrumentation & Measurement

Filter-Bank-Enabled Leaky-Wave Antenna Array Technique for Full-Band-Locked Radar System in Stitched Frequency-Space Domain

Authors: Dongze Zheng, Ke Wu

Abstract: Inspired by the filter-bank (FB) concept that is normally used for multi-rate signal processing, an FB-enabled array technique of leaky-wave antennas (LWAs) is proposed and studied for creating full-band-locked frequency-scanning radar (FSR) systems in a stitched frequency-space domain. This is mainly for addressing the coupling dilemma between the range and angle resolutions, which is historicall… ▽ More Inspired by the filter-bank (FB) concept that is normally used for multi-rate signal processing, an FB-enabled array technique of leaky-wave antennas (LWAs) is proposed and studied for creating full-band-locked frequency-scanning radar (FSR) systems in a stitched frequency-space domain. This is mainly for addressing the coupling dilemma between the range and angle resolutions, which is historically and naturally inherited from a conventional FSR. First of all, the frequency-modulated continuous-wave system architecture is selected to exemplify and recall the characteristics of a conventional FSR with an emphasis on that resolution coupling. Then, a radar solution featuring a stitched frequency-space domain is introduced for the resolution decoupling, depending on an array of FB-enabled LWA channels. With the radar equation, FB-related conditions for realizing the critical frequency-space stitching are analytically derived and equivalently converted into several design specifications of an array of LWAs, i.e., engineered beam-scanning functions, beam-crossovers, and phase alignments. To facilitate the practical implementation of such an array technique, a detailed and generalized design flow is developed. Finally, for a simple proof of concept, an FB-enabled two-channel LWA array is modeled, fabricated, and measured. Simulated and measured results are in a reasonable agreement, and both demonstrate the desired frequency-space stitching behavior, i.e., enhanced spectrum bandwidth and widened radiation beamwidth. The proposed array solution may be potentially deployed for FSR systems with a full-band-locked beam illumination and a decoupled range-angle resolution, which may present a competitor against the phased array technique. △ Less

Submitted 15 June, 2022; originally announced June 2022.

Comments: This work will be submitted to the IEEE for possible publication. Copyright may be transferred without notice, after which this version may no longer be accessible

KSSOLV 2.0: An efficient MATLAB toolbox for solving the Kohn-Sham equations with plane-wave basis set

Authors: Shizhe Jiao, Zhenlin Zhang, Kai Wu, Lingyun Wan, Huanhuan Ma, Jielan Li, Sheng Chen, Xinming Qin, Jie Liu, Zijing Ding, Jinlong Yang, Yingzhou Li, Wei Hu, Lin Lin, Chao Yang

Abstract: KSSOLV (Kohn-Sham Solver) is a MATLAB toolbox for performing Kohn-Sham density functional theory (DFT) calculations with a plane-wave basis set. KSSOLV 2.0 preserves the design features of the original KSSOLV software to allow users and developers to easily set up a problem and perform ground-state calculations as well as to prototype and test new algorithms. Furthermore, it includes new functiona… ▽ More KSSOLV (Kohn-Sham Solver) is a MATLAB toolbox for performing Kohn-Sham density functional theory (DFT) calculations with a plane-wave basis set. KSSOLV 2.0 preserves the design features of the original KSSOLV software to allow users and developers to easily set up a problem and perform ground-state calculations as well as to prototype and test new algorithms. Furthermore, it includes new functionalities such as new iterative diagonalization algorithms, k-point sampling for electron band structures, geometry optimization and advanced algorithms for performing DFT calculations with local, semi-local, and hybrid exchange-correlation functionals. It can be used to study the electronic structures of both molecules and solids. We describe these new capabilities in this work through a few use cases. We also demonstrate the numerical accuracy and computational efficiency of KSSOLV on a variety of examples. △ Less

Submitted 6 August, 2022; v1 submitted 4 June, 2022; originally announced June 2022.

Comments: 25 pages, 11 figures, software release

Journal ref: Comp. Phys. Comm. 279, 108424 (2022)

Self-mixing in microtubule-kinesin active fluid from nonuniform to uniform distribution of activity

Authors: Teagan E Bate, Megan E Varney, Ezra H Taylor, Joshua H Dickie, Chih-Che Chueh, Michael M Norton, Kun-Ta Wu

Abstract: Active fluids have applications in micromixing, but little is known about the mixing kinematics of systems with spatiotemporally-varying activity. To investigate, UV-activated caged ATP was used to activate controlled regions of microtubule-kinesin active fluid and the mixing process was observed with fluorescent tracers and molecular dyes. At low Péclet numbers (diffusive transport), the active-i… ▽ More Active fluids have applications in micromixing, but little is known about the mixing kinematics of systems with spatiotemporally-varying activity. To investigate, UV-activated caged ATP was used to activate controlled regions of microtubule-kinesin active fluid and the mixing process was observed with fluorescent tracers and molecular dyes. At low Péclet numbers (diffusive transport), the active-inactive interface progressed toward the inactive area in a diffusion-like manner that was described by a simple model combining diffusion with Michaelis-Menten kinetics. At high Péclet numbers (convective transport), the active-inactive interface progressed in a superdiffusion-like manner that was qualitatively captured by an active-fluid hydrodynamic model coupled to ATP transport. Results showed that active fluid mixing involves complex coupling between distribution of active stress and active transport of ATP and reduces mixing time for suspended components with decreased impact of initial component distribution. This work will inform application of active fluids to promote micromixing in microfluidic devices. △ Less

Submitted 16 October, 2022; v1 submitted 24 May, 2022; originally announced May 2022.

Comments: 46 pages, 19 figures

Design and testing of LGAD sensor with shallow carbon implantation

Authors: Kewei Wu, Xuewei Jia, Tao Yang, Mengzhao Li, Wei Wang, Mei Zhao, Zhijun Liang, Joao Guimaraes da Costa, Yunyun Fan, Han Cui, Alissa Howard, Gregor Kramberger, Xin Shi, Yuekun Heng, Yuhang Tan, Bo Liu, Yuan Feng, Shuqi Li, Mengran Li, Chengjun Yu, Xuan Yang, Mingjie Zhai, Gaobo Xu, Gangping Yan, Qionghua Zhai , et al. (4 additional authors not shown)

Abstract: The low gain avalanche detectors (LGADs) are thin sensors with fast charge collection which in combination with internal gain deliver an outstanding time resolution of about 30 ps. High collision rates and consequent large particle rates crossing the detectors at the upgraded Large Hadron Collider (LHC) in 2028 will lead to radiation damage and deteriorated performance of the LGADs. The main conse… ▽ More The low gain avalanche detectors (LGADs) are thin sensors with fast charge collection which in combination with internal gain deliver an outstanding time resolution of about 30 ps. High collision rates and consequent large particle rates crossing the detectors at the upgraded Large Hadron Collider (LHC) in 2028 will lead to radiation damage and deteriorated performance of the LGADs. The main consequence of radiation damage is loss of gain layer doping (acceptor removal) which requires an increase of bias voltage to compensate for the loss of charge collection efficiency and consequently time resolution. The Institute of High Energy Physics (IHEP), Chinese Academy of Sciences (CAS) has developed a process based on the Institute of Microelectronics (IME), CAS capability to enrich the gain layer with carbon to reduce the acceptor removal effect by radiation. After 1 MeV neutron equivalent fluence of 2.5$\times$10$^{15}$ n$_{eq}$/cm$^{2}$, which is the maximum fluence to which sensors will be exposed at ATLAS High Granularity Timing Detector (HGTD), the IHEP-IME second version (IHEP-IMEv2) 50 $μ$m LGAD sensors already deliver adequate charge collection > 4 fC and time resolution < 50 ps at voltages < 400 V. The operation voltages of these 50 $μ$m devices are well below those at which single event burnout may occur. △ Less

Submitted 31 May, 2022; v1 submitted 10 May, 2022; originally announced May 2022.

Front Propagation from Radiative Sources

Authors: Theodore Steele, Kinwah Wu

Abstract: Fronts are regions of transition from one state to another in a medium. They are present in many areas of science and applied mathematics, and modelling them and their evolution is often an effective way of treating the underlying phenomena responsible for them. In this paper, we propose a new approach to modelling front propagation, which characterises the evolution of structures surrounding radi… ▽ More Fronts are regions of transition from one state to another in a medium. They are present in many areas of science and applied mathematics, and modelling them and their evolution is often an effective way of treating the underlying phenomena responsible for them. In this paper, we propose a new approach to modelling front propagation, which characterises the evolution of structures surrounding radiative sources. This approach is generic and has a wide range of applications, particularly when dealing with the propagation of phase or state transitions in media surrounding radiation emitting objects. As an illustration, we show an application in modelling the propagation of ionisation fronts around early stars during the cosmological Epoch of Reionisation (EoR) and show that the results are consistent with those of existing equations but provide much richer sources of information. △ Less

Submitted 9 May, 2022; originally announced May 2022.

Comments: 10 pages, 4 figures

Observation of biradical spin coupling through hydrogen bonds

Authors: Yang He, Na Li, Ivano E. Castelli, Ruoning Li, Yajie Zhang, Xue Zhang, Chao Li, Bingwu Wang, Song Gao, Lianmao Peng, Shimin Hou, Ziyong Shen, Jing-Tao Lü, Kai Wu, Per Hedegård, Yongfeng Wang

Abstract: Investigation of intermolecular electron spin interaction is of fundamental importance in both science and technology.Here, radical pairs of all-trans retinoic acid molecules on Au(111) are created using an ultra-low temperature scanning tunneling microscope. Antiferromagnetic coupling between two radicals is identified by magnetic-field dependent spectroscopy.The measured exchange energies are fr… ▽ More Investigation of intermolecular electron spin interaction is of fundamental importance in both science and technology.Here, radical pairs of all-trans retinoic acid molecules on Au(111) are created using an ultra-low temperature scanning tunneling microscope. Antiferromagnetic coupling between two radicals is identified by magnetic-field dependent spectroscopy.The measured exchange energies are from 0.1 to 1.0 meV. The biradical spin coupling is mediated through O-H$\cdots$O hydrogen bonds, as elucidated from analysis combining density functional theory calculation and a modern version of valence bond theory. △ Less

Submitted 4 May, 2022; originally announced May 2022.

Comments: 10 pages with supplemental materials, PRL accepted