Showing 1–50 of 111 results for author: Lopez, A

Design, Construction, and Test of Compact, Distributed-Charge, X-Band Accelerator Systems that Enable Image-Guided, VHEE FLASH Radiotherapy

Authors: Christopher P. J. Barty, J. Martin Algots, Alexander J. Amador, James C. R. Barty, Shawn M. Betts, Marcelo A. Casteñada, Matthew M. Chu, Michael E. Daley, Ricardo A. De Luna Lopez, Derek A. Diviak, Haytham H. Effarah, Roberto Feliciano, Adan Garcia, Keith J. Grabiel, Alex S. Griffin, Frederic V. Hartemann, Leslie Heid, Yoonwoo Hwang, Gennady Imeshev, Michael Jentschel, Christopher A. Johnson, Kenneth W. Kinosian, Agnese Lagzda, Russell J. Lochrie, Michael W. May , et al. (18 additional authors not shown)

Abstract: The design and optimization of laser-Compton x-ray systems based on compact distributed charge accelerator structures can enable micron-scale imaging of disease and the concomitant production of beams of Very High Energy Electrons (VHEEs) capable of producing FLASH-relevant dose rates. The physics of laser-Compton x-ray scattering ensures that the scattered x-rays follow exactly the trajectory of… ▽ More The design and optimization of laser-Compton x-ray systems based on compact distributed charge accelerator structures can enable micron-scale imaging of disease and the concomitant production of beams of Very High Energy Electrons (VHEEs) capable of producing FLASH-relevant dose rates. The physics of laser-Compton x-ray scattering ensures that the scattered x-rays follow exactly the trajectory of the incident electrons, thus providing a route to image-guided, VHEE FLASH radiotherapy. The keys to a compact architecture capable of producing both laser-Compton x-rays and VHEEs are the use of X-band RF accelerator structures which have been demonstrated to operate with over 100 MeV/m acceleration gradients. The operation of these structures in a distributed charge mode in which each radiofrequency (RF) cycle of the drive RF pulse is filled with a low-charge, high-brightness electron bunch is enabled by the illumination of a high-brightness photogun with a train of UV laser pulses synchronized to the frequency of the underlying accelerator system. The UV pulse trains are created by a patented pulse synthesis approach which utilizes the RF clock of the accelerator to phase and amplitude modulate a narrow band continuous wave (CW) seed laser. In this way it is possible to produce up to 10 μみゅーA of average beam current from the accelerator. Such high current from a compact accelerator enables production of sufficient x-rays via laser-Compton scattering for clinical imaging and does so from a machine of "clinical" footprint. At the same time, the production of 1000 or greater individual micro-bunches per RF pulse enables > 10 nC of charge to be produced in a macrobunch of < 100 ns. The design, construction, and test of the 100-MeV class prototype system in Irvine, CA is also presented. △ Less

Submitted 7 August, 2024; originally announced August 2024.

Comments: 28 pages, 13 figures, 2 tables

Exact boundary-value solution for an electromagnetic wave propagating in a linearly-varying index of refraction

Authors: N. A. Lopez

Abstract: The propagation of electromagnetic waves in a linearly-varying index of refraction is a fundamental problem in wave physics, being relevant in fusion science for describing certain wave-based heating and diagnostic schemes. Here, an exact solution is obtained for a given incoming wavefield specified on the boundary transverse to the direction of inhomogeneity by performing a spectral, rather than… ▽ More The propagation of electromagnetic waves in a linearly-varying index of refraction is a fundamental problem in wave physics, being relevant in fusion science for describing certain wave-based heating and diagnostic schemes. Here, an exact solution is obtained for a given incoming wavefield specified on the boundary transverse to the direction of inhomogeneity by performing a spectral, rather than asymptotic, matching. Two case studies are then presented: a Gaussian beam at oblique incidence and a speckled wavefield at normal incidence. For the Gaussian beam, it is shown that when the waist $W$ is sufficiently large, oblique incidence manifests simply as rigid translation and focal shift of the corresponding diffraction pattern at normal incidence. The destruction of the hyperbolic umbilic caustic (corresponding to a critically focused beam) as $W$ is reduced is then demonstrated. The caustic disappears once $W \lesssim δでるた_a \sqrt{L}$ ($L$ being the medium lengthscale normalized by the Airy skin depth $δでるた_a$), at which point the wave behavior is increasingly described by Airy functions, but experiences less focusing as a result. To maximize the intensity of a launched Gaussian beam at a turning point, one should therefore minimize the imaginary part of the launched complex beam parameter while having the real part satisfy critical focusing. For the speckled wavefield, it is shown that the transverse speckle pattern only couples to the Airy longitudinal pattern when the coupling parameter $ηいーた= \sqrt{L}/f_\#$ is large, with $f_\#$ being the f-number of the launching aperture. When $ηいーた\ll 1$, a reduced description of the total wavefield can be obtained by simply multiplying the incoming speckle pattern with the Airy swelling. △ Less

Submitted 30 July, 2024; originally announced July 2024.

Comments: 16 pages, 8 figures, 2 appendices. Submission to special issue on contributions to EFTC 2023

Reconfigurable classifier based on spin torque driven magnetization switching in electrically connected magnetic tunnel junctions

Authors: A. Lopez, D. Costa, T. Bohnert, P. P. Freitas, R. Ferreira, I. Barbero, J. Camarero, C. Leon, J. Grollier, M. Romera

Abstract: A promising branch of neuromorphic computing aims to perform cognitive operations in hardware leveraging the physics of efficient and well-established nano-devices. In this work, we present a reconfigurable classifier based on a network of electrically connected magnetic tunnel junctions that categorizes information encoded in the amplitude of input currents through the spin torque driven magnetiz… ▽ More A promising branch of neuromorphic computing aims to perform cognitive operations in hardware leveraging the physics of efficient and well-established nano-devices. In this work, we present a reconfigurable classifier based on a network of electrically connected magnetic tunnel junctions that categorizes information encoded in the amplitude of input currents through the spin torque driven magnetization switching output configuration. The network can be trained to classify new data by adjusting additional programming currents applied selectively to the junctions. We experimentally demonstrate that a network composed of three magnetic tunnel junctions can learn to classify spoken vowels with a recognition rate that surpasses the performance of software multilayer neural networks with the same number of trained parameters in this task. These results, obtained with the same nano-devices and working principle employed in industrial spin-transfer torque magnetic random-access memories (STT-MRAM), constitute an important step towards the development of large-scale neuromorphic networks based on well-established technology. △ Less

Submitted 9 July, 2024; originally announced July 2024.

Comments: To be published in Physical Review Applied. 13 pages, 4 figures

The multi$-$physics analysis and design of CUSP, a two CubeSat constellation for Space Weather and Solar flares X-ray polarimetry

Authors: Giovanni Lombardi, Sergio Fabiani, Ettore Del Monte, Enrico Costa, Paolo Soffitta, Fabio Muleri, Ilaria Baffo, Marco E. Biancolini, Sergio Bonomo, Daniele Brienza, Riccardo Campana, Mauro Centrone, Gessica Contini, Giovanni Cucinella, Andrea Curatolo, Nicolas De Angelis, Giovanni De Cesare, Andrea Del Re, Sergio Di Cosimo, Simone Di Filippo, Alessandro Di Marco, Emanuele Di Meo, Giuseppe Di Persio, Immacolata Donnarumma, Pierluigi Fanelli , et al. (16 additional authors not shown)

Abstract: The CUbesat Solar Polarimeter (CUSP) project aims to develop a constellation of two CubeSats orbiting the Earth to measure the linear polarization of solar flares in the hard X-ray band by means of a Compton scattering polarimeter on board of each satellite. CUSP will allow to study the magnetic reconnection and particle acceleration in the flaring magnetic structures. CUSP is a project approved f… ▽ More The CUbesat Solar Polarimeter (CUSP) project aims to develop a constellation of two CubeSats orbiting the Earth to measure the linear polarization of solar flares in the hard X-ray band by means of a Compton scattering polarimeter on board of each satellite. CUSP will allow to study the magnetic reconnection and particle acceleration in the flaring magnetic structures. CUSP is a project approved for a Phase B study by the Italian Space Agency in the framework of the Alcor program aimed to develop CubeSat technologies and missions. In this paper we describe the a method for a multi-physical simulation analysis while analyzing some possible design optimization of the payload design solutions adopted. In particular, we report the mechanical design for each structural component, the results of static and dynamic finite element analysis, the preliminary thermo-mechanical analysis for two specific thermal cases (hot and cold orbit) and a topological optimization of the interface between the platform and the payload. △ Less

Submitted 4 July, 2024; originally announced July 2024.

Comments: Proceeding of SPIE Conference "Astronomical Telescopes+ Instrumentation", Yokohama (Japan), 16-21 June 2024

Report number: Paper no. 13093-307

The role of the thermal properties of electrons on the dispersion properties of Alfvén waves in space plasmas

Authors: Nicolás Villarroel-Sepúlveda, Pablo S. Moya, Rodrigo A. López, Daniel Verscharen

Abstract: Context. The transition from left-hand to right-hand polarized Alfvén waves depends on the wavenumber, the ratio of kinetic to magnetic pressure $βべーた$, temperature anisotropy, and ion composition of the plasma. Along with the temperature anisotropy, the electron-to-proton temperature ratio $T_e/T_p$ is of great relevance for the characterization of the thermal properties of a plasma. This ratio vari… ▽ More Context. The transition from left-hand to right-hand polarized Alfvén waves depends on the wavenumber, the ratio of kinetic to magnetic pressure $βべーた$, temperature anisotropy, and ion composition of the plasma. Along with the temperature anisotropy, the electron-to-proton temperature ratio $T_e/T_p$ is of great relevance for the characterization of the thermal properties of a plasma. This ratio varies significantly between different space plasma environments. Thus, studying how variations on this ratio affect the polarisation properties of electromagnetic waves becomes highly relevant for our understanding of the dynamics of space plasmas. Aim. We present an extensive study on the effect of the thermal properties of electrons on the behaviour and characteristics of Alfvénic waves in fully kinetic linear theory, as well as on the transition from electromagnetic ion-cyclotron (EMIC) to kinetic Alfvén waves (KAW). Method. We solve the fully kinetic dispersion relation for oblique electromagnetic waves of the Alfvén branch in a homogenous Maxwellian electron-proton plasma. We quantify the effect of the thermal properties of electrons by varying the electron-to-proton temperature ratio for different configurations of the propagation angle, $βべーた_p=8πぱいnkT_p/B^2$, and wavenumber. Results. We show that the temperature ratio $T_e/T_p$ has strong and non-trivial effects on the polarisation of the Alfvénic modes, especially at kinetic scales and $βべーた_e+βべーた_p>0.5$. We conclude that electron inertia plays an important role in the kinetic scale physics of the KAW in the warm plasma regime, and thus cannot be excluded in hybrid models for computer simulations. △ Less

Submitted 4 July, 2024; originally announced July 2024.

Comments: 13 pages, 7 figures

Report number: AA49471-24

Collisional whistler instability and electron temperature staircase in inhomogeneous plasma

Authors: N. A. Lopez, A. F. A. Bott, A. A. Schekochihin

