Search

Measurements of the galaxy density and weak-lensing profiles of galaxy clusters typically rely on an assumed cluster center, which is taken to be the brightest cluster galaxy or other proxies for the true halo center. Departure of the assumed cluster center from the true halo center bias the resultant profile measurements, an effect known as miscentering bias. Currently, miscentering is typically modeled in stacked profiles of clusters with a two parameter model. We use an alternate approach in which the profiles of individual clusters are used with the corresponding likelihood computed using a Gaussian mixture model. We test the approach using halos and the corresponding subhalo profiles from the IllustrisTNG hydrodynamic simulations. We obtain significantly improved estimates of the miscentering parameters for both 3D and projected 2D profiles relevant for imaging surveys. We discuss applications to upcoming cosmological surveys., Comment: 7 pages, 6 figures, to be submitted to A&A The Python package used in this paper is publicly available at https://github.com/KyleMiller1/Halo-Miscentering-Mixture-Model

The main challenge of 21 cm cosmology experiments is astrophysical foregrounds which are difficult to separate from the signal due to telescope systematics. An earlier study has shown that foreground residuals induced by antenna gain errors can be estimated and subtracted using the hybrid foreground residual subtraction (HyFoReS) technique which relies on cross-correlating linearly filtered data. In this paper, we apply a similar technique to the CHIME stacking analysis to subtract beam-induced foreground contamination. Using a linear high-pass delay filter for foreground suppression, the CHIME collaboration reported a $11.1\sigma$ detection in the 21 cm signal stacked on eBOSS quasar locations, despite foreground residual contamination mostly due to the instrument chromatic transfer function. We cross-correlate the foreground-dominated data at low delay with the contaminated signal at high delay to estimate residual foregrounds and subtract them from the signal. We find foreground residual subtraction can improve the signal-to-noise ratio of the stacked 21 cm signal by $ 10 - 20\%$ after the delay foreground filter, although some of the improvement can also be achieved with an alternative flagging technique. We have shown that it is possible to use HyFoReS to reduce beam-induced foreground contamination, benefiting the analysis of the HI auto power spectrum with CHIME and enabling the recovery of large scale modes.

A laser pulse composed of a fundamental and an appropriately phased second harmonic can drive a time-dependent current of photoionized electrons that generates broadband THz radiation. Over the propagation distances relevant to many experiments, dispersion causes the relative phase between the harmonics to evolve. This "dephasing" slows the accumulation of THz energy and results in a multi-cycle THz pulse with significant angular dispersion. Here, we introduce a novel optical configuration that compensates the relative phase evolution, allowing for the formation of a half-cycle THz pulse with almost no angular dispersion. The configuration uses the spherical aberration of an axilens to map a prescribed radial phase variation in the near field to a desired longitudinal phase variation in the far field. Simulations that combine this configuration with an ultrashort flying focus demonstrate the formation of a half-cycle THz pulse with a controlled emission angle and 1/4 the angular divergence of the multi-cycle pulse created by a conventional optical configuration.

Collision between relativistic electron sheets and counter-propagating laser pulses is recognized as a promising way to produce intense attosecond X-rays through coherent Thomson backscattering (TBS). In a double-layer scheme, the electrons in an ultrathin solid foil are first pushed out by an intense laser driver and then interact with the laser reflected off a second foil to form a high-density relativistic electron sheet with vanishing transverse momentum. However, the repulsion between these concentrated electrons can increase the thickness of the layer, reducing both its density and subsequently the coherent TBS. Here, we present a systematic study on the evolution of the flying electron layer and find that its resulting thickness is determined by the interplay between the intrinsic space-charge expansion and the velocity compression induced by the drive laser. How the laser driver, the target areal density, the reflector and the collision laser intensity affect the properties of the produced X-rays is explored. Multi-dimensional particle-in-cell simulations indicate that employing this scheme in the nonlinear regime has the potential to stably produce soft X-rays with several GW peak power in hundreds of TW ultrafast laser facilities. The pulse duration can be tuned to tens of attoseconds. This compact and intense attosecond X-ray source may have broad applications in attosecond science., Comment: 10 pages, 6 figures

