Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “sparse”

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.

At least 937 records · Page 52

Universal topological quantum computation with strongly correlated Majorana edge modes

Abstract Majorana-based quantum gates are not complete for performing universal topological quantum computation while Fibonacci-based gates are difficult to be realized electronically and hardly coincide with the conventional quantum circuit models. In reference Hu and Kane (2018 Phys. Rev. Lett. 120 066801), it has been shown that a strongly correlated Majorana edge mode in a chiral topological superconductor can be decomposed into a Fibonacci anyon τ and a thermal operator anyon ɛ in the tricritical Ising model. The deconfinement of τ and ɛ via the interaction between the fermion modes yields the anyon collisions and gives the braiding of either τ or ɛ . With these braidings, the complete members of a set of universal gates, the Pauli gates, the Hadamard gate and extra phase gates for one-qubit as well as controlled-NOT (CNOT) gate for two-qubits, are topologically assembled. Encoding quantum information and reading out the computation results can be carried out through electric signals. With the sparse-dense mixed encodings, we set up the quantum circuit where the CNOT gate turns out to be a probabilistic gate and design the corresponding devices with thin films of the chiral topological superconductor. As an example of the universal topological quantum computing, we show the application to Shor’s integer factorization algorithm.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Shape-driven, emergent behavior in active particle mixtures

Abstract Active particle systems can vary greatly from one-component systems of spheres to mixtures of particle shapes at different composition ratios. We investigate computationally the combined effect of anisotropy and stoichiometry on the collective behavior of two-dimensional active colloidal mixtures of polygons. We uncover three emergent phenomena not yet reported in active Brownian particle systems. First, we find that mixtures containing hexagons exhibit micro-phase separation with large grains of hexagonal symmetry. We quantify a measurable, implicit ‘steric attraction’ between the active particles as a result of shape anisotropy and activity. This calculation provides further evidence that implicit interactions in active systems, even without explicit attraction, can lead to an effective preferential attraction between particles. Next, we report stable fluid clusters in mixtures containing one triangle or square component. We attribute the fluidization of the dense cluster to the interplay of cluster destabilizing particles, which introduce grain boundaries and slip planes into the system, causing solid-like clusters to break up into fluid clusters. Third, we show that fluid clusters can coexist with solid clusters within a sparse gas of particles in a steady state of three coexisting phases. Our results highlight the potential for a wide variety of behavior to be accessible to active matter systems and establish a route to control active colloidal systems through simple parameter designs.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Diffusion energy barrier of Au on Bi 2 Se 3 : theory and experiment

The stability and diffusion of ultra-thin thermally deposited Au films on Bi 2 Se 3 was studied using scanning tunneling microscopy and density functional theory calculations. The Au/Bi 2 Se 3 interface is of interest as gold is predicted to provide excellent electrical contact while maintaining the spin-polarized characteristics of the electronic states in Bi 2 Se 3 that make the material attractive for spintronic applications. When deposited at room temperature, Au 10 covers the surface with tightly packed islands of nanometer scale dimension. The surface morphology is stable up to 400K. At this annealing temperature, Au atoms have sufficient energy to diffuse across the surface and aggregate into larger nanostructures. At 550K, the Bi 2 Se 3 surface is only sparsely covered, and the Au has formed clusters with length scales 5-10 times larger than the original islands formed at room temperature. Comparison of the experiment and first principle calculation lead to the conclusion that the diffusion energy barrier for Au on Bi 2 Se 3 is as high as 0.47 eV, 15 which is much larger than diffusion barriers on other van der Waals materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Analyzing the cosmic web environment in the vicinity of grand-design and flocculent spirals with local geometric index

We explore the environment of a combined set of 367 grand-design and 619 flocculent spiral galaxies. We introduce a novel estimator called the local geometric index to quantify the morphology of the local environment of these 986 spirals. The local geometric index allows us to classify the environment of galaxies into voids, sheets, filaments, and clusters. We find that grand-designs are mostly located in dense environments like clusters and filaments (∼ 78%), whereas the fraction of the flocculents lying in sparse environments like voids and sheets is significantly higher (> 10%) than that of the grand-designs. A p-value < 10 -3 from a Kolmogorov-Smirnov test indicates that our results are statistically significant at 99.9% confidence level. Further, we note that dense environments with large tidal flows are dominated by the grand-designs. On the other hand, low-density environments such as sheets and voids favor the growth of flocculents.

79 ASTRONOMY AND ASTROPHYSICS↗

Measurement of the small-scale 3D Lyman- α forest power spectrum