Abstract: High-beta magnetized plasmas often exhibit anomalously structured temperature profiles, as seen from galaxy cluster observations and recent experiments. It is well known that when such plasmas are collisionless, temperature gradients along the magnetic field can excite whistler waves that efficiently scatter electrons to limit their heat transport. Only recently has it been shown that parallel tem… ▽ More High-beta magnetized plasmas often exhibit anomalously structured temperature profiles, as seen from galaxy cluster observations and recent experiments. It is well known that when such plasmas are collisionless, temperature gradients along the magnetic field can excite whistler waves that efficiently scatter electrons to limit their heat transport. Only recently has it been shown that parallel temperature gradients can excite whistler waves also in collisional plasmas. Here we develop a Wigner--Moyal theory for the collisional whistler instability starting from Braginskii-like fluid equations in a slab geometry. This formalism is necessary because, for a large region in parameter space, the fastest-growing whistler waves have wavelengths comparable to the background temperature gradients. We find additional damping terms in the expression for the instability growth rate involving inhomogeneous Nernst advection and resistivity. They (i) enable whistler waves to re-arrange the electron temperature profile via growth, propagation, and subsequent dissipation, and (ii) allow non-constant temperature profiles to exist stably. For high-beta plasmas, the marginally stable solutions take the form of a temperature staircase along the magnetic field lines. The electron heat flux can also be suppressed by the Ettingshausen effect when the whistler intensity profile is sufficiently peaked and oriented opposite the background temperature gradient. This mechanism allows cold fronts without magnetic draping, might reduce parallel heat losses in inertial fusion experiments, and generally demonstrates that whistler waves can regulate transport even in the collisional limit. △ Less

Submitted 16 June, 2024; originally announced June 2024.

Comments: 40 pages (31 pages main article), 15 figures, 7 appendices

Regarding the extension of metaplectic geometrical optics to modelling evanescent waves in ray-tracing codes

Authors: N. A. Lopez, R. Højlund, M. G. Senstius

Abstract: Metaplectic geometrical optics (MGO) is a recently developed ray-tracing framework to accurately compute the wavefield behavior near a caustic (turning point or focal point), where traditional ray-tracing breaks down. However, MGO has thus far been restricted to having real-valued wavevectors. This is disadvantageous because often upon crossing a caustic from the `illuminated' region to the `shado… ▽ More Metaplectic geometrical optics (MGO) is a recently developed ray-tracing framework to accurately compute the wavefield behavior near a caustic (turning point or focal point), where traditional ray-tracing breaks down. However, MGO has thus far been restricted to having real-valued wavevectors. This is disadvantageous because often upon crossing a caustic from the `illuminated' region to the `shadow' region, two real-valued rays coalesce into one complex-valued ray corresponding to the transition from propagating to evanescent behavior. One can distinguish caustics as having either `illuminated shadows' or `proper shadows' -- the former corresponds to when the shadow still contains real-valued rays (albeit in a fewer quantity than in the illuminated region), while the latter corresponds to when the shadow contains no real-valued rays. Here, by means of examples, we show how MGO can be used to model both types of shadows. First, for illuminated shadows we show that MGO can actually be used `as is', provided a corrected quadrature-angle bias is used compared to that proposed in the original references. This is then implemented and demonstrated in a recently developed MGO ray-tracing code. Second, we show that for proper shadows, the MGO formalism can still be used if the symplectic rotation matrix that removes caustics along rays is allowed to be complex-valued. In both cases, strong agreement is seen between the MGO and the exact solution, demonstrating the potential of MGO for improving the predictive capability of ray-tracing codes and laying the foundations for modeling more complicated evanescent phenomena such as tunnelling with MGO. △ Less

Submitted 6 June, 2024; v1 submitted 3 June, 2024; originally announced June 2024.

Comments: 10 pages, 5 figures, 2 appendices. Updated funding acknowledgement

Experimental observations of bifurcated power decay lengths in the near Scrape-Off Layer of tokamak plasmas

Authors: X Zhang, C Marsden, M Moscheni, E Maartensson, A Rengle, M Robinson, T OGorman, H F Lowe, E Vekshina, S Janhunen, P F Buxton, M Sertoli, M Romanelli, A Scarabosio, T K Gray, N A Lopez, the ST40 Team

Abstract: The scrape-off layer parallel heat flux decay lengths measured at ST40, a high field, low aspect ratio spherical tokamak, have been observed to bifurcate into two groups. The wide group matches closely with the scale of ion poloidal Larmour radius and follows existing H-mode scalings, while the narrow group falls up to 10 times below scalings, on the scale of ion total Larmour radius. The onset of… ▽ More The scrape-off layer parallel heat flux decay lengths measured at ST40, a high field, low aspect ratio spherical tokamak, have been observed to bifurcate into two groups. The wide group matches closely with the scale of ion poloidal Larmour radius and follows existing H-mode scalings, while the narrow group falls up to 10 times below scalings, on the scale of ion total Larmour radius. The onset of the narrow scrape-off layer width is observed to be associated with suppressed magnetic fluctuations, suggesting reduced electromagnetic turbulence levels in the SOL. △ Less

Submitted 23 May, 2024; originally announced May 2024.

Comments: 6 pages, 4 figures

Fully Kinetic Simulations of Proton-Beam-Driven Instabilities from Parker Solar Probe Observations

Authors: Luca Pezzini, Andrei N. Zhukov, Fabio Bacchini, Giuseppe Arrò, Rodrigo A. López, Alfredo Micera, Maria Elena Innocenti, Giovanni Lapenta

Abstract: The expanding solar wind plasma ubiquitously exhibits anisotropic non-thermal particle velocity distributions. Typically, proton Velocity Distribution Functions (VDFs) show the presence of a core and a field-aligned beam. Novel observations made by Parker Solar Probe (PSP) in the innermost heliosphere have revealed new complex features in the proton VDFs, namely anisotropic beams that sometimes ex… ▽ More The expanding solar wind plasma ubiquitously exhibits anisotropic non-thermal particle velocity distributions. Typically, proton Velocity Distribution Functions (VDFs) show the presence of a core and a field-aligned beam. Novel observations made by Parker Solar Probe (PSP) in the innermost heliosphere have revealed new complex features in the proton VDFs, namely anisotropic beams that sometimes experience perpendicular diffusion. This phenomenon gives rise to VDFs that resemble a "hammerhead". In this study, we use a 2.5D fully kinetic simulation to investigate the stability of proton VDFs with anisotropic beams observed by PSP. Our setup consists of a core and an anisotropic beam populations that drift with respect to each other. This configuration triggers a proton-beam instability from which nearly parallel fast magnetosonic modes develop. Our results demonstrate that before this instability reaches saturation, the waves resonantly interact with the beam protons, causing significant perpendicular heating at the expense of the parallel temperature. Furthermore, the proton perpendicular heating induces a hammerhead-like shape in the resulting VDF. Our results suggest that this mechanism probably contribute to producing the observed hammerhead distributions. △ Less

Submitted 13 May, 2024; originally announced May 2024.

Electro-phononic and magneto-phononic frequency conversion

Authors: Daniel A. Bustamante Lopez, Wanzheng Hu, Dominik M. Juraschek

Abstract: Nonlinear frequency conversion by optical rectification, as well as difference- and sum-frequency generation are fundamental processes for producing electromagnetic radiation at different frequencies. Here, we demonstrate that coherently excited infrared-active phonons can be used as transducers for generating nonlinear electric polarizations and magnetizations via phonon-phonon and phonon-magnon… ▽ More Nonlinear frequency conversion by optical rectification, as well as difference- and sum-frequency generation are fundamental processes for producing electromagnetic radiation at different frequencies. Here, we demonstrate that coherently excited infrared-active phonons can be used as transducers for generating nonlinear electric polarizations and magnetizations via phonon-phonon and phonon-magnon interactions, in a way similar to nonlinear optical frequency conversion. We derive analytical solutions for the time-dependent polarizations and magnetizations for the second-order response to the electric field component of an ultrashort laser pulse. These allow us to define second-order nonlinear electric and magneto-electric susceptibilities that capture the rectification, as well as the impulsive and sum-frequency excitation of coherent phonons and magnons. Our theoretical framework naturally incorporates existing mechanisms and further leads to the prediction of a hybrid magneto-opto-phononic inverse Faraday effect involving photon-phonon-magnon scattering. Our work demonstrates nonlinear phononics as a pathway to controlling the electric polarization and magnetization in solids. △ Less

Submitted 30 April, 2024; originally announced April 2024.

Nano-Electromagnetic Super-dephasing in Collective Atom-Atom Interactions

Authors: Wenbo Sun, Adrian E. Rubio López, Zubin Jacob

Abstract: Pure dephasing and spontaneous emission are two non-unitary processes of atoms or spins interacting with fluctuating electromagnetic (EM) modes. Collective spontaneous emission (e.g., superradiance) originates from interactions with EM modes in resonance with atoms and has received considerable attention. Meanwhile, the analogous collective dephasing phenomena remain poorly understood. Here, we in… ▽ More Pure dephasing and spontaneous emission are two non-unitary processes of atoms or spins interacting with fluctuating electromagnetic (EM) modes. Collective spontaneous emission (e.g., superradiance) originates from interactions with EM modes in resonance with atoms and has received considerable attention. Meanwhile, the analogous collective dephasing phenomena remain poorly understood. Here, we introduce the nano-EM super-dephasing phenomenon arising in the photonic environment near lossy material interfaces. We show that this effect is enhanced by over 10 orders of magnitude compared to free space or photonic cavities due to the presence of long-range correlations in low-frequency evanescent EM fluctuations. We unravel the universality of nano-EM super-dephasing behaviors near ferrimagnets, metals, and superconductors and their dependence on low-frequency material properties. We demonstrate that the scaling of nano-EM super-dephasing is independent of EM modes' wavelengths and differs from the conventional $N^2$ scaling of superradiance by analyzing the decoherence of entangled states, including GHZ states. Finally, we show how to experimentally isolate and control super-dephasing to open interesting frontiers for scalable quantum systems. △ Less

Submitted 28 February, 2024; originally announced February 2024.

Comments: 6+17 pages, 3+3 figures

Critical mobility in policy making for epidemic containment

Authors: Jesús A. Moreno López, Sandro Meloni, Jose J. Ramasco

Abstract: When considering airborne epidemic spreading in social systems, a natural connection arises between mobility and epidemic contacts. As individuals travel, possibilities to encounter new people either at the final destination or during the transportation process appear. Such contacts can lead to new contagion events. In fact, mobility has been a crucial target for early non-pharmaceutical containme… ▽ More When considering airborne epidemic spreading in social systems, a natural connection arises between mobility and epidemic contacts. As individuals travel, possibilities to encounter new people either at the final destination or during the transportation process appear. Such contacts can lead to new contagion events. In fact, mobility has been a crucial target for early non-pharmaceutical containment measures against the recent COVID-19 pandemic, with a degree of intensity ranging from public transportation line closures to regional, city or even home confinements. Nonetheless, quantitative knowledge on the relationship between mobility-contagions and, consequently, on the efficiency of containment measures remains elusive. Here we introduce an agent-based model with a simple interaction between mobility and contacts. Despite its simplicity our model shows the emergence of a critical mobility level, inducing major outbreaks when surpassed. We explore the interplay between mobility restrictions and the infection in recent intervention policies seen across many countries, and how interventions in the form of closures triggered by incidence rates can guide the epidemic into an oscillatory regime with recurrent waves. We consider how the different interventions impact societal well-being, the economy and the population. Finally, we propose a mitigation framework based on the critical nature of mobility in an epidemic, able to suppress incidence and oscillations at will, preventing extreme incidence peaks with potential to saturate health care resources. △ Less

Submitted 9 May, 2024; v1 submitted 8 February, 2024; originally announced February 2024.

Comments: 13 pages, 5 figures

Modified Mach Zehnder interferometer for spatial coherence measurement

Authors: FJ Torcal-Milla, J Lobera, EM Roche, AM Lopez, V Palero, N Andres, P Arroyo