Laser wakefield accelerators (LWFAs) have electric fields that are orders of magnitude larger than those of conventional accelerators, promising an attractive, small-scale alternative for next-generation light sources and lepton colliders. The maximum energy gain in a single-stage LWFA is limited by dephasing, which occurs when the trapped particles outrun the accelerating phase of the wakefield. Here, we demonstrate that a single space-time structured laser pulse can be used for ionization injection and electron acceleration over many dephasing lengths in the bubble regime. Simulations of a dephasingless laser wakefield accelerator driven by a 6.2-J laser pulse show 25 pC of injected charge accelerated over 20 dephasing lengths (1.3 cm) to a maximum energy of 2.1 GeV. The space-time structured laser pulse features an ultrashort, programmable-trajectory focus. Accelerating the focus, reducing the focused spot-size variation, and mitigating unwanted self-focusing stabilize the electron acceleration, which improves beam quality and leads to projected energy gains of 125 GeV in a single, sub-meter stage driven by a 500-J pulse., Comment: 18 pages, 4 figures

We introduce machine learning (ML) models that predict the electronic structure of materials across a wide temperature range. Our models employ neural networks and are trained on density functional theory (DFT) data. Unlike most other ML models that use DFT data, our models directly predict the local density of states (LDOS) of the electronic structure. This provides several advantages, including access to multiple observables such as the electronic density and electronic total free energy. Moreover, our models account for both the electronic and ionic temperatures independently, making them ideal for applications like laser-heating of matter. We validate the efficacy of our LDOS-based models on a metallic test system. They accurately capture energetic effects induced by variations in ionic and electronic temperatures over a broad temperature range, even when trained on a subset of these temperatures. These findings open up exciting opportunities for investigating the electronic structure of materials under both ambient and extreme conditions.

A laser pulse composed of a fundamental and properly phased second harmonic exhibits an asymmetric electric field that can drive a time-dependent current of photoionized electrons. The current produces an ultrashort burst of terahertz (THz) radiation. When driven by a conventional laser pulse, the THz radiation is emitted into a cone with an angle determined by the dispersion of the medium. Here we demonstrate that the programmable-velocity intensity peak of a spatiotemporally structured, two-color laser pulse can be used to control the emission angle, focal spot, and spectrum of the THz radiation. Of particular interest for applications, a structured pulse with a subluminal intensity peak can drive highly focusable, on-axis THz radiation.

The quality of electron beams produced from plasma-based accelerators, i.e., normalized brightness and energy spread, has made transformative progress in the past several decades in both simulation and experiment. Recently, full-scale particle-in-cell (PIC) simulations have shown that electron beams with unprecedented brightness ($10^{20}\sim10^{21}~\mathrm{A}/\mathrm{m}^2/\mathrm{rad}^2$) and $0.1\sim 1$ MeV energy spread can be produced through controlled injection in a slowly expanding bubble that arises when a particle beam or laser pulse propagates in density gradient, or when a particle beam self-focuses in uniform plasma or has a superluminal flying focus. However, in previous simulations of work on self-injection triggered by an evolving laser driver in a uniform plasma, the resulting beams did not exhibit comparable brightnesses and energy spreads. Here, we demonstrate through the use of large-scale high-fidelity PIC simulations that a slowly expanding bubble driven by a laser pulse in a uniform plasma can indeed produce self-injected electron beams with similar brightnesses and energy spreads as for an evolving bubble driven by an electron beam driver. We consider laser spot sizes roughly equal to the matched spot sizes in a uniform plasma and find that the evolution of the bubble occurs naturally through the evolution of the laser. The effects of the electron beam quality on the choice of physical as well as numerical parameters, e.g. grid sizes and field solvers used in the PIC simulations are presented. It is found that this original and simplest injection scheme can produce electron beams with beam quality exceeding that of the more recent concepts.

In a laser wakefield accelerator (LWFA), an intense laser pulse excites a plasma wave that traps and accelerates electrons to relativistic energies. When the pulse overlaps the accelerated electrons, it can enhance the energy gain through direct laser acceleration (DLA) by resonantly driving the betatron oscillations of the electrons in the plasma wave. The particle-in-cell (PIC) algorithm, although often the tool of choice to study DLA, contains inherent errors due to numerical dispersion and the time staggering of the electric and magnetic fields. Further, conventional PIC implementations cannot reliably disentangle the fields of the plasma wave and laser pulse, which obscures interpretation of the dominant acceleration mechanism. Here, a customized field solver that reduces errors from both numerical dispersion and time staggering is used in conjunction with a field decomposition into azimuthal modes to perform PIC simulations of DLA in an LWFA. Comparisons with traditional PIC methods, model equations, and experimental data show improved accuracy with the customized solver and convergence with an order-of-magnitude fewer cells. The azimuthal-mode decomposition reveals that the most energetic electrons receive comparable energy from DLA and LWFA., Comment: 10 pages, 5 figures, to submit to Physics of Plasmas

