SEARCH · Search NASA

Results for “distributed quantum computing”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,027 records · Page 57

All sources NASA NTRS DOE OSTI

DOE OSTI2020

ATOMIC Simulations and Experimental Data for Basalt-like Compounds

This data consists of simulations and experimental measurements of laser-induced breakdown spectroscopy (LIBS). The simulations are produced by ATOMIC, a general purpose plasma modeling and kinetics code that has been designed to compute emission (or absorption) spectra from plasmas [2]. The makeup of the plasma was considered to be divided into some proportion water, some proportion Martian atmosphere (CO2), and some proportion target (from the rock or object impacted by the laser), where these proportions add to 1. Based on expert knowledge, the proportion of water was kept in the range [0.0,0.5] and the proportion of atmosphere was kept in the range [0.02, 0.9]. Our overall suite of simulations contains six sets of simulations that differ in which elements were considered to make up the target. Within each set, we used uniformly drawn temperatures and log mass densities within pre-specified ranges. The temperature range was [0.5,1.5] eV and the log (base 10) mass density range was [-7,-4]. The proportion of water, atmosphere, and target were drawn from a symmetric Dirichlet distribution, but draws in which the propor- tion of water or atmosphere exceeded the pre-specified limits were rejected from the design. Up to eleven constituent elements (Si, Al, Fe, Mg, Ca, O, Ti, Mn, Na, K, P) were considered for the target, as they are the most common elements found in basalt compounds and were used in [1]. For each run, the proportions of the constituent elements making up the target were drawn from a symmetric Dirichlet distribution. We ran 1,350 simulations that included nonzero proportions of all eleven elements. We also ran simulations which excluded some of these elements. In particular, we ran 1,000 simulations that only included nonzero proportions for the six most common elements (Si, Al, Fe, Mg, Ca, O). We also ran five sets, each with 500 simulations, that only included nonzero proportions for five of the six most common elements (but where all sets included O). Thus, we generated a total of 4,850 spectra representing basalt-like compounds in which the target was comprised of oxygen and between four and ten other elements. The ATOMIC code produced spectra over a range of 240nm - 880nm that roughly mimics the range collected by the ChemCam instrument on the Mars rover Curiosity. Each spectra had 32,000 wavelengths split across three spectrometer ranges (to mimic ChemCam). The experimental data, described in [1], measures a prepared basalt sample. All files are kept in directories whose names indicate the set of elements considered for the target with file names numbered to indicate the line in the design files used to produce the simulation. The designs are provided as text files with names indicating their purpose. The experimental data is provided as a CSV file which contains a header with measurement information, followed by a collection of 50 shots across a collection of wavelengths, along with the computed median and mean across shots.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2014

Migration of finite sized particles in a laminar square channel flow from low to high Reynolds numbers

The migration of neutrally buoyant finite sized particles in a Newtonian square channel flow is investigated in the limit of very low solid volumetric concentration, within a wide range of channel Reynolds numbers Re = [0.07-120]. In situ microscope measurements of particle distributions, taken far from the channel inlet (at a distance several thousand times the channel height), revealed that particles are preferentially located near the channel walls at Re > 10 and near the channel center at Re < 1. Whereas the cross-streamline particle motion is governed by inertia-induced lift forces at high inertia, it seems to be controlled by shear-induced particle interactions at low (but finite) Reynolds numbers, despite the low solid volume fraction (<1%). The transition between both regimes is observed in the range Re = [1-10]. In order to exclude the effect of multi-body interactions, the trajectories of single freely moving particles are calculated thanks to numerical simulations based on the force coupling method. With the deployed numerical tool, the complete particle trajectories are accessible within a reasonable computational time only in the inertial regime (Re > 10). In this regime, we show that (i) the particle undergoes cross-streamline migration followed by a cross-lateral migration (parallel to the wall) in agreement with previous observations, and (ii) the stable equilibrium positions are located at the midline of the channel faces while the diagonal equilibrium positions are unstable. At low flow inertia, the first instants of the numerical simulations (carried at Re = O(1)) reveal that the cross-streamline migration of a single particle is oriented towards the channel wall, suggesting that the particle preferential positions around the channel center, observed in the experiments, are rather due to multi-body interactions.

42 ENGINEERING↗

DOE OSTI2025

Orientation microscopy–assisted grain boundary analysis for protonic ceramic cell electrolytes

Abstract Grain boundaries in protonic ceramic cell (PCC) electrolytes hinder proton transport, reducing interfacial conductivity. In multicomponent PCC electrolytes, the inclusion of sintering aids further accentuates the complexity of grain boundaries. In this study, we synthesize nanocrystalline BaCe 0.4 Zr 0.4 Y 0.1 Yb 0.1 O 3− δ thin films via pulsed laser deposition and analyze their grain boundary character distributions using orientation data collected by precession electron diffraction technique. The results reveal an anisotropic distribution of grain boundary characters, with notably high populations of 180°‐tilt and twist grain boundaries. These findings provide critical insights into identifying the predominant grain boundaries in this PCC electrolyte material, assessing the vast five‐dimensional grain boundary space.