Abstract: Spatial coherence of light sources is usually obtained by using the classical Youngs interferometer. Despite the original experiment has been improved in successive works, some drawbacks still remain. In addition, several pairs of points must be used to characterize the totally the 2D complex-coherence degree of the source. In this work, an alternative based on a modified Mach-Zehnder interferomet… ▽ More Spatial coherence of light sources is usually obtained by using the classical Youngs interferometer. Despite the original experiment has been improved in successive works, some drawbacks still remain. In addition, several pairs of points must be used to characterize the totally the 2D complex-coherence degree of the source. In this work, an alternative based on a modified Mach-Zehnder interferometer which includes a pair of lenses is presented. With this modified Mach-Zehnder interferometer, we are able to measure the spatial coherence length between any two points of the laser beam section simultaneously. The set-up does not have any movable part, which makes it robust and portable. To test it, the two-dimensional spatial coherence of a high-speed laser with two cavities has been measured for different pulse energy values. We observe from the experimental measurements that the complex degree of coherence changes with the selected output energy. Both laser cavities seem to have similar complex coherence degree for maximum energy although it is not symmetrical. Thus, this analysis will allow us to determine the best configuration of the double cavity laser for interferometric applications. Furthermore, the proposed approach can be applied to any other light sources. △ Less

Submitted 16 January, 2024; originally announced January 2024.

Quantum-like behavior of an active particle in a double-well potential

Authors: Rahil N. Valani, Álvaro G. López

Abstract: A macroscopic, self-propelled wave-particle entity (WPE) that emerges as a walking droplet on the surface of a vibrating liquid bath exhibits several hydrodynamic quantum analogs. We explore the rich dynamical and quantum-like features emerging in a model of an idealized one-dimensional wave-particle entity in a double-well potential. The integro-differential equation of motion for the WPE transfo… ▽ More A macroscopic, self-propelled wave-particle entity (WPE) that emerges as a walking droplet on the surface of a vibrating liquid bath exhibits several hydrodynamic quantum analogs. We explore the rich dynamical and quantum-like features emerging in a model of an idealized one-dimensional wave-particle entity in a double-well potential. The integro-differential equation of motion for the WPE transforms to a Lorenz-like system, which we explore in detail. We observe the analog of quantized eigenstates as discrete limit cycles that arise by varying the width of the double-well potential, and also in the form of multistability with coexisting limit cycles. These states show narrow as well as wide energy level splitting. Tunneling-like behavior is also observed where the WPE erratically transitions between the two wells of the double-well potential. We rationalize this phenomena in terms of crisis-induced intermittency. Further, we discover a fractal structure in the escape time distribution of the particle from a well based on initial conditions, indicating unpredictability of this tunneling-like intermittent behavior at all scales. The chaotic intermittent dynamics lead to wave-like emergent features in the probability distribution of particle's position that show qualitative similarity with its quantum counterpart. Lastly, rich dynamical features are also observed such as a period doubling route to chaos as well as self-similar periodic islands in the chaotic parameter set. △ Less

Submitted 17 June, 2024; v1 submitted 10 January, 2024; originally announced January 2024.

Comments: 15 pages, 12 figures

Microfluidics for Hydrodynamics Investigations of Sand Dollar Larvae

Authors: Wesley A. Chen, Bryant A. Lopez, Haley B. Obenshain, Moses Villeda, Brian T. Le, Brenda AAB. Ametepe, Ariana Lee, Douglas A. Pace, Siavash Ahrar

Abstract: The life cycle of most marine invertebrates includes a planktonic larval stage before metamorphosis to bottom-dwelling adulthood. During larval stage, ciliary-mediated activity enables feeding (capture unicellular algae) and transport of materials (oxygen) required for the larva's growth, development, and successful metamorphosis. Investigating the underlying hydrodynamics of these behaviors is va… ▽ More The life cycle of most marine invertebrates includes a planktonic larval stage before metamorphosis to bottom-dwelling adulthood. During larval stage, ciliary-mediated activity enables feeding (capture unicellular algae) and transport of materials (oxygen) required for the larva's growth, development, and successful metamorphosis. Investigating the underlying hydrodynamics of these behaviors is valuable for addressing fundamental biological questions (e.g., phenotypic plasticity) and advancing engineering applications. In this work, we combined microfluidics and fluorescence microscopy as a miniaturized PIV (mPIV) to study ciliary-medicated hydrodynamics during suspension feeding in sand dollar larvae (Dendraster excentricus). First, we confirmed the approach's feasibility by examining the underlying hydrodynamics (vortex patterns) for low- and high-fed larvae. Next, ciliary hydrodynamics were tracked from 11 days post-fertilization (DPF) to 20 DPF for 21 low-fed larvae. Microfluidics enabled the examination of baseline activities (without external flow) and behaviors in the presence of environmental cues (external flow). A library of qualitative vortex patterns and quantitative hydrodynamics was generated and shared as a stand alone repository. Results from mPIV (velocities) were used to examine the role of ciliary activity in transporting materials (oxygen). Given the laminar flow and the viscosity-dominated environments surrounding the larvae, overcoming the diffusive boundary layer is critical for the organism's survival. Peclet number analysis for oxygen transport suggested that ciliary velocities help overcome the diffusion dominated transport (max Pe numbers between 30-60). Microfluidics serving as mPIV provided a scalable and accessible approach for investigating the ciliary hydrodynamics of marine organisms. △ Less

Submitted 29 December, 2023; originally announced January 2024.

Comments: 21 pages and 11 figures (videos not included)

Optimizing carbon dioxide trapping for geological storage

Authors: Jaione Tirapu Azpiroz, Ronaldo Giro, Rodrigo Neumann Barros Ferreira, Marcio Nogueira Pereira da Silva, Manuela Fernandes Blanco Rodriguez, Adolfo E. Correa Lopez, David A. Lazo Vasquez, Matheus Esteves Ferreira, Mariana Del Grande, Ademir Ferreira Da Silva, Mathias B. Steiner

Abstract: Carbon dioxide (CO2) trapping in capillary networks of reservoir rocks is a pathway to long-term geological storage. At pore scale, the CO2 trapping potential depends on injection pressure, temperature, and the rock's interaction with the surrounding fluids. Modeling this interaction requires adequate representations of both capillary volume and surface. For the lack of scalable representations, h… ▽ More Carbon dioxide (CO2) trapping in capillary networks of reservoir rocks is a pathway to long-term geological storage. At pore scale, the CO2 trapping potential depends on injection pressure, temperature, and the rock's interaction with the surrounding fluids. Modeling this interaction requires adequate representations of both capillary volume and surface. For the lack of scalable representations, however, the prediction of a rock's CO2 storage potential has been challenging. Here, we report how to represent a rock's pore space by statistically sampled capillary networks (ssCN) that preserve morphological rock characteristics. We have used the ssCN method to simulate CO2 drainage within a representative sandstone sample at reservoir pressures and temperatures, exploring intermediate- and CO2-wet conditions. This wetting regime is often neglected, despite evidence of plausibility. By raising pressure and temperature we observe increasing CO2 penetration within the capillary network. For contact angles approaching 90 degrees, the CO2 saturation exhibits a pronounced maximum reaching 80 percent of the accessible pore volume. This is about twice as high as the saturation values reported previously. For enabling validation of our results and a broader application of our methodology, we have made available the rock tomography data, the digital rock computational workflows, and the ssCN models used in this study. △ Less

Submitted 20 December, 2023; originally announced December 2023.

Comments: Main article: 16 pages and 3 figures. Supplementary Information: 17 pages, 3 equations, 13 figures and 3 tables. Includes DOI for accompanying digital rock data, URL of github repository for flow simulator code, and of accompanying digital rock processing python code and automation of the scientific workflows

Accurate and Fast Fischer-Tropsch Reaction Microkinetics using PINNs

Authors: Harshil Patel, Aniruddha Panda, Tymofii Nikolaienko, Stanislav Jaso, Alejandro Lopez, Kaushic Kalyanaraman

Abstract: Microkinetics allows detailed modelling of chemical transformations occurring in many industrially relevant reactions. Traditional way of solving the microkinetics model for Fischer-Tropsch synthesis (FTS) becomes inefficient when it comes to more advanced real-time applications. In this work, we address these challenges by using physics-informed neural networks(PINNs) for modelling FTS microkinet… ▽ More Microkinetics allows detailed modelling of chemical transformations occurring in many industrially relevant reactions. Traditional way of solving the microkinetics model for Fischer-Tropsch synthesis (FTS) becomes inefficient when it comes to more advanced real-time applications. In this work, we address these challenges by using physics-informed neural networks(PINNs) for modelling FTS microkinetics. We propose a computationally efficient and accurate method, enabling the ultra-fast solution of the existing microkinetics models in realistic process conditions. The proposed PINN model computes the fraction of vacant catalytic sites, a key quantity in FTS microkinetics, with median relative error (MRE) of 0.03%, and the FTS product formation rates with MRE of 0.1%. Compared to conventional equation solvers, the model achieves up to 1E+06 times speed-up when running on GPUs, thus being fast enough for multi-scale and multi-physics reactor modelling and enabling its applications in real-time process control and optimization. △ Less

Submitted 17 November, 2023; originally announced November 2023.

Neural ODEs as a discovery tool to characterize the structure of the hot galactic wind of M82

Authors: Dustin D. Nguyen, Yuan-Sen Ting, Todd A. Thompson, Sebastian Lopez, Laura A. Lopez