The arrow polynomial is an invariant of framed oriented virtual links that generalizes the virtual Kauffman bracket. In this paper we define the homological arrow polynomial, which generalizes the arrow polynomial to framed oriented virtual links with labeled components. The key observation is that, given a link in a thickened surface, the homology class of the link defines a functional on the surface's skein module, and by applying it to the image of the link in the skein module this gives a virtual link invariant. We give a graphical calculus for the homological arrow polynomial by taking the usual diagrams for the Kauffman bracket and including labeled "whiskers" that record intersection numbers with each labeled component of the link. We use the homological arrow polynomial to study $(\mathbb{Z}/n\mathbb{Z})$-nullhomologous virtual links and checkerboard colorability, giving a new way to complete Imabeppu's characterization of checkerboard colorability of virtual links with up to four crossings. We also prove a version of the Kauffman-Murasugi-Thistlethwaite theorem that the breadth of an evaluation of the homological arrow polynomial for an "h-reduced" diagram $D$ is $4(c(D)-g(D)+1)$., Comment: 32 pages, 16 figures, 5 tables

Bremsstrahlung x-rays generated in laser-solid interactions can be used as light sources for high-energy-density science. We present electron and x-ray spectra from particle-in-cell and Monte Carlo simulations, varying laser pulse intensity and duration at fixed energy of 200$\,$J. Superponderomotive electron temperatures are observed at low intensity; a new temperature scaling is given that depends on pulse duration and density scale length. Short, high-intensity pulses create low-divergence electron beams before self-generated magnetic fields evolve, resulting in more forward-going MeV x-rays., Comment: 6 pages, 5 figures

This paper analyses images from 43 to 340 GHz to trace the structure of the Source I disk in Orion-KL with $\sim$12 AU resolution. The data reveal an almost edge-on disk with an outside diameter $\sim$ 100 AU which is heated from the inside. The high opacity at 220-340 GHz hides the internal structure and presents a surface temperature $\sim$500 K. Images at 43, 86 and 99 GHz reveal structure within the disk. At 43 GHz there is bright compact emission with brightness temperature $\sim$1300 K. Another feature, most prominent at 99 GHz, is a warped ridge of emission. The data can be explained by a simple model with a hot inner structure, seen through cooler material. A wide angle outflow mapped in SiO emission ablates material from the interior of the disk, and extends in a bipolar outflow over 1000 AU along the rotation axis of the disk. SiO $v=0$ $J=5-4$ emission appears to have a localized footprint in the warped ridge. These observations suggest that the ridge is the working surface of the disk, and heated by accretion and the outflow. The disk structure may be evolving, with multiple accretion and outflow events. We discuss two sources of variability: 1) variable accretion onto the disk as Source I travels through the filamentary debris from the BN-Source I encounter $\sim$550 yr ago; and 2) episodic accretion from the disk onto the protostar which may trigger multiple outflows. The warped inner disk structure is direct evidence that SrcI could be a binary experiencing episodic accretion., Comment: Accepted in the Astrophysical Journal

An extended thermal particle boundary condition is devised to more efficiently and accurately model laser-plasma interactions in overdense plasmas. Particle-in-cell simulations of such interactions require many particles per cell, and a large region of background plasma is often necessary to correctly mimic a semi-infinite plasma and avoid electron refluxing from a truncated plasma. For long-pulse lasers of many picoseconds, such constraints can become prohibitively expensive. Here, an extended particle boundary condition (absorber) is designed that instantaneously stops and re-emits energetic particles streaming toward the simulation boundary over a defined region, allowing sufficient time and space for a suitably cool return current to develop in the background plasma. Tunable parameters of the absorber are explained, and simulations using the absorber with a 3-ps laser are shown to accurately reproduce those of a causally separated boundary while requiring only 20% the number of particles., Comment: 12 pages, 9 figures

