Probing the scale of grand unification with gravitational waves
Not Available
SEARCH · Engineering Papers
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Not Available
Not Available
The tetragonal intermetallic compound EuAl 4 hosts an exciting variety of low-temperature phases. In addition to a charge density wave below 140 K, four ordered magnetic phases are observed below 15.4 K. Recently, a skyrmion phase was proposed based on Hall effect measurements under a c-axis magnetic field. Here, we present a detailed investigation of the phase transitions in EuAl 4 under c-axis magnetic field. Our dilatometry, heat-capacity, DC magnetometry, AC magnetic susceptibility, and resonant ultrasound spectroscopy measurements reveal three magnetic phase transitions not previously reported. The first key result is a detailed H ∥ [001] magnetic phase diagram mapping the seven phases we observe. Second, we identify a high-field phase, phase VII, which directly corresponds to the region were skyrmions have been previously suggested. Our results provide guidance for future studies exploring the complex magnetic interactions and spin structures in EuAl 4 .
In this paper, we describe the development of a reactive force field for C/O systems under extreme temperatures and pressures, based on the many-body Chebyshev Interaction Model for Efficient Simulation (ChIMES). The resulting model, which targets carbon condensation under thermodynamic conditions of 6500 K and 2.5 g cm –3 , affords a balance between model accuracy, complexity, and training set generation expense. We show that the model recovers much of the accuracy of density functional theory for the prediction of structure, dynamics, and chemistry when applied to dissociative condensed phase systems at 1:1 and 1:2 C:O ratios, as well as molten carbon. Our C/O modeling approach exhibits a 104 increase in efficiency for the same system size (i.e., 128 atoms) and a linear system size scalability over standard quantum molecular dynamics methods, allowing the simulation of significantly larger systems than previously possible. We find that the model captures the condensed-phase reaction-coupled formation of carbon clusters implied by recent experiments, and that this process is susceptible to strong finite size effects. Overall, we find the present ChIMES model to be well suited for studying chemical processes and cluster formation at pressures and temperatures typical of shock waves. We expect that the present C/O modeling paradigm can serve as a template for the development of a broader high pressure–high temperature force-field for condensed phase chemistry in organic materials.
This final technical report summarizes the work of the Michigan State University (MSU) team supported by grant # DE-SC0021423 awarded by the Office of Nuclear Physics, Department of Energy. Objectives: The successful fulfillment of the FRIB science mission hinges on ensuring the availability of fast, stopped, and reaccelerated beams consisting of rare isotopes. This project's research and development focus was dedicated to supporting the advancement and creation of a cutting-edge linear gas stopper. The primary aim is to efficiently convert the high-intensity fast beams of rare isotopes provided by FRIB into high-quality, low-energy beams. These beams are essential for conducting stopped beam experiments or for subsequent reacceleration. The overarching goal is to advance technology, aiming to increase the beam rate capability of the linear gas stopper for medium-to-heavy-mass rare isotopes by more than tenfold compared to the currently most effective gas stopper in operation, and to improve the purity of the delivered beams. Project Description: The existing technology employed in gas stopping devices designed for low-energy, rare-isotope beams presents limitations in both the purity of extracted beams and the intensities of injected beams. These limitations are incompatible with the requirements of the recently commissioned rare isotope beam facility, FRIB. Our research and development efforts, aligned with the previously outlined objectives, focused on addressing the most critical aspects for enhancing beam-rate capability and purity. Specifically, advanced particle-in-cell simulations were developed and integrated into a simulation pipeline to explore the efficacy of multi-layer RF carpets on increasing ion transport efficiency with high incoming beam rates that generate space charge fields which can limit it. We also explored the possibility of using a collision-induced-dissociation (CID) gas cell to break up molecular contaminant ions that are generated during the stopping process. A prototype CID gas cell was constructed and tested with beams from an offline ion source, validating the concept with the successful demonstration of breaking of molecular ions. The outcome of this research enabled the formulation a conceptual design for a next-generation linear gas stopping device specifically tailored for FRIB. This device is envisioned to deliver rare-isotope-ion beams at a rate of 10 8 particles per second or higher, accompanied by advancements in purity. Methods employed: This project leverages advancements in technologies initially designed for the Advanced Cryogenic Gas Stopper (ACGS), the current state-of-the-art linear gas stopper, through the use of new simulations and beam purification via collision-induced-dissociation. The methods include: 1. Development of a prototype low-energy, low-pressure CID gas stopper. This prototype features a thin, approximately 20 nm, Si 3 N 4 entrance window designed for dissociating stable and rare-isotope molecular ions. The goal is to enhance beam purification and overall efficiency. 2. Creation of Particle-in-cell (PIC) simulations to assess the advantages of multi-layer RF carpets and multi-point extraction for ion transport efficiency. These simulations rely on the 3DCylPIC package, specifically designed for studying devices of this nature. The goal is to quantify and mitigate ion transport losses due to space charge generated in the stopping process of large numbers of ions. 3. Perform ion transport simulations across an RF carpet using an 8-phase travelling wave and evaluate its performance. Compared to the 4-phase RF carpets used in ACGS, the 8-phase carpets will double the wavelength of the generated traveling wave allowing for larger maximum RF amplitudes that could result in improved ion transport efficiency for high-intensity incoming beams when large space charge fields are present. Impact: Tackling the primary challenges associated with transforming high-energy projectile fragment beams into low-energy beams—specifically, addressing efficiency and purity—holds significant promise for advancing FRIB science. This advancement will particularly impact precision mass measurements, laser spectroscopy of short-lived nuclei, and studies in astrophysics and nuclear reactions using reaccelerated beams. These domains play a crucial role in addressing key questions outlined in the 2023 NSAC long-range plan, spanning nuclear structure, nuclear astrophysics, and fundamental symmetries. Additionally, they contribute to addressing 10 out of the 17 benchmarks identified by the NRC RISAC. The development of a next-generation gas stopping device capable of delivering low-energy, rare-isotope beams at a rate of 10 8 particles per second, or more, with high purity holds the potential to unlock experiments that would otherwise be unfeasible. Furthermore, it is expected to reduce the time required for experiments at FRIB, thereby maximizing scientific output. The research and development activities performed as part of this project bolstered essential competencies at FRIB in beam physics and ion source technologies, provided valuable training opportunities for junior scientists.
We revisit the problem of small Bjorken-x evolution of the gluon and flavor-singlet quark helicity distributions in the shock wave (s-channel) formalism. Earlier works on the subject in the same framework resulted in an evolution equation for the gluon field-strength F 12 and quark “axial current” $\overline{\psi}\gamma$ + γ 5 ψ operators (sandwiched between the appropriate light-cone Wilson lines) in the double-logarithmic approximation (summing powers of α s ln 2 (1/x) with α s the strong coupling constant). In this work, we observe that an important mixing of the above operators with another gluon operator, ${}_D{}^{\leftarrow i} {}_D{}^{i}$, also sandwiched between the light-cone Wilson lines (with the repeated transverse index i = 1, 2 summed over), was missing in the previous works. This operator has the physical meaning of the sub-eikonal (covariant) phase: its contribution to helicity evolution is shown to be proportional to another sub-eikonal operator, ${}_D{}^{i} {}_D{}^{\leftarrow i}$, which is related to the Jaffe-Manohar polarized gluon distribution. In this work we include this new operator into small-x helicity evolution, and construct novel evolution equations mixing all three operators (${}_D{}^{i} {}_D{}^{\leftarrow i}$, F 12 , and $\overline{\psi}\gamma$ + γ 5 ψ), generalizing the results of . We also construct closed double-logarithmic evolution equations in the large-N c and large-N c &N f limits, with N c and N f the numbers of quark colors and flavors, respectively. Solving the large-N c equations numerically we obtain the following small-x asymptotics of the quark and gluon helicity distributions ΔΣ and ΔG, along with the g 1 structure function, $$ \Delta \Sigma \left(x,{Q}^2\right)\sim \Delta G\left(x,{Q}^2\right)\sim {g}_1\left(x,{Q}^2\right)\sim {\left(\frac{1}{x}\right)}^{3.66\sqrt{\frac{\alpha_s{N}_c}{2\pi }}} $$ in complete agreement with the earlier work by Bartels, Ermolaev and Ryskin.
The computation of two-electron repulsion integrals (ERIs) is often the most expensive step of integral-direct self-consistent field methods. Formally it scales as O(N 4 ), where N is the number of Gaussian basis functions used to represent the molecular wave function. In practice, this scaling can be reduced to O(N 2 ) or less by neglecting small integrals with screening methods. The contributions of the ERIs to the Fock matrix are of Coulomb (J) and exchange (K) type and require separate algorithms to compute matrix elements efficiently. We previously implemented highly efficient GPU-accelerated J-matrix and K-matrix algorithms in the electronic structure code TeraChem. Although these implementations supported the use of multiple GPUs on a node, they did not support the use of multiple nodes. This presents a key bottleneck to cutting-edge ab initio simulations of large systems, e.g., excited state dynamics of photoactive proteins. We present our implementation of multinode multi-GPU J- and K-matrix algorithms in TeraChem using the Regent programming language. Regent directly supports distributed computation in a task-based model and can generate code for a variety of architectures, including NVIDIA GPUs. We demonstrate multinode scaling up to 45 GPUs (3 nodes) and benchmark against hand-coded TeraChem integral code. Finally, we also outline our metaprogrammed Regent implementation, which enables flexible code generation for integrals of different angular momenta.
A three-dimensional breakdown behavior of the magnetic Bloch-point domain wall (BP-DW) propagation was investigated in a ferromagnetic nanowire with the variation of the external magnetic field by means of micromagnetic simulation. As magnetic field strength increases up to a threshold value, the BP-DW velocity approaches a critical phase velocity of the spontaneously emitting spin wave (SW), where a Cherenkov-type DW breakdown phenomenon is observed originating from an interaction between the spontaneously emitting SW and the BP-DW. It is found that the velocity of the BP-DW approaches a maximum value (~2000 m/s) due to the intrinsic reversal time of a BP spin texture. This suggests that although cylindrical ferromagnetic nanowires might be free from the two-dimensional Walker breakdown phenomenon, there exists a Cherenkov-type three-dimensional breakdown behavior.
Abstract The interaction of a relativistic laser pulse ( > 10 18 W cm −2 µ m 2 ) with foil targets can accelerate ions to energies of tens of MeV u −1 with optimized laser and target parameters. We report the results on simulations of the interaction of a 6.0 × 10 20 W cm −2 laser pulse incident on ultrathin (10–500 nm) titanium foils to investigate the roles of laser pulse duration and ionization mechanisms in the acceleration of titanium ions. While holding peak intensity constant, two laser pulse durations were investigated, 140 and 650 fs. The optimum thickness is dependent on pulse duration, as it requires the concurrence of target transparency with the incidence of the peak laser intensity. The collisional processes do not play a significant role in Ti ion beam generation from the 140 fs laser pulse duration at the optimum thickness (30 nm). However, for the 650 fs laser optimum, collisions improve the conversion efficiency of highly energetic ( > 10 MeV u −1 ), high charge titanium (Ti 21 − 22 + ) by a factor of 20, and the titanium ion cutoff energy by ∼15%. This improvement is due to the fact that collisional ionization increases the electron density of the plasma, which delays the time of relativistic transparency, causing collisions to decrease the optimum foil thickness from 150 to 100 nm. Additionally, collisional ionization increases the charge-to-mass ratio of the titanium, and injects more electrons into the accelerating sheath field. At the optimum thickness, target normal sheath acceleration is the dominant mechanism of acceleration, with additional contributions from radiation pressure and shock wave acceleration.
Abstract Scattering amplitudes are both a wonderful playground to discover novel ideas in quantum field theory and simultaneously of immense phenomenological importance to make precision predictions for e.g. particle collider observables and more recently also for gravitational wave signals. In this review chapter, we give an overview of some of the exciting recent progress on reformulating QFT in terms of mathematical, geometric quantities, such as polytopes, associahedra, Grassmanians, and the amplituhedron. In this novel approach, standard notions of locality and unitarity are derived concepts rather than fundamental ingredients in the construction which might give us a handle on a number of open questions in QFT that have evaded an answer for decades. We first give a basic summary of positive geometry before discussing the associahedron—one of the simplest physically relevant geometric examples—and its relation to tree-level scattering amplitudes in bi-adjoint ϕ 3 theory. Our second example is the amplituhedron construction for scattering amplitudes in planar maximally supersymmetric Yang–Mills theory.
Axions are some of the best motivated particles beyond the Standard Model. We show how the attractive self-interactions of dark matter (DM) axions over a broad range of masses, from 10 –22 eV to 10 7 GeV, can lead to nongravitational growth of density fluctuations and the formation of bound objects. This structure formation enhancement is driven by parametric resonance when the initial field misalignment is large, and it affects axion density perturbations on length scales of order the Hubble horizon when the axion field starts oscillating, deep inside the radiation-dominated era. This effect can turn an otherwise nearly scale-invariant spectrum of adiabatic perturbations into one that has a spike at the aforementioned scales, producing objects ranging from dense DM halos to scalar-field configurations such as solitons and oscillons. We call this class of cosmological scenarios for axion DM production “the large-misalignment mechanism.” In this work, we explore observational consequences of this mechanism for axions with masses up to 10 eV. For axions heavier than 10 –5 eV, the compact axion halos are numerous enough to significantly impact Earth-bound direct detection experiments, yielding intermittent but coherent signals with repetition rates exceeding one per decade and crossing times less than a day. These episodic increases in the axion density and kinematic coherence suggest new approaches for axion DM searches, including for the QCD axion. Dense structures made up of axions from 10 –22 eV to 10 –5 eV are detectable through gravitational lensing searches, and their gravitational interactions can also perturb baryonic structures and alter star formation. At very high misalignment amplitudes, the axion field can undergo self-interaction-induced implosions long before matter-radiation equality, producing potentially-detectable low-frequency stochastic gravitational waves.
We present universal relations between entanglement entropy, which quantifies the quantum correlation between subsystems, and the cross section, which is the primary observable for high-energy particle scattering, by employing a careful formulation of wave packets for the incoming particles. For 2-to-2 elastic scattering with no initial entanglement and subdividing the system along particle labels, we show that both the Rényi and Tsallis entropies in the final states are directly proportional to the elastic cross section in units of the transverse size for the initial wave packets, which is then interpreted as the elastic scattering probability. The relations do not depend on the underlying dynamics of the quantum field theory and are valid to all orders in coupling strengths. Furthermore, computing quantum correlations between momentum and nonkinematic data leads to entanglement entropies expressed as various semi-inclusive elastic cross sections. Our result gives rise to a novel “area law” for entanglement entropy in a two-body system. Published by the American Physical Society 2025
A multi-mode approach for real-time inspection of structural components may be applied to rapid, wide area measurement of thickness of thick plate-like structures using full-field multi-mode response measurement and analysis. The approach may allow estimation of thickness from full-field multi-mode response to single-tone ultrasonic excitation in thick plates. The approach may utilize wavenumber information across all available wave modes in order to make wavenumber spectroscopy sensitive to changes in thickness for a broader range of nominal initial thicknesses.
In addition to a plethora of emergent phenomena, the spatial topology of optical vortices enables an array of applications in optical communications and quantum information science. Multibeam nonlinear optical processes, augmented by optical vortices, are essential in this context, providing robust access to an infinitely large set of quantum states associated with the orbital angular momentum of light. Here, we push the boundaries of vortex nonlinear optics to the ultimate limits of material dimensionality. By exploiting multipulse difference frequency, sum frequency, and four-wave mixing in monolayer quantum materials, we demonstrate their ability to independently control the orbital angular momentum and radial distribution of vortex light-fields in addition to their wavelength. Due to the atomically thin nature of the host crystal, this control spans a broad spectral bandwidth in a highly integrable platform that is unconstrained by the traditional limits of bulk nonlinear optical materials. Our work heralds an innovative path for ultracompact and scalable hybrid nanophotonic technologies empowered by twisted nonlinear light–matter interactions in van der Waals nanomaterials.
There is controversy surrounding the origin and evolution of our universe's largest supermassive black holes (SMBHs). In this study, we consider the possibility that some of these black holes formed from the direct collapse of primordial density perturbations. Since the mass of a primordial black hole is limited by the size of the cosmological horizon at the time of collapse, these SMBHs must form rather late, and are naively in conflict with constraints from CMB spectral distortions. These limits can be avoided, however, if the distribution of primordial curvature perturbations is highly non-Gaussian. After quantifying the departure from Gaussianity needed to evade these bounds, we explore a model of multi-field inflation — a non-minimal, self-interacting curvaton model — which has all the necessary ingredients to yield such dramatic non-Gaussianities. We leave the detailed model building and numerics to a future study, however, as our goal is to highlight the challenges associated with forming SMBHs from direct collapse and to identify features that a successful model would need to have. This study is particularly timely in light of recent observations of high-redshift massive galaxy candidates by the James Webb Space Telescope as well as evidence from the NANOGrav experiment for a stochastic gravitational wave background consistent with SMBH mergers.
There is a deep connection between the ground states of transverse-field spin systems and the late-time distributions of evolving viral populations—within simple models, both are obtained from the principal eigenvector of the same matrix. However, that vector is the wave-function amplitude in the quantum spin model, whereas it is the probability itself in the population model. We show that this seemingly minor difference has significant consequences: Phase transitions that are discontinuous in the spin system become continuous when viewed through the population perspective, and transitions that are continuous become governed by new critical exponents. We introduce a more general class of models that encompasses both cases and that can be solved exactly in a mean-field limit. Numerical results are also presented for a number of one-dimensional chains with power-law interactions. We see that well-worn spin models of quantum statistical mechanics can contain unexpected new physics and insights when treated as population-dynamical models and beyond, motivating further studies.
Here, we combine hydrodynamic and field theoretic methods to develop a general theory of phonons as Goldstone bosons in crystals, glasses, and liquids based on nonaffine displacements and the consequent Goldstone phase relaxation. We relate the conservation, or lack thereof, of specific higher-form currents with properties of the underlying deformation field—nonaffinity—which dictates how molecules move under an applied stress or deformation. In particular, the single-valuedness of the deformation field is associated with conservation of higher-form charges that count the number of topological defects. Our formalism predicts, from first principles, the presence of propagating shear waves above a critical wave vector in liquids, thus giving a formal derivation of the phenomenon in terms of fundamental symmetries. The same picture provides also a theoretical explanation of the corresponding “positive sound dispersion” phenomenon for longitudinal sound. Importantly, accordingly to our theory, the main collective relaxation timescale of a liquid or a glass (known as the α relaxation for the latter) is given by the phase relaxation time, which is not necessarily related to the Maxwell time. Finally, we build a nonequilibrium effective action using the in-in formalism defined on the Schwinger-Keldysh contour, that further supports the emerging picture. In summary, our work suggests that the fundamental difference between solids, fluids, and glasses has to be identified with the associated generalized higher-form global symmetries and their topological structure, and that the Burgers vector for the displacement fields serves as a suitable topological order parameter distinguishing the solid (ordered) phase and the amorphous ones (fluids, glasses).
Power electronic converters are an enabling technology for the emerging marine energy applications, such as using ocean waves to produce electricity. This paper outlines the power take-off system and its key components used in a wave energy converter offering modularity and scalability to generate power efficiently. The proposed power take-off system was implemented based on a modular multilevel converter and could be deployed to convert any alternating current electrical energy to a different alternating current for interconnection to grid or non-grid applications. Examples of widespread deployment are supplying electricity to coastal communities or producing clean drinking water. The analysis using both the simulation tests and laboratory experiments verified the design objectives and basic functionality of the developed power take-off system. An acceptable response using a field programmable gate array-based controlled laboratory testbench was achieved, complying with guidelines specified in the prevalent industry standards. Seamless operation during steady-state and transients for the studied wave energy converter was achieved as supported by the obtained results. The key findings of this work were experimentally examined under different load conditions, direct current bus voltage fluctuations, and generator speed–torque regulation. The ability of the power take-off system to generate high-power quality of the waveforms, e.g., against adhering to the IEEE 519-2022 standard for total harmonic distortion limits, is also confirmed.