Showing 1–50 of 92 results for author: Lai, C

SubsurfaceGen: Procedural Generation of Field-Scale Earth Models and Seismic Data

Authors: Joseph Stitt, Pratik Rathore, Madeleine Udell, Ching-Yao Lai

Abstract: Full waveform inversion (FWI) is the gold standard for subsurface imaging, with applications from carbon sequestration to energy and mineral exploration to earthquake hazard assessment. Machine learning approaches to FWI need field-scale, geologically diverse, and physically realistic training data, but existing resources such as Marmousi, SEAM, and OpenFWI fall short on spatial extent, temporal e… ▽ More Full waveform inversion (FWI) is the gold standard for subsurface imaging, with applications from carbon sequestration to energy and mineral exploration to earthquake hazard assessment. Machine learning approaches to FWI need field-scale, geologically diverse, and physically realistic training data, but existing resources such as Marmousi, SEAM, and OpenFWI fall short on spatial extent, temporal extent, geological diversity, and physical realism. We address these limitations with SubsurfaceGen, a GPU-accelerated generator for 3D velocity models and seismic data. Along with SubsurfaceGen, we release a paired dataset of 4,276 2D velocity slices, 5 s wavefields, and 8 s shot gathers drawn from 42 realistic, field-scale 3D velocity models, each spanning 10 km x 10 km laterally and 6.19 km deep at 10 m resolution. The dataset spans six geological settings -- four built with SubsurfaceGen and two drawn from prior sources -- relevant for carbon sequestration and hydrocarbon exploration. We use this dataset to evaluate neural operators on wavefield prediction and encoder-decoders on end-to-end velocity inversion, holding out one geological setting for out-of-distribution testing. These experiments surface failure modes at field-scale and demonstrate how SubsurfaceGen and the associated dataset can impact ML-based FWI. △ Less

Submitted 28 May, 2026; originally announced May 2026.

Comments: 38 pages

Generative AI impacts on intra-urban inequality and skill premium in Beijing

Authors: Xiliu He, Haoxiang Zhao, Mingyi Ma, Edward Wen Chuan Lai, Koei Enomoto, Anni Hu, Jiatong Li, Lingyun Chu, Yuan Lai

Abstract: Generative artificial intelligence (GenAI) is the first automation wave to reach high-cognitive tasks at scale, yet its effects on intra-urban inequality remain largely unknown. Using 5 million job postings from Beijing (2018--2024), we construct a neighborhood-level GenAI Exposure Index by aggregating task-level assessments from five leading large language models. We examine the spatial, structur… ▽ More Generative artificial intelligence (GenAI) is the first automation wave to reach high-cognitive tasks at scale, yet its effects on intra-urban inequality remain largely unknown. Using 5 million job postings from Beijing (2018--2024), we construct a neighborhood-level GenAI Exposure Index by aggregating task-level assessments from five leading large language models. We examine the spatial, structural and causal mechanisms of this shock. We find that GenAI exposure is highly concentrated in the city's core districts, deepening the intra-urban AI divide. Since 2023, high-exposure neighborhoods have experienced wage stagnation even as they continue to attract high-skilled workers -- a "high-skill trap." This wage penalty is driven by task de-skilling and intensified labor-market crowding. A difference-in-differences design centered on ChatGPT's release supports a causal interpretation. These findings challenge the prevailing theory of skill-biased technological change and provide a basis for inclusive AI governance in global technology hubs. △ Less

Submitted 25 May, 2026; originally announced May 2026.

Comments: 21 pages, 8 figures

Alpha Background in Multi-Grid Neutron Detectors

Authors: A. Backis, C. -C. Lai, J. R. M. Annand, K. G. Fissum, G. Zuzel, M. Czubak, K. Livingston

Abstract: Alpha emission from actinide impurities in Al is a source of background counting rate in Multi-Grid type detectors of thermal neutrons. The alpha emission rates from samples of radio-purity Al and \mathrm{Al/B_{4}C} composite, used in grid construction, were measured on a large-area, low background spectrometer. Although the alpha emission rate from the composite was a factor \sim280 higher than r… ▽ More Alpha emission from actinide impurities in Al is a source of background counting rate in Multi-Grid type detectors of thermal neutrons. The alpha emission rates from samples of radio-purity Al and \mathrm{Al/B_{4}C} composite, used in grid construction, were measured on a large-area, low background spectrometer. Although the alpha emission rate from the composite was a factor \sim280 higher than radio-pure Al, \mathrm{25\:μm} Ni plating of the composite reduced the rate by a factor \sim1170. Background counting rates in two Multi-Grid prototypes were compared. They used identical configurations of \mathrm{B_{4}C}-coated, radio-pure Al normal blades for the grids, but the first employed radio-purity Al for the radial blades, while the second used Ni-plated \mathrm{Al/B_{4}C} on the radial blades. The background rate from the second prototype was around 20% of that from the first. △ Less

Submitted 21 May, 2026; originally announced May 2026.

Comments: 12 figures and text

Graduate Training in Quantum Information Science and Engineering: Lessons, Challenges, and a Roadmap from the NSF Research Traineeship Programs

Authors: Yohannes Abate, Victor Acosta, Alessandro Alabastri, Mehmet Aydeniz, Viktoriia E. Babicheva, Lincoln D. Carr, I-Tung Chen, Wandi Ding, Tara Drake, Mattias Fitzpatrick, Kai-Mei C. Fu, Jay Gupta, Kaden R. A. Hazzard, Sophia E. Hayes, Jin Hu, Hilary M. Hurst, Sohrab Ismail-Beigi, Ehsan Khatami, Junichiro Kono, Cheng-Yu Lai, Xiuling Li, Yingmei Liu, Sara Mouradian, Kater Murch, Borja Peropadre , et al. (9 additional authors not shown)

Abstract: Since 2019, eighteen NSF Research Traineeship (NRT) awards in quantum information science and engineering (QISE) and adjacent fields have been funded, constituting the largest NSF-coordinated investment in graduate QISE training in the United States. Synthesizing lessons from our programs, we work through the central tensions that every QISE graduate program must negotiate: between depth in a home… ▽ More Since 2019, eighteen NSF Research Traineeship (NRT) awards in quantum information science and engineering (QISE) and adjacent fields have been funded, constituting the largest NSF-coordinated investment in graduate QISE training in the United States. Synthesizing lessons from our programs, we work through the central tensions that every QISE graduate program must negotiate: between depth in a home discipline and breadth across the field, between structured instruction and open-ended experiential and hands-on learning, and between training individual specialists and cultivating teams that collectively cover all areas of QISE. We describe the structural and pedagogical innovations the NRT programs have developed in response, assess what is working and what remains unresolved, and sketch 12 open problems the community will need to address as QISE graduate education scales beyond the well-resourced research universities where it has up till now been mainly concentrated. Eight concrete recommendations follow: (1) adopt the startup model of team-based training as an organizing philosophy; (2) invest immediately in sensing and communication curriculum development; (3) build student agency into program governance, not just activities; (4) establish structural mechanisms for industrial engagement rather than depending on goodwill; (5) design for sustainability from year one; (6) develop graduate-level textbooks spanning all three QISE pillars: computing, sensing, and communications; (7) establish shared outcome assessment instruments across programs; and (8) develop structured mechanisms for faculty professional development in QISE. △ Less

Submitted 8 May, 2026; originally announced May 2026.

Comments: 49 pages, 3 figures, 1 table

Transferable Physics-Informed Representations via Closed-Form Head Adaptation

Authors: Jian Cheng Wong, Isaac Yin Chung Lai, Pao-Hsiung Chiu, Chin Chun Ooi, Abhishek Gupta, Yew-Soon Ong

Abstract: Physics-informed neural networks (PINNs) have garnered significant interest for their potential in solving partial differential equations (PDEs) that govern a wide range of physical phenomena. By incorporating physical laws into the learning process, PINN models have demonstrated the ability to learn physical outcomes reasonably well. However, current PINN approaches struggle to predict or solve n… ▽ More Physics-informed neural networks (PINNs) have garnered significant interest for their potential in solving partial differential equations (PDEs) that govern a wide range of physical phenomena. By incorporating physical laws into the learning process, PINN models have demonstrated the ability to learn physical outcomes reasonably well. However, current PINN approaches struggle to predict or solve new PDEs effectively when there is a lack of training examples, indicating they do not generalize well to unseen problem instances. In this paper, we present a transferable learning approach for PINNs premised on a fast Pseudoinverse PINN framework (Pi-PINN). Pi-PINN learns a transferable physics-informed representation in a shared embedding space and enables rapid solving of both known and unknown PDE instances via closed-form head adaptation using a least-squares-optimal pseudoinverse under PDE constraints. We further investigate the synergies between data-driven multi-task learning loss and physics-informed loss, providing insights into the design of more performant PINNs. We demonstrate the effectiveness of Pi-PINN on various PDE problems, including Poisson's equation, Helmholtz equation, and Burgers' equation, achieving fast and accurate physics-informed solutions without requiring any data for unseen instances. Pi-PINN can produce predictions 100-1000 times faster than a typical PINN, while producing predictions with 10-100 times lower relative error than a typical data-driven model even with only two training samples. Overall, our findings highlight the potential of transferable representations with closed-form head adaptation to enhance the efficiency and generalization of PINNs across PDE families and scientific and engineering applications. △ Less

Submitted 23 April, 2026; originally announced April 2026.

Comments: Accepted at IJCNN 2026

Brillouin-Enhanced Photonic Stepped-Frequency Radar

Authors: Ziqian Zhang, Ryan L. Russell, Choon Kong Lai, Benjamin J. Eggleton

Abstract: Photonic stepped-frequency (SF) radar offers high range resolution and only requires low-speed driving electronics, but existing architectures face challenges in achieving low phase noise and uniform frequency steps simultaneously. Here, we demonstrate a photonic SF radar system that exploits dual Brillouin lasers in a shared fiber cavity to simultaneously suppress phase noise and ensure uniform f… ▽ More Photonic stepped-frequency (SF) radar offers high range resolution and only requires low-speed driving electronics, but existing architectures face challenges in achieving low phase noise and uniform frequency steps simultaneously. Here, we demonstrate a photonic SF radar system that exploits dual Brillouin lasers in a shared fiber cavity to simultaneously suppress phase noise and ensure uniform frequency stepping. Phase noise is reduced through Brillouin optomechanical suppression and common-mode noise rejection upon photomixing. Frequency-step uniformity is enforced via lasing at a series of uniformly spaced cavity resonances. The system generates an X-band SF waveform spanning 1.31 GHz, achieving >23 dB of phase-noise improvement at a 100 kHz offset relative to a low-cost driving voltage-controlled oscillator. The demonstrated system reduces the dependence of the output waveform quality on noise in the driving electronics, offering a path towards high-performance radar sensing. △ Less

Submitted 21 April, 2026; originally announced April 2026.

The SIREN Program: A Scalable Model for Short-Term Undergraduate Research Experiences at Community Colleges