Small-scale correlations measured in the Lyman-α (Lyα) forest encode information about the intergalactic medium and the primordial matter power spectrum. In this article, we present and implement a simple method to measure the 3-dimensional power spectrum, P 3D , of the Lyα forest at wavenumbers k corresponding to small, ~ Mpc scales. In order to estimate P 3D from sparsely and unevenly distributed data samples, we rely on averaging 1-dimensional Fourier Transforms, as previously carried out to estimate the 1-dimensional power spectrum of the Lyα forest, P 1D . Further, this methodology exhibits a very low computational cost. We confirm the validity of this approach through its application to Nyx cosmological hydrodynamical simulations. Subsequently, we apply our method to the eBOSS DR16 Lyα forest sample, providing as a proof of principle, a first P 3D measurement averaged over two redshift bins z = 2.2 and z = 2.4. This work highlights the potential for forthcoming P 3D measurements, from upcoming large spectroscopic surveys, to untangle degeneracies in the cosmological interpretation of P 1D .

79 ASTRONOMY AND ASTROPHYSICS↗

Trigonometric Parallaxes of Two T Dwarfs With Keck and ShaneAO Astrometry

We present trigonometric parallax and proper motion measurements for two T-type brown dwarfs. We derive our measurements from infrared laser guide star adaptive optics observations spanning five years from the ShaneAO/SHARCS and NIRC2/medium-cam instruments on the Shane and Keck telescopes, respectively. To improve our astrometric precision, we measure and apply a distortion correction to our fields for both instruments. We also transform the Keck and ShaneAO astrometric reference frames onto the ICRS using five-parameter parallax and proper motion solutions for background reference stars from Gaia DR2. Fitting for parallax and proper motion, we measure parallaxes of 73.5 ± 9.2 mas and 70.1 ± 6.7 mas for WISEJ19010703+47181688 (WISE1901) and WISEJ21543294+59421370 (WISE2154), respectively. We utilize Monte Carlo methods to estimate the error in our sparse field methods, taking into account overfitting and differential atmospheric refraction. Comparing to previous measurements in the literature, all of our parallax and proper motion values fall within 2σ of the published measurements, and 4 of 6 measurements are within 1σ. These data are among the first parallax measurements of these T dwarfs and serve as precise measurements for calibrating stellar formation models. These two objects are the first results of an ongoing survey of T dwarfs with Keck/NIRC2 and the Shane Adaptive Optics system at Lick Observatory.

79 ASTRONOMY AND ASTROPHYSICS↗

Time-dependent SOLPS-ITER simulations of the tokamak plasma boundary for model predictive control using SINDy *

Abstract Time-dependent SOLPS-ITER simulations have been used to identify reduced models with the sparse identification of nonlinear dynamics (SINDy) method and develop model-predictive control of the boundary plasma state using main ion gas puff actuation. A series of gas actuation sequences are input into SOLPS-ITER to produce a dynamic response in upstream and divertor plasma quantities. The SINDy method is applied to identify reduced linear and nonlinear models for the electron density at the outboard midplane n e , s e p O M P and the electron temperature at the outer divertor T e , s e p d i v . Note that T e , s e p d i v is not necessarily the peak value of T e along the divertor. The identified reduced models are interpretable by construction (i.e. not black box), and have the form of coupled ordinary differential equations. Despite significant noise in T e , s e p d i v , the reduced models can be used to predict the response over a range of actuation levels to a maximum deviation of 0.5% in n e , s e p O M P and 5%–10% in T e , s e p d i v for the cases considered. Model retraining using time history data triggered by a preset error threshold is also demonstrated. A model predictive control strategy for nonlinear models is developed and used to perform feedback control of a SOLPS-ITER simulation to produce a setpoint trajectory in n e , s e p O M P using the integrated plasma simulator framework. The developed techniques are general and can be applied to time-dependent data from other boundary simulations or experimental data. Ongoing work is extending the approach to model identification and control for divertor detachment, which will present transient nonlinear behavior from impurity seeding, including realistic latency and synthetic diagnostic signals derived from the full SOLPS-ITER output.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A framework for discrete optimization of stellarator coils

Designing magnets for three-dimensional plasma confinement is a key task for advancing the stellarator as a fusion reactor concept. Stellarator magnets must produce an accurate field while leaving adequate room for other components and being reasonably simple to construct and assemble. In this paper, a framework for coil design and optimization is introduced that enables the attainment of sparse magnet solutions with arbitrary restrictions on where coils may be located. The solution space is formulated as a 'wireframe' consisting of a mesh of interconnected wire segments enclosing the plasma. Two methods are developed for optimizing the current distribution on a wireframe: Regularized Constrained Least Squares, which uses a linear least-squares approach to optimize the currents in each segment, and Greedy Stellarator Coil Optimization, a fully discrete procedure in which loops of current are added to the mesh one by one to achieve the desired magnetic field on the plasma boundary. Examples are presented of solutions obtainable with each method, some of which achieve high field accuracy while obeying spatial constraints that permit easy assembly.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Evidence of Kramer extrapolation inaccuracy for predicting high field Nb 3 Sn properties