We formalize Hall's Marriage Theorem in the Lean theorem prover for inclusion in mathlib, which is a community-driven effort to build a unified mathematics library for Lean. One goal of the mathlib project is to contain all of the topics of a complete undergraduate mathematics education. We provide three presentations of the main theorem statement: in terms of indexed families of finite sets, of relations on types, and of matchings in bipartite graphs. We also formalize a version of K\H{o}nig's lemma (in terms of inverse limits) to boost the theorem to the case of countably infinite index sets. We give a description of the design of the recent mathlib library for simple graphs, and we also give a necessary and sufficient condition for a simple graph to carry a function., Comment: 15 pages

We show using particle-in-cell (PIC) simulations and theoretical analysis that a high-quality electron beam whose density is modulated at angstrom scales can be generated directly using density downramp injection in a periodically modulated density in nonlinear plasma wave wakefields. The density modulation turns on and off the injection of electrons at the period of the modulation. Due to the unique longitudinal mapping between the electrons' initial positions and their final trapped positions inside the wake, this results in an electron beam with density modulation at a wavelength orders of magnitude shorter than the plasma density modulation. The ponderomotive force of two counter propagating lasers of the same frequency can generate a density modulation at half the laser wavelength. Assuming a laser wavelength of $0.8\micro\meter$, fully self-consistent OSIRIS PIC simulations show that this scheme can generate high quality beams modulated at wavelengths between 10s and 100 angstroms. Such beams could produce fully coherent, stable, hundreds of GW X-rays by going through a resonant undulator., Comment: 6 pages, 4 figures

Metal-insulator transition (MIT) compounds are materials that may exhibit insulating or metallic behavior, depending on the physical conditions, and are of immense fundamental interest owing to their potential applications in emerging microelectronics. There is a dearth of thermally-driven MIT materials, however, which makes delineating these compounds from those that are exclusively insulating or metallic challenging. Here we report a material database comprising temperature-controlled MITs (and metals and insulators with similar chemical composition and stoichiometries to the MIT compounds) from high quality experimental literature, built through a combination of materials-domain knowledge and natural language processing. We featurize the dataset using compositional, structural, and energetic descriptors, including two MIT relevant energy scales, an estimated Hubbard interaction and the charge transfer energy, as well as the structure-bond-stress metric referred to as the global-instability index (GII). We then perform supervised classification, constructing three electronic-state classifiers: metal vs non-metal (M), insulator vs non-insulator (I), and MIT vs non-MIT (T). We identify two important descriptors that separate metals, insulators, and MIT materials in a 2D feature space: the average deviation of the covalent radius and the range of the Mendeleev number. We further elaborate on other important features (GII and Ewald energy), and examine how they affect classification of binary vanadium and titanium oxides. We discuss the relationship of these atomic features to the physical interactions underlying MITs in the rare-earth nickelate family. Last, we implement an online version of the classifiers, enabling quick probabilistic class predictions by uploading a crystallographic structure file.

Next-generation high-power lasers that can be focused to intensities exceeding 10^23 W/cm^2 are enabling new physics and applications. The physics of how these lasers interact with matter is highly nonlinear, relativistic, and can involve lowest-order quantum effects. The current tool of choice for modeling these interactions is the particle-in-cell (PIC) method. In strong fields, the motion of charged particles and their spin is affected by radiation reaction. Standard PIC codes usually use Boris or its variants to advance the particles, which requires very small time steps in the strong-field regime to obtain accurate results. In addition, some problems require tracking the spin of particles, which creates a 9D particle phase space (x, u, s). Therefore, numerical algorithms that enable high-fidelity modeling of the 9D phase space in the strong-field regime are desired. We present a new 9D phase space particle pusher based on analytical solutions to the position, momentum and spin advance from the Lorentz force, together with the semi-classical form of RR in the Landau-Lifshitz equation and spin evolution given by the Bargmann-Michel-Telegdi equation. These analytical solutions are obtained by assuming a locally uniform and constant electromagnetic field during a time step. The solutions provide the 9D phase space advance in terms of a particle's proper time, and a mapping is used to determine the proper time step for each particle from the simulation time step. Due to the analytical integration, the constraint on the time step needed to resolve trajectories in ultra-high fields can be greatly reduced. We present single-particle simulations and full PIC simulations to show that the proposed particle pusher can greatly improve the accuracy of particle trajectories in 9D phase space for given laser fields. A discussion on the numerical efficiency of the proposed pusher is also provided.