Authors: Emilie Hein, Polin Yadak, Denise Hum, Jessica Hurless, Luis Jibaja Prado, Susanne Schubert, Pia Walawalkar, Marco Wehrfritz, Daria Baitazarova, Yuting Duan, Tin Htoo, Codie Lai, Aryanna Mendoza, Eslin Villalta, Valeria Zarco

Abstract: Providing meaningful research experiences for undergraduate students is a well-recognized challenge, particularly at community colleges and teaching-focused institutions where resources are limited and faculty time is dedicated to instruction. To address this, the Summer Introduction to Research and Experimentation in Nuclear physics (SIREN) was developed as a three-week, intensive summer program… ▽ More Providing meaningful research experiences for undergraduate students is a well-recognized challenge, particularly at community colleges and teaching-focused institutions where resources are limited and faculty time is dedicated to instruction. To address this, the Summer Introduction to Research and Experimentation in Nuclear physics (SIREN) was developed as a three-week, intensive summer program that engaged 24 students working in six teams on original research projects centered around cosmic detection. Each team was supported by an advisor and a peer mentor, combining structured guidance with near-peer learning. After developing a range of technical skills, students built cosmic detectors and used them to explore a variety of projects. The program also included multiple workshops on research skills and seminars led by guest speakers, giving students access to broader STEM pathways. The evaluation of the program was based on advisor and student surveys, as well as peer mentor observations. The findings suggest that the short-term team-based model promoted engagement, collaboration, and skill development while leading to the completion of advanced projects. The advisors highlighted the effectiveness of the program structure and the excitement generated by the direct application of newly acquired skills to research. Data showed that students were able to significantly develop transferable skills through hands-on activities. Post-program surveys indicated that participants experienced an increased sense of belonging in their majors and greater confidence in pursuing careers in STEM. △ Less

Submitted 15 December, 2025; originally announced December 2025.

Resolving Sharp Gradients of Unstable Singularities to Machine Precision via Neural Networks

Authors: Yongji Wang, Tristan Léger, Ching-Yao Lai, Tristan Buckmaster

Abstract: Recent work introduced a robust computational framework combining embedded mathematical structures, advanced optimization, and neural network architecture, leading to the discovery of multiple unstable self-similar solutions for key fluid dynamics equations, including the Incompressible Porous Media (IPM) and 2D Boussinesq systems. While this framework confirmed the existence of these singularitie… ▽ More Recent work introduced a robust computational framework combining embedded mathematical structures, advanced optimization, and neural network architecture, leading to the discovery of multiple unstable self-similar solutions for key fluid dynamics equations, including the Incompressible Porous Media (IPM) and 2D Boussinesq systems. While this framework confirmed the existence of these singularities, an accuracy level approaching double-float machine precision was only achieved for stable and 1st unstable solutions of the 1D Córdoba-Córdoba-Fontelos model. For highly unstable solutions characterized by extreme gradients, the accuracy remained insufficient for validation. The primary obstacle is the presence of sharp solution gradients. Those gradients tend to induce large, localized PDE residuals during training, which not only hinder convergence, but also obscure the subtle signals near the origin required to identify the correct self-similar scaling parameter lambda of the solutions. In this work, we introduce a gradient-normalized PDE residual re-weighting scheme to resolve the high-gradient challenge while amplifying the critical residual signals at the origin for lambda identification. Coupled with the multi-stage neural network architecture, the PDE residuals are reduced to the level of round-off error across a wide spectrum of unstable self-similar singularities previously discovered. Furthermore, our method enables the discovery of new highly unstable singularities, i.e. the 4th unstable solution for IPM equations and a novel family of highly unstable solitons for the Nonlinear Schrödinger equations. This results in achieving high-gradient solutions with high precision, providing an important ingredient for bridging the gap between numerical discovery and computer-assisted proofs for unstable phenomena in nonlinear PDEs. △ Less

Submitted 27 November, 2025; originally announced November 2025.

Comments: 27 pages, 12 figures

Neutron Interaction Properties of Structural Materials for Multi-Grid Neutron Detectors

Authors: A. Backis, C. -C. Lai, M. Aouane, P. P. Deen, K. G. Fissum, J. R. M. Annand, K. Livingston, D. Raspino

Abstract: The T-REX neutron time-of-flight spectrometer at the European Spallation Source will use Multi-Grid Technology, which relies on thin B4C coatings on the Al blades of the grids to detect scattered thermal neutrons. Following a Monte Carlo study of internal shielding to suppress neutron multiple scattering in T-REX, the neutron transmission and scattering properties of 12 shielding-material samples… ▽ More The T-REX neutron time-of-flight spectrometer at the European Spallation Source will use Multi-Grid Technology, which relies on thin B4C coatings on the Al blades of the grids to detect scattered thermal neutrons. Following a Monte Carlo study of internal shielding to suppress neutron multiple scattering in T-REX, the neutron transmission and scattering properties of 12 shielding-material samples have been measured at the ISIS spallation neutron source. Neutron transmission was measured on the EMMA beam line at wavelengths 0.5-4.7 A, using a 2D-position-sensitive, neutron GEM detector, while neutron scattering was measured for 6 of the samples at the Merlin spectrometer, at wavelengths 0.72, 1.28, 1.85 and 2.41 A. The present tests show that a B4C/Al composite material, plated with Ni to stop intrinsic alpha background, is an effective neutron absorber, suitable for incorporation in the Multi-Grid structures which detect the neutrons in inelastic neutron spectrometers . △ Less

Submitted 5 November, 2025; v1 submitted 4 November, 2025; originally announced November 2025.

Discovery of Unstable Singularities

Authors: Yongji Wang, Mehdi Bennani, James Martens, Sébastien Racanière, Sam Blackwell, Alex Matthews, Stanislav Nikolov, Gonzalo Cao-Labora, Daniel S. Park, Martin Arjovsky, Daniel Worrall, Chongli Qin, Ferran Alet, Borislav Kozlovskii, Nenad Tomašev, Alex Davies, Pushmeet Kohli, Tristan Buckmaster, Bogdan Georgiev, Javier Gómez-Serrano, Ray Jiang, Ching-Yao Lai

Abstract: Whether singularities can form in fluids remains a foundational unanswered question in mathematics. This phenomenon occurs when solutions to governing equations, such as the 3D Euler equations, develop infinite gradients from smooth initial conditions. Historically, numerical approaches have primarily identified stable singularities. However, these are not expected to exist for key open problems,… ▽ More Whether singularities can form in fluids remains a foundational unanswered question in mathematics. This phenomenon occurs when solutions to governing equations, such as the 3D Euler equations, develop infinite gradients from smooth initial conditions. Historically, numerical approaches have primarily identified stable singularities. However, these are not expected to exist for key open problems, such as the boundary-free Euler and Navier-Stokes cases, where unstable singularities are hypothesized to play a crucial role. Here, we present the first systematic discovery of new families of unstable singularities. A stable singularity is a robust outcome, forming even if the initial state is slightly perturbed. In contrast, unstable singularities are exceptionally elusive; they require initial conditions tuned with infinite precision, being in a state of instability whereby infinitesimal perturbations immediately divert the solution from its blow-up trajectory. In particular, we present multiple new, unstable self-similar solutions for the incompressible porous media equation and the 3D Euler equation with boundary, revealing a simple empirical asymptotic formula relating the blow-up rate to the order of instability. Our approach combines curated machine learning architectures and training schemes with a high-precision Gauss-Newton optimizer, achieving accuracies that significantly surpass previous work across all discovered solutions. For specific solutions, we reach near double-float machine precision, attaining a level of accuracy constrained only by the round-off errors of the GPU hardware. This level of precision meets the requirements for rigorous mathematical validation via computer-assisted proofs. This work provides a new playbook for exploring the complex landscape of nonlinear partial differential equations (PDEs) and tackling long-standing challenges in mathematical physics. △ Less

Submitted 17 September, 2025; originally announced September 2025.

Comments: 20 pages, 6 figures. Supplementary information will be uploaded in a forthcoming version of the manuscript

Moment Estimates and DeepRitz Methods on Learning Diffusion Systems with Non-gradient Drifts

Authors: Fanze Kong, Chen-Chih Lai, Yubin Lu

Abstract: Conservative-dissipative dynamics are ubiquitous across a variety of complex open systems. We propose a data-driven two-phase method, the Moment-DeepRitz Method, for learning drift decompositions in generalized diffusion systems involving conservative-dissipative dynamics. The method is robust to noisy data, adaptable to rough potentials and oscillatory rotations. We demonstrate its effectiveness… ▽ More Conservative-dissipative dynamics are ubiquitous across a variety of complex open systems. We propose a data-driven two-phase method, the Moment-DeepRitz Method, for learning drift decompositions in generalized diffusion systems involving conservative-dissipative dynamics. The method is robust to noisy data, adaptable to rough potentials and oscillatory rotations. We demonstrate its effectiveness through several numerical experiments. △ Less

Submitted 31 August, 2025; originally announced September 2025.

Moment Estimate and Variational Approach for Learning Generalized Diffusion with Non-gradient Structures

Authors: Fanze Kong, Chen-Chih Lai, Yubin Lu

Abstract: This paper proposes a data-driven learning framework for identifying governing laws of generalized diffusions with non-gradient components. By combining energy dissipation laws with a physically consistent penalty and first-moment evolution, we design a two-stage method to recover the pseudo-potential and rotation in the pointwise orthogonal decomposition of a class of non-gradient drifts in gener… ▽ More This paper proposes a data-driven learning framework for identifying governing laws of generalized diffusions with non-gradient components. By combining energy dissipation laws with a physically consistent penalty and first-moment evolution, we design a two-stage method to recover the pseudo-potential and rotation in the pointwise orthogonal decomposition of a class of non-gradient drifts in generalized diffusions. Our two-stage method is applied to complex generalized diffusion processes including dissipation-rotation dynamics, rough pseudo-potentials and noisy data. Representative numerical experiments demonstrate the effectiveness of our approach for learning physical laws in non-gradient generalized diffusions. △ Less

Submitted 5 August, 2025; v1 submitted 3 August, 2025; originally announced August 2025.

Unraveling the Molecular Structure of Lipid Nanoparticles through in-silico Self-Assembly for Rational Delivery Design

Authors: Xuan Bai, Yu Lu, Tianhao Yu, Kangjie Lv, Cai Yao, Feng Shi, Andong Liu, Kai Wang, Wenshou Wang, Chris Lai

Abstract: Lipid nanoparticles (LNPs) are a leading platform in the delivery of RNA-based therapeutics, playing a pivotal role in the clinical success of mRNA vaccines and other nucleic acid drugs. Their performance in RNA encapsulation and delivery is critically governed by the molecular structure of ionizable lipids and the overall formulation composition. However, mechanistic insight into how these factor… ▽ More Lipid nanoparticles (LNPs) are a leading platform in the delivery of RNA-based therapeutics, playing a pivotal role in the clinical success of mRNA vaccines and other nucleic acid drugs. Their performance in RNA encapsulation and delivery is critically governed by the molecular structure of ionizable lipids and the overall formulation composition. However, mechanistic insight into how these factors govern LNP architecture and function remains limited, primarily owing to the challenges of capturing nanoscale assembly and organization using experimental techniques. Here, we employ coarse-grained molecular dynamics simulations to systematically investigate how ionizable lipid chemistry influences LNP self-assembly, internal organization, and surface properties. We further explore the effects of formulation ratios and pH-dependent deprotonation on both the internal structure and surface morphology of LNPs. Leveraging these insights, we demonstrate how in silico structural characteristics can inform the rational design of novel ionizable lipids and optimization of formulation ratios, supported with experimental validations. Our findings offer a molecular-level understanding of LNP assembly dynamics and architecture, thereby establishing a computational framework linking lipid chemistry and LNP formulation to the structure and performance of LNP, to advance the rational design of novel LNP delivery systems. △ Less