Future applications requiring high magnetic fields, such as the proposed Future Circular Collider, demand a substantially higher critical current density, J c , at fields ≥16 T than is presently available in any commercial strand, so there is a strong effort to develop new routes to higher J c Nb 3 Sn. As a consequence, evaluating the irreversibility field (H irr ) of any new conductor to ensure reliable performance at these higher magnetic fields becomes essential. To predict the irreversibility field for Nb 3 Sn wires, critical current measurements, I c , are commonly performed in the 12-15 T range and the Kramer extrapolation is used to predict higher field properties. Here, the Kramer extrapolation typically models the contribution only for sparse grain boundary pinning, yet Nb 3 Sn wires rely on a high density of grain boundaries to provide the flux pinning that enables their high critical current density. However, whole-field range VSM measurements up to 30 T recently showed for Nb 3 Sn RRP® wires that the field dependence of the pinning force curve significantly deviates from the typical grain boundary shape, leading to a 1-2 T overestimation of Hirr when extrapolated from the typical mid-field data taken only up to about 15 T. In this work we characterized a variety of both RRP® and PIT Nb 3 Sn wires by transport measurements up to 29 T at the Laboratoire National des Champs Magnétiques Intenses (LNCMI), part of the European Magnetic Field Laboratory in Grenoble, to verify whether or not such overestimation is related to the measurement technique and whether or not it is a common feature across different designs. Indeed we also found that when measured in transport the 12-15 T Kramer extrapolation overestimates the actual H irr in both types of conductor with an inaccuracy of up to 1.6 T, confirming that high field characterization is a necessary tool to evaluate the actual high field performance of each Nb 3 Sn wire.

36 MATERIALS SCIENCE↗

The SparkPix-S ASIC for the sparsified readout of 1 MHz frame-rate X-ray cameras at LCLS-II: pixel design and simulation results

Exploiting the “sparse” nature of the information in XPCS (X-ray Photon Correlation Spectroscopy) and XSVS (Speckle Visibility Spectroscopy) experiments, we present the SparkPix-S, a 3-sides buttable Application Specific Integrated Circuit (ASIC) based on a sparsified readout strategy for large-format hybrid detectors. The SparkPix-S architecture, based on the successful ePix family, will be composed as follows: a front-end 2-D matrix of 384×352 square pixels with 50 µm pitch is arranged to match the dimensions of a PIN Si-sensor matrix; charge readout, signal shaping and amplitude discrimination is performed at pixel-level, by means of a low-power (<18 µW) analog processor, which, in case of an event, negotiates access to an analog bus placed every other column; on the chip periphery (balcony), the information on each bus is digitized by an array of successive approximation analog-to-digital converters (SAR-ADCs) running at 10 Msps; on the digital back-end the global logic will generate the output data stream using low-voltage differential signalling (LVDS). A first prototype of the SparkPix-S, with a reduced matrix size of 96×96 pixels, is currently under production on a 130 nm CMOS technology. Simulated performance results show an equivalent noise charge <60 el. r.m.s. at 1 MHz repetition rate, with a maximum input energy of 60 keV and capability to discriminate charge signals with equivalent energy as low as 900 eV.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Comparison of point cloud and image-based models for calorimeter fast simulation

Score based generative models are a new class of generative models that have been shown to accurately generate high dimensional calorimeter datasets. Recent advances in generative models have used images with 3D voxels to represent and model complex calorimeter showers. Point clouds, however, are likely a more natural representation of calorimeter showers, particularly in calorimeters with high granularity. Furthermore, point clouds preserve all of the information of the original simulation, more naturally deal with sparse datasets, and can be implemented with more compact models and data files. In this work, two state-of-the-art score based models are trained on the same set of calorimeter simulation and directly compared.

47 OTHER INSTRUMENTATION↗

Measurement of the mean excitation energy of liquid argon

The mean excitation energy (I-value) of liquid argon is a critical input for energy estimation in neutrino oscillation experiments. It is measured to be 199 +4 -5 eV using the range of 402.2 MeV protons from the Fermilab Linac. This compares to the author's recent evaluation of (197 ± 7) eV based on a combination of an oscillator strength distribution analysis, gaseous argon range measurements, sparse stopping power data on solid argon, and an extrapolation of data on the effect of phase from other substances. Using all sources of information, we recommend a value of (198 ± 4) eV for liquid argon, which is significantly higher than 188 eV, from ICRU-37's gaseous argon evaluation, commonly used in Monte Carlo codes such as Geant4.