This paper reviews current literature in the field of predictive maintenance from the system point of view. We differentiate the existing capabilities of condition estimation and failure risk forecasting as currently applied to simple components, from the capabilities needed to solve the same tasks for complex assets. System-level analysis faces more complex latent degradation states, it has to comprehensively account for active maintenance programs at each component level and consider coupling between different maintenance actions, while reflecting increased monetary and safety costs for system failures. As a result, methods that are effective for forecasting risk and informing maintenance decisions regarding individual components do not readily scale to provide reliable sub-system or system level insights. A novel holistic modeling approach is needed to incorporate available structural and physical knowledge and naturally handle the complexities of actively fielded and maintained assets., Comment: 24 pages, 3 figures

A customized finite-difference field solver for the particle-in-cell (PIC) algorithm that provides higher fidelity for wave-particle interactions in intense electromagnetic waves is presented. In many problems of interest, particles with relativistic energies interact with intense electromagnetic fields that have phase velocities near the speed of light. Numerical errors can arise due to (1) dispersion errors in the phase velocity of the wave, (2) the staggering in time between the electric and magnetic fields and between particle velocity and position and (3) errors in the time derivative in the momentum advance. Errors of the first two kinds are analyzed in detail. It is shown that by using field solvers with different $\mathbf{k}$-space operators in Faraday's and Ampere's law, the dispersion errors and magnetic field time-staggering errors in the particle pusher can be simultaneously removed for electromagnetic waves moving primarily in a specific direction. The new algorithm was implemented into OSIRIS by using customized higher-order finite-difference operators. Schemes using the proposed solver in combination with different particle pushers are compared through PIC simulation. It is shown that the use of the new algorithm, together with an analytic particle pusher (assuming constant fields over a time step), can lead to accurate modeling of the motion of a single electron in an intense laser field with normalized vector potentials, $eA/mc^2$, exceeding $10^4$ for typical cell sizes and time steps.

Furthering our understanding of many of today's interesting problems in plasma physics---including plasma based acceleration and magnetic reconnection with pair production due to quantum electrodynamic effects---requires large-scale kinetic simulations using particle-in-cell (PIC) codes. However, these simulations are extremely demanding, requiring that contemporary PIC codes be designed to efficiently use a new fleet of exascale computing architectures. To this end, the key issue of parallel load balance across computational nodes must be addressed. We discuss the implementation of dynamic load balancing by dividing the simulation space into many small, self-contained regions or "tiles," along with shared-memory (e.g., OpenMP) parallelism both over many tiles and within single tiles. The load balancing algorithm can be used with three different topologies, including two space-filling curves. We tested this implementation in the code OSIRIS and show low overhead and improved scalability with OpenMP thread number on simulations with both uniform load and severe load imbalance. Compared to other load-balancing techniques, our algorithm gives order-of-magnitude improvement in parallel scalability for simulations with severe load imbalance issues., Comment: 13 pages, 7 figures, submitted to Computer Physics Communications

In Dunfield's catalog of the hyperbolic manifolds in the SnapPy census which are complements of L-space knots in $S^3$, we determine that $22$ have tunnel number $2$ while the remaining all have tunnel number $1$. Notably, these $22$ manifolds contain $9$ asymmetric L-space knot complements. Furthermore, using SnapPy and KLO we find presentations of these $22$ knots as closures of positive braids that realize the Morton-Franks-Williams bound on braid index. The smallest of these has genus $12$ and braid index $4$., Comment: This is a significant expansion/rewrite of the original version. The title and abstract have both changed to reflect this. Relevant code is included in the body of the article itself. Originally, this was a note that collected the results of an "office hour" session on the use of SnapPy and KLO held at the ICERM workshop Perspectives on Dehn Surgery, July 15--19, 2019

In order to apply quantum topology methods to nonplanar graphs, we define a planar diagram category that describes the local topology of embeddings of graphs into surfaces. These \emph{virtual graphs} are a categorical interpretation of ribbon graphs. We describe an extension of the flow polynomial to virtual graphs, the $S$-polynomial, and formulate the $\mathfrak{sl}(N)$ Penrose polynomial for non-cubic graphs, giving contraction-deletion relations. The $S$-polynomial is used to define an extension of the Yamada polynomial to virtual spatial graphs, and with it we obtain a sufficient condition for non-classicality of virtual spatial graphs. We conjecture the existence of local relations for the $S$-polynomial at squares of integers., Comment: 37 pages, 27 figures