Patel, Sooraj [School of Aerospace and Mechanical ↗

DOE OSTI2022

Path integral Monte Carlo approach to the structural properties and collective excitations of liquid $$^3{\text {He}}$$ without fixed nodes

Abstract Due to its nature as a strongly correlated quantum liquid, ultracold helium is characterized by the nontrivial interplay of different physical effects. Bosonic $$^4{\text {He}}$$ 4 He exhibits superfluidity and Bose-Einstein condensation. Its physical properties have been accurately determined on the basis of ab initio path integral Monte Carlo (PIMC) simulations. In contrast, the corresponding theoretical description of fermionic $$^3{\text {He}}$$ 3 He is severely hampered by the notorious fermion sign problem, and previous PIMC results have been derived by introducing the uncontrolled fixed-node approximation. In this work, we present extensive new PIMC simulations of normal liquid $$^3{\text {He}}$$ 3 He without any nodal constraints. This allows us to to unambiguously quantify the impact of Fermi statistics and to study the effects of temperature on different physical properties like the static structure factor $$S({\mathbf {q}})$$ S ( q ) , the momentum distribution $$n({\mathbf {q}})$$ n ( q ) , and the static density response function $$\chi ({\mathbf {q}})$$ χ ( q ) . In addition, the dynamic structure factor $$S({\mathbf {q}},\omega )$$ S ( q , ω ) is rigorously reconstructed from imaginary-time PIMC data. From simulations of $$^3{\text {He}}$$ 3 He , we derived the familiar phonon–maxon–roton dispersion function that is well-known for $$^4{\text {He}}$$ 4 He and has been reported previously for two-dimensional $$^3{\text {He}}$$ 3 He films (Nature 483:576–579 (2012)). The comparison of our new results for both $$S({\mathbf {q}})$$ S ( q ) and $$S({\mathbf {q}},\omega )$$ S ( q , ω ) with neutron scattering measurements reveals an excellent agreement between theory and experiment.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2016

Scalable posterior approximations for large-scale Bayesian inverse problems via likelihood-informed parameter and state reduction

Two prime bottlenecks to the solution of large-scale Bayesian inverse problems are the scaling of posterior sampling algorithms to high-dimensional parameter spaces and the computational cost of forward model evaluations. However, incomplete or noisy data, the state variation and parameter dependence of the forward model, and correlations in the prior collectively provide useful structure that can be exploited for dimension reduction in this setting—both in the parameter space of the inverse problem and in the state space of the forward model. To this end, we show how to jointly construct low-dimensional subspaces of the parameter space and the state space in order to accelerate the Bayesian solution of the inverse problem. As a byproduct of state dimension reduction, we also show how to identify low-dimensional subspaces of the data in problems with high-dimensional observations. These subspaces enable approximation of the posterior as a product of two factors: (i) a projection of the posterior onto a low-dimensional parameter subspace, wherein the original likelihood is replaced by an approximation involving a reduced model; and (ii) the marginal prior distribution on the high-dimensional complement of the parameter subspace. We introduce and compare several strategies for constructing these subspaces using only a limited number of forward and adjoint model simulations. The resulting posterior approximations can rapidly be characterized using standard sampling techniques, e.g., Markov chain Monte Carlo. Two numerical examples demonstrate the accuracy and efficiency of our approach: inversion of an integral equation in atmospheric remote sensing, where the data dimension is very high; and the inference of a heterogeneous transmissivity field in a groundwater system, which involves a partial differential equation forward model with high dimensional state and parameters.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2021

Security proof of practical quantum key distribution with detection-efficiency mismatch

Quantum key distribution (QKD) protocols with threshold detectors are driving high-performance QKD demonstrations. The corresponding security proofs usually assume that all physical detectors have the same detection efficiency. However, the efficiencies of the detectors used in practice might show a mismatch depending on the manufacturing and setup of these detectors. A mismatch can also be induced as the different spatial-temporal modes of an incoming signal might couple differently to a detector. Here we develop a method that allows to provide security proofs without the usual assumption. Our method can take the detection-efficiency mismatch into account without having to restrict the attack strategy of the adversary. Especially, we do not rely on any photon-number cutoff of incoming signals such that our security proof is directly applicable to practical situations. We illustrate our method for a receiver that is designed for polarization encoding and is sensitive to a number of spatial-temporal modes. In our detector model, the absence of quantum interference between any pair of spatial-temporal modes is assumed. For a QKD protocol with this detector model, we can perform a security proof with characterized efficiency mismatch and without photon-number cutoff assumption. Our method also shows that in the absence of efficiency mismatch in our detector model, the key rate increases if the loss due to detection inefficiency is assumed to be outside of the adversary's control, as compared to the view where for a security proof this loss is attributed to the action of the adversary.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2019

Next generation of extreme-resolution electron beam lithography

In top down nanofabrication research facilities around the world, the direct-write high-resolution patterning tool of choice is overwhelmingly electron beam lithography. Remarkably small features can be written in a variety of polymeric resists [V. R. Manfrinato et al., Nano Lett. 14, 4406 (2014); V. R. Manfrinato, A. Stein, L. Zhang, Y. Nam, K. G. Yager, E. A. Stach, and C. T. Black, Nano Lett. 17, 4562 (2017)]. However, this technology, which in this article the authors will refer to as conventional electron beam lithography (CEBL), is reaching its practical resolution and precision limits [V. R. Manfrinato et al., Nano Lett. 14, 4406 (2014)]. Hydrogen depassivation lithography (HDL) [J. N. Randall, J. W. Lyding, S. Schmucker, J. R. Von Ehr, J. Ballard, R. Saini, and Y. Ding, J. Vac. Sci. Technol. B 27, 2764 (2009); J. N. Randall, J. B. Ballard, J. W. Lyding, S. Schmucker, J. R. Von Ehr, R. Saini, H. Xu, and Y. Ding, Microelectron. Eng. 87, 955 (2010)] is a different version of electron beam lithography that is not limited in resolution and precision in the way that CEBL is. It uses a cold field emitter, a scanning tunneling microscope (STM) tip, to deliver a small spot of electrons on a Si (100) 2 × 1 H-passivated surface to expose a self-developing resist that is a monolayer of H adsorbed to the Si surface. Subnanometer features [S. Chen, H. Xu, K. E. J. Goh, L. Liu, and J. N. Randall, Nanotechnology 23, 275301 (2012)], and even the removal of single H atoms can be routinely accomplished [M. A. Walsh and M. C. Hersam, Annu. Rev. Phys. Chem. 60, 193 (2009)]. It is known that the H desorption process at low biases is a multielectron process [E. Foley, A. Kam, J. Lyding, and P. Avouris, Phys. Rev. Lett. 80, 1336 (1998)], but the tunneling distribution of the electrons from the STM tip to the Si surface lattice is not known. The authors have developed a simple model that demonstrates that the combination of two highly nonlinear processes creates a much higher contrast exposure mechanism than CEBL. Currently, HDL has been used almost exclusively on the Si (100) surface and has a limited number of pattern transfer techniques including Si and Ge patterned epitaxy, selective atomic layer deposition of TiO2 followed by reactive ion etching [J. B. Ballard, T. W. Sisson, J. H. G. Owen, W. R. Owen, E. Fuchs, J. Alexander, J. N. Randall, and J. R. Von Ehr, J. Vac. Sci. Technol. B 31, 06FC01 (2013)], and selective deposition of dopant atoms for quantum devices and materials [Workshop on 2D Quantum MetaMaterials held at NIST, Gaithersburg, MD, April 25–26, 2018, edited by J. Owen and W. P. Kirk]. While the throughput of HDL is very low, going parallel in a big way appears promising [J. N. Randall, J. H. G. Owen, J. Lake, R. Saini, E. Fuchs, M. Mahdavi, S. O. R. Moheimani, and B. C. Schaefer, J. Vac. Sci. Technol. B 36, 6 (2018)]. However, the most exciting aspect of HDL is its atomic-scale resolution and precision, which is key to nanoscale research. The authors see HDL emerging as the ultimate high-resolution patterning tool in top down nanofabrication research facilities.

42 ENGINEERING↗

DOE OSTI2015

Effects of forcing time scale on the simulated turbulent flows and turbulent collision statistics of inertial particles

In this paper, we study systematically the effects of forcing time scale in the large-scale stochastic forcing scheme of Eswaran and Pope [“An examination of forcing in direct numerical simulations of turbulence,” Comput. Fluids 16, 257 (1988)] on the simulated flow structures and statistics of forced turbulence. Using direct numerical simulations, we find that the forcing time scale affects the flow dissipation rate and flow Reynolds number. Other flow statistics can be predicted using the altered flow dissipation rate and flow Reynolds number, except when the forcing time scale is made unrealistically large to yield a Taylor microscale flow Reynolds number of 30 and less. We then study the effects of forcing time scale on the kinematic collision statistics of inertial particles. We show that the radial distribution function and the radial relative velocity may depend on the forcing time scale when it becomes comparable to the eddy turnover time. This dependence, however, can be largely explained in terms of altered flow Reynolds number and the changing range of flow length scales present in the turbulent flow. We argue that removing this dependence is important when studying the Reynolds number dependence of the turbulent collision statistics. The results are also compared to those based on a deterministic forcing scheme to better understand the role of large-scale forcing, relative to that of the small-scale turbulence, on turbulent collision of inertial particles. To further elucidate the correlation between the altered flow structures and dynamics of inertial particles, a conditional analysis has been performed, showing that the regions of higher collision rate of inertial particles are well correlated with the regions of lower vorticity. Regions of higher concentration of pairs at contact are found to be highly correlated with the region of high energy dissipation rate.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2013

Development of a fast and feasible spectrum modeling technique for flattening filter free beams

Purpose: To develop a fast and robust technique for the determination of optimized photon spectra for flattening filter free (FFF) beams to be applied in convolution/superposition dose calculations. Methods: A two-step optimization method was developed to derive optimal photon spectra for FFF beams. In the first step, a simple functional form of the photon spectra proposed by Ali ['Functional forms for photon spectra of clinical linacs,' Phys. Med. Biol. 57, 31-50 (2011)] is used to determine generalized shapes of the photon spectra. In this method, the photon spectra were defined for the ranges of field sizes to consider the variations of the contributions of scattered photons with field size. Percent depth doses (PDDs) for each field size were measured and calculated to define a cost function, and a collapsed cone convolution (CCC) algorithm was used to calculate the PDDs. In the second step, the generalized functional form of the photon spectra was fine-tuned in a process whereby the weights of photon fluence became the optimizing free parameters. A line search method was used for the optimization and first order derivatives with respect to the optimizing parameters were derived from the CCC algorithm to enhance the speed of the optimization. The derived photon spectra were evaluated, and the dose distributions using the optimized spectra were validated. Results: The optimal spectra demonstrate small variations with field size for the 6 MV FFF beam and relatively large variations for the 10 MV FFF beam. The mean energies of the optimized 6 MV FFF spectra were decreased from 1.31 MeV for a 3 Multiplication-Sign 3 cm{sup 2} field to 1.21 MeV for a 40 Multiplication-Sign 40 cm{sup 2} field, and from 2.33 MeV at 3 Multiplication-Sign 3 cm{sup 2} to 2.18 MeV at 40 Multiplication-Sign 40 cm{sup 2} for the 10 MV FFF beam. The developed method could significantly improve the agreement between the calculated and measured PDDs. Root mean square differences on the optimized PDDs were observed to be 0.41% (3 Multiplication-Sign 3 cm{sup 2}) down to 0.21% (40 Multiplication-Sign 40 cm{sup 2}) for the 6 MV FFF beam, and 0.35% (3 Multiplication-Sign 3 cm{sup 2}) down to 0.29% (40 Multiplication-Sign 40 cm{sup 2}) for the 10 MV FFF beam. The first order derivatives from the functional form were found to improve the speed of computational time up to 20 times compared to the other techniques. Conclusions: The derived photon spectra resulted in good agreements with measured PDDs over the range of field sizes investigated. The suggested method is easily applicable to commercial radiation treatment planning systems since it only requires measured PDDs as input.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2024

First-Principles Simulation of Beam-Induced Processes Underlying Atomic Manipulation in Electron Microscopes

The development of experimental methods and apparatuses capable of promoting atomically precise material manipulations holds great promise for realizing the ultimate limit of feature miniaturization in materials and devices. The ability to modify materials atom by atom is anticipated to usher in new technologies in areas as diverse as separation science, medicine, and quantum information science. Historically, scanning probe-based techniques have been the most prominent approaches in this space. However, these methods are best suited for the manipulation of surface-exposed regions of materials, as the strong perturbations required for bond scission are delivered most effectively to atoms in the near-proximity to the scanning probe. In contrast, convergent electron beams with energies tuned slightly below the threshold for inducing irreversible knock-on damage have recently been employed (within scanning transmission electron microscopy) to promote atomic-scale bond rearrangements in various beam-stable solids. Currently, however, the efficiency and selectivity of beam-induced atomic manipulation processes with focused electron beams are such that long irradiation times are required to induce a desired atomic rearrangement. With a better understanding of the underlying physics dictating the outcome of a given irradiation event, methods can be devised to improve the efficiency of these techniques so that their promise can be fully realized through widespread adoption.To this end, this Account details our recent efforts to develop and apply tractable first-principles simulation approaches for studying the response of materials to electric beam-like external electric potentials applied in real space. We briefly review the concepts and capabilities in the area of atomically precise materials manipulation and review the early demonstrations of accomplishments in this area, focusing on studies using scanned convergent electron beam probes in particular. We expound upon the depth of the challenge and identify critical shortcomings of theoretical methods that have previously been employed in the simulation of beam-induced processes. We then describe the computational methods that we have generalized from the concepts and tools most commonly applied to the study of molecular photochemistry and how our adaptations of these methods can be employed to capture the relevant dynamical phenomena for beam-induced processes ranging from the initial electron scattering to the ensuing multistate reactions. Here, we contextualize these methods within the current state of the art in this area, which has historically focused primarily on the simulation of inelastic image formation in the electron microscope for the purpose of interpreting the results of quantitative electron microscopy experiments. We demonstrate that the spatial distribution of state-specific excitation rates due to the presence of an external (probe) electric charge is inhomogeneous, such that irradiation at particular locations in materials can favor specific electronic transitions (and disallow others). In addition to the potential for excited-state reaction pathways to be accessed through the initial inelastic scattering of the tightly focused electron beam from the targeted atoms, we also identify favorable conditions for the electronically nonadiabatic evolution of the highly vibrationally excited system to open complex multistate reaction pathways. Implications of the early results for understanding the mechanisms and potential routes to improved efficiency and selectivity in beam-induced reactions are discussed. We conclude with a summary of the current state of theory and modeling capabilities in this area and provide our perspective on future directions for theoretical and experimental developments that we view as crucial to advancing the use of convergent electron beams in mode-specific, atomically precise platforms for direct-write materials modifications.

36 MATERIALS SCIENCE↗

DOE OSTI2016

Topics in Theoretical Physics

This award supported a broadly based research effort in theoretical particle physics, including research aimed at uncovering the laws of nature at short (subatomic) and long (cosmological) distances. These theoretical developments apply to experiments in laboratories such as CERN, the facility that operates the Large Hadron Collider outside Geneva, as well as to cosmological investigations done using telescopes and satellites. The results reported here apply to physics beyond the so-called Standard Model of particle physics; physics of high energy collisions such as those observed at the Large Hadron Collider; theoretical and mathematical tools and frameworks for describing the laws of nature at short distances; cosmology and astrophysics; and analytic and computational methods to solve theories of short distance physics. Some specific research accomplishments include + Theories of the electroweak interactions, the forces that give rise to many forms of radioactive decay; + Physics of the recently discovered Higgs boson. + Models and phenomenology of dark matter, the mysterious component of the universe, that has so far been detected only by its gravitational effects. + High energy particles in astrophysics and cosmology. + Algorithmic research and Computational methods for physics of and beyond the Standard Model. + Theory and applications of relativity and its possible limitations. + Topological effects in field theory and cosmology. + Conformally invariant systems and AdS/CFT. This award also supported significant training of students and postdoctoral fellows to lead the research effort in particle theory for the coming decades. These students and fellows worked closely with other members of the group as well as theoretical and experimental colleagues throughout the physics community. Many of the research projects funded by this grant arose in response to recently obtained experimental results in the areas of particle physics and cosmology. We describe a few of these below. Relativity is founded on a symmetry property of nature called "Lorentz Invariance". Like all symmetry properties, it is essential to determine precisely how symmetric nature actually is; that is, do the laws of nature fully respect the symmetry or is there room for tiny symmetry violating effects? An important consequence of Lorentz invariance is the existence of a universal limiting velocity for all physical particles. Light travels at this limiting velocity so it is frequently referred to as simply "the speed of light", but relativity requires that ALL particles travel more slowly than this speed. Once the Higgs particle was discovered in 2012 a natural question was whether or not this particle's speed was consistent with relativity. Although the speed of the Higgs particle is not measurable directly, Cohen has shown that, if the maximal speed of the Higgs particle was not precisely the same as the speed of light, then the Higgs would have some unusual properties. In some cases the Higgs would be unstable to some unusual decay modes; in other cases the interactions of the Higgs with other particles would change the properties of these other particles in ways that could be observed in so-called cosmic rays, very energetic particles (such as photons, protons and other atomic nuclei) coming from space. Once these particles hit the upper atmosphere they produce a "shower" of particles that can be seen by ground-based instruments. If the Higgs has a maximal speed that differs even a tiny bit from the speed of light these showers would look quite different from what is observed. In this way Cohen was able to establish that the Higgs travels with a maximal speed that cannot differ from the speed of light by more than one part in a thousand-trillion. This is by far the most precisely determined property of the Higgs particle. Cohen and Schmaltz reviewed evidence from the Large Hadron Collider (LHC), a particle physics experiment operating at the CERN laboratory near Geneva, for a new particle sometimes called a W'. This evidence included certain unexpected by-products in collisions of protons at very high energy. While the evidence was not significant enough to claim a discovery, it was sufficiently intriguing that many particle theorists worked to construct explanations for this signal. Cohen and Schmaltz were able to determine that such explanations are highly constrained by previous experiments involving collisions of very energetic particles. Nevertheless they were able to construct a theory that adequately explains the LHC data and remain consistent with prior experiments. Their explanation predicts the existence of yet another new particle, called a Z', with a mass slightly greater than that of the W'. This additional particle, if it exists, should be seen as more data is collected from the LHC. Amusingly, there is one collision by-product that has already been seen by the CMS experiment at the LHC that supports the existence of this new particle; however, it is not unlikely that this single event is a so-called "background" event, that is a somewhat atypical by-product of a conventional Standard Model process. This theory for the anomalous LHC data will either be confirmed or excluded with further data-taking at the LHC. The ratio of the number of electrons produced in bottom quark decays over the number of muons produced has been measured at the LHC. This ratio is interesting because it can be predicted very precisely from a basic property of the Standard Model: lepton universality. If lepton universality is correct, the ratio of electrons to muons is predicted to be equal to 1. The first measurements of this ratio find a value different from 1 with a statistical significance of about 3 standard deviations. Schmaltz and collaborator proposed a new extension of the Standard Model which can explain the new data. In addition, Schmaltz and collaborators proposed several new measurements of ratios of decay rates which can confirm or rule out the surprising results from the earlier LHC data. The most recent and precise measurements of the cosmic microwave background from the Planck satellite, from a combination of measurements of the dark matter distribution in the universe, and from a measurement of the expansion rate of the universe today show some disagreement when interpreted in terms of the so-called LambdaCDM model. Schmaltz and collaborators proposed an alternative model to LambdaCDM in which the usual cold dark matter is replaced by a new ``dark sector". This sector consists of a cold dark matter particle which interacts with a newly postulated dark radiation component of the universe. The dark radiation can help explain the discrepancy in measurements of the expansion rate, and the dark matter interactions subtly modify the clumping of dark matter at large scales, thus potentially explaining both kinds of tensions in the data. In two publications Schmaltz described the new model and then performed a precision comparison of the predictions of the model with all currently available cosmological data. The results favor the new model at the level of three standard deviations with current data. Quantum Field Theory (QFT) is the language we use to describe quantum systems which are consistent with Einstein’s theory of Special Relativity. In particular, the requirement of Einstein’s theory that signals not travel faster than the speed of light constrains the types of interactions which particles can engage in. One consequence of relativity is that these interactions cannot preserve particle number. The stronger the interactions, the more severe the particle number violation in a given Relativistic QFT. When particle number violation is strong, it becomes very difficult to adequately parameterize the quantum wave function (which characterizes the state of a quantum system). For example, though we can formulate the QFT which describes the strong force as a set of interactions between quark and gluon particles, we have no clear idea how to express the proton state in terms of these quarks and gluons. This is because the proton, though a bound state of quarks and gluons, is not a state of a fixed number of particles due to strong interactions. Yet, understanding the proton state is very important in order to theoretically predict the reaction rates observed at the LHC in Geneva, which is a proton-proton collider. Katz has formulated a new approach to QFT, which among other things offers a way to adequately approximate the quantum wave function of a bound state at strong coupling. The approximation scheme is related to the fact that any sensible QFT (including that of the strong interactions) is at short distances approximately self-similar upon rescaling of space and time. It turns out that keeping track of the response upon this rescaling is important in efficiently parameterizing the state. Katz and collaborators have used this observation to approximate the state of the proton in toy versions of the strong force. In the late 60s Sheldon Glashow, Abdus Salam and Steven Weinberg (1979 Nobel Prize awardees) proposed a theory unifying weak and electromagnetic interaction which assumed the existence of new particles, the W and Z bosons. The W and Z bosons were eventually detected in high-energy collision in a particle accelerator at CERN, and the recent discovery of the Higgs meson at the Large Hadron Collider (LHC), always at CERN, completed the picture. However, deep theoretical considerations indicate that the theory by Glashow, Weinberg and Salam, often referred to as "the standard model" cannot be the whole story: the existence of new particles and new interactions at yet higher energies is widely anticipated. The experiments at the LHC are looking for these, while theorists, like Brower, Rebbi and collaborators, are investigating models for these new interactions. Working in a large national collaboration with access to the most powerful DOE computers Brower, Rebbi and colleagues have been using calculational techniques, similar to those successfully employed for many years to investigate the interactions among quarks in nucleons, to study theories that can describe the expected "beyond the standard model" (BSM) interactions. Their results, which include also a model for dark matter, have been published in several refereed papers in prestigious journals. Various ideas in topologically interesting field theories predict hypothetical objects such as fractional charges and Majorana excitations. However, such fascinating objects have not been seen in particle physics. Nevertheless, these objects demonstrate possible phenomena that quantum field theory can support. Pi used condensed matter physics as a laboratory to study possible realizations and observable effects of these objects predicted by quantum field theory. In recent times there has developed considerable interest among condensed matter field theorists in precisely the same geometrical and topological structures, which were first discovered in particle physics field theories. From particle physicists' point of view, this is an interesting development, since condensed matter provides an arena in which one can concretely realize particle physics ideas. Moreover, particle physicists can learn new ideas from condensed matter physics. Higgs phenomenon is precisely an important particle physics realization of condensed matter ideas. In contrast to the small distance characterizing condensed matter systems, field theory also describes large distance physics characterizing cosmology. Pi worked on various geometrical effects in the standard theory of cosmology, viz general relativity.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2019

Scrambling via braiding of nonabelions

We study how quantum states are scrambled via braiding in systems of non-Abelian anyons through the lens of entanglement spectrum statistics. In particular, we focus on the degree of scrambling, defined as the randomness produced by braiding, at the same amount of entanglement entropy. To quantify the degree of randomness, we define a distance between the entanglement spectrum level spacing distribution of a state evolved under random braids and that of a Haar-random state, using the Kullback-Leibler divergence D KL . We study D KL numerically for random braids of Majorana fermions (supplemented with random local four-body interactions) and Fibonacci anyons. For comparison, we also obtain D KL for the Sachdev-Ye-Kitaev model of Majorana fermions with all-to-all interactions, random unitary circuits built out of (a) Hadamard (H), π/8 (T), and CNOT gates, and (b) random unitary circuits built out of two-qubit Haar random unitaries. To compare the degree of randomness that different systems produce beyond entanglement entropy, we look at D KL as a function of the Page limit-normalized entanglement entropy S/S max . Furthermore, our results reveal a hierarchy of scrambling among various models—even for the same amount of entanglement entropy—at intermediate times, whereas all models exhibit the same late-time behavior. In particular, we find that braiding of Fibonacci anyons randomizes initial product states more efficiently than the universal H+T+CNOT set.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2014

Numerical analysis of flows of rarefied gases in long channels with octagonal cross section shapes

Isothermal, pressure driven rarefied gas flows through long channels with octagonal cross section shapes are analyzed computationally. The capillary is between inlet and outlet reservoirs. The cross section is constant along the axial direction. The boundary condition at the solid-gas interface is assumed to be diffuse reflection. Since the channel is long, the gaseous velocity is small compared to the average molecular speed. Consequently, a linearized description can be used. The flow is described by the linearized Bhatnagar-Gross-Krook kinetic model. The solution of the problem is divided into two stages. First, the local flow field is determined by assuming the local pressure gradient. Secondly, the global flow behavior is deduced by the consideration of the conservation of the mass along the axis of the capillary. The kinetic equation is solved by the discrete velocity method on the cross section. Both spatial and velocity spaces are discretized. A body fitted rectangular grid is used for the spatial space. Near the boundary, first-order, while in the interior part of the flow domain, second-order finite-differences are applied to approximate the spatial derivatives. This combination results into an efficient and straightforward numerical treatment. The velocity space is represented by a Gauss-Legendre quadrature. The kinetic equation is solved in an iterative manner. The local dimensionless flow rate is calculated and tabulated for a wide range of the gaseous rarefaction for octagonal cross sections with various geometrical parameters. It exhibits the Knudsen minimum phenomenon. The flow rates in the octagonal channel are compared to those through capillaries with circular and square cross sections. Typical velocity profiles are also shown. The mass flow rate and the distribution of the pressure are determined and presented for global pressure driven flows.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2014

Stochastic generation of explicit pore structures by thresholding Gaussian random fields

We provide a description and computational investigation of an efficient method to stochastically generate realistic pore structures. Smolarkiewicz and Winter introduced this specific method in pores resolving simulation of Darcy flows (Smolarkiewicz and Winter, 2010 [1]) without giving a complete formal description or analysis of the method, or indicating how to control the parameterization of the ensemble. We address both issues in this paper. The method consists of two steps. First, a realization of a correlated Gaussian field, or topography, is produced by convolving a prescribed kernel with an initial field of independent, identically distributed random variables. The intrinsic length scales of the kernel determine the correlation structure of the topography. Next, a sample pore space is generated by applying a level threshold to the Gaussian field realization: points are assigned to the void phase or the solid phase depending on whether the topography over them is above or below the threshold. Hence, the topology and geometry of the pore space depend on the form of the kernel and the level threshold. Manipulating these two user prescribed quantities allows good control of pore space observables, in particular the Minkowski functionals. Extensions of the method to generate media with multiple pore structures and preferential flow directions are also discussed. To demonstrate its usefulness, the method is used to generate a pore space with physical and hydrological properties similar to a sample of Berea sandstone. -- Graphical abstract: -- Highlights: •An efficient method to stochastically generate realistic pore structures is provided. •Samples are generated by applying a level threshold to a Gaussian field realization. •Two user prescribed quantities determine the topology and geometry of the pore space. •Multiple pore structures and preferential flow directions can be produced. •A pore space based on Berea sandstone is generated.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2013

Spatial frequency spectrum of the x-ray scatter distribution in CBCT projections

Purpose: X-ray scatter is a source of significant image quality loss in cone-beam computed tomography (CBCT). The use of Monte Carlo (MC) simulations separating primary and scattered photons has allowed the structure and nature of the scatter distribution in CBCT to become better elucidated. This work seeks to quantify the structure and determine a suitable basis function for the scatter distribution by examining its spectral components using Fourier analysis.Methods: The scatter distribution projection data were simulated using a CBCT MC model based on the EGSnrc code. CBCT projection data, with separated primary and scatter signal, were generated for a 30.6 cm diameter water cylinder [single angle projection with varying axis-to-detector distance (ADD) and bowtie filters] and two anthropomorphic phantoms (head and pelvis, 360 projections sampled every 1°, with and without a compensator). The Fourier transform of the resulting scatter distributions was computed and analyzed both qualitatively and quantitatively. A novel metric called the scatter frequency width (SFW) is introduced to determine the scatter distribution's frequency content. The frequency content results are used to determine a set basis functions, consisting of low-frequency sine and cosine functions, to fit and denoise the scatter distribution generated from MC simulations using a reduced number of photons and projections. The signal recovery is implemented using Fourier filtering (low-pass Butterworth filter) and interpolation. Estimates of the scatter distribution are used to correct and reconstruct simulated projections.Results: The spatial and angular frequencies are contained within a maximum frequency of 0.1 cm{sup −1} and 7/(2π) rad{sup −1} for the imaging scenarios examined, with these values varying depending on the object and imaging setup (e.g., ADD and compensator). These data indicate spatial and angular sampling every 5 cm and π/7 rad (∼25°) can be used to properly capture the scatter distribution, with reduced sampling possible depending on the imaging scenario. Using a low-pass Butterworth filter, tuned with the SFW values, to denoise the scatter projection data generated from MC simulations using 10{sup 6} photons resulted in an error reduction of greater than 85% for the estimating scatter in single and multiple projections. Analysis showed that the use of a compensator helped reduce the error in estimating the scatter distribution from limited photon simulations by more than 37% when compared to the case without a compensator for the head and pelvis phantoms. Reconstructions of simulated head phantom projections corrected by the filtered and interpolated scatter estimates showed improvements in overall image quality.Conclusions: The spatial frequency content of the scatter distribution in CBCT is found to be contained within the low frequency domain. The frequency content is modulated both by object and imaging parameters (ADD and compensator). The low-frequency nature of the scatter distribution allows for a limited set of sine and cosine basis functions to be used to accurately represent the scatter signal in the presence of noise and reduced data sampling decreasing MC based scatter estimation time. Compensator induced modulation of the scatter distribution reduces the frequency content and improves the fitting results.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

DOE OSTI2018

Self-organized Segregation on the Grid

We consider an agent-based model with exponentially distributed waiting times in which two types of agents interact locally over a graph, and based on this interaction and on the value of a common intolerance threshold τ, decide whether to change their types. This is equivalent to a zero-temperature ising model with Glauber dynamics, an asynchronous cellular automaton with extended Moore neighborhoods, or a Schelling model of self-organized segregation in an open system, and has applications in the analysis of social and biological networks, and spin glasses systems. Some rigorous results were recently obtained in the theoretical computer science literature, and this work provides several extensions. We enlarge the intolerance interval leading to the expected formation of large segregated regions of agents of a single type from the known size ϵ>0 to size ≈0.134. Namely, we show that for 0.433<τ<1/2 (and by symmetry 1/2<τ<0.567), the expected size of the largest segregated region containing an arbitrary agent is exponential in the size of the neighborhood. We further extend the interval leading to expected large segregated regions to size ≈0.312 considering “almost segregated” regions, namely regions where the ratio of the number of agents of one type and the number of agents of the other type vanishes quickly as the size of the neighborhood grows. In this case, we show that for 0.344<τ≤0.433 (and by symmetry for 0.567≤τ<0.656) the expected size of the largest almost segregated region containing an arbitrary agent is exponential in the size of the neighborhood. This behavior is reminiscent of supercritical percolation, where small clusters of empty sites can be observed within any sufficiently large region of the occupied percolation cluster. The exponential bounds that we provide also imply that complete segregation, where agents of a single type cover the whole grid, does not occur with high probability for p=1/2 and the range of intolerance considered.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2022

Learning functional priors and posteriors from data and physics

In this work, we develop a new Bayesian framework based on deep neural networks to be able to extrapolate in space-time using historical data and to quantify uncertainties arising from both noisy and gappy data in physical problems. Specifically, the proposed approach has two stages: (1) prior learning and (2) posterior estimation. At the first stage, we employ the physics-informed Generative Adversarial Networks (PI-GAN) to learn a functional prior either from a prescribed function distribution, e.g., Gaussian process, or from historical data and physics. At the second stage, we employ the Hamiltonian Monte Carlo (HMC) method to estimate the posterior in the latent space of PI-GANs. In addition, we use two different approaches to encode the physics: (1) automatic differentiation, used in the physicsinformed neural networks (PINNs) for scenarios with explicitly known partial differential equations (PDEs), and (2) operator regression using the deep operator network (DeepONet) for PDE-agnostic scenarios. We then test the proposed method for (1) meta-learning for one-dimensional regression, and forward/inverse PDE problems (combined with PINNs); (2) PDE-agnostic physical problems (combined with DeepONet), e.g., fractional diffusion as well as saturated stochastic (100-dimensional) flows in heterogeneous porous media; and (3) spatial-temporal regression problems, i.e., inference of a marine riser displacement field using experimental data from the Norwegian Deepwater Programme (NDP). The results demonstrate that the proposed approach can provide accurate predictions as well as uncertainty quantification given very limited scattered and noisy data, since historical data could be available to provide informative priors. In summary, the proposed method is capable of learning flexible functional priors, e.g., both Gaussian and non-Gaussian process, and can be readily extended to big data problems by enabling mini-batch training using stochastic HMC or normalizing flows since the latent space is generally characterized as low dimensional.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2013

Optical photon transport in powdered-phosphor scintillators. Part II. Calculation of single-scattering transport parameters

Purpose: Monte Carlo methods based on the Boltzmann transport equation (BTE) have previously been used to model light transport in powdered-phosphor scintillator screens. Physically motivated guesses or, alternatively, the complexities of Mie theory have been used by some authors to provide the necessary inputs of transport parameters. The purpose of Part II of this work is to: (i) validate predictions of modulation transform function (MTF) using the BTE and calculated values of transport parameters, against experimental data published for two Gd{sub 2}O{sub 2}S:Tb screens; (ii) investigate the impact of size-distribution and emission spectrum on Mie predictions of transport parameters; (iii) suggest simpler and novel geometrical optics-based models for these parameters and compare to the predictions of Mie theory. A computer code package called phsphr is made available that allows the MTF predictions for the screens modeled to be reproduced and novel screens to be simulated. Methods: The transport parameters of interest are the scattering efficiency (Q{sub sct}), absorption efficiency (Q{sub abs}), and the scatter anisotropy (g). Calculations of these parameters are made using the analytic method of Mie theory, for spherical grains of radii 0.1-5.0 {mu}m. The sensitivity of the transport parameters to emission wavelength is investigated using an emission spectrum representative of that of Gd{sub 2}O{sub 2}S:Tb. The impact of a grain-size distribution in the screen on the parameters is investigated using a Gaussian size-distribution ({sigma}= 1%, 5%, or 10% of mean radius). Two simple and novel alternative models to Mie theory are suggested: a geometrical optics and diffraction model (GODM) and an extension of this (GODM+). Comparisons to measured MTF are made for two commercial screens: Lanex Fast Back and Lanex Fast Front (Eastman Kodak Company, Inc.). Results: The Mie theory predictions of transport parameters were shown to be highly sensitive to both grain size and emission wavelength. For a phosphor screen structure with a distribution in grain sizes and a spectrum of emission, only the average trend of Mie theory is likely to be important. This average behavior is well predicted by the more sophisticated of the geometrical optics models (GODM+) and in approximate agreement for the simplest (GODM). The root-mean-square differences obtained between predicted MTF and experimental measurements, using all three models (GODM, GODM+, Mie), were within 0.03 for both Lanex screens in all cases. This is excellent agreement in view of the uncertainties in screen composition and optical properties. Conclusions: If Mie theory is used for calculating transport parameters for light scattering and absorption in powdered-phosphor screens, care should be taken to average out the fine-structure in the parameter predictions. However, for visible emission wavelengths ({lambda} < 1.0 {mu}m) and grain radii (a > 0.5 {mu}m), geometrical optics models for transport parameters are an alternative to Mie theory. These geometrical optics models are simpler and lead to no substantial loss in accuracy.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