47 OTHER INSTRUMENTATION↗

Dynamic global vegetation models underestimate net CO 2 flux mean and inter-annual variability in dryland ecosystems

Despite their sparse vegetation, dryland regions exert a huge influence over global biogeochemical cycles because they cover more than 40% of the world surface (Schimel 2010 Science 327 418–9). It is thought that drylands dominate the inter-annual variability (IAV) and long-term trend in the global carbon (C) cycle (Poulter et al 2014 Nature 509 600–3, Ahlstrom et al 2015 Science 348 895–9, Zhang et al 2018 Glob. Change Biol. 24 3954–68). Projections of the global land C sink therefore rely on accurate representation of dryland C cycle processes; however, the dynamic global vegetation models (DGVMs) used in future projections have rarely been evaluated against dryland C flux data. Here, we carried out an evaluation of 14 DGVMs (TRENDY v7) against net ecosystem exchange (NEE) data from 12 dryland flux sites in the southwestern US encompassing a range of ecosystem types (forests, shrub- and grasslands). We find that all the models underestimate both mean annual C uptake/release as well as the magnitude of NEE IAV, suggesting that improvements in representing dryland regions may improve global C cycle projections. Across all models, the sensitivity and timing of ecosystem C uptake to plant available moisture was at fault. Spring biases in gross primary production (GPP) dominate the underestimate of mean annual NEE, whereas models' lack of GPP response to water availability in both spring and summer monsoon are responsible for inability to capture NEE IAV. Errors in GPP moisture sensitivity at high elevation forested sites were more prominent during the spring, while errors at the low elevation shrub and grass-dominated sites were more important during the monsoon. We propose a range of hypotheses for why model GPP does not respond sufficiently to changing water availability that can serve as a guide for future dryland DGVM developments. Our analysis suggests that improvements in modeling C cycle processes across more than a quarter of the Earth's land surface could be achieved by addressing the moisture sensitivity of dryland C uptake.

54 ENVIRONMENTAL SCIENCES↗

Increasing freshwater and dissolved organic carbon flows to Northwest Alaska’s Elson lagoon

Manifestations of climate change in the Arctic are numerous and include hydrological cycle intensification and permafrost thaw, both expected as a result of atmospheric and surface warming. Across the terrestrial Arctic dissolved organic carbon (DOC) entrained in arctic rivers may be providing a carbon subsidy to coastal food webs. Yet, data from field sampling is too often of limited duration to confidently ascertain impacts of climate change on freshwater and DOC flows to coastal waters. This study applies numerical modeling to investigate trends in freshwater and DOC exports from land to Elson Lagoon in Northwest Alaska over the period 1981-2020. While the modeling approach has limitations, the results point to significant increases in freshwater and DOC exports to the lagoon over the past four decades. The model simulation reveals significant increases in surface, subsurface (suprapermafrost), and total freshwater exports. Significant increases are also noted in surface and subsurface DOC production and export, influenced by warming soils and associated active-layer thickening. The largest changes in subsurface components are noted in September, which has experienced a ~50% increase in DOC export emanating from suprapermafrost flow. Direct coastal suprapermafrost freshwater and DOC exports in late summer more than doubled between the first and last five years of the simulation period, with a large anomaly in September 2019 representing a more than fourfold increase over September direct coastal export during the early 1980s. These trends highlight the need for dedicated measurement programs that will enable improved understanding of climate change impacts on coastal zone processes in this data sparse region of Northwest Alaska.

54 ENVIRONMENTAL SCIENCES↗

Overlooked cooling effects of albedo in terrestrial ecosystems

Radiative forcing (RF) resulting from changes in surface albedo is increasingly recognized as a significant driver of global climate change but has not been adequately estimated, including by Intergovernmental Panel on Climate Change (IPCC) assessment reports, compared with other warming agents. Here, we first present the physical foundation for modeling albedo-induced RF and the consequent global warming impact (GWI Δα ). We then highlight the shortcomings of available current databases and methodologies for calculating GWI Δα at multiple temporal scales. There is a clear lack of comprehensive in situ measurements of albedo due to sparse geographic coverage of ground-based stations, whereas estimates from satellites suffer from biases due to the limited frequency of image collection, and estimates from earth system models (ESMs) suffer from very coarse spatial resolution land cover maps and associated albedo values in pre-determined lookup tables. Field measurements of albedo show large differences by ecosystem type and large diurnal and seasonal changes. As indicated from our findings in southwest Michigan, GWI Δα is substantial, exceeding the RF Δα values of IPCC reports. Inclusion of GWI Δα to landowners and carbon credit markets for specific management practices are needed in future policies. We further identify four pressing research priorities: developing a comprehensive albedo database, pinpointing accurate reference sites within managed landscapes, refining algorithms for remote sensing of albedo by integrating geostationary and other orbital satellites, and integrating the GWI Δα component into future ESMs.