Abstract: Dynamic astrophysical phenomena are predominantly described by differential equations, yet our understanding of these systems is constrained by our incomplete grasp of non-linear physics and scarcity of comprehensive datasets. As such, advancing techniques in solving non-linear inverse problems becomes pivotal to addressing numerous outstanding questions in the field. In particular, modeling hot g… ▽ More Dynamic astrophysical phenomena are predominantly described by differential equations, yet our understanding of these systems is constrained by our incomplete grasp of non-linear physics and scarcity of comprehensive datasets. As such, advancing techniques in solving non-linear inverse problems becomes pivotal to addressing numerous outstanding questions in the field. In particular, modeling hot galactic winds is difficult because of unknown structure for various physical terms, and the lack of \textit{any} kinematic observational data. Additionally, the flow equations contain singularities that lead to numerical instability, making parameter sweeps non-trivial. We leverage differentiable programming, which enables neural networks to be embedded as individual terms within the governing coupled ordinary differential equations (ODEs), and show that this method can adeptly learn hidden physics. We robustly discern the structure of a mass-loading function which captures the physical effects of cloud destruction and entrainment into the hot superwind. Within a supervised learning framework, we formulate our loss function anchored on the astrophysical entropy ($K \propto P/ρろー^{5/3}$). Our results demonstrate the efficacy of this approach, even in the absence of kinematic data $v$. We then apply these models to real Chandra X-Ray observations of starburst galaxy M82, providing the first systematic description of mass-loading within the superwind. This work further highlights neural ODEs as a useful discovery tool with mechanistic interpretability in non-linear inverse problems. We make our code public at this GitHub repository (https://github.com/dustindnguyen/2023_NeurIPS_NeuralODEs_M82). △ Less

Submitted 28 November, 2023; v1 submitted 3 November, 2023; originally announced November 2023.

Comments: 9 Pages, 2 Figures, Accepted at the NeurIPS 2023 workshop on Machine Learning and the Physical Sciences

New solution to Airy's equation for modeling beams near turning points

Authors: N. A. Lopez

Abstract: Here I present a summarised description of the general behavior of waves near turning points. I show that when asymptotic matching to a prescribed incident wavefield is performed, the solution changes from being the Airy function to being more generally the hyperbolic umbilic function. This solution accounts for (de-)focusing of the incident wavefield due to gradient-index lensing in addition to t… ▽ More Here I present a summarised description of the general behavior of waves near turning points. I show that when asymptotic matching to a prescribed incident wavefield is performed, the solution changes from being the Airy function to being more generally the hyperbolic umbilic function. This solution accounts for (de-)focusing of the incident wavefield due to gradient-index lensing in addition to the familiar Airy physics, which importantly do not superimpose. Possible applications to fusion research are discussed. △ Less

Submitted 26 September, 2023; originally announced September 2023.

Comments: Proceeding for 2023 EPS plasma physics conference

Photon discerner: Adaptive quantum optical sensing near the shot noise limit

Authors: F. Bao, L. Bauer, A. E. Rubio Lopez, Z. Jacob

Abstract: Photon statistics of an optical field can be used for quantum optical sensing in low light level scenarios free of bulky optical components. However, photon-number-resolving detection to unravel the photon statistics is challenging. Here, we propose a novel detection approach, that we call `photon discerning', which uses adaptive photon thresholding for photon statistical estimation without record… ▽ More Photon statistics of an optical field can be used for quantum optical sensing in low light level scenarios free of bulky optical components. However, photon-number-resolving detection to unravel the photon statistics is challenging. Here, we propose a novel detection approach, that we call `photon discerning', which uses adaptive photon thresholding for photon statistical estimation without recording exact photon numbers. Our photon discerner is motivated by the field of neural networks where tunable thresholds have proven efficient for isolating optimal decision boundaries in machine learning tasks. The photon discerner maximizes Fisher information per photon by iteratively choosing the optimal threshold in real-time to approach the shot noise limit. Our proposed scheme of adaptive photon thresholding leads to unique remote-sensing applications of quantum DoLP (degree of linear polarization) camera and quantum LiDAR. We investigate optimal thresholds and show that the optimal photon threshold can be counter-intuitive (not equal to 1) even for weak signals (mean photon number much less than 1), due to the photon bunching effect. We also put forth a superconducting nanowire realization of the photon discerner which can be experimentally implemented in the near-term. We show that the adaptivity of our photon discerner enables it to beat realistic photon-number-resolving detectors with limited photon-number resolution. Our work suggests a new class of detectors for information-theory driven, compact, and learning-based quantum optical sensing. △ Less

Submitted 27 July, 2023; originally announced July 2023.

Comments: 12 pages, 9 figures

The effect of heavy ions on the dispersion properties of kinetic Alfvén waves in astrophysical plasmas

Authors: Nicolás Villarroel-Sepúlveda, Rodrigo A. López, Pablo S. Moya

Abstract: Context. Spacecraft measurements have shown Kinetic Alfvén Waves propagating in the terrestrial magnetosphere at lower wavenormal angles than predicted by linear Vlasov theory of electron-proton plasmas. To explain these observations, it has been suggested that the abundant heavy ion populations in this region may have strong, non-trivial effects that allow Alfvénic waves to acquire right-handed p… ▽ More Context. Spacecraft measurements have shown Kinetic Alfvén Waves propagating in the terrestrial magnetosphere at lower wavenormal angles than predicted by linear Vlasov theory of electron-proton plasmas. To explain these observations, it has been suggested that the abundant heavy ion populations in this region may have strong, non-trivial effects that allow Alfvénic waves to acquire right-handed polarization at lower angles with respect to the background magnetic field, as in the case of typical electron-proton plasma. Aims. We study the dispersion properties of Alfvénic waves in plasmas with stationary phase-space distribution functions with different heavy ion populations. Our extensive numerical analysis has allowed us to quantify the role of the heavy ion components on the transition from the left-hand polarized electromagnetic ion-cyclotron (EMIC) mode to the right-hand polarized kinetic Alfvén wave (KAW) mode. Methods. We used linear Vlasov-Maxwell theory to obtain the dispersion relation for oblique electromagnetic waves. The dispersion relation of Alfvén waves was obtained numerically by considering four different oxygen ion concentrations ranging between 0.0 and 0.2 for all propagation angles, as a function of both the wavenumber and the plasma beta parameter. Results. The inclusion of the heavy O+ ions is found to considerably reduce the transition angle from EMIC to KAW both as a function of the wave number and plasma beta. With increasing O+ concentrations, waves become more damped in specific wavenumber regions. However, the inclusion of oxygen ions may allow weakly damped KAW to effectively propagate at smaller wave-normal angles than in the electron-proton case, as suggested by observations. △ Less

Submitted 2 June, 2023; originally announced June 2023.

Comments: 11 pages, 8 figures

Journal ref: A&A 675, A84 (2023)

The LHCb upgrade I

Authors: LHCb collaboration, R. Aaij, A. S. W. Abdelmotteleb, C. Abellan Beteta, F. Abudinén, C. Achard, T. Ackernley, B. Adeva, M. Adinolfi, P. Adlarson, H. Afsharnia, C. Agapopoulou, C. A. Aidala, Z. Ajaltouni, S. Akar, K. Akiba, P. Albicocco, J. Albrecht, F. Alessio, M. Alexander, A. Alfonso Albero, Z. Aliouche, P. Alvarez Cartelle, R. Amalric, S. Amato , et al. (1298 additional authors not shown)

Abstract: The LHCb upgrade represents a major change of the experiment. The detectors have been almost completely renewed to allow running at an instantaneous luminosity five times larger than that of the previous running periods. Readout of all detectors into an all-software trigger is central to the new design, facilitating the reconstruction of events at the maximum LHC interaction rate, and their select… ▽ More The LHCb upgrade represents a major change of the experiment. The detectors have been almost completely renewed to allow running at an instantaneous luminosity five times larger than that of the previous running periods. Readout of all detectors into an all-software trigger is central to the new design, facilitating the reconstruction of events at the maximum LHC interaction rate, and their selection in real time. The experiment's tracking system has been completely upgraded with a new pixel vertex detector, a silicon tracker upstream of the dipole magnet and three scintillating fibre tracking stations downstream of the magnet. The whole photon detection system of the RICH detectors has been renewed and the readout electronics of the calorimeter and muon systems have been fully overhauled. The first stage of the all-software trigger is implemented on a GPU farm. The output of the trigger provides a combination of totally reconstructed physics objects, such as tracks and vertices, ready for final analysis, and of entire events which need further offline reprocessing. This scheme required a complete revision of the computing model and rewriting of the experiment's software. △ Less

Submitted 17 May, 2023; originally announced May 2023.

Comments: All figures and tables, along with any supplementary material and additional information, are available at http://lhcbproject.web.cern.ch/lhcbproject/Publications/LHCbProjectPublic/LHCb-DP-2022-002.html (LHCb public pages)

Report number: LHCb-DP-2022-002

Salt-rejecting continuous passive solar thermal desalination via convective flow and thin-film condensation

Authors: Patrick I. Babb, S. Farzad Ahmadi, Forrest Brent, Ruby Gans, Mabel Aceves Lopez, Jiuxu Song, Qixian Wang, Brandon Zou, Xiangying Zuo, Amanda Strom, Jaya Nolt, Tyler Susko, Kirk Fields, Yangying Zhu

Abstract: Passive solar desalination is an emerging low-cost technology for fresh water production. State of the art desalinators typically evaporate water using wicking structures to achieve high solar-to-vapor efficiency by minimizing heat loss. However, wicking structures cannot reject salt continuously which limits the operating duration of the desalinators to several hours before the devices are turned… ▽ More Passive solar desalination is an emerging low-cost technology for fresh water production. State of the art desalinators typically evaporate water using wicking structures to achieve high solar-to-vapor efficiency by minimizing heat loss. However, wicking structures cannot reject salt continuously which limits the operating duration of the desalinators to several hours before the devices are turned off to reject salt. While significant research has focused on developing efficient evaporators to achieve high solar-to-vapor efficiency, inefficient condensers have become the bottleneck for the overall solar-to-water efficiency. To overcome these challenges, we designed a passive inverted single stage solar membrane desalinator that achieves continuous desalination and salt rejection. By flowing salt water on a radiative absorbing, porous, hydrophobic evaporator membrane using gravity, salt continuously diffuses away from the membrane while allowing heated water vapor to transport to and condense on a cooler microporous membrane below. Our design utilizes thin-film condensation on a microporous membrane which offers ample three-phase contact region to enhance condensation phase change heat transfer. By condensing within the microporous membrane, we reduce the gap distance between the condenser and evaporator membranes, which reduces the vapor transport resistance. We experimentally demonstrated a record-high continuous desalination and salt rejection test duration of 7 days under one-sun. Despite an increased convection heat loss necessary for salt rejection on the evaporator, our desalinator still achieved a water-collection rate of 0.487 $kg$ $m^{-2}h^{-1}$, which corresponds to a 32.2% solar-to-water efficiency. This work signifies an improvement in the robustness of current state of the art desalinators and presents a new architecture to further optimize passive solar desalinators. △ Less

Submitted 16 May, 2023; originally announced May 2023.

Comments: 24 pages, 9 figures

On the intensity of focused waves near turning points

Authors: N. A. Lopez, E. Kur, D. J. Strozzi

Abstract: A wave near an isolated turning point is typically assumed to have an Airy function profile with respect to the separation distance. This description is incomplete, however, and is insufficient to describe the behavior of more realistic wavefields that are not simple plane waves. Asymptotic matching to a prescribed incoming wavefield generically introduces a phasefront curvature term that changes… ▽ More A wave near an isolated turning point is typically assumed to have an Airy function profile with respect to the separation distance. This description is incomplete, however, and is insufficient to describe the behavior of more realistic wavefields that are not simple plane waves. Asymptotic matching to a prescribed incoming wavefield generically introduces a phasefront curvature term that changes the characteristic wave behavior from the Airy function to that of the hyperbolic umbilic function. This function, which is one of the seven classic 'elementary' functions from catastrophe theory along with the Airy function, can be understood intuitively as the solution for a linearly focused Gaussian beam propagating in a linearly varying density profile, as we show. The morphology of the caustic lines that govern the intensity maxima of the diffraction pattern as one alters the density lengthscale of the plasma, the focal length of the incident beam, and also the injection angle of the incident beam are presented in detail. This morphology includes a Goos-Hänchen shift and focal shift at oblique incidence that do not appear in a reduced ray-based description of the caustic. The enhancement of the intensity swelling factor for a focused wave compared to the typical Airy solution is highlighted, and the impact of finite lens aperture is discussed. Collisional damping and finite beam waist are included in the model and appear as complex components to the arguments of the hyperbolic umbilic function. The observations presented here on the behavior of waves near turning points should aid the development of improved reduced wave models to be used, for example, in designing modern nuclear fusion experiments. △ Less

Submitted 30 January, 2023; originally announced January 2023.

Comments: 12 pages, 7 figures, 3 appendices

Journal ref: Phys. Rev. E 107, 055204 (2023)

The aperiodic firehose instability of counter-beaming electrons in space plasmas

Authors: M. Lazar, R. A. López, P. S. Moya, S. Poedts, S. M. Shaaban

Abstract: Recent studies have revealed new unstable regimes of the counter-beaming electrons specific to hot and dilute plasmas from astrophysical scenarios. The (counter-)beaming electron firehose instability (BEFI) is induced for highly oblique angles of propagation relative to the magnetic field, resembling the fast growing and aperiodic mode triggered by the temperature anisotropy. It is investigated he… ▽ More Recent studies have revealed new unstable regimes of the counter-beaming electrons specific to hot and dilute plasmas from astrophysical scenarios. The (counter-)beaming electron firehose instability (BEFI) is induced for highly oblique angles of propagation relative to the magnetic field, resembling the fast growing and aperiodic mode triggered by the temperature anisotropy. It is investigated here for space plasma conditions that includes the influence of an embedding background plasma of electrons and protons. Kinetic theory is applied to prescribe the unstable regimes, and differentiate from the regimes of interplay with other instabilities. Linear theory predicts a systematic inhibition of the BEFI, by reducing the growth rates and the range of unstable wave-number with increasing the relative density of the background electrons. To obtain finite growth rates, the beam speed does not need to be high (just comparable to thermal speed), but beams must be dense enough, with a relative density at least 15-20\% of the total density. The plasma conditions favorable to this instability are reduced under the influence of background electrons. PIC simulations confirm not only that BEFI can be excited in the presence of background electrons, but also the inhibiting effect of this population. In the regimes of transition to electrostatic (ES) instabilities, BEFI is still robust enough to develop as a secondary instability, after the relaxation of beams under a quick interaction with ES fluctuations. BEFI resembles the properties of firehose heat-flux instability triggered by the electron strahl. However, BEFI is driven by a double (counter-beaming) strahl, and develops at oblique angles, which makes it effective in the regularization of the electron counter-beams observed in closed magnetic field topologies and interplanetary shocks. △ Less

Submitted 13 December, 2022; originally announced December 2022.

Comments: Paper accepted for publication in Astronomy and Astrophysics (12.12.2022)

Journal ref: A&A 670, A85 (2023)

Field-assisted birefringent Compton scattering

Authors: N. Ahmadiniaz, T. E. Cowan, M. Ding, M. A. Lopez Lopez, R. Sauerbrey, R. Shaisultanov, R. Schützhold

Abstract: Motivated by experimental initiatives such as the Helmholtz International Beamline for Extreme Fields (HIBEF), we study Compton scattering of x-rays at electrons in a strong external field (e.g., a strong optical laser) with special emphasis on the polarization-changing (i.e., birefringent) contribution on the amplitude level. Apart from being a potential background process for the planned vacuum… ▽ More Motivated by experimental initiatives such as the Helmholtz International Beamline for Extreme Fields (HIBEF), we study Compton scattering of x-rays at electrons in a strong external field (e.g., a strong optical laser) with special emphasis on the polarization-changing (i.e., birefringent) contribution on the amplitude level. Apart from being a potential background process for the planned vacuum birefringence experiments, this effect could be used for diagnostic purposes. Since the birefringent signal from free electrons (i.e., without the external field) vanishes in forward direction, the ratio of the birefringent and the normal (polarization conserving) contribution yields information about the field strength at the interaction point. △ Less

Submitted 6 December, 2022; originally announced December 2022.

Comments: 6 pages, 2 figures

Metaplectic Geometrical Optics

Authors: Nicolas A. Lopez

Abstract: Ray optics is an intuitive and computationally efficient model for wave propagation through nonuniform media. However, the underlying geometrical-optics (GO) approximation of ray optics breaks down at caustics, erroneously predicting the wave intensity to be infinite and thereby limiting the predictive capabilities of GO-based codes. Full-wave modeling can be used instead, but the added computatio… ▽ More Ray optics is an intuitive and computationally efficient model for wave propagation through nonuniform media. However, the underlying geometrical-optics (GO) approximation of ray optics breaks down at caustics, erroneously predicting the wave intensity to be infinite and thereby limiting the predictive capabilities of GO-based codes. Full-wave modeling can be used instead, but the added computational cost brings its own set of tradeoffs. Developing more efficient caustic remedies has therefore been an active area of research for the past few decades. In this thesis, I present a new ray-based approach called 'metaplectic geometrical optics' (MGO) that can be applied to any linear wave equation. Instead of evolving waves in the usual x (coordinate) or k (spectral) representation, MGO uses a mixed X=Ax+B representation. By continuously adjusting the coefficients A and B along the rays via sequenced metaplectic transforms (MTs) of the wavefield, corresponding to symplectic transformations of the ray phase space, one can ensure that GO remains valid in the X coordinates without caustic singularities. The caustic-free result is then mapped back onto the original x space using metaplectic transforms, as demonstrated on a number of examples. Besides outlining the basic theory, this thesis also presents specialized fast algorithms for MGO. These algorithms focus on the MT, which is a unitary integral mapping that can be considered a generalization of the Fourier transform. First, a discrete representation of the MT is developed that can be computed in linear time when evaluated in the near-identity limit; finite MTs can then be implemented by iterating $K\gg 1$ near-identity MTs. Second, an algorithm based on Gauss--Freud quadrature is developed for efficiently computing finite MTs along their steepest-descent curves, which may be useful in catastrophe-optics applications beyond MGO. △ Less

Submitted 6 October, 2022; originally announced October 2022.

Comments: PhD thesis

Comparing the counter-beaming and temperature anisotropy driven aperiodic electron firehose instabilities in collisionless plasma environments

Authors: Pablo S. Moya, Rodrigo A Lopez, Marian Lazar, Stefaan Poedts, Shaaban M Shaaban

Abstract: The electron firehose instabilities are among the most studied kinetic instabilities, especially in the context of space plasmas, whose dynamics is mainly controlled by collisionless wave-particle interactions. This paper undertakes a comparative analysis of the aperiodic electron firehose instabilities excited either by the anisotropic temperature or by the electron counter-beaming populations. T… ▽ More The electron firehose instabilities are among the most studied kinetic instabilities, especially in the context of space plasmas, whose dynamics is mainly controlled by collisionless wave-particle interactions. This paper undertakes a comparative analysis of the aperiodic electron firehose instabilities excited either by the anisotropic temperature or by the electron counter-beaming populations. Two symmetric counter-beams provide an effective kinetic anisotropy similar to the temperature anisotropy of a single (non-drifting) population, with temperature along the magnetic field direction larger than that in perpendicular direction. Therefore, the counter-beaming plasma is susceptible to firehose-like instabilities (FIs), parallel and oblique branches. Here we focus on the oblique beaming FI, which is also aperiodic when the free energy is provided by symmetric counter-beams. Our results show that, for relative small drifts or beaming speeds ($U$), not exceeding the thermal speed ($αあるふぁ$), the aperiodic FIs exist in the same interval of wave-numbers and the same range of oblique angles (with respect to the magnetic field direction), but the growth rates of counter-beaming FI (CBFI) are always higher than those of temperature anisotropy FI (TAFI). For $U/αあるふぁ> 1$, however, another electrostatic two-stream instability (ETSI) is also predicted, which may have growth rates higher than those of CBFI, and may dominate in that case the dynamics. △ Less

Submitted 1 September, 2022; originally announced September 2022.

Comments: 19 pages, 6 figures, 2 tables, accepted for publication in The Astrophysical Journal

Final Result of the MAJORANA DEMONSTRATOR's Search for Neutrinoless Double-$βべーた$ Decay in $^{76}$Ge

Authors: I. J. Arnquist, F. T. Avignone III, A. S. Barabash, C. J. Barton, P. J. Barton, K. H. Bhimani, E. Blalock, B. Bos, M. Busch, M. Buuck, T. S. Caldwell, Y-D. Chan, C. D. Christofferson, P. -H. Chu, M. L. Clark, C. Cuesta, J. A. Detwiler, Yu. Efremenko, H. Ejiri, S. R. Elliott, G. K. Giovanetti, M. P. Green, J. Gruszko, I. S. Guinn, V. E. Guiseppe , et al. (35 additional authors not shown)

Abstract: The MAJORANA DEMONSTRATOR searched for neutrinoless double-$βべーた$ decay ($0νββ$) of $^{76}$Ge using modular arrays of high-purity Ge detectors operated in vacuum cryostats in a low-background shield. The arrays operated with up to 40.4 kg of detectors (27.2 kg enriched to $\sim$88\% in $^{76}$Ge). From these measurements, the DEMONSTRATOR has accumulated 64.5 kg yr of enriched active exposure. With a… ▽ More The MAJORANA DEMONSTRATOR searched for neutrinoless double-$βべーた$ decay ($0νββ$) of $^{76}$Ge using modular arrays of high-purity Ge detectors operated in vacuum cryostats in a low-background shield. The arrays operated with up to 40.4 kg of detectors (27.2 kg enriched to $\sim$88\% in $^{76}$Ge). From these measurements, the DEMONSTRATOR has accumulated 64.5 kg yr of enriched active exposure. With a world-leading energy resolution of 2.52 keV FWHM at the 2039 keV $Q_{βべーたβべーた}$ (0.12\%), we set a half-life limit of $0νββ$ in $^{76}$Ge at $T_{1/2}>8.3\times10^{25}$ yr (90\% C.L.). This provides a range of upper limits on $m_{βべーたβべーた}$ of $(113-269)$ meV (90\% C.L.), depending on the choice of nuclear matrix elements. △ Less

Submitted 10 February, 2023; v1 submitted 15 July, 2022; originally announced July 2022.

Comments: 8 pages, 2 figures

Journal ref: Phys. Rev. Lett. 130, 062501 (2022)

Is the Observable Universe Consistent with the Cosmological Principle?

Authors: Pavan Kumar Aluri, Paolo Cea, Pravabati Chingangbam, Ming-Chung Chu, Roger G. Clowes, Damien Hutsemékers, Joby P. Kochappan, Alexia M. Lopez, Lang Liu, Niels C. M. Martens, C. J. A. P. Martins, Konstantinos Migkas, Eoin Ó Colgáin, Pratyush Pranav, Lior Shamir, Ashok K. Singal, M. M. Sheikh-Jabbari, Jenny Wagner, Shao-Jiang Wang, David L. Wiltshire, Shek Yeung, Lu Yin, Wen Zhao

Abstract: The Cosmological Principle (CP) -- the notion that the Universe is spatially isotropic and homogeneous on large scales -- underlies a century of progress in cosmology. It is conventionally formulated through the Friedmann-Lemaître-Robertson-Walker (FLRW) cosmologies as the spacetime metric, and culminates in the successful and highly predictive $Λらむだ$-Cold-Dark-Matter ($Λらむだ$CDM) model. Yet, tensions ha… ▽ More The Cosmological Principle (CP) -- the notion that the Universe is spatially isotropic and homogeneous on large scales -- underlies a century of progress in cosmology. It is conventionally formulated through the Friedmann-Lemaître-Robertson-Walker (FLRW) cosmologies as the spacetime metric, and culminates in the successful and highly predictive $Λらむだ$-Cold-Dark-Matter ($Λらむだ$CDM) model. Yet, tensions have emerged within the $Λらむだ$CDM model, most notably a statistically significant discrepancy in the value of the Hubble constant, $H_0$. Since the notion of cosmic expansion determined by a single parameter is intimately tied to the CP, implications of the $H_0$ tension may extend beyond $Λらむだ$CDM to the CP itself. This review surveys current observational hints for deviations from the expectations of the CP, highlighting synergies and disagreements that warrant further study. Setting aside the debate about individual large structures, potential deviations from the CP include variations of cosmological parameters on the sky, discrepancies in the cosmic dipoles, and mysterious alignments in quasar polarizations and galaxy spins. While it is possible that a host of observational systematics are impacting results, it is equally plausible that precision cosmology may have outgrown the FLRW paradigm, an extremely pragmatic but non-fundamental symmetry assumption. △ Less

Submitted 27 February, 2023; v1 submitted 12 July, 2022; originally announced July 2022.

Comments: extended contents and references, 73 pages (excluding references), 30 figures, version accepted for publication in Class. Quant. Grav. "Focus issue on the Hubble constant tension"

Journal ref: Classical and Quantum Gravity, Vol. 40, Issue No. 9, Page No. 094001 (2023)

Multivalent optical cycling centers in polyatomic molecules

Authors: Phelan Yu, Adrian Lopez, William A. Goddard III, Nicholas R. Hutzler

Abstract: Optical control of polyatomic molecules promises new opportunities in precision metrology, fundamental chemistry, quantum information, and many-body science. Contemporary experimental and theoretical efforts have mostly focused on cycling photons via excitation of a single electron localized to an alkaline earth (group 2)-like metal center. In this manuscript, we consider pathways towards optical… ▽ More Optical control of polyatomic molecules promises new opportunities in precision metrology, fundamental chemistry, quantum information, and many-body science. Contemporary experimental and theoretical efforts have mostly focused on cycling photons via excitation of a single electron localized to an alkaline earth (group 2)-like metal center. In this manuscript, we consider pathways towards optical cycling in polyatomic molecules with multi-electron degrees of freedom, which arise from two or more cycling electrons localized to $p$-block post-transition metal and metalloid (group 13, 14, and 15) centers. We characterize the electronic structure and rovibrational branching of several prototypical candidates using ab initio quantum chemical methods. Despite increased internal complexity and challenging design parameters, we find several molecules possessing quasi-closed photon cycling schemes with highly diagonal, visible and near-infrared transitions. Furthermore, we identify new heuristics for engineering optically controllable and laser-coolable polyatomic molecules with multi-electron cycling centers. Our results help elucidate the interplay between hybridization, repulsion, and ionicity in optically active species and provide a first step towards using polyatomic molecules with complex electronic structure as a resource for quantum science and measurement. △ Less

Submitted 24 May, 2022; originally announced May 2022.

Comments: 29 pages, 11 figures, 5 tables, and 6 supplemental data tables

Journal ref: Phys. Chem. Chem. Phys. 25, 154 (2023)

Mixing the solar wind proton and electron scales. Theory and 2D-PIC simulations of firehose instability

Authors: R. A. López, A. Micera, M. Lazar, S. Poedts, G. Lapenta, A. N. Zhukov, E. Boella, S. M. Shaaban

Abstract: Firehose-like instabilities (FIs) are cited in multiple astrophysical applications. Of particular interest are the kinetic manifestations in weakly-collisional or even collisionless plasmas, where these instabilities are expected to contribute to the evolution of macroscopic parameters. Relatively recent studies have initiated a realistic description of FIs, as induced by the interplay of both spe… ▽ More Firehose-like instabilities (FIs) are cited in multiple astrophysical applications. Of particular interest are the kinetic manifestations in weakly-collisional or even collisionless plasmas, where these instabilities are expected to contribute to the evolution of macroscopic parameters. Relatively recent studies have initiated a realistic description of FIs, as induced by the interplay of both species, electrons and protons, dominant in the solar wind plasma. This work complements the current knowledge with new insights from linear theory and the first disclosures from 2D PIC simulations, identifying the fastest growing modes near the instability thresholds and their long-run consequences on the anisotropic distributions. Thus, unlike previous setups, these conditions are favorable to those aperiodic branches that propagate obliquely to the uniform magnetic field, with (maximum) growth rates higher than periodic, quasi-parallel modes. Theoretical predictions are, in general, confirmed by the simulations. The aperiodic electron FI (a-EFI) remains unaffected by the proton anisotropy, and saturates rapidly at low-level fluctuations. Regarding the firehose instability at proton scales, we see a stronger competition between the periodic and aperiodic branches. For the parameters chosen in our analysis, the a-PFI is excited before than the p-PFI, with the latter reaching a significantly higher fluctuation power. However, both branches are significantly enhanced by the presence of anisotropic electrons. The interplay between EFIs and PFIs also produces a more pronounced proton isotropization. △ Less

Submitted 4 May, 2022; originally announced May 2022.

Comments: Accepted in ApJ

Experimental study of 13C(αあるふぁ,n)16O reactions in the Majorana Demonstrator calibration data

Authors: MAJORANA Collaboration, I. J. Arnquist, F. T. Avignone III, A. S. Barabash, C. J. Barton, K. H. Bhimani, E. Blalock, B. Bos, M. Busch, M. Buuck, T. S. Caldwell, Y-D. Chan, C. D. Christofferson, P. -H. Chu, M. L. Clark, C. Cuesta, J. A. Detwiler, Yu. Efremenko, H. Ejiri, S. R. Elliott, G. K. Giovanetti, M. P. Green, J. Gruszko, I. S. Guinn, V. E. Guiseppe , et al. (33 additional authors not shown)

Abstract: Neutron captures and delayed decays of reaction products are common sources of backgrounds in ultra-rare event searches. In this work, we studied $^{13}$C($αあるふぁ,n)^{16}$O reactions induced by $αあるふぁ$-particles emitted within the calibration sources of the \textsc{Majorana Demonstrator}. These sources are thorium-based calibration standards enclosed in carbon-rich materials. The reaction rate was estimate… ▽ More Neutron captures and delayed decays of reaction products are common sources of backgrounds in ultra-rare event searches. In this work, we studied $^{13}$C($αあるふぁ,n)^{16}$O reactions induced by $αあるふぁ$-particles emitted within the calibration sources of the \textsc{Majorana Demonstrator}. These sources are thorium-based calibration standards enclosed in carbon-rich materials. The reaction rate was estimated by using the 6129-keV $γがんま$-rays emitted from the excited $^{16}$O states that are populated when the incoming $αあるふぁ$-particles exceed the reaction Q-value. Thanks to the excellent energy performance of the \textsc{Demonstrator}'s germanium detectors, these characteristic photons can be clearly observed in the calibration data. Facilitated by \textsc{Geant4} simulations, a comparison between the observed 6129-keV photon rates and predictions by a TALYS-based software was performed. The measurements and predictions were found to be consistent, albeit with large statistical uncertainties. This agreement provides support for background projections from ($αあるふぁ,n$)-reactions in future double-beta decay search efforts. △ Less

Submitted 11 July, 2022; v1 submitted 27 March, 2022; originally announced March 2022.

Comments: Published version

Journal ref: Physical Review C 105, 064610 (2022)

Gaussian Processes for Finite Size Extrapolation of Many-Body Simulations

Authors: Edgar Josué Landinez Borda, Kenneth O. Berard, Annette Lopez, Brenda Rubenstein

Abstract: Key to being able to accurately model the properties of realistic materials is being able to predict their properties in the thermodynamic limit. Nevertheless, because most many-body electronic structure methods scale as a high-order polynomial, or even exponentially, with system size, directly simulating large systems in their thermodynamic limit rapidly becomes computationally intractable. As a… ▽ More Key to being able to accurately model the properties of realistic materials is being able to predict their properties in the thermodynamic limit. Nevertheless, because most many-body electronic structure methods scale as a high-order polynomial, or even exponentially, with system size, directly simulating large systems in their thermodynamic limit rapidly becomes computationally intractable. As a result, researchers typically estimate the properties of large systems that approach the thermodynamic limit by extrapolating the properties of smaller, computationally-accessible systems based on relatively simple scaling expressions. In this work, we employ Gaussian processes to more accurately and efficiently extrapolate many-body simulations to their thermodynamic limit. We train our Gaussian processes on Smooth Overlap of Atomic Positions (SOAP) descriptors to extrapolate the energies of one-dimensional hydrogen chains obtained using two high-accuracy many-body methods: Coupled Cluster theory and Auxiliary Field Quantum Monte Carlo (AFQMC). In so doing, we show that Gaussian processes trained on relatively short, 10-30-atom chains can predict the energies of both homogeneous and inhomogeneous hydrogen chains in their thermodynamic limit with sub-milliHartree accuracy. Unlike standard scaling expressions, our GPR-based approach is highly generalizable given representative training data and is not dependent on systems' geometries or dimensionality. This work highlights the potential for machine learning to correct for the finite size effects that routinely complicate the interpretation of finite size many-body simulations. △ Less

Submitted 3 April, 2024; v1 submitted 19 December, 2021; originally announced December 2021.

Comments: 14 pages, 8 figures

Metaplectic geometrical optics for ray-based modeling of caustics: Theory and algorithms

Authors: N. A. Lopez, I. Y. Dodin

Abstract: The optimization of radiofrequency-wave (RF) systems for fusion experiments is often performed using ray-tracing codes, which rely on the geometrical-optics (GO) approximation. However, GO fails at caustics such as cutoffs and focal points, erroneously predicting the wave intensity to be infinite. This is a critical shortcoming of GO, since the caustic wave intensity is often the quantity of inter… ▽ More The optimization of radiofrequency-wave (RF) systems for fusion experiments is often performed using ray-tracing codes, which rely on the geometrical-optics (GO) approximation. However, GO fails at caustics such as cutoffs and focal points, erroneously predicting the wave intensity to be infinite. This is a critical shortcoming of GO, since the caustic wave intensity is often the quantity of interest, e.g. RF heating. Full-wave modeling can be used instead, but the computational cost limits the speed at which such optimizations can be performed. We have developed a less expensive alternative called metaplectic geometrical optics (MGO). Instead of evolving waves in the usual $\textbf{x}$ (coordinate) or $\text{k}$ (spectral) representation, MGO uses a mixed $\textbf{X} \equiv \textsf{A}\textbf{x} + \textsf{B}\textbf{k}$ representation. By continuously adjusting the matrix coefficients $\textsf{A}$ and $\textsf{B}$ along the rays, one can ensure that GO remains valid in the $\textbf{X}$ coordinates without caustic singularities. The caustic-free result is then mapped back onto the original $\textbf{x}$ space using metaplectic transforms. Here, we overview the MGO theory and review algorithms that will aid the development of an MGO-based ray-tracing code. We show how using orthosymplectic transformations leads to considerable simplifications compared to previously published MGO formulas. We also prove explicitly that MGO exactly reproduces standard GO when evaluated far from caustics (an important property which until now has only been inferred from numerical simulations), and we relate MGO to other semiclassical caustic-removal schemes published in the literature. This discussion is then augmented by an explicit comparison of the computed spectrum for a wave bounded between two cutoffs. △ Less

Submitted 10 May, 2022; v1 submitted 14 December, 2021; originally announced December 2021.

Comments: Invited paper for APS DPP 2021. 21 pages, 6 figures, 5 appendices

Journal ref: Phys.Plasmas 29, 052111 (2022)

The MAJORANA DEMONSTRATOR Readout Electronics System

Authors: N. Abgrall, M. Amman, I. J. Arnquist, F. T. Avignone III, A. S. Barabash, C. J. Barton, P. J. Barton, F. E. Bertrand, K. H. Bhimani, B. Bos, A. W. Bradley, T. H. Burritt, M. Busch, M. Buuck, T. S. Caldwell, Y-D. Chan, C. D. Christofferson, P. -H. Chu, M. L. Clark, R. J. Cooper, C. Cuesta, J. A. Detwiler, A. Drobizhev, D. W. Edwins, Yu. Efremenko , et al. (54 additional authors not shown)

Abstract: The MAJORANA DEMONSTRATOR comprises two arrays of high-purity germanium detectors constructed to search for neutrinoless double-beta decay in 76-Ge and other physics beyond the Standard Model. Its readout electronics were designed to have low electronic noise, and radioactive backgrounds were minimized by using low-mass components and low-radioactivity materials near the detectors. This paper prov… ▽ More The MAJORANA DEMONSTRATOR comprises two arrays of high-purity germanium detectors constructed to search for neutrinoless double-beta decay in 76-Ge and other physics beyond the Standard Model. Its readout electronics were designed to have low electronic noise, and radioactive backgrounds were minimized by using low-mass components and low-radioactivity materials near the detectors. This paper provides a description of all components of the MAJORANA DEMONSTRATOR readout electronics, spanning the front-end electronics and internal cabling, back-end electronics, digitizer, and power supplies, along with the grounding scheme. The spectroscopic performance achieved with these readout electronics is also demonstrated. △ Less

Submitted 23 February, 2022; v1 submitted 17 November, 2021; originally announced November 2021.

Comments: For submission to JINST, 17 figures. v2: revised version

Signatures of muonic activation in the Majorana Demonstrator

Authors: I. J. Arnquist, F. T. Avignone III, A. S. Barabash, C. J. Barton, F. E. Bertrand, E. Blalock, B. Bos, M. Busch, M. Buuck, T. S. Caldwell, Y-D. Chan, C. D. Christofferson, P. -H. Chu, M. L. Clark, C. Cuesta, J. A. Detwiler, T. R. Edwards, Yu. Efremenko, H. Ejiri, S. R. Elliott, G. K. Giovanetti, M. P. Green, J. Gruszko, I. S. Guinn, V. E. Guiseppe , et al. (33 additional authors not shown)

Abstract: Experiments searching for very rare processes such as neutrinoless double-beta decay require a detailed understanding of all sources of background. Signals from radioactive impurities present in construction and detector materials can be suppressed using a number of well-understood techniques. Background from in-situ cosmogenic interactions can be reduced by siting an experiment deep underground.… ▽ More Experiments searching for very rare processes such as neutrinoless double-beta decay require a detailed understanding of all sources of background. Signals from radioactive impurities present in construction and detector materials can be suppressed using a number of well-understood techniques. Background from in-situ cosmogenic interactions can be reduced by siting an experiment deep underground. However, the next generation of such experiments have unprecedented sensitivity goals of 10$^{28}$ years half-life with background rates of 10$^{-5}$cts/(keV kg yr) in the region of interest. To achieve these goals, the remaining cosmogenic background must be well understood. In the work presented here, Majorana Demonstrator data is used to search for decay signatures of meta-stable germanium isotopes. Contributions to the region of interest in energy and time are estimated using simulations, and compared to Demonstrator data. Correlated time-delayed signals are used to identify decay signatures of isotopes produced in the germanium detectors. A good agreement between expected and measured rate is found and different simulation frameworks are used to estimate the uncertainties of the predictions. The simulation campaign is then extended to characterize the background for the LEGEND experiment, a proposed tonne-scale effort searching for neutrinoless double-beta decay in $^{76}$Ge. △ Less

Submitted 27 October, 2021; originally announced October 2021.

Comments: 14 pages, 12 figures

On the role of solar wind expansion as a source of whistler waves: scattering of suprathermal electrons and heat flux regulation in the inner heliosphere

Authors: A. Micera, A. N. Zhukov, R. A. López, E. Boella, A. Tenerani, M. Velli, G. Lapenta, M. E. Innocenti

Abstract: The role of solar wind expansion in generating whistler waves is investigated using the EB-iPic3D code, which models solar wind expansion self-consistently within a fully kinetic semi-implicit approach. The simulation is initialized with an electron velocity distribution function modeled after Parker Solar Probe observations during its first perihelion at 0.166 au, consisting of a dense core and a… ▽ More The role of solar wind expansion in generating whistler waves is investigated using the EB-iPic3D code, which models solar wind expansion self-consistently within a fully kinetic semi-implicit approach. The simulation is initialized with an electron velocity distribution function modeled after Parker Solar Probe observations during its first perihelion at 0.166 au, consisting of a dense core and an anti-sunward strahl. This distribution function is initially stable with respect to kinetic instabilities. Expansion drives the solar wind into successive regimes where whistler heat flux instabilities are triggered. These instabilities produce sunward whistler waves initially characterized by predominantly oblique propagation with respect to the interplanetary magnetic field. The excited waves interact with the electrons via resonant scattering processes. As a consequence, the strahl pitch angle distribution broadens and its drift velocity reduces. Strahl electrons are scattered in the direction perpendicular to the magnetic field, and an electron halo is formed. At a later stage, resonant electron firehose instability is triggered and further affects the electron temperature anisotropy as the solar wind expands. Wave-particle interaction processes are accompanied by a substantial reduction of the solar wind heat flux. The simulated whistler waves are in qualitative agreement with observations in terms of wave frequencies, amplitudes and propagation angles. Our work proposes an explanation for the observations of oblique and parallel whistler waves in the solar wind. We conclude that solar wind expansion has to be factored in when trying to explain kinetic processes at different heliocentric distances. △ Less

Submitted 1 July, 2021; v1 submitted 30 June, 2021; originally announced June 2021.

Steepest-descent algorithm for simulating plasma-wave caustics via metaplectic geometrical optics

Authors: Sean M. Donnelly, Nicolas A. Lopez, I. Y. Dodin

Abstract: The design and optimization of radiofrequency-wave systems for fusion applications is often performed using ray-tracing codes, which rely on the geometrical-optics (GO) approximation. However, GO fails at wave cutoffs and caustics. To accurately model the wave behavior in these regions, more advanced and computationally expensive "full-wave" simulations are typically used, but this is not strictly… ▽ More The design and optimization of radiofrequency-wave systems for fusion applications is often performed using ray-tracing codes, which rely on the geometrical-optics (GO) approximation. However, GO fails at wave cutoffs and caustics. To accurately model the wave behavior in these regions, more advanced and computationally expensive "full-wave" simulations are typically used, but this is not strictly necessary. A new generalized formulation called metaplectic geometrical optics (MGO) has been proposed that reinstates GO near caustics. The MGO framework yields an integral representation of the wavefield that must be evaluated numerically in general. We present an algorithm for computing these integrals using Gauss-Freud quadrature along the steepest-descent contours. Benchmarking is performed on the standard Airy problem, for which the exact solution is known analytically. The numerical MGO solution provided by the new algorithm agrees remarkably well with the exact solution and significantly improves upon previously derived analytical approximations of the MGO integral. △ Less

Submitted 13 August, 2021; v1 submitted 27 April, 2021; originally announced April 2021.

Comments: 13 pages, 5 figures, 2 tables, 1 appendix

Journal ref: Phys. Rev. E 104, 025304 (2021)

Benchmarking the Scrape-Off-Layer Fast Ion (SOLFI) particle tracer code in a collisionless magnetic mirror with electrostatic potential drop

Authors: Xin Zhang, Nicolas A Lopez, Francesca M Poli

Abstract: Optimizing the confinement and transport of fast ions is an important consideration in the design of modern fusion reactors. For spherical tokamaks in particular, fast ions can significantly influence global plasma behavior because their large drift orbits often sample both core and scrape-off-layer (SOL) plasma conditions. Their Larmor radii are also comparable to the SOL width, rendering the com… ▽ More Optimizing the confinement and transport of fast ions is an important consideration in the design of modern fusion reactors. For spherical tokamaks in particular, fast ions can significantly influence global plasma behavior because their large drift orbits often sample both core and scrape-off-layer (SOL) plasma conditions. Their Larmor radii are also comparable to the SOL width, rendering the commonly chosen guiding center approximations inappropriate. Accurately modeling the behavior of fast ions therefore requires retaining a complete description of the fast ion orbit including its Larmor motion. Here, we introduce the Scrape-Off-Layer Fast Ion (SOLFI) code, which is a new and versatile full-orbit Monte Carlo particle tracer being developed to follow fast ion orbits inside and outside the separatrix. We benchmark SOLFI in a simple straight mirror geometry and show that the code (i) conserves particle energy and magnetic moment, (ii) obtains the correct passing boundary for particles moving in magnetic mirror field with an imposed electrostatic field, and (iii) correctly observes equal ion and electron current at the ambipolar potential predicted from analytical theory. △ Less

Submitted 29 June, 2021; v1 submitted 15 April, 2021; originally announced April 2021.

Comments: 11 pages, 9 figures

Characterization of Nb22O54 microrods grown from niobium oxide powders recovered from mine tailings

Authors: B. Sotillo, F. A. Lopez, L. Alcaraz, P. Fernandez

Abstract: In this work, the possibility of using niobium oxide recovered from tailings from the Penouta Sn-Nb-Ta deposit (located at the North of Spain) as starting material for growing microstructures is demonstrated. The properties of the starting material have been studied to understand its crystal structure, quality and purity. Recovered niobium oxide powders are mainly of TT-Nb2O5. These powders have b… ▽ More In this work, the possibility of using niobium oxide recovered from tailings from the Penouta Sn-Nb-Ta deposit (located at the North of Spain) as starting material for growing microstructures is demonstrated. The properties of the starting material have been studied to understand its crystal structure, quality and purity. Recovered niobium oxide powders are mainly of TT-Nb2O5. These powders have been used to grow Nb22O54 microrods by an evaporation method in an argon atmosphere. Different characterization techniques (X-Ray diffraction, scanning electron microscopy, micro-Raman spectroscopy, luminescence) have been used to determine the properties of Nb22O54 microrods, mainly focusing on the crystal quality and refractive index. The present study opens the way to the transformation of waste (mine tailings) into a material of high technological value as niobium oxide, and its reintroduction into the value chain for a wide range of applications, from coatings to batteries and supercapacitors. △ Less

Submitted 4 April, 2021; originally announced April 2021.

Journal ref: Ceramics International 47, 13859-13864 (2021)

General dispersion properties of magnetized plasmas with drifting bi-Kappa distributions. DIS-K: DIspersion Solver for Kappa plasmas

Authors: R. A. López, S. M. Shaaban, M. Lazar

Abstract: Space plasmas are known to be out of (local) thermodynamic equilibrium, as observations show direct or indirect evidences of non-thermal velocity distributions of plasma particles. Prominent are the anisotropies relative to the magnetic field, anisotropic temperatures, field-aligned beams or drifting populations, but also, the suprathermal populations enhancing the high-energy tails of the observe… ▽ More Space plasmas are known to be out of (local) thermodynamic equilibrium, as observations show direct or indirect evidences of non-thermal velocity distributions of plasma particles. Prominent are the anisotropies relative to the magnetic field, anisotropic temperatures, field-aligned beams or drifting populations, but also, the suprathermal populations enhancing the high-energy tails of the observed distributions. Drifting bi-Kappa distribution functions can provide a good representation of these features and enable for a kinetic fundamental description of the dispersion and stability of these collision-poor plasmas, where particle-particle collisions are rare but wave-particle interactions appears to play a dominant role in the dynamic. In the present paper we derive the full set of components of the dispersion tensor for magnetized plasma populations modeled by drifting bi-Kappa distributions. A new solver called DIS-K (DIspersion Solver for Kappa plasmas) is proposed to solve numerically the dispersion relations of high complexity. The solver is validated by comparing to the damped and unstable wave solutions obtained with other codes, operating in the limits of drifting Maxwellian and non-drifting Kappa models. These new theoretical tools enable more realistic characterizations, both analytical and numerical, of wave fluctuations and instabilities in complex kinetic configurations measured in-situ in space plasmas. △ Less

Submitted 6 May, 2021; v1 submitted 24 February, 2021; originally announced February 2021.

Comments: Accepted for publication in Journal of Plasma Physics

Journal ref: Journal of Plasma Physics, Volume 87, Issue 3, June 2021, 905870310

On the interplay of solar wind proton and electron instabilities: Linear and quasi-linear approaches

Authors: S. M. Shaaban, M. Lazar, R. A. López, R. F. Wimmer-Schweingruber

Abstract: Important efforts are currently made for understanding the so-called kinetic instabilities, driven by the anisotropy of different species of plasma particles present in the solar wind and terrestrial magnetosphere. These instabilities are fast enough to efficiently convert the free energy of plasma particles into enhanced (small-scale) fluctuations with multiple implications, regulating the anisot… ▽ More Important efforts are currently made for understanding the so-called kinetic instabilities, driven by the anisotropy of different species of plasma particles present in the solar wind and terrestrial magnetosphere. These instabilities are fast enough to efficiently convert the free energy of plasma particles into enhanced (small-scale) fluctuations with multiple implications, regulating the anisotropy of plasma particles. In this paper we use both linear and quasilinear (QL) frameworks to describe complex unstable regimes, which realistically combine different temperature anisotropies of electrons and ions (protons). Thus parameterized are various instabilities, e.g., proton and electron firehose, electromagnetic ion cyclotron, and whistler instability, showing that their main linear properties are markedly altered by the interplay of anisotropic electrons and protons. Linear theory may predict a strong competition of two instabilities of different nature when their growth rates are comparable. In the QL phase wave fluctuations grow and saturate at different levels and temporal scales, by comparison to the individual excitation of the proton or electron instabilities. In addition, cumulative effects of the combined proton and electron induced fluctuations can markedly stimulate the relaxations of their temperature anisotropies. Only whistler fluctuations inhibit the efficiency of proton firehose fluctuations in the relaxation of anisotropic protons. These results offer valuable premises for further investigations in numerical simulations, to decode the full spectrum of kinetic instabilities resulting from the interplay of anisotropic electrons and protons in space plasmas. △ Less

Submitted 11 January, 2021; originally announced January 2021.

Comments: Accepted for publication at MNRAS

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 503, Issue 3, May 2021, Page 3134

Average conservative chaos in quantum dusty plasmas

Authors: Álvaro G. López, Rustam Ali, Laxmikanta Mandi, Prasanta Chatterjee

Abstract: We consider a hydrodynamic model of a quantum dusty plasma. We prove mathematically that the resulting dust ion acoustic plasma waves present the property of being conservative on average. Furthermore, we test this property numerically, confirming its validity. Using standard techniques from the study of dynamical systems, as for example the Lyapunov characteristic exponents, we investigate the ch… ▽ More We consider a hydrodynamic model of a quantum dusty plasma. We prove mathematically that the resulting dust ion acoustic plasma waves present the property of being conservative on average. Furthermore, we test this property numerically, confirming its validity. Using standard techniques from the study of dynamical systems, as for example the Lyapunov characteristic exponents, we investigate the chaotic dynamics of the plasma and show numerically its existence for a wide range of parameter values. Finally, we illustrate how chaotic dynamics organizes in the parameter space for fixed values of the initial conditions, as the Mach number and the quantum diffraction parameter are continuously varied. △ Less

Submitted 21 December, 2020; originally announced January 2021.

Exactly unitary discrete representations of the metaplectic transform for linear-time algorithms

Authors: N. A. Lopez, I. Y. Dodin

Abstract: The metaplectic transform (MT), a generalization of the Fourier transform sometimes called the linear canonical transform, is a tool used ubiquitously in modern optics, for example, when calculating the transformations of light beams in paraxial optical systems. The MT is also an essential ingredient of the geometrical-optics modeling of caustics that was recently proposed by the authors. In parti… ▽ More The metaplectic transform (MT), a generalization of the Fourier transform sometimes called the linear canonical transform, is a tool used ubiquitously in modern optics, for example, when calculating the transformations of light beams in paraxial optical systems. The MT is also an essential ingredient of the geometrical-optics modeling of caustics that was recently proposed by the authors. In particular, this application relies on the near-identity MT (NIMT); however, the NIMT approximation used so far is not exactly unitary and leads to numerical instability. Here, we develop a discrete MT that is exactly unitary, and approximate it to obtain a discrete NIMT that is also unitary and can be computed in linear time. We prove that the discrete NIMT converges to the discrete MT when iterated, thereby allowing the NIMT to compute MTs that are not necessarily near-identity. We then demonstrate the new algorithms with a series of examples. △ Less

Submitted 8 March, 2021; v1 submitted 10 December, 2020; originally announced December 2020.

Comments: 10 pages, 3 figures, 1 appendix

Metaplectic geometrical optics for modeling caustics in uniform and nonuniform media

Authors: N. A. Lopez, I. Y. Dodin

Abstract: As an approximate theory that is highly regarded for its computational efficiency, geometrical optics (GO) is widely used for modeling waves in various areas of physics. However, GO fails at caustics, which significantly limits its applicability. A new framework, called metaplectic geometrical optics (MGO), has recently been developed that allows caustics of certain types to be modeled accurately… ▽ More As an approximate theory that is highly regarded for its computational efficiency, geometrical optics (GO) is widely used for modeling waves in various areas of physics. However, GO fails at caustics, which significantly limits its applicability. A new framework, called metaplectic geometrical optics (MGO), has recently been developed that allows caustics of certain types to be modeled accurately within the GO framework. Here, we extend MGO to the most general case. To illustrate our new theory, we also apply it to several sample problems, including calculations of two-dimensional wavefields near fold and cusp caustics. In contrast with traditional-GO solutions, the corresponding MGO solutions are finite everywhere and approximate well the true wavefield near these caustics. △ Less

Submitted 8 March, 2021; v1 submitted 7 September, 2020; originally announced September 2020.

Comments: 14 pages, 3 figures, 3 appendices. Complete derivations in supplementary pdf

Journal ref: J. Opt. 23, 025601 (2021)

$αあるふぁ$-event Characterization and Rejection in Point-Contact HPGe Detectors

Authors: The MAJORANA Collaboration, I. J. Arnquist, F. T. Avignone III, A. S. Barabash, C. J. Barton, F. E. Bertrand, E. Blalock, B. Bos, M. Busch, M. Buuck, T. S. Caldwell, Y-D. Chan, C. D. Christofferson, P. -H. Chu, M. L. Clark, C. Cuesta, J. A. Detwiler, A. Drobizhev, T. R. Edwards, D. W. Edwins, Yu. Efremenko, S. R. Elliott, T. Gilliss, G. K. Giovanetti, M. P. Green , et al. (40 additional authors not shown)

Abstract: P-type point contact (PPC) HPGe detectors are a leading technology for rare event searches due to their excellent energy resolution, low thresholds, and multi-site event rejection capabilities. We have characterized a PPC detector's response to $αあるふぁ$ particles incident on the sensitive passivated and p+ surfaces, a previously poorly-understood source of background. The detector studied is identical… ▽ More P-type point contact (PPC) HPGe detectors are a leading technology for rare event searches due to their excellent energy resolution, low thresholds, and multi-site event rejection capabilities. We have characterized a PPC detector's response to $αあるふぁ$ particles incident on the sensitive passivated and p+ surfaces, a previously poorly-understood source of background. The detector studied is identical to those in the MAJORANA DEMONSTRATOR experiment, a search for neutrinoless double-beta decay ($0νββ$) in $^{76}$Ge. $αあるふぁ$ decays on most of the passivated surface exhibit significant energy loss due to charge trapping, with waveforms exhibiting a delayed charge recovery (DCR) signature caused by the slow collection of a fraction of the trapped charge. The DCR is found to be complementary to existing methods of $αあるふぁ$ identification, reliably identifying $αあるふぁ$ background events on the passivated surface of the detector. We demonstrate effective rejection of all surface $αあるふぁ$ events (to within statistical uncertainty) with a loss of only 0.2% of bulk events by combining the DCR discriminator with previously-used methods. The DCR discriminator has been used to reduce the background rate in the $0νββ$ region of interest window by an order of magnitude in the MAJORANA DEMONSTRATOR, and will be used in the upcoming LEGEND-200 experiment. △ Less

Submitted 14 March, 2022; v1 submitted 23 June, 2020; originally announced June 2020.

Comments: Accepted for publication in European Journal of Physics C

Generalized Sturmian Functions in prolate spheroidal coordinates

Authors: D. M. Mitnik, F. A. Lopez, L. U. Ancarani

Abstract: With the aim of describing bound and continuum states for diatomic molecules, we develop and implement a spectral method that makes use of Generalized Sturmian Functions (GSF) in prolate spheroidal coordinates. In order to master all computational issues, we apply here the method to one--electron molecular ions and compare it with benchmark data for both ground and excited states. We actually prop… ▽ More With the aim of describing bound and continuum states for diatomic molecules, we develop and implement a spectral method that makes use of Generalized Sturmian Functions (GSF) in prolate spheroidal coordinates. In order to master all computational issues, we apply here the method to one--electron molecular ions and compare it with benchmark data for both ground and excited states. We actually propose two different computational schemes to solve the two coupled differential equations. The first one is an iterative 1d procedure in which one solves alternately the angular and the radial equations, the latter yielding the state energy. The second, named direct $2d$ method, consists in representing the Hamiltonian matrix in a two--dimensional GSF basis set, and its further diagonalization. Both spectral schemes are timewise computationally efficient since the basis elements are such that no derivatives have to be calculated numerically. Moreover, very accurate results are obtained with minimal basis sets. This is related on one side to the use of the natural coordinate system and, on the other, to the intrinsic good property of all GSF basis elements that are constructed as to obey appropriate physical boundary conditions. The present implementation for bound states paves the way for the study of continuum states involved in ionization of one or two--electron diatomic targets. △ Less

Submitted 3 December, 2020; v1 submitted 11 June, 2020; originally announced June 2020.

Electromagnetic Ion-Ion Instabilities in Space Plasmas: Effects of Suprathermal Populations

Authors: S. M. Shaaban, M. Lazar, R. A. López, S. Poedts

Abstract: In collision-poor plasmas from space, three distinct ion-ion instabilities can be driven by the proton beams streaming along the background magnetic field: left-hand resonant, non-resonant, and right-hand resonant instabilities. These instabilities are in general investigated considering only idealized proton beams with Maxwellian velocity distributions, and ignoring the implications of supratherm… ▽ More In collision-poor plasmas from space, three distinct ion-ion instabilities can be driven by the proton beams streaming along the background magnetic field: left-hand resonant, non-resonant, and right-hand resonant instabilities. These instabilities are in general investigated considering only idealized proton beams with Maxwellian velocity distributions, and ignoring the implications of suprathermal populations, usually reproduced by the Kappa power-laws. Moreover, the existing theories minimize the kinetic effects of electrons, assuming them isotropic and Maxwellian distributed. In an attempt to overcome these limitations, in the present paper we present the results of an extended investigation of ion-ion instabilities, which show that their dispersion and stability properties (e.g. growth rates, wave frequencies, and the unstable wave numbers) are highly sensitive to the influence of suprathermal populations and anisotropic electrons. These results offer valuable explanations for the origin of the enhanced low-frequency fluctuations, frequently observed in space plasmas and associated with proton beams. △ Less

Submitted 13 June, 2020; v1 submitted 10 June, 2020; originally announced June 2020.

Comments: Accepted for publication at ApJ

Alternative high plasma beta regimes of electron heat-flux instabilities in the solar wind

Authors: R. A. López, M. Lazar, S. M. Shaaban, S. Poedts, P. S. Moya

Abstract: The heat transport in the solar wind is dominated by the suprathermal electron populations, i.e., a tenuous halo and a field-aligned beam/strahl, with high energies and antisunward drifts along the magnetic field. Their evolution may offer plausible explanations for the rapid decrease of the heat flux with the solar wind expansion, typically invoked being the self-generated instabilities, or the s… ▽ More The heat transport in the solar wind is dominated by the suprathermal electron populations, i.e., a tenuous halo and a field-aligned beam/strahl, with high energies and antisunward drifts along the magnetic field. Their evolution may offer plausible explanations for the rapid decrease of the heat flux with the solar wind expansion, typically invoked being the self-generated instabilities, or the so-called heat flux instabilities (HFIs). The present paper provides a unified description of the full spectrum of HFIs, as prescribed by the linear kinetic theory for high beta conditions ($βべーた_e \gg 0.1$) and different relative drifts ($U$) of the suprathermals. HFIs of different nature are distinguished, i.e., electromagnetic, electrostatic or hybrid, propagating parallel or obliquely to the magnetic field, etc., as well as their regimes of interplay (co-existence) or dominance. These alternative regimes of HFIs complement each other and may be characteristic to different relative drifts of suprathermal electrons and various conditions in the solar wind, e.g., in the slow or fast winds, streaming interaction regions and interplanetary shocks. Moreover, these results strongly suggest that heat flux regulation may not involve only one but several HFIs, concomitantly or successively in time. Conditions for a single, well defined instability with major effects on the suprathermal electrons and, implicitly, the heat flux, seem to be very limited. Whistler HFIs are more likely to occur but only for minor drifts (as also reported by recent observations), which may explain a modest implication in their regulation, shown already in quasilinear studies and numerical simulations. △ Less

Submitted 17 August, 2020; v1 submitted 7 June, 2020; originally announced June 2020.

Comments: Submitted to ApJ