Submitted 3 August, 2025; originally announced August 2025.

Nonspherically symmetric equilibrium bubbles in a steadily rotating incompressible fluid

Authors: Chen-Chih Lai, Michael I. Weinstein

Abstract: This note presents two nontrivial, rotational equilibrium solutions to the spatial uniform gas pressure (isobaric) approximate model of Prosperetti in the inviscid case. Building on Gavrilov's work [GAFA 2019], we first establish the existence of equilibrium solutions with nontrivial (rotational) liquid flow. Second, we construct a nonspherically symmetric, horn-torus-shaped equilibrium bubble und… ▽ More This note presents two nontrivial, rotational equilibrium solutions to the spatial uniform gas pressure (isobaric) approximate model of Prosperetti in the inviscid case. Building on Gavrilov's work [GAFA 2019], we first establish the existence of equilibrium solutions with nontrivial (rotational) liquid flow. Second, we construct a nonspherically symmetric, horn-torus-shaped equilibrium bubble under mild spatial decay conditions of the liquid flow. In addition, we extend earlier results on the characterization of spherical equilibrium bubbles to the axisymmetric, purely azimuthal setting. Finally, we implement a numerical simulation of the equilibrium bubble shape using the Physics-Informed Neural Network (PINN) approximation. △ Less

Submitted 23 July, 2025; originally announced July 2025.

Comments: 12 pages, 1 figure

Brillouin lasers in Bragg grating microresonators

Authors: Ryan L. Russell, Moritz Merklein, Choon Kong Lai, Cong Tinh Bui, Alvaro Casas-Bedoya, Duk-Yong Choi, Stephen J. Madden, Benjamin J. Eggleton

Abstract: Chip-scale coherent light sources are required in applications spanning metrology and sensing to telecommunications. Brillouin lasers (BLs) offer a route to ultra-coherent optical sources in compact microresonators with free spectral range (FSR) matched to the Brillouin frequency shift (BFS). However, BFS - FSR matching typically facilitates cascaded Brillouin scattering, constraining achievable B… ▽ More Chip-scale coherent light sources are required in applications spanning metrology and sensing to telecommunications. Brillouin lasers (BLs) offer a route to ultra-coherent optical sources in compact microresonators with free spectral range (FSR) matched to the Brillouin frequency shift (BFS). However, BFS - FSR matching typically facilitates cascaded Brillouin scattering, constraining achievable BL output power and coherence. Here, we demonstrate inhibition of cascading in a planar-integrated chalcogenide microresonator by exploiting the photonic bandgap (PBG) associated with a post-fabrication inscribed, reconfigurable intracavity Bragg grating. The PBG inhibits energy transfer within the target Brillouin scattering pathway, such as from pump to first-order Stokes wave. As a quantitative measure of Brillouin scattering inhibition, we report at least six-fold increase in threshold for onset of BL oscillation, which is ultimately limited by thermorefraction. For on-chip pump power of 399 mW, sufficient for a tenth-order Brillouin cascade, complete inhibition was achieved. Our work positions Bragg grating microresonators as an enabling platform for high performance on-chip BL sources, with reconfigurable modes of operation. △ Less

Submitted 4 June, 2025; originally announced June 2025.

Multiple Descents in Deep Learning as a Sequence of Order-Chaos Transitions in LSTM Networks

Authors: Wenbo Wei, Fan Xu, Nicholas Chong Jia Le, Choy Heng Lai, Ling Feng