54 ENVIRONMENTAL SCIENCES↗

Adaptive pruning-based optimization of parameterized quantum circuits

Abstract Variational hybrid quantum–classical algorithms are powerful tools to maximize the use of noisy intermediate-scale quantum devices. While past studies have developed powerful and expressive ansatze, their near-term applications have been limited by the difficulty of optimizing in the vast parameter space. In this work, we propose a heuristic optimization strategy for such ansatze used in variational quantum algorithms, which we call ‘parameter-efficient circuit training (PECT)’. Instead of optimizing all of the ansatz parameters at once, PECT launches a sequence of variational algorithms, in which each iteration of the algorithm activates and optimizes a subset of the total parameter set. To update the parameter subset between iterations, we adapt the Dynamic Sparse Reparameterization scheme which was originally proposed for training deep convolutional neural networks. We demonstrate PECT for the Variational Quantum Eigensolver, in which we benchmark unitary coupled-cluster ansatze including UCCSD and k -UpCCGSD, as well as the Low-Depth Circuit Ansatz (LDCA), to estimate ground state energies of molecular systems. We additionally use a layerwise variant of PECT to optimize a hardware-efficient circuit for the Sycamore processor to estimate the ground state energy densities of the one-dimensional Fermi-Hubbard model. From our numerical data, we find that PECT can enable optimizations of certain ansatze that were previously difficult to converge and more generally can improve the performance of variational algorithms by reducing the optimization runtime and/or the depth of circuits that encode the solution candidate(s).

Physics↗

Noise dynamics of quantum annealers: estimating the effective noise using idle qubits

Quantum annealing is a type of analog computation that aims to use quantum mechanical fluctuations in search of optimal solutions of QUBO (quadratic unconstrained binary optimization) or, equivalently, Ising problems. Since NP-hard problems can in general be mapped to Ising and QUBO formulations, the quantum annealing paradigm has the potential to help solve various NP-hard problems. Current quantum annealers, such as those manufactured by D-Wave Systems, Inc. have various practical limitations including the size (number of qubits) of the problem that can be solved, the qubit connectivity, and error due to the environment or system calibration, which can reduce the quality of the solutions. Typically, for an arbitrary problem instance, the corresponding QUBO (or Ising) structure will not natively embed onto the available qubit architecture on the quantum chip. Thus, in these cases, a minor embedding of the problem structure onto the device is necessary. However, minor embeddings on these devices do not always make use of the full sparse chip hardware graph, and a large portion of the available qubits stay unused during quantum annealing. In this work, we embed a disjoint random QUBO on the unused parts of the chip alongside the QUBO to be solved, which acts as an indicator of the solution quality of the device over time. Using experiments on three different D-Wave quantum annealers, we demonstrate that (i) long term trends in solution quality exist on the D-Wave device, and (ii) the unused qubits can be used to measure the current level of noise of the quantum system.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A subpolar-focused stratospheric aerosol injection deployment scenario

Stratospheric aerosol injection (SAI) is a prospective climate intervention technology that would seek to abate climate change by deflecting back into space a small fraction of the incoming solar radiation. While most consideration given to SAI assumes a global intervention, this paper considers an alternative scenario whereby SAI might be deployed only in the subpolar regions. Subpolar deployment would quickly envelope the poles as well and could arrest or reverse ice and permafrost melt at high latitudes. This would yield global benefit by retarding sea level rise. Given that effective SAI deployment could be achieved at much lower altitudes in these regions than would be required in the tropics, it is commonly assumed that subpolar deployment would present fewer aeronautical challenges. An SAI deployment intended to reduce average surface temperatures in both the Arctic and Antarctic regions by 2 °C is deemed here to be feasible at relatively low cost with conventional technologies. However, we do not find that such a deployment could be undertaken with a small fleet of pre-existing aircraft, nor that relegating such a program to these sparsely populated regions would obviate the myriad governance challenges that would confront any such deployment. Nevertheless, given its feasibility and potential global benefit, the prospect of subpolar-focused SAI warrants greater attention.

54 ENVIRONMENTAL SCIENCES↗