Abstract: We observe a novel `multiple-descent' phenomenon during the learning process of a recurrent neural network called long-short-term memory (LSTM) networks during its training on real-world task, in which the performance goes through long cycles of up and down trends multiple times after the model is overtrained. By carrying out asymptotic stability analysis of the models, we found that the cycles in… ▽ More We observe a novel `multiple-descent' phenomenon during the learning process of a recurrent neural network called long-short-term memory (LSTM) networks during its training on real-world task, in which the performance goes through long cycles of up and down trends multiple times after the model is overtrained. By carrying out asymptotic stability analysis of the models, we found that the cycles in performance -- indicated by loss function in test data -- are closely associated with the phase transition process between order and chaos of the model, and the local optimal training step are consistently at the critical transition point between the two phases. More importantly, the most optimal point of the model usually occurs at the first transition from order to chaos, where the `width' of the `edge of chaos' is often the widest, allowing the best exploration of weight configurations for learning. △ Less

Submitted 15 June, 2026; v1 submitted 26 May, 2025; originally announced May 2025.

The First Hardware Demonstration of a Universal Programmable RRAM-based Probabilistic Computer for Molecular Docking

Authors: Yihan He, Ming-Chun Hong, Qiming Ding, Chih-Sheng Lin, Chih-Ming Lai, Chao Fang, Xiao Gong, Tuo-Hung Hou, Gengchiau Liang

Abstract: Molecular docking is a critical computational strategy in drug design and discovery, but the complex diversity of biomolecular structures and flexible binding conformations create an enormous search space that challenges conventional computing methods. Although quantum computing holds promise for these challenges, it remains constrained by scalability, hardware limitations, and precision issues. H… ▽ More Molecular docking is a critical computational strategy in drug design and discovery, but the complex diversity of biomolecular structures and flexible binding conformations create an enormous search space that challenges conventional computing methods. Although quantum computing holds promise for these challenges, it remains constrained by scalability, hardware limitations, and precision issues. Here, we report a prototype of a probabilistic computer (p-computer) that efficiently and accurately solves complex molecular docking for the first time, overcoming previously encountered challenges. At the core of the system is a p-computing chip based upon our artificial tunable probabilistic bits (p-bits), which are compatible with computing-in-memory schemes, based upon 180 nm CMOS technology and BEOL HfO2 RRAM. We successfully demonstrated the superior performance of the p-computer in practical ligand-protein docking scenarios. A 42-node molecular docking problem of lipoprotein with LolA-LolCDE complex-a key point in developing antibiotics against Gram-negative bacteria, was successfully solved. Our results align well with the Protein-Ligand Interaction Profiler tool. This work marks the first application of p-computing in molecular docking-based computational biology, which has great potential to overcome the limitations in success rate and efficiency of current technologies in addressing complex bioinformatics problems. △ Less

Submitted 26 March, 2025; originally announced March 2025.

Comments: The main text comprises 24 pages with 5 figures. The supplementary information includes 9 pages, containing Supplementary Figures S1 to S3 and Supplementary Tables S1 to S9

Journal ref: Nat Commun 17, 702 (2026)

Scaffold-Assisted Window Junctions for Superconducting Qubit Fabrication

Authors: Chung-Ting Ke, Jun-Yi Tsai, Yen-Chun Chen, Zhen-Wei Xu, Elam Blackwell, Matthew A. Snyder, Spencer Weeden, Peng-Sheng Chen, Chih-Ming Lai, Shyh-Shyuan Sheu, Zihao Yang, Cen-Shawn Wu, Alan Ho, R. McDermott, John Martinis, Chii-Dong Chen

Abstract: The superconducting qubit is one of the promising directions in realizing fault-tolerant quantum computing (FTQC), which requires many high-quality qubits. To achieve this, it is desirable to leverage modern semiconductor industry technology to ensure quality, uniformity, and reproducibility. However, conventional Josephson junction fabrication relies mainly on resist-assistant double-angle evapor… ▽ More The superconducting qubit is one of the promising directions in realizing fault-tolerant quantum computing (FTQC), which requires many high-quality qubits. To achieve this, it is desirable to leverage modern semiconductor industry technology to ensure quality, uniformity, and reproducibility. However, conventional Josephson junction fabrication relies mainly on resist-assistant double-angle evaporation, posing integration challenges. Here, we demonstrate a lift-off-free qubit fabrication that integrates seamlessly with existing industrial technologies. This method employs a silicon oxide (SiO$_2$) scaffold to define an etched window with a well-controlled size to form a Josephson junction. The SiO$_2$, which has a large dielectric loss, is etched away in the final step using vapor HF leaving little residue. This Window junction (WJ) process mitigates the degradation of qubit quality during fabrication and allows clean removal of the scaffold. The WJ process is validated by inspection and Josephson junction measurement. The scaffold removal process is verified by measuring the quality factor of the resonators. Furthermore, compared to scaffolds fabricated by plasma-enhanced chemical vapor deposition (PECVD), qubits made by WJ through physical vapor deposition (PVD) achieve relaxation time up to $57\,μ\text{s}$. Our results pave the way for a lift-off-free qubit fabrication process, designed to be compatible with modern foundry tools and capable of minimizing damage to the substrate and material surfaces. △ Less

Submitted 13 March, 2025; originally announced March 2025.

Isotropic Metamaterial Stiffness Beyond Hashin-Shtrikman Upper Bound

Authors: Manish Kumar Singh, Chang Quan Lai

Abstract: Since its introduction more than 60 years ago, the Hashin-Shtrikman upper bound has stood as the theoretical limit for the stiffness of isotropic composites and porous solids, acting as an important reference against which the moduli of heterogeneous structural materials are assessed. Here, we show through first-principles calculations, supported by finite element simulations, that the Hashin-Shtr… ▽ More Since its introduction more than 60 years ago, the Hashin-Shtrikman upper bound has stood as the theoretical limit for the stiffness of isotropic composites and porous solids, acting as an important reference against which the moduli of heterogeneous structural materials are assessed. Here, we show through first-principles calculations, supported by finite element simulations, that the Hashin-Shtrikman upper bound can be exceeded by the isotropic elastic response of an anisotropic structure constructed from an anisotropic material. The material and structural anisotropies mutually reinforce each other to realize the overall isotropic response, without incurring the mass penalty faced by the hybridization of geometries with complementary anisotropies. 3 designs were investigated (plate BCC, plate FCC and plate SC) but only plate SC yielded a solution for the anisotropic properties of the material, which are remarkably similar to that of single crystal nickel and single crystal ferrite. △ Less

Submitted 18 November, 2024; originally announced November 2024.

Comments: Supplementary Information available at the end of the article

LipidBERT: A Lipid Language Model Pre-trained on METiS de novo Lipid Library

Authors: Tianhao Yu, Cai Yao, Zhuorui Sun, Feng Shi, Lin Zhang, Kangjie Lyu, Xuan Bai, Andong Liu, Xicheng Zhang, Jiali Zou, Wenshou Wang, Chris Lai, Kai Wang

Abstract: In this study, we generate and maintain a database of 10 million virtual lipids through METiS's in-house de novo lipid generation algorithms and lipid virtual screening techniques. These virtual lipids serve as a corpus for pre-training, lipid representation learning, and downstream task knowledge transfer, culminating in state-of-the-art LNP property prediction performance. We propose LipidBERT,… ▽ More In this study, we generate and maintain a database of 10 million virtual lipids through METiS's in-house de novo lipid generation algorithms and lipid virtual screening techniques. These virtual lipids serve as a corpus for pre-training, lipid representation learning, and downstream task knowledge transfer, culminating in state-of-the-art LNP property prediction performance. We propose LipidBERT, a BERT-like model pre-trained with the Masked Language Model (MLM) and various secondary tasks. Additionally, we compare the performance of embeddings generated by LipidBERT and PhatGPT, our GPT-like lipid generation model, on downstream tasks. The proposed bilingual LipidBERT model operates in two languages: the language of ionizable lipid pre-training, using in-house dry-lab lipid structures, and the language of LNP fine-tuning, utilizing in-house LNP wet-lab data. This dual capability positions LipidBERT as a key AI-based filter for future screening tasks, including new versions of METiS de novo lipid libraries and, more importantly, candidates for in vivo testing for orgran-targeting LNPs. To the best of our knowledge, this is the first successful demonstration of the capability of a pre-trained language model on virtual lipids and its effectiveness in downstream tasks using web-lab data. This work showcases the clever utilization of METiS's in-house de novo lipid library as well as the power of dry-wet lab integration. △ Less

Submitted 3 May, 2025; v1 submitted 12 August, 2024; originally announced August 2024.

Asymmetric deformations of a perturbed spherical bubble in an incompressible fluid

Authors: Chen-Chih Lai, Michael I. Weinstein

Abstract: We study the dynamics of a gas bubble in a fluid with surface tension, initially near a spherical equilibrium. While there are many studies and applications of radial bubble dynamics, the theory of general deformations from a spherical equilibrium is less developed. We aim to understand how asymmetrically perturbed equilibrium bubbles evolve toward spherical equilibrium due to thermal or viscous d… ▽ More We study the dynamics of a gas bubble in a fluid with surface tension, initially near a spherical equilibrium. While there are many studies and applications of radial bubble dynamics, the theory of general deformations from a spherical equilibrium is less developed. We aim to understand how asymmetrically perturbed equilibrium bubbles evolve toward spherical equilibrium due to thermal or viscous dissipation in an incompressible liquid. We focus on the isobaric approximation [Prosperetti, JFM, 1991], under which the gas pressure within the bubble is spatially uniform and obeys the ideal gas law. The liquid outside the bubble is incompressible, irrotational, and has surface tension. We prove that any equilibrium gas bubble must be spherical by showing that the bubble boundary is a closed surface of constant mean curvature. We then study the initial value problem (IVP) for the coupled PDEs, constitutive laws and interface conditions of the isobaric approximation for general (asymmetric) small initial perturbations of the spherical bubble in the linearized approximation. Our first result, considering thermal damping without viscosity, proves that the linearized IVP is globally well-posed. The monopole (radial) component of the perturbation decays exponentially over time, while the multipole (non-radial) components undergo undamped oscillations. This indicates a limitation of the isobaric model for non-spherical dynamics. Our second result, incorporating viscous dissipation, shows that the IVP is linear and nonlinearly ill-posed due to an incompatibility of normal stress boundary conditions, for non-spherical solutions, and the irrotationality assumption. Our study concludes that to accurately capture the dynamics of general deformations of a gas bubble, the model must account for either vorticity generated at the bubble-fluid boundary, spatial non-uniformities in the gas pressure, or both. △ Less

Submitted 7 August, 2024; originally announced August 2024.

Gravity-driven viscous flow over partially lubricated bed

Authors: Joshua H. Rines, Ching-Yao Lai, Yongji Wang

Abstract: We present an investigation into the response of a viscous fluid flowing over a sloped bed across a spatially finite patch of basal lubrication. We present a simple analytical model that captures the fundamental structure of such lubrication-induced stress and velocity perturbations in Newtonian fluids, as well as scaling arguments and numerical experiments that extend our analysis to power-law fl… ▽ More We present an investigation into the response of a viscous fluid flowing over a sloped bed across a spatially finite patch of basal lubrication. We present a simple analytical model that captures the fundamental structure of such lubrication-induced stress and velocity perturbations in Newtonian fluids, as well as scaling arguments and numerical experiments that extend our analysis to power-law fluids. These analyses concisely reveal the underlying relationships between the system parameters (fluid thickness, $h$, slope, $α$, slippery patch length, $\ell$, and sliding condition outside of the slippery patch, $γ$) and the magnitude and spatial extent of the resulting perturbed stresses, $τ_{xx}$, and velocities, $u_p$. From these results, we conclude that the induced stresses are exponentially decaying functions of distance away from the patch location, and show that the amplitude of the perturbations scales linearly with surface slope and patch length while the decay length scales with thickness and patch length, and is critically dependent on the basal boundary condition outside of the slippery patch. These fundamental relationships can be incorporated into more complex models to investigate whether rapid lake drainages on ice sheets, which create a partially lubricated bed, can generate sufficient stress and velocity perturbations in the overlying ice flow to trigger lake drainage cascades. △ Less

Submitted 14 October, 2025; v1 submitted 30 July, 2024; originally announced July 2024.

Photocathode characterisation for robust PICOSEC Micromegas precise-timing detectors

Authors: M. Lisowska, R. Aleksan, Y. Angelis, S. Aune, J. Bortfeldt, F. Brunbauer, M. Brunoldi, E. Chatzianagnostou, J. Datta, K. Dehmelt, G. Fanourakis, S. Ferry, D. Fiorina, K. J. Floethner, M. Gallinaro, F. Garcia, I. Giomataris, K. Gnanvo, F. J. Iguaz, D. Janssens, A. Kallitsopoulou, M. Kovacic, B. Kross, C. C. Lai, P. Legou , et al. (33 additional authors not shown)

Abstract: The PICOSEC Micromegas detector is a~precise-timing gaseous detector based on a~Cherenkov radiator coupled with a~semi-transparent photocathode and a~Micromegas amplifying structure, targeting a~time resolution of tens of picoseconds for minimum ionising particles. Initial single-pad prototypes have demonstrated a~time resolution below 25 ps, prompting ongoing developments to adapt the concept for… ▽ More The PICOSEC Micromegas detector is a~precise-timing gaseous detector based on a~Cherenkov radiator coupled with a~semi-transparent photocathode and a~Micromegas amplifying structure, targeting a~time resolution of tens of picoseconds for minimum ionising particles. Initial single-pad prototypes have demonstrated a~time resolution below 25 ps, prompting ongoing developments to adapt the concept for High Energy Physics applications, where sub-nanosecond precision is essential for event separation, improved track reconstruction and particle identification. The achieved performance is being transferred to robust multi-channel detector modules suitable for large-area detection systems requiring excellent timing precision. To enhance the robustness and stability of the PICOSEC Micromegas detector, research on robust carbon-based photocathodes, including Diamond-Like Carbon (DLC) and Boron Carbide (B4C), is pursued. Results from prototypes equipped with DLC and B4C photocathodes exhibited a~time resolution of approximately 32 ps and 34.5 ps, respectively. Efforts dedicated to improve detector robustness and stability enhance the feasibility of the PICOSEC Micromegas concept for large experiments, ensuring sustained performance while maintaining excellent timing precision. △ Less

Submitted 9 December, 2024; v1 submitted 13 July, 2024; originally announced July 2024.

$\mathcal{PT}$-Symmetry induced Bi-Stability in Non-Hermitian Cavity Magnomechanics

Authors: Chaoyi Lai, Shah Fahad, Kashif Ammar Yasir

Abstract: We study the steady-state non-Hermitian magnomechanical system driven by a transverse magnetic field directly interacting with YIG sphere and excites cavity magnons and photons. To make the system non-Hermitian, we use a traveling field directly interacting with magnons generating gain to the system. We start by illustrating PT-configuration of the system, which contains two PT broken region aroun… ▽ More We study the steady-state non-Hermitian magnomechanical system driven by a transverse magnetic field directly interacting with YIG sphere and excites cavity magnons and photons. To make the system non-Hermitian, we use a traveling field directly interacting with magnons generating gain to the system. We start by illustrating PT-configuration of the system, which contains two PT broken region around exceptional point and PT protected region along the axis of exceptional point. Late, we discover that the numbers of cavity photons and magnons show bistable behavior depending upon the PT configuration, which becomes more significant as the values of the magnon-photon coupling and traveling field strength increases. We illustrate that steady-state photon only shows bistable behavior when the system in in lossy PT broken configuration, means strength of traveling field is less than the magnon-photon coupling. Otherwise, it will just contain a single stable state because of bistability suppression with gain in the system, which is unlike with any other investigation in this direction. Further, a larger magnon-photon coupling increases photon intensity and decreases magnon intensity, because of photon and magnon energy exchange, leading to enhanced photon bistablity and decreased magnon bistability. However, in case of increasing strength of traveling field, both photon as well as magnon bistability is appeared to be decreasing. We also study the steady-state effective potential of the system and illustrate the occurrence of bistability with nonlinear interactions between contour trajectories, which similarly depends on the PT broken configuration of the system. △ Less

Submitted 1 July, 2024; originally announced July 2024.

Comments: 10 pages, 6 figures

Machine learning for climate physics and simulations

Authors: Ching-Yao Lai, Pedram Hassanzadeh, Aditi Sheshadri, Maike Sonnewald, Raffaele Ferrari, Venkatramani Balaji

Abstract: We discuss the emerging advances and opportunities at the intersection of machine learning (ML) and climate physics, highlighting the use of ML techniques, including supervised, unsupervised, and equation discovery, to accelerate climate knowledge discoveries and simulations. We delineate two distinct yet complementary aspects: (1) ML for climate physics and (2) ML for climate simulations. While p… ▽ More We discuss the emerging advances and opportunities at the intersection of machine learning (ML) and climate physics, highlighting the use of ML techniques, including supervised, unsupervised, and equation discovery, to accelerate climate knowledge discoveries and simulations. We delineate two distinct yet complementary aspects: (1) ML for climate physics and (2) ML for climate simulations. While physics-free ML-based models, such as ML-based weather forecasting, have demonstrated success when data is abundant and stationary, the physics knowledge and interpretability of ML models become crucial in the small-data/non-stationary regime to ensure generalizability. Given the absence of observations, the long-term future climate falls into the small-data regime. Therefore, ML for climate physics holds a critical role in addressing the challenges of ML for climate simulations. We emphasize the need for collaboration among climate physics, ML theory, and numerical analysis to achieve reliable ML-based models for climate applications. △ Less

Submitted 17 August, 2024; v1 submitted 19 April, 2024; originally announced April 2024.

In-beam test results of an RPC-based module for position-sensitive neutron detectors with timing readout

Authors: G. Canezin, L. M. S. Margato, A. Morozov, A. Blanco, J. Saraiva, L. Lopes, P. Fonte, Chung Chuan Lai, Per-Olof Svensson, G. Markaj, Florian M. Piegsa

Abstract: Recently we have proposed a new concept of a thermal neutron detector based on resistive plate chambers and 10B4C solid neutron converters, enabling to readout with high resolution in both the 3D position of neutron capture and the neutron time of flight (ToF). In this paper, we report the results of the first beam tests conducted with a new neutron RPC detection module, coupled to the position re… ▽ More Recently we have proposed a new concept of a thermal neutron detector based on resistive plate chambers and 10B4C solid neutron converters, enabling to readout with high resolution in both the 3D position of neutron capture and the neutron time of flight (ToF). In this paper, we report the results of the first beam tests conducted with a new neutron RPC detection module, coupled to the position readout units of a new design. The main focus is on the measurements of the neutron ToF and identification of the converter layer where the neutron is captured, giving the position along the beam direction. △ Less

Submitted 23 February, 2024; originally announced February 2024.

Multi-Blade detector with VMM3a-ASIC-based readout: installation and commissioning at the reflectometer Amor at PSI

Authors: F. Piscitelli, F. Ghazi Moradi, F. S. Alves, M. J. Christensen, J. Hrivnak, A. Johansson, K. Fissum, C. C. Lai, A. Monera Martinez, D. Pfeiffer, E. Shahu, J. Stahn, P. O. Svensson

Abstract: The Multi-Blade (MB) Boron-10-based neutron detector is the chosen technology for three instruments at the European Spallation Source (ESS): the two ESS reflectometers, ESTIA and FREIA, and the Test Beam Line. A fourth MB detector has been built, installed and commissioned for the user operation of the reflectometer Amor at PSI (Switzerland). Amor can be considered a downscaled version of the ESS… ▽ More The Multi-Blade (MB) Boron-10-based neutron detector is the chosen technology for three instruments at the European Spallation Source (ESS): the two ESS reflectometers, ESTIA and FREIA, and the Test Beam Line. A fourth MB detector has been built, installed and commissioned for the user operation of the reflectometer Amor at PSI (Switzerland). Amor can be considered a downscaled version of the ESS reflectometer ESTIA. They are based on the same Selene guide concept, optimized for performing focusing reflectometry on small samples. The experience gained at Amor is invaluable for the future deployment of the MB detector at the ESS. This manuscript describes the MB detector construction and installation at Amor along with the readout electronics chain based on the VMM3a ASIC. The readout chain deployed at Amor is equivalent of that of the ESS, including the readout master module (RMM), event-formation-units (EFUs), Kafka, FileWriter and live visualisation tools. △ Less

Submitted 18 March, 2024; v1 submitted 13 February, 2024; originally announced February 2024.

Comments: 16 pages, 12 figures

Journal ref: 2024 JINST 19 P05010

On-Chip Stimulated Brillouin Scattering via Surface Acoustic Waves

Authors: Govert Neijts, Choon Kong Lai, Maren Kramer Riseng, Duk-Yong Choi, Kunlun Yan, David Marpaung, Stephen J. Madden, Benjamin J. Eggleton, Moritz Merklein

Abstract: Surface acoustic wave (SAW) devices are ubiquitously used for signal processing and filtering, as well as mechanical, chemical and biological sensing, and show promise as quantum transducers. However, nowadays most SAWs are excited and driven via electromechanical coupling and interdigital transducers (IDTs), limiting operation bandwidth and flexibility. Novel ways to coherently excite and detect… ▽ More Surface acoustic wave (SAW) devices are ubiquitously used for signal processing and filtering, as well as mechanical, chemical and biological sensing, and show promise as quantum transducers. However, nowadays most SAWs are excited and driven via electromechanical coupling and interdigital transducers (IDTs), limiting operation bandwidth and flexibility. Novel ways to coherently excite and detect SAWs all-optically interfaced with photonic integrated circuits are yet elusive. Backward Stimulated Brillouin scattering (SBS) provides strong coherent interactions between optical and acoustic waves in chip-scale waveguides, however, demonstrations have been limited to single longitudinal waves in the waveguide core. Here, we numerically model and experimentally demonstrate surface acoustic wave stimulated Brillouin scattering (SAW-SBS) on a photonic chip. We designed and fabricated tailored waveguides made out of GeAsSe glass that show good overlap between SAWs at 3.81 GHz and guided optical modes, without requiring a top cladding. We measure a 225 W$^{-1}$m$^{-1}$ Brillouin gain coefficient of the surface acoustic resonance and linewidth narrowing to 40 MHz. Experimentally accessing this new regime of stimulated Brillouin scattering opens the door for novel on-chip sensing and signal processing applications, strong Brillouin interactions in materials that do not provide sufficient acoustic guidance in the waveguide core as well as excitation of surface acoustic waves in non-piezoelectric materials. △ Less

Submitted 2 October, 2023; originally announced October 2023.

Comments: 7 pages, 4 figures

Journal ref: APL Photonics 1 October 2024; 9 (10): 106114

Soft matter physics of the ground beneath our feet

Authors: Anne Voigtländer, Morgane Houssais, Karol A. Bacik, Ian C. Bourg, Justin C. Burton, Karen E. Daniels, Sujit S. Datta, Emanuela Del Gado, Nakul S. Deshpande, Olivier Devauchelle, Behrooz Ferdowsi, Rachel Glade, Lucas Goehring, Ian J. Hewitt, Douglas Jerolmack, Ruben Juanes, Arshad Kudrolli, Ching-Yao Lai, Wei Li, Claire Masteller, Kavinda Nissanka, Allan M. Rubin, Howard A. Stone, Jenny Suckale, Nathalie M. Vriend , et al. (2 additional authors not shown)

Abstract: Inspired by presentations by the authors during a workshop organized at the Princeton Center for Theoretical Science (PCTS) in January 2022, we present a perspective on some of the outstanding questions related to the "physics of the ground beneath our feet." These identified challenges are intrinsically shared with the field of Soft Matter but also have unique aspects when the natural environment… ▽ More Inspired by presentations by the authors during a workshop organized at the Princeton Center for Theoretical Science (PCTS) in January 2022, we present a perspective on some of the outstanding questions related to the "physics of the ground beneath our feet." These identified challenges are intrinsically shared with the field of Soft Matter but also have unique aspects when the natural environment is studied. △ Less

Submitted 31 July, 2023; originally announced August 2023.

Comments: Perspective Paper, 30 pages, 15 figures

A Fuzzy Classification Framework to Identify Equivalent Atoms in Complex Materials and Molecules

Authors: King Chun Lai, Sebastian Matera, Christoph Scheurer, Karsten Reuter

Abstract: The nature of an atom in a bonded structure -- such as in molecules, in nanoparticles or solids, at surfaces or interfaces -- depends on its local atomic environment. In atomic-scale modeling and simulation, identifying groups of atoms with equivalent environments is a frequent task, to gain an understanding of the material function, to interpret experimental results or to simply restrict demandin… ▽ More The nature of an atom in a bonded structure -- such as in molecules, in nanoparticles or solids, at surfaces or interfaces -- depends on its local atomic environment. In atomic-scale modeling and simulation, identifying groups of atoms with equivalent environments is a frequent task, to gain an understanding of the material function, to interpret experimental results or to simply restrict demanding first-principles calculations. While routine, this task can often be challenging for complex molecules or non-ideal materials with breaks of symmetries or long-range order. To automatize this task, we here present a general machine-learning framework to identify groups of (nearly) equivalent atoms. The initial classification rests on the representation of the local atomic environment through a high-dimensional smooth overlap of atomic positions (SOAP) vector. Recognizing that not least thermal vibrations may lead to deviations from ideal positions, we then achieve a fuzzy classification by mean-shift clustering within a low-dimensional embedded representation of the SOAP points as obtained through multidimensional scaling. The performance of this classification framework is demonstrated for simple aromatic molecules and crystalline Pd surface examples. △ Less

Submitted 28 June, 2023; originally announced June 2023.

Comments: Accepted manuscript in Journal of Chemical Physics. Repositories of the package (DECAF): DOI:10.17617/3.U7VKBM or https://gitlab.mpcdf.mpg.de/klai/decaf

Provable Routing Analysis of Programmable Photonics

Authors: Zhengqi Gao, Xiangfeng Chen, Zhengxing Zhang, Chih-Yu Lai, Uttara Chakraborty, Wim Bogaerts, Duane S. Boning

Abstract: Programmable photonic integrated circuits (PPICs) are an emerging technology recently proposed as an alternative to custom-designed application-specific integrated photonics. Light routing is one of the most important functions that need to be realized on a PPIC. Previous literature has investigated the light routing problem from an algorithmic or experimental perspective, e.g., adopting graph the… ▽ More Programmable photonic integrated circuits (PPICs) are an emerging technology recently proposed as an alternative to custom-designed application-specific integrated photonics. Light routing is one of the most important functions that need to be realized on a PPIC. Previous literature has investigated the light routing problem from an algorithmic or experimental perspective, e.g., adopting graph theory to route an optical signal. In this paper, we also focus on the light routing problem, but from a complementary and theoretical perspective, to answer questions about what is possible to be routed. Specifically, we demonstrate that not all path lengths (defined as the number of tunable basic units that an optical signal traverses) can be realized on a square-mesh PPIC, and a rigorous realizability condition is proposed and proved. We further consider multi-path routing, where we provide an analytical expression on path length sum, upper bounds on path length mean/variance, and the maximum number of realizable paths. All of our conclusions are proven mathematically. Illustrative potential optical applications using our observations are also presented. △ Less

Submitted 21 June, 2023; originally announced June 2023.

Optimizing performance for on-chip SBS-based isolator

Authors: Choon Kong Lai, Moritz Merklein, Alvaro Casas Bedoya, Yang Liu, Stephen J. Madden, Christopher G. Poulton, Michael J. Steel, Benjamin J. Eggleton

Abstract: Non-reciprocal optical components such as isolators and circulators are crucial for preventing catastrophic back-reflection and controlling optical crosstalk in photonic systems. While non-reciprocal devices based on Brillouin intermodal transitions have been experimentally demonstrated in chip-scale platforms, harnessing such interactions has required a suspended waveguide structure, which is cha… ▽ More Non-reciprocal optical components such as isolators and circulators are crucial for preventing catastrophic back-reflection and controlling optical crosstalk in photonic systems. While non-reciprocal devices based on Brillouin intermodal transitions have been experimentally demonstrated in chip-scale platforms, harnessing such interactions has required a suspended waveguide structure, which is challenging to fabricate and is potentially less robust than a non-suspended structure, thereby limiting the design flexibility. In this paper, we numerically investigate the performance of a Brillouin-based isolation scheme in which a dual-pump-driven optoacoustic interaction is used to excite confined acoustic waves in a traditional ridge waveguide. We find that acoustic confinement, and therefore the amount of Brillouin-driven mode conversion, can be enhanced by selecting an appropriate optical mode pair and waveguide geometry of two arsenic based chalcogenide platforms. Further, we optimize the isolator design in its entirety, including the input couplers, mode filters, the Brillouin-active waveguide as well as the device fabrication tolerances. We predict such a device can achieve 30 dB isolation over a 38 nm bandwidth when 500 mW pump power is used; in the presence of a +/- 10 nm fabrication-induced width error, such isolation can be maintained over a 5-10 nm bandwidth. △ Less

Submitted 25 December, 2022; originally announced December 2022.

Double-GEM based thermal neutron detector prototype

Authors: L. A. Serra Filho, R. Felix dos Santos, G. G. A. de Souza, M. M. M. Paulino, F. A. Souza, M. Moralles, H. Natal da Luz, M. Bregant, M. G. Munhoz, Chung-Chuan Lai, Carina Höglund, Per-Olof Svensson, Linda Robinson, Richard Hall-Wilton

Abstract: The Helium-3 shortage and the growing interest in neutron science constitute a driving factor in developing new neutron detection technologies. In this work, we report the development of a double-GEM detector prototype that uses a $^{10}$B$_4$C layer as a neutron converter material. GEANT4 simulations were performed predicting an efficiency of 3.14(10) %, agreeing within 2.7 $σ$ with the experimen… ▽ More The Helium-3 shortage and the growing interest in neutron science constitute a driving factor in developing new neutron detection technologies. In this work, we report the development of a double-GEM detector prototype that uses a $^{10}$B$_4$C layer as a neutron converter material. GEANT4 simulations were performed predicting an efficiency of 3.14(10) %, agreeing within 2.7 $σ$ with the experimental and analytic detection efficiencies obtained by the detector when tested in a 41.8 meV thermal neutron beam. The detector is position sensitive, equipped with a 256+256 strip readout connected to resistive chains, and achieves a spatial resolution better than 3 mm. The gain stability over time was also measured with a fluctuation of about 0.2 %h$^{-1}$ of the signal amplitude. A simple data acquisition with only 5 electronic channels is sufficient to operate this detector. △ Less

Submitted 19 July, 2022; v1 submitted 14 May, 2022; originally announced May 2022.

Comments: Prepared for submission to JINST

Journal ref: JINST 17 P09018 (2022)

Sub-40nm Nanogratings Self-Organized in PVP-based Polymer Composite Film by Photoexcitation and Two Sequent Splitting under Femtosecond Laser Irradiation

Authors: Li-Yun Chen, Cheng-Cheng Guo, Ming-Ming Pan, Chen Lai, Yun-Xia Wang, Guo-Cai Liao, Zi-Wei Ma, Fan-Wei Zhang, Jagadeesh Suriyaprakash, Lijing Guo, Eser Akinoglu, Qiang Li, Li-Jun Wu

Abstract: Laser-induced periodic surface structures (LIPSSs) on various materials have been extensively investigated because of their wide applications. The combination of different materials allows for greater freedom in tailoring their functions and achieving responses not possible in a homogeneous material. By utilizing a femtosecond (fs) laser to irradiate the Fe-doped Polyvinyl Pyrrolidone (PVP) compos… ▽ More Laser-induced periodic surface structures (LIPSSs) on various materials have been extensively investigated because of their wide applications. The combination of different materials allows for greater freedom in tailoring their functions and achieving responses not possible in a homogeneous material. By utilizing a femtosecond (fs) laser to irradiate the Fe-doped Polyvinyl Pyrrolidone (PVP) composite film, highly regular ultrafine nanogratings (U-nanogratings) with a period as small as 35.0 ($\pm$ 2.0) nm can be self-organized on the surface with extremely high efficiency. The period of the U-nanogratings can be controlled by varying the scanning speed of the laser beam (deposited energy) and the thickness of the composite film. Based on the experimental, theoretical, and simulation results, we propose a two-step formation mechanism: composite film excitation and two sequent grating-splitting. The high photosensitivity and low glass transition temperature of the composite film facilitate the fabrication of the ultrafine nanostructures. The proposed design method for the composite material and fabrication process could not only provide a strategy for obtaining highly regular U-nanogratings, but also offer a platform to explore the interaction physics between ultra-short pulses and matter under extreme conditions. △ Less

Submitted 15 March, 2022; originally announced March 2022.

Multitube monitors: a new-generation of neutron beam monitors

Authors: F. Lafont, D. Barkats, J-C. Buffet, S. Cuccaro, B. Guerard, C-C Lai, J. Marchal, J. Pentenero, N. Sartor, R. Hall-Wilton, K. Kanaki, L. Robinson, P O. Svensson

Abstract: With the renewal of many neutron science instruments and the commissioning of new neutron facilities, there is a rising demand for improved neutron beam monitoring systems with reduced beam perturbations and higher counting rate capability. Fission chambers are the most popular beam monitors; however, their use on some instruments may be prevented by the background generated by fast neutrons emitt… ▽ More With the renewal of many neutron science instruments and the commissioning of new neutron facilities, there is a rising demand for improved neutron beam monitoring systems with reduced beam perturbations and higher counting rate capability. Fission chambers are the most popular beam monitors; however, their use on some instruments may be prevented by the background generated by fast neutrons emitted during neutron captures in 235U and by neutrons scattered in the material of the fission chamber. Multitube detectors, on the other hand, offer a good alternative with minimum beam perturbations. The purpose of this paper is to report and analyse the results of the measurements performed with several Multitubes used for beam monitoring. We show that the transparency of Multitube beam monitors is 97.6 +/-0.4 %, and that their detection efficiency is uniform, with a deviation from the mean value < 0.7%. A counting rate reduction of 10% due to pile-up effects is measured at a rate of 550 kHz. In addition to neutron beam intensity monitoring, the Multitube can be configured for 1-dimensional or 2-dimensional localisation. We present the preliminary results of these additional functionalities. △ Less

Submitted 19 May, 2022; v1 submitted 3 February, 2022; originally announced February 2022.

Comments: 15 pages, 16 figures

Effect of substrate roughness and material selection on the microstructure of sputtering deposited boron carbide thin films

Authors: Chung-Chuan Lai, Robert Boyd, Per-Olof Svensson, Carina Höglund, Linda Robinson, Jens Birch, Richard Hall-Wilton

Abstract: Amorphous boron carbide (B4C) thin films are by far the most popular form for the neutron converting layers in the 10B-based neutron detectors, which are a rising trend in detector technologies in response to the increasing scarcity and price of 3He, the standard material for neutron detection. The microstructure of the B4C films is closely related to the important properties, e.g. density and adh… ▽ More Amorphous boron carbide (B4C) thin films are by far the most popular form for the neutron converting layers in the 10B-based neutron detectors, which are a rising trend in detector technologies in response to the increasing scarcity and price of 3He, the standard material for neutron detection. The microstructure of the B4C films is closely related to the important properties, e.g. density and adhesion, for the converting layers, which eventually affect the detection efficiency and the long-term stability of the detectors. To study the influence from substrates of different roughness and materials, the B4C films were deposited on polished Si substrates with Al, Ti, and Cu buffer layers and unpolished Si, Al, Ti, and Cu substrates by direct current magnetron sputtering at a substrate temperature of 623 K. The tapered columnar grains and nodular defects, generally observed in SEM images, indicated a strong shadowing effect where voids were introduced around the grains. The change in the grain size did not show a direct dependence to the substrate roughness, acquired from the surface profile, nor to the mass density of the films, obtained from reflectivity patterns. However, films with non-uniform size of columnar grains were deposited on substrates with high skewness, leading to a drop of mass density from ~95 % down to ~70 % of tabulated bulk density. On the other hand, similar microstructures and mass density were obtained from the films deposited on Al, Ti, and Cu of different roughness and good adhesion were observed from cross-cut adhesion tests, showing the reliability of sputtering deposited B4C films on common structural materials in neutron detectors. △ Less

Submitted 29 January, 2022; originally announced January 2022.

Comments: Available online at Surface & Coatings Technology

Asymptotic self-similar blow-up profile for three-dimensional axisymmetric Euler equations using neural networks

Authors: Yongji Wang, Ching-Yao Lai, Javier Gómez-Serrano, Tristan Buckmaster

Abstract: Whether there exist finite time blow-up solutions for the 2-D Boussinesq and the 3-D Euler equations are of fundamental importance to the field of fluid mechanics. We develop a new numerical framework, employing physics-informed neural networks (PINNs), that discover, for the first time, a smooth self-similar blow-up profile for both equations. The solution itself could form the basis of a future… ▽ More Whether there exist finite time blow-up solutions for the 2-D Boussinesq and the 3-D Euler equations are of fundamental importance to the field of fluid mechanics. We develop a new numerical framework, employing physics-informed neural networks (PINNs), that discover, for the first time, a smooth self-similar blow-up profile for both equations. The solution itself could form the basis of a future computer-assisted proof of blow-up for both equations. In addition, we demonstrate PINNs could be successfully applied to find unstable self-similar solutions to fluid equations by constructing the first example of an unstable self-similar solution to the Córdoba-Córdoba-Fontelos equation. We show that our numerical framework is both robust and adaptable to various other equations. △ Less

Submitted 7 May, 2023; v1 submitted 18 January, 2022; originally announced January 2022.

Comments: Main paper: 6 pages, 3 figures. Supplementary material: 15 pages, 12 figures. To appear in Physical Review Letters

uRANIA-V: an innovative solution for neutron detection in homeland security

Authors: R. Farinelli, I. Balossino, G. Bencivenni, G. Cibinetto, G. Felici, S. Fiore, I. Garzia, M. Gatta, M. Giovannetti, R. Hall-Wilton, C. C. Lai, L. Lavezzi, G. Mezzadri, G. Morello, E. Paoletti, G. Papalino, A. Pietropaolo, M. Pillon, M. Poli Lener, L. Robinson, M. Scodeggio, P. O. Svensson

Abstract: Detection of neutrons is becoming of the utmost importance, especially in the studies of radioactive waste and in homeland security applications. The crisis of 3He availability has required the development of innovative techniques. One solution is to develop light gas detectors for neutron counting to be used as portals for ports and airports. The neutron is converted on the Boron-coated cathode,… ▽ More Detection of neutrons is becoming of the utmost importance, especially in the studies of radioactive waste and in homeland security applications. The crisis of 3He availability has required the development of innovative techniques. One solution is to develop light gas detectors for neutron counting to be used as portals for ports and airports. The neutron is converted on the Boron-coated cathode, releasing a charged particle, whose passage can be identified by the gas detector. While several technologies have been deployed in the past, the project μRANIA-V ( μRwell Advanced Neutron Identification Apparatus) aims to detect thermal neutrons by means of the μRwell technology, an innovative gas detector. The goal is to produce tiles to operate as portals in homeland security or for radioactive waste management. The technological transfer towards the industry has started, thus the production can be cost-effective also owing to a construction process relatively easier compared to similar apparatus. By reading directly the signals from the amplification stage, the neutrons can be counted with simplified electronics further reducing the total cost. In this paper, the project will be described, with details on the μRwell technology and on the neutron counting, on the test beam performed, and on the future plans. △ Less

Submitted 2 September, 2021; originally announced September 2021.

Comments: ANIMMA 2021 conference proceeding

Edge of chaos as a guiding principle for modern neural network training

Authors: Lin Zhang, Ling Feng, Kan Chen, Choy Heng Lai

Abstract: The success of deep neural networks in real-world problems has prompted many attempts to explain their training dynamics and generalization performance, but more guiding principles for the training of neural networks are still needed. Motivated by the edge of chaos principle behind the optimal performance of neural networks, we study the role of various hyperparameters in modern neural network tra… ▽ More The success of deep neural networks in real-world problems has prompted many attempts to explain their training dynamics and generalization performance, but more guiding principles for the training of neural networks are still needed. Motivated by the edge of chaos principle behind the optimal performance of neural networks, we study the role of various hyperparameters in modern neural network training algorithms in terms of the order-chaos phase diagram. In particular, we study a fully analytical feedforward neural network trained on the widely adopted Fashion-MNIST dataset, and study the dynamics associated with the hyperparameters in back-propagation during the training process. We find that for the basic algorithm of stochastic gradient descent with momentum, in the range around the commonly used hyperparameter values, clear scaling relations are present with respect to the training time during the ordered phase in the phase diagram, and the model's optimal generalization power at the edge of chaos is similar across different training parameter combinations. In the chaotic phase, the same scaling no longer exists. The scaling allows us to choose the training parameters to achieve faster training without sacrificing performance. In addition, we find that the commonly used model regularization method - weight decay - effectively pushes the model towards the ordered phase to achieve better performance. Leveraging on this fact and the scaling relations in the other hyperparameters, we derived a principled guideline for hyperparameter determination, such that the model can achieve optimal performance by saturating it at the edge of chaos. Demonstrated on this simple neural network model and training algorithm, our work improves the understanding of neural network training dynamics, and can potentially be extended to guiding principles of more complex model architectures and algorithms. △ Less

Submitted 20 July, 2021; originally announced July 2021.

CO$_2$-driven diffusiophoresis for removal of bacteria

Authors: Suin Shim, Sepideh Khodaparast, Ching-Yao Lai, Jing Yan, Jesse T. Ault, Bhargav Rallabandi, Orest Shardt, Howard A. Stone

Abstract: We investigate CO$_2$-driven diffusiophoresis of colloidal particles and bacterial cells in a Hele-Shaw geometry. Combining experiments and a model, we understand the characteristic length and time scales of CO$_2$-driven diffusiophoresis in relation to system dimensions and CO$_2$ diffusivity. Directional migration of wild-type V. cholerae and a mutant lacking flagella, as well as S. aureus and P… ▽ More We investigate CO$_2$-driven diffusiophoresis of colloidal particles and bacterial cells in a Hele-Shaw geometry. Combining experiments and a model, we understand the characteristic length and time scales of CO$_2$-driven diffusiophoresis in relation to system dimensions and CO$_2$ diffusivity. Directional migration of wild-type V. cholerae and a mutant lacking flagella, as well as S. aureus and P. aeruginosa, near a dissolving CO$_2$ source shows that diffusiophoresis of bacteria is achieved independent of cell shape and Gram stain. Long-time experiments suggest possible applications for bacterial diffusiophoresis to cleaning systems or anti-biofouling surfaces. △ Less

Submitted 14 September, 2020; originally announced September 2020.

Time- and energy-resolved effects in the boron-10 based Multi-Grid and helium-3 based thermal neutron detectors

Authors: A. Backis, A. Khaplanov, R. Al Jebali, R. Ammer, I. Apostolidis, J. Birch, C. -C. Lai, P. P. Deen, M. Etxegarai, N. de Ruette, J. Freita Ramos, D. F. Förster, E. Haettner, R. Hall-Wilton, D. Hamilton, C. Höglund, P. M. Kadletz, K. Kanaki, E. Karnickis, O. Kirstein, S. Kolya, Z. Kraujalyte, A. Laloni, K. Livingston, O. Löhman , et al. (11 additional authors not shown)

Abstract: The boron-10 based Multi-Grid detector is being developed as an alternative to helium-3 based neutron detectors. At the European Spallation Source, the detector will be used for time-of-flight neutron spectroscopy at cold to thermal neutron energies. The objective of this work is to investigate fine time- and energy-resolved effects of the Multi-Grid detector, down to a few $μ$eV, while comparing… ▽ More The boron-10 based Multi-Grid detector is being developed as an alternative to helium-3 based neutron detectors. At the European Spallation Source, the detector will be used for time-of-flight neutron spectroscopy at cold to thermal neutron energies. The objective of this work is to investigate fine time- and energy-resolved effects of the Multi-Grid detector, down to a few $μ$eV, while comparing it to the performance of a typical helium-3 tube. Furthermore, it is to characterize differences between the detector technologies in terms of internal scattering, as well as the time reconstruction of ~ $μ$s short neutron pulses. The data were taken at the Helmholtz Zentrum Berlin, where the Multi-Grid detector and a helium-3 tube were installed at the ESS test beamline, V20. Using a Fermi-chopper, the neutron beam of the reactor was chopped into a few tens of $μ$s wide pulses before reaching the detector, located a few tens of cm downstream. The data of the measurements show an agreement between the derived and calculated neutron detection efficiency curve. The data also provide fine details on the effect of internal scattering, and how it can be reduced. For the first time, the chopper resolution was comparable to the timing resolution of the Multi-Grid detector. This allowed a detailed study of time- and energy resolved effects, as well as a comparison with a typical helium-3 tube. △ Less

Submitted 14 January, 2021; v1 submitted 2 June, 2020; originally announced June 2020.

Comments: 37 pages, 21 figures

Journal ref: Meas. Sci. Technol. 32 035903 (2021)

u-RANIA: a neutron detector based on μ-RWELL technology

Authors: I. Balossino, G. Bencivenni, P. Bielowka, G. Cibinetto, R. Farinelli, G. Felici, I. Garzia, M. Gatta, P. Giacomelli, M. Giovannetti, R. Hall Wilton, C. -C. Lai, L. Lavezzi, F. Messi, G. Mezzadri, G. Morello, M. Pinamonti, M. Poli Lener, L. Robinson, M. Scodeggio, P. -O. Svensson

Abstract: In the framework of the ATTRACT-uRANIA project, funded by the European Community, we are developing an innovative neutron imaging detector based on micro-Resistive WELL ($μ$ -RWELL) technology. The $μ$ -RWELL, based on the resistive detector concept, ensuring an efficient spark quenching mechanism, is a highly reliable device. It is composed by two main elements: a readout-PCB and a cathode. The a… ▽ More In the framework of the ATTRACT-uRANIA project, funded by the European Community, we are developing an innovative neutron imaging detector based on micro-Resistive WELL ($μ$ -RWELL) technology. The $μ$ -RWELL, based on the resistive detector concept, ensuring an efficient spark quenching mechanism, is a highly reliable device. It is composed by two main elements: a readout-PCB and a cathode. The amplification stage for this device is embedded in the readout board through a resistive layer realized by means of an industrial process with DLC (Diamond-Like Carbon). A thin layer of B$_4$C on the copper surface of the cathode allows the thermal neutrons detection through the release of $^7$Li and $α$ particles in the active volume. This technology has been developed to be an efficient and convenient alternative to the $^3$He shortage. The goal of the project is to prove the feasibility of such a novel neutron detector by developing and testing small planar prototypes with readout boards suitably segmented with strip or pad read out, equipped with existing electronics or readout in current mode. Preliminary results from the test with different prototypes, showing a good agreement with the simulation, will be presented together with construction details of the prototypes and the future steps of the project. △ Less

Submitted 17 August, 2020; v1 submitted 13 May, 2020; originally announced May 2020.

Comments: Prepared for the INSTR20 Conference Proceeding for JINST

Evidence for spontaneous arrangement of two-way flow in water bridge via particle image velocimetry

Authors: Ping-Rui Tsai, Hong-Yue Huang, Cheng-Wei Lai, Yu-Ting Cheng, Chih-Yung Huang, Cheng-En Tsai, Yi-Chun Lee, Hong Hao, Tzay-Ming Hong

Abstract: By revisiting the century-old problem of water bridge, we demonstrate that it is in fact dynamic and comprises of two coaxial water currents that carry different charges and flow in opposite directions. This spontaneous separation is triggered by the different stages to construct the water bridge. Initially, a flow is facilitated by the cone jet that is powered by H+ and flows out of the positive-… ▽ More By revisiting the century-old problem of water bridge, we demonstrate that it is in fact dynamic and comprises of two coaxial water currents that carry different charges and flow in opposite directions. This spontaneous separation is triggered by the different stages to construct the water bridge. Initially, a flow is facilitated by the cone jet that is powered by H+ and flows out of the positive-electrode beaker. An opposing cone-jet from negative beaker is established later and forced to take the outer route. This spontaneous arrangement of two-way flow is revealed by using fluorescein and carbon powder as tracers, and the Particle Image Velocimetry, These two opposing flows are found to carry non-equal flux that results in a net transport of water to the negative beaker. We manage to estimate the flow speed and cross-sectional area of these co-axial flows as a function of time and applied voltage. Note that the water on the outer layer functions as a millimeter tube that confines and interacts strongly with the flow inside. This provides a first natural and yet counter example to the recently reported near-frictionless flow in an equally miniatureized soft wall made from ferrofluid. △ Less

Submitted 25 July, 2023; v1 submitted 10 April, 2020; originally announced April 2020.

Comments: 6 pages, 5 figures

Multi-Functional Variable Thickness Structure for Broadband and Omnidirectional Focusing and Collimation

Authors: Liuxian Zhao, Changquan Lai, Miao Yu

Abstract: Luneburg lens is a symmetric gradient-index lens with a refractive index that increases from the outer surface to the center in a radial manner. It has the ability to focus and collimate waves, which makes it useful for energy harvesting, waveguiding and as a component in transducers. An ideal Luneburg lens should be easy to fabricate, has broadband and omnidirectional characteristics, as well as… ▽ More Luneburg lens is a symmetric gradient-index lens with a refractive index that increases from the outer surface to the center in a radial manner. It has the ability to focus and collimate waves, which makes it useful for energy harvesting, waveguiding and as a component in transducers. An ideal Luneburg lens should be easy to fabricate, has broadband and omnidirectional characteristics, as well as a focal length that can be easily tuned. However, existing structural Luneburg lenses based on phononic crystals can hardly achieve these requirements. Here, we propose an alternative structural Luneburg lens which has a refractive index that varies smoothly with its radial distance as a result of a changing thickness. Theoretical calculations, numerical simulations and experimental measurements of flexural wave propagation through the lens showed that flexural wave focusing can be obtained inside, at the edge and outside of the variable thickness lens for different frequencies and propagation directions. Flexural wave collimation was also demonstrated when a point source was placed at the respective focal points for each lens. Furthermore, it was shown that flexural waves that were focused onto a piezoelectric energy harvester by the Luneburg lens can lead to a significant increase in the harvested voltage compared to that obtained without focusing. △ Less

Submitted 20 February, 2020; originally announced February 2020.

The Multi-Blade Boron-10-based neutron detector performance using a focusing reflectometer

Authors: G. Mauri, I. Apostolidis, M. J. Christensen, A. Glavic, C. C. Lai, A. Laloni, F. Messi, A. Lindh Olsson, L. Robinson, J. Stahn, P. O. Svensson, R. Hall-Wilton, F. Piscitelli

Abstract: The Multi-Blade is a Boron-10-based neutron detector designed for neutron reflectometers and developed for the two instruments (Estia and FREIA) planned for the European Spallation Source in Sweden. A reflectometry demonstrator has been installed at the AMOR reflectometer at the Paul Scherrer Institut (PSI - Switzerland). The setup exploits the Selene guide concept and it can be considered a scale… ▽ More The Multi-Blade is a Boron-10-based neutron detector designed for neutron reflectometers and developed for the two instruments (Estia and FREIA) planned for the European Spallation Source in Sweden. A reflectometry demonstrator has been installed at the AMOR reflectometer at the Paul Scherrer Institut (PSI - Switzerland). The setup exploits the Selene guide concept and it can be considered a scaled-down demonstrator of Estia. The results of these tests are discussed. It will be shown how the characteristics of the Multi-Blade detector are features that allow the focusing reflectometry operation mode. Additionally the performance of the Multi-Blade, in terms of rate capability, exceeds current state-of-the-art technology. The improvements with respect to the previous prototypes are also highlighted; from background considerations to the linear and angular uniformity response of the detector. △ Less

Submitted 9 January, 2020; originally announced January 2020.

Comments: 32 pages, 28 figures

Journal ref: Journal of Instrumentation 15, P03010, (2020)

Zero-field spin-orbit-torque switching driven by magnetic spin Hall effect

Authors: Po-Hung Lin, Po-Wei Lee, Yu-Hsuan Lin, Bo-Yuan Yang, Vinod Kumar, Hsiu-Hau Lin, Chih-Huang Lai

Abstract: Spin Hall effect plays an essential role in generating spin current from the injected charge current, following the Dyakonov-Perel rule that the directions of charge flow, spin flow and spin polarization are mutually perpendicular to each other. Recently, its generalization from an antiferromagnet, so-called magnetic spin Hall effect, has been studied and verified by measuring anomalous spin accum… ▽ More Spin Hall effect plays an essential role in generating spin current from the injected charge current, following the Dyakonov-Perel rule that the directions of charge flow, spin flow and spin polarization are mutually perpendicular to each other. Recently, its generalization from an antiferromagnet, so-called magnetic spin Hall effect, has been studied and verified by measuring anomalous spin accumulations. Here, we investigate the magnetic spin Hall effect in bilayer materials made of a heavy metal and an antiferromagnet. The spin current generated by the magnetic spin Hall effect accomplishes spin-orbit-torque switching for ferromagnetic magnetization and exchange bias concurrently without any external magnetic field. The switching mechanism crucially relies on the non-collinear spin texture in the antiferromagnet, capable of generating symmetry-breaking components in the spin-current tensor so that the external magnetic field is no longer necessary. The zero-field concurrent switching of magnetization and exchange bias is a significant technological breakthrough. Furthermore, our findings pave the way to explore the magnetic spin Hall effects in various spin textures through spin-orbit-torque switching. △ Less

Submitted 5 November, 2019; originally announced November 2019.

Propagation of Uncertainty in a Rotating Pipe Mechanism to Generate an Impinging Swirling Jet Flow for Heat Transfer from a Flat Plate

Authors: F. -J. Granados-Ortiz, J. Ortega-Casanova, C. -H. Lai

Abstract: In Computational Fluid Dynamics (CFD) studies composed of the coupling of different simulations, the uncertainty in one stage may be propagated to the following stage and affect the accuracy of the prediction. In this paper, a framework for uncertainty quantification is applied to the two-step simulation of the mechanical design of a swirling jet flow generated by a rotating pipe (\textit{Simulati… ▽ More In Computational Fluid Dynamics (CFD) studies composed of the coupling of different simulations, the uncertainty in one stage may be propagated to the following stage and affect the accuracy of the prediction. In this paper, a framework for uncertainty quantification is applied to the two-step simulation of the mechanical design of a swirling jet flow generated by a rotating pipe (\textit{Simulation 1}) impinging on a flat plate to provide convective heat transfer (\textit{Simulation 2}). The first approach is the Stochastic Collocation Method (SCM) with Clenshaw-Curtis sparse grids. The conclusion drawn from the analysis is that the simulated system does not exhibit a significant sensitivity to stochastic variations of model input parameters, over the tested uncertainty ranges. Additionally, a set of non-linear regression models for the stochastic velocity and turbulent profiles for the pipe nozzle are created and tested, since impinging jets at Reynolds number of $Re=23000$ are very frequent in the literature, but stochastic inlet conditions have never been provided. Numerical results demonstrate a negligible difference in the predicted convective heat transfer with respect to the use of the profiles simulated via CFD. These suggested surrogate models can be directly embedded onto other CFD applications (e.g arrays of jets or jet flows impinging on plates with different shapes) in which a realistic swirling flow under uncertainty can be of interest. △ Less

Submitted 27 October, 2019; originally announced October 2019.

Turbulence statistics in a negatively buoyant multiphase plume

Authors: Ankur D. Bordoloi, Chris C. K. Lai, Laura K. Clark, Gerardo Veliz, Evan Variano

Abstract: We investigate the turbulence statistics in a {multiphase plume made of heavy particles (particle Reynolds number at terminal velocity is 450)}. Using refractive-index-matched stereoscopic particle image velocimetry, we measure the locations of particles {whose buoyancy drives the formation of a multiphase plume,} {together with the local velocity of the induced flow in the ambient salt-water}. {M… ▽ More We investigate the turbulence statistics in a {multiphase plume made of heavy particles (particle Reynolds number at terminal velocity is 450)}. Using refractive-index-matched stereoscopic particle image velocimetry, we measure the locations of particles {whose buoyancy drives the formation of a multiphase plume,} {together with the local velocity of the induced flow in the ambient salt-water}. {Measurements in the plume centerplane exhibit self-similarity in mean flow characteristics consistent with classic integral plume theories.} The turbulence characteristics resemble those measured in a bubble plume, {including strong anisotropy in the normal Reynolds stresses. However, we observe structural differences between the two multiphase plumes. First, the skewness of the probability density function (PDF) of the axial velocity fluctuations is not that which would be predicted by simply reversing the direction of a bubble plume. Second, in contrast to a bubble plume, the particle plume has a non-negligible fluid-shear production term in the turbulent kinetic energy (TKE) budget. Third, the radial decay of all measured terms in the TKE budget is slower than those in a bubble plume.} Despite these dissimilarities, a bigger picture emerges that applies to both flows. The TKE production by particles (or bubbles) roughly balances the viscous dissipation, except near the plume centerline. The one-dimensional power-spectra of the velocity fluctuations show a -3 power-law that puts both the particle and bubble plume in a category different from single-phase shear-flow turbulence. △ Less

Submitted 18 July, 2019; originally announced July 2019.

Journal ref: J. Fluid Mech. 896 (2020) A19

An immersed boundary method for fluid--structure--acoustics interactions involving large deformations and complex geometries

Authors: Li Wang, Fang-Bao Tian, Joseph C. S. Lai

Abstract: This paper presents an immersed boundary (IB) method for fluid--structure--acoustics interactions involving large deformations and complex geometries. In this method, the fluid dynamics is solved by a finite difference method where the temporal, viscous and convective terms are respectively discretized by the third-order Runge-Kutta scheme, the fourth-order central difference scheme and a fifth-or… ▽ More This paper presents an immersed boundary (IB) method for fluid--structure--acoustics interactions involving large deformations and complex geometries. In this method, the fluid dynamics is solved by a finite difference method where the temporal, viscous and convective terms are respectively discretized by the third-order Runge-Kutta scheme, the fourth-order central difference scheme and a fifth-order W/TENO (Weighted/Targeted Essentially Non-oscillation) scheme. Without loss of generality, a nonlinear flexible plate is considered here, and is solved by a finite element method based on the absolute nodal coordinate formulation. The no-slip boundary condition at the fluid--structure interface is achieved by using a diffusion-interface penalty IB method. With the above proposed method, the aeroacoustics field generated by the moving boundaries and the associated flows are inherently solved. In order to validate and verify the current method, several benchmark cases are conducted: acoustic waves scattered from a stationary cylinder in a quiescent flow, sound generation by a stationary and a rotating cylinder in a uniform flow, sound generation by an insect in hovering flight, deformation of a red blood cell induced by acoustic waves and acoustic waves scattered by a stationary sphere. The comparison of the sound scattered by a cylinder shows that the present IB--WENO scheme, a simple approach, has an excellent performance which is even better than the implicit IB--lattice Boltzmann method. For the sound scattered by a sphere, the IB--TENO scheme has a lower dissipation compared with the IB--WENO scheme. Applications of this technique to model fluid-structure-acoustics interactions of flapping foils mimicking an insect wing section during forward flight and flapping foil energy harvester are also presented, considering the effects of foil shape and flexibility. △ Less

Submitted 1 May, 2019; originally announced May 2019.

Production and Integration of the ATLAS Insertable B-Layer

Authors: B. Abbott, J. Albert, F. Alberti, M. Alex, G. Alimonti, S. Alkire, P. Allport, S. Altenheiner, L. Ancu, E. Anderssen, A. Andreani, A. Andreazza, B. Axen, J. Arguin, M. Backhaus, G. Balbi, J. Ballansat, M. Barbero, G. Barbier, A. Bassalat, R. Bates, P. Baudin, M. Battaglia, T. Beau, R. Beccherle , et al. (352 additional authors not shown)

Abstract: During the shutdown of the CERN Large Hadron Collider in 2013-2014, an additional pixel layer was installed between the existing Pixel detector of the ATLAS experiment and a new, smaller radius beam pipe. The motivation for this new pixel layer, the Insertable B-Layer (IBL), was to maintain or improve the robustness and performance of the ATLAS tracking system, given the higher instantaneous and i… ▽ More During the shutdown of the CERN Large Hadron Collider in 2013-2014, an additional pixel layer was installed between the existing Pixel detector of the ATLAS experiment and a new, smaller radius beam pipe. The motivation for this new pixel layer, the Insertable B-Layer (IBL), was to maintain or improve the robustness and performance of the ATLAS tracking system, given the higher instantaneous and integrated luminosities realised following the shutdown. Because of the extreme radiation and collision rate environment, several new radiation-tolerant sensor and electronic technologies were utilised for this layer. This paper reports on the IBL construction and integration prior to its operation in the ATLAS detector. △ Less

Submitted 6 June, 2018; v1 submitted 2 March, 2018; originally announced March 2018.

Comments: 90 pages in total. Author list: ATLAS IBL Collaboration, starting page 2. 69 figures, 20 tables. Published in Journal of Instrumentation. All figures available at: https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PLOTS/PIX-2018-001

Journal ref: Journal of Instrumentation JINST 13 T05008 (2018)