Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “state distribution”

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 613 records · Page 34

Parameter estimation in nonlinear distributed systems - Approximation theory and convergence results

An abstract approximation framework and convergence theory is described for Galerkin approximations applied to inverse problems involving nonlinear distributed parameter systems. Parameter estimation problems are considered and formulated as the minimization of a least-squares-like performance index over a compact admissible parameter set subject to state constraints given by an inhomogeneous nonlinear distributed system. The theory applies to systems whose dynamics can be described by either time-independent or nonstationary strongly maximal monotonic operators defined on a reflexive Banach space which is densely and continuously embedded in a Hilbert space. It is demonstrated that if readily verifiable conditions on the system's dependence on the unknown parameters are satisfied, and the usual Galerkin approximation assumption holds, then solutions to the approximating problems exist and approximate a solution to the original infinite-dimensional identification problem.

Banks, H. T.↗

Resource Prospector, the Decadal Survey and the Scientific Context for the Exploration of the Moon

The Inner Planets Panel of the Planetary Exploration Decadal Survey defined several science questions related to the origins, emplacement, and sequestration of lunar polar volatiles: 1. What is the lateral and vertical distribution of the volatile deposits? 2. What is the chemical composition and variability of polar volatiles? 3. What is the isotopic composition of the volatiles? 4. What is the physical form of the volatiles? 5. What is the rate of the current volatile deposition? A mission concept study, the Lunar Polar Volatiles Explorer (LPVE), defined a approximately $1B New Frontiers mission to address these questions. The NAS/NRC report, 'Scientific Context for the Exploration of the Moon' identified he lunar poles as special environments with important implications. It put forth the following goals: Science Goal 4a-Determine the compositional state (elemental, isotopic, mineralogic) and compositional distribution (lateral and depth) of the volatile component in lunar polar regions. Science Goal 4b-Determine the source(s) for lunar polar volatiles. Science Goal 4c-Understand the transport, retention, alteration, and loss processes that operate on volatile materials at permanently shaded lunar regions. Science Goal 4d-Understand the physical properties of the extremely cold (and possibly volatile rich) polar regolith. Science Goal 4e-Determine what the cold polar regolith reveals about the ancient solar environment.

Moon↗

Investigating Low-Altitude Constellations of Ad-Hoc Lunar PNT System for Distributed Spacecraft Autonomy

In this study, we examine a low-altitude Lunar Position, Navigation, and Timing (LPNT) constellations and the localization performance of Centralized Extended Kalman Filter (CEKF) and Decentralized Extended Kalman Filter (DEKF) algorithms. The primary investigation involves a 100-node swarm operating at a 100 km altitude, in contrast to previous studies that examined a 21-node asset in a frozen-orbit at 5,500 km. The autonomous operation of large-scale swarm is based on two-way Inter-Satellite Link (ISL) measurements, which involve pseudoranges and relative velocities among swarm nodes. We perform a numerical assessment of the two filtering approaches, utilizing ‘fully sampled’ measurements from all available assets as well as ‘two ISL’ measurements where each spacecraft is restricted to only two antennas. This research includes an analysis of CEKF under 2-ISL constraints and evaluates the performance of DEKF in a 100-node swarm, which has not been explored in previous studies. In addition, we examine the impact of increasing the sampling frequency for DEKF, showing that the update cycle can be shortened from a 10-minute interval. A novel approach for ‘2-ISL limited’ DEKF will also be introduced, using a matching formulation that exhaustively enumerates all potential matches. This study provides valuable insights into large-scale distributed swarm operations, considering various filter configurations, sampling frequencies, matching strategies, and scalability of CEKF and DEKF for low-altitude LPNT applications. The Lunar PNT technology plays a key role in providing reliable and robust navigation services on the Moon's surface and the South pole, where the primary Lunar missions are planned. To support upcoming Lunar missions, including small satellites from NASA's Commercial Lunar Payload Services program, the Lunar PNT system must be adaptable to smaller platforms like CubeSats. Driven by the growing involvement of public and private exploration partnerships, the traditional low Earth orbit missions are shifting to beyond geosynchronous orbit [1]. These upcoming missions aim to foster a sustainable and innovative exploration program, in collaboration with commercial and international partners, to facilitate human expansion throughout the solar system and return new knowledge and opportunities to Earth [2]. As part of this trend, there are increasing efforts to utilize science missions in Lunar orbit to develop a non-dedicated and ad-hoc PNT network system. Two traditional approaches, the Deep Space Network (DSN) and the weak signal Global Positioning System (GPS), are established deep-space navigation technologies for missions beyond the geosynchronous orbit. Beginning in 1958, the DSN was developed to communicate with the Explorer 1 spacecraft based on the use of radiometric tracking in spacecraft navigation [3]. The DSN is capable of providing nearly unfettered coverage to spacecraft beyond low-Earth orbit (LEO), however, increased space mission volume has created concerns about future expectations of DSN usage for spacecraft navigation [4]. For cislunar mission applications, the position accuracy using DSN achieves 100 m (3σ) with at least three geometrically diverse ground stations when using radiometric tracking alone [5]. The DSN's dependence on Earth-based ground stations restricts its operational capabilities to periods of Earth visibility. This limitation, coupled with its poor localization performance, renders the DSN unsuitable for future lunar missions that demand continuous tracking and precise positioning. To satisfy the increasing requirements of DSN in Lunar applications, spacecrafts are also required to improve their onboard antenna power and efficiency of the transmission. However, there is an important aggregate cost trade between adding capabilities to every spacecraft and adding to a capacity on the ground that serves multiple spacecraft [6]. A weak GPS system can provide PNT service while the user spacecraft is bound to the Moon, leveraging a single, steerable high gain antenna with the relatively narrow beam which includes all the sources in its field of view [7]. However, the higher the altitude the receiver is above the GPS constellations, the poorer and the weaker are the relative geometry and the received signal powers, respectively, leading to a significant navigation accuracy reduction [8]. The transmitted power becomes weaker with increasing distance from the Earth as well as signals tracked from one of the side lobes of the GPS antenna pattern. As a results, the number of visible satellites and relative geometric condition of the GPS satellites at very high altitude drops dramatically and reduces the navigation solution accuracy. Therefore, the weak GPS system is also not an ideal way to provide PNT service to upcoming Lunar missions when considering its limited geometric condition and the recued navigation accuracy. Another navigation approach on the Moon is being developed, similar to the Global Navigation Satellite System (GNSS) on Earth, aiming to offer navigation service with continuous 24/7 coverage across the entire Lunar surface. For example, lunar communications relay and navigation systems (LCRNS) by NASA and Lunar navigation satellite systems (LNSS) by JAXA are designed to serve as dedicated Position, Navigation, and Timing (PNT) systems for the Moon. However, designing a dedicated LNSS and PNT service involves additional challenges, which are unique to the lunar environment, including limited payload capacity for the CubeSat platform, i.e., the size, weight, and power (SWaP) of the onboard clock, limited lunar ground monitoring stations, and limited financial investment as compared to the legacy Earth-GPS [9]. NASA’s focus on utilizing CubeSat platforms on the Moon leads to an alternative Lunar navigation platform that leverages the existing Lunar science and exploration assets. The small satellites used in Lunar missions can be used to create a low-cost, autonomous, ad-hoc, and on-demand mission-centric Lunar PNT swarm capable of providing PNT services to these low-cost lunar missions [10]. As upcoming Lunar missions will often operate at low-altitude about 30 km to 100 km for scientific observations and mapping purposes, the low-altitude orbital constellations could be employed to create an ad-hoc Lunar PNT system. However, several issues must be addressed, such as the instability of these orbits, which often require maintenance or are only suitable for short-duration missions, operating for fewer than 90 days. Additionally, at an altitude of 100 km, the satellites have a limited period during which they are above the horizon and capable of providing PNT service to users. The implementation of a non-dedicated, ad-hoc Lunar navigation constellation facilitates on-demand PNT services. A preliminary study of ad-hoc Lunar PNT system was conducted using 21 spacecraft in 5,5000 km altitude frozen orbits to test its feasibility and a basic performance of orbital asset localization among ad-hoc Lunar constellations in small satellites format [10]. These swarm assets are designed for autonomous localization with minimal Earth interaction, reducing dependency on bandwidth and ground resources. The design in [10] demonstrated the feasibility of a decentralized PNT approach, specifically employing a DEKF approach for state estimation, which helps minimize onboard operating costs. The DEKF method distributes computation across individual satellites, which lightens the computational load while maintaining accuracy in orbit ephemeris and clock offsets, similar to centralized systems [11]. In a follow-on study [12], each spacecraft was limited to 2 communications antennae, forcing the selection of measurements and scheduling spacecraft activities to perform the measurements. A matching algorithm is implemented to select the best measurements and schedule position estimation updates. The decentralized localization performance is also investigated with increasing levels of network degradation for swarm assets considering the impact of intermittent and permanent communication failure, to demonstrate the robustness and fidelity of the decentralized Lunar PNT service [13]. This study confirmed that the ad-hoc PNT constellations in frozen orbit are highly robust and resilient to communication failures. However, unlike frozen orbit swarm assets, the low-altitude satellites have a limited ground view at an altitude of 100 km, where the ad-hoc Lunar constellation consists of 98 low-altitude satellites, evenly distributed across seven circular polar orbital planes, alongside two satellites in a frozen orbit at an altitude of 5,500 km (Figure 1). Therefore, the number of satellites visible to ground users is significantly limited in low-altitude orbit constellations. As each visibility of a spacecraft remains intact for only a few ticks before it moves out of the field of view, the ground user encounters challenges in maintaining continuous navigation service, resulting in sparse availability and provision of Lunar PNT system. Consequently, service availability is primarily restricted to the Lunar South Pole region (Figure 2). Given these limitations and concerns, the localization performance of low-altitude swarm assets will be assessed in this study. We focus on the investigation of the localization performance of low-altitude swarm assets and ground users near the Lunar South Pole. The overall flow of the Lunar PNT simulation incorporates the DEKF approach of asset localization and the weighted least-squares approach in user localization (Figure 3). The autonomous Lunar PNT simulation is primarily implemented in MATLAB, where the DEKF based on the matching scheduler is implemented with Google’s OR-tools as a model builder and Gurobi optimization tool as a backend solver. The General Mission Analysis Tool (GMAT) is utilized to generate ephemeris data for swarm assets, and accounts for satellite orbital details, mass, and perturbations like solar radiation pressure and drag coefficients. Each ephemeris dataset is produced in the Moon International Celestial Reference Frame (ICRF) inertial coordinate system. For state estimation, the distributed swarm assets rely on two-way Inter-Satellite Link (ISL) measurements, which involve tracking pseudoranges and relative velocities between visible satellites and anchor nodes during each observation. Numerical evaluations of the decentralized localization process are conducted to demonstrate the feasibility of the low-altitude PNT system in providing reliable navigation services. The main approach involves using DEKF and CEKF to localize 100 satellites in low-altitude constellations, where the CEKF is implemented to serve as a baseline for comparing the performance of distributed algorithms. In both cases, we evaluate ‘fully sampled’ measurements from all available assets, and ‘two ISL’ measurements when spacecraft are constrained to have only two antennas. We test four estimation techniques: CEKF fully sampled, CEKF two ISL, DEKF fully sampled, and DEKF two ISL filters. As the DEKF update cycle is comprised of network setup, communication, and computations, a global broadcast network and 2-way ISL network setup will take from 4 to 6 minutes as maximum [12]. In this simulation, the DEKF update cycle is set to 10 minutes, including a 4-minute latency for obtaining and computing the actual measurement updates. We experiment an increased update cycle to demonstrate the feasibility and evaluate the impact on localization performance using various tuning values for measurement noise covariances (Figures 4 and 5). By comparing centralized and decentralized approaches using a matching algorithm, we analyze the influence of cross-correlation factors in the covariance matrix, assuming 100% reliability of all assets and measurements. The increased frequency and the adjustments of tuning parameters reveal distinct error patterns between the two scenarios. The localization accuracy of the swarm assets and ground users is assessed by taking the median error across 100 assets and one ground user (84.9°S, 137.5°E) over 7-day simulation period (Table 1). Since the user localization accuracy is significantly affected by the performance of the swarm assets, it is crucial to maintain high localization accuracy within the swarm. This study will continue to explore decentralized filtering for autonomous LPNT operations, with further investigation of an 'iterative' matching approach which enumerates every valid matching pair, planned for the following month.

Yeji Kim↗

CFDP Performance over Weather-Dependent Ka-Band Channel

This study presents an analysis of the delay performance of the CCSDS File Delivery Protocol (CFDP) over weather-dependent Ka-band channel. The Ka-band channel condition is determined by the strength of the atmospheric noise temperature, which is weather dependent. Noise temperature data collected from the Deep Space Network (DSN) Madrid site is used to characterize the correlations between good and bad channel states in a two-state Markov model. Specifically, the probability distribution of file delivery latency using the CFDP deferred Negative Acknowledgement (NAK) mode is derived and quantified. Deep space communication scenarios with different file sizes and bit error rates (BERs) are studied and compared. Furthermore, we also examine the sensitivity of our analysis with respect to different data sampling methods. Our analysis shows that while the weather-dependent channel only results in fairly small increases in the average number of CFDP retransmissions required, the maximum number of transmissions required to complete 99 percentile, on the other hand, is significantly larger for the weather-dependent channel due to the significant correlation of poor weather states.

weather↗

CFDP Performance over Weather-dependent Ka-band Channel

This study presents an analysis of the delay performance of the CCSDS File Delivery Protocol (CFDP) over weather-dependent Ka-band channel. The Ka-band channel condition is determined by the strength of the atmospheric noise temperature, which is weather dependent. Noise temperature data collected from the Deep Space Network (DSN) Madrid site is used to characterize the correlations between good and bad channel states in a two-state Markov model. Specifically, the probability distribution of file delivery latency using the CFDP deferred Negative Acknowledgement (NAK) mode is derived and quantified. Deep space communication scenarios with different file sizes and bit error rates (BERs) are studied and compared. Furthermore, we also examine the sensitivity of our analysis with respect to different data sampling methods. Our analysis shows that while the weather-dependent channel only results in fairly small increases in the average number of CFDP retransmissions required, the maximum number of transmissions required to complete 99 percentile, on the other hand, is significantly larger for the weather-dependent channel due to the significant correlation of poor weather states.

deep space communications↗

Solid state remote power controllers for 120 VDC power systems

Remote Power controllers (RPCs) are devices that combine in one unit the capability to perform all the needed functions of load switching and provide total system protection of equipment and wires. The unique developments of solid-state RPCs for 120 Vdc power distribution systems are reviewed. The discussion covers design guidelines, power switch design concepts, and performance effectiveness. An NPN transistor is used as the the basic switch element. Since the ultimate goal of the 120 Vdc RPC program is to demonstrate technology readiness, the final phase is directed to the design, fabrication, and testing of multi-chip hybrid prototypes in hermetically sealed packages. RPCs have potential application in spacecraft and aircraft electrical systems, in transportation systems, and various industrial applications. The upper voltage limitation on the RPC design is related to the capability and availability of suitable high-voltage power transistors. The merits of solid-state RPCs are noted.

Sundberg, G. R.↗

Global Interoperability of High Definition Video Streams Via ACTS and Intelsat

In 1993, a proposal at the Japan.-U.S. Cooperation in Space Program Workshop lead to a subsequent series of satellite communications experiments and demonstrations, under the title of Trans-Pacific High Data Rate Satellite Communications Experiments. The first of which is a joint collaboration between government and industry teams in the United States and Japan that successfully demonstrated distributed high definition video (HDV) post-production on a global scale using a combination of high data rate satellites and terrestrial fiber optic asynchronous transfer mode (ATM) networks. The HDV experiment is the first GIBN experiment to establish a dual-hop broadband satellite link for the transmission of digital HDV over ATM. This paper describes the team's effort in using the NASA Advanced Communications Technology Satellite (ACTS) at rates up to OC-3 (155 Mbps) between Los Angeles and Honolulu, and using Intelsat at rates up to DS-3 (45 Mbps) between Kapolei and Tokyo, with which HDV source material was transmitted between Sony Pictures High Definition Center (SPHDC) in Los Angeles and Sony Visual Communication Center (VCC) in Shinagawa, Tokyo. The global-scale connection also used terrestrial networks in Japan, the States of Hawaii and California. The 1.2 Gbps digital HDV stream was compressed down to 22.5 Mbps using a proprietary Mitsubishi MPEG-2 codec that was ATM AAL-5 compatible. The codec: employed four-way parallel processing. Improved versions of the codec are now commercially available. The successful post-production activity performed in Tokyo with a HDV clip transmitted from Los Angeles was predicated on the seamless interoperation of all the equipment between the sites, and was an exciting example in deploying a global-scale information infrastructure involving a combination of broadband satellites and terrestrial fiber optic networks. Correlation of atmospheric effects with cell loss, codec drop-out, and picture quality were made. Current efforts in the Trans-Pacific series plan to examine the use of Internet Protocol (IP)-related technologies over such an infrastructure. The use of IP allows the general public to be an integral part of the exciting activities, helps to examine issues in constructing the solar-system internet, and affords an opportunity to tap the research results from the (reliable) multicast and distributed systems communities. The current Trans- Pacific projects, including remote astronomy and digital library (visible human) are briefly described.

Hsu, Eddie↗

Energy content of stormtime ring current from phase space mapping simulations

We perform a phase space mapping study to estimate the enhancement in energy content that results from stormtime particle transport in the equatorial magnetosphere. Our pre-storm phase space distribution is based on a steady-state transport model. Using results from guiding-center simulations of ion transport during model storms having main phases of 3 hr, 6 hr, and 12 hr, we map phase space distributions of ring current protons from the pre-storm distribution in accordance with Liouville's theorem. We find that transport can account for the entire ten to twenty-fold increase in magnetospheric particle energy content typical of a major storm if a realistic stormtime enhancement of the phase space density f is imposed at the nightside tail plasma sheet (represented by an enhancement of f at the neutral line in our model).

Chen, Margaret W.↗

Global simulation of auroral arcs

Numerical simulation of global formation of auroral arcs is carried out for a three-dimensional coupled ionosphere-magnetosphere system. With the presence of the stationary convection electric field, the field-aligned current driven by the polarization and drift of the ionospheric conductivity enhancement is shown to be responsible for the feedback instability of the shear Alfven wave. The simulation result demonstrates the dynamic process of auroral arc formation, where an ionospheric perturbation nearly resonant with the magnetospheric Alfven wave (toroidal mode) can eventually grow into auroral arcs. Many of the important characteristics of quiet auroral arcs are demonstrated, and it is shown that the overall distribution of auroras depends critically on the stationary state of the auroral oval characterized by the large-scale Birkeland current, the electric field, and the electron density distribution.

Miura, A.↗

Distributed Damage Estimation for Prognostics based on Structural Model Decomposition

Model-based prognostics approaches capture system knowledge in the form of physics-based models of components, and how they fail. These methods consist of a damage estimation phase, in which the health state of a component is estimated, and a prediction phase, in which the health state is projected forward in time to determine end of life. However, the damage estimation problem is often multi-dimensional and computationally intensive. We propose a model decomposition approach adapted from the diagnosis community, called possible conflicts, in order to both improve the computational efficiency of damage estimation, and formulate a damage estimation approach that is inherently distributed. Local state estimates are combined into a global state estimate from which prediction is performed. Using a centrifugal pump as a case study, we perform a number of simulation-based experiments to demonstrate the approach.

Daigle, Matthew↗

Study of GaAs-oxide interface by transient capacitance spectroscopy - Discrete energy interface states

Interface states and bulk GaAs energy levels were simultaneously investigated in GaAs MOS structures prepared by anodic oxidation. These two types of energy levels were successfully distinguished by carrying out a comparative analysis of deep level transient capacitance spectra of the MOS structures and MS structures prepared on the same samples of epitaxially grown GaAs. The identification and study of the interface states and bulk levels was also performed by investigating the transient capacitance spectra as a function of the filling pulse magnitude. It was found that in the GaAs-anodic oxide interface there are states present with a discrete energy rather than with a continuous energy distribution. The value of the capture cross section of the interface states was found to be 10 to the 14th to 10 to the 15th/sq cm, which is more accurate than the extremely large values of 10 to the -8th to 10 to the -9th/sq cm reported on the basis of conductance measurements.

Kamieniecki, E.↗

Transient Numerical Modeling of Catalytic Channels

This paper presents a transient model of catalytic combustion suitable for isolated channels and monolith reactors. The model is a lumped two-phase (gas and solid) model where the gas phase is quasi-steady relative to the transient solid. Axial diffusion is neglected in the gas phase; lateral diffusion, however, is accounted for using transfer coefficients. The solid phase includes axial heat conduction and external heat loss due to convection and radiation. The combustion process utilizes detailed gas and surface reaction models. The gas-phase model becomes a system of stiff ordinary differential equations while the solid phase reduces, after discretization, into a system of stiff ordinary differential-algebraic equations. The time evolution of the system came from alternating integrations of the quasi-steady gas and transient solid. This work outlines the numerical model and presents some sensitivity studies on important parameters including internal transfer coefficients, catalytic surface site density, and external heat-loss (if applicable). The model is compared to two experiments using CO fuel: (1) steady-state conversion through an isothermal platinum (Pt) tube and (2) transient propagation of a catalytic reaction inside a small Pt tube. The model requires internal mass-transfer resistance to match the experiments at lower residence times. Under mass-transport limited conditions, the model reasonably predicted exit conversion using global mass-transfer coefficients. Near light-off, the model results did not match the experiment precisely even after adjustment of mass-transfer coefficients. Agreement improved for the first case after adjusting the surface kinetics such that the net rate of CO adsorption increased compared to O2. The CO / O2 surface mechanism came from a sub-set of reactions in a popular CH4 / O2 mechanism. For the second case, predictions improved for lean conditions with increased external heat loss or adjustment of the kinetics as in the first case. Finally, the results show that different initial surface-species distribution leads to different steady-states under certain conditions. These results demonstrate the utility of a lumped two-phase model of a transient catalytic combustor with detailed chemistry.

Struk, Peter M.↗

Orbiting Debris: a Space Environmental Problem. Background Paper

Man-made debris, now circulating in a multitude of orbits about Earth as the result of the exploration and use of the space environment, poses a growing hazard to future space operations. The 6,000 or so debris objects large enough to be cataloged by the U.S. Space Surveillance Network are only a small percentage of the total debris capable of damaging spacecraft. Unless nations reduce the amount of orbital debris they produce, future space activities could suffer loss of capability, destruction of spacecraft, and perhaps even loss of life as a result of collisions between spacecraft and debris. Better understanding of the extent and character of "space junk" will be crucial for planning future near-Earth missions, especially those projects involving humans in space. This OTA background paper summarizes the current state of knowledge about the causes and distribution of orbiting debris, and examines R&D needs for reducing the problem. As this background paper notes, addressing the problem will require the involvement of all nations active in space. The United States has taken the lead to increase international understanding of the issue but much work lies ahead. In undertaking this background paper, OTA sought the contributions of a broad spectrum of knowledgeable individuals and organizations. Some provided information, others reviewed drafts. OTA gratefully acknowledges their contributions of time and intellectual effort. As with all OTA studies, the content of this background paper is the sole responsibility .of the Office of Technology Assessment and does not necessarily represent the views of our advisors or reviewers.

Lionel S. Johns↗

Feedback control of two beam, two joint systems with distributed flexibility

The control of the flexible motion in a plane of two pinned beams is addressed with application to remote manipulators. Three types of linear feedback control schemes are considered: joint angle and velocity feedback with (GRC) and without (IJC) cross joint feedback, and feedback of flexible state variables (FFC). Two models of the distributed flexibility are presented along with some results obtained from them. The relative merit of the three control schemes is discussed.

Book, W. J.↗

Coupling between fluid dynamics and energy addition in arcjet and microwave thrusters

A new approach to numerically solving the problem of the constricted electric arcjet is presented. An Euler Implicit finite difference scheme is used to solve the full compressible Navier Stokes equations in two dimensions. The boundary and initial conditions represent the constrictor section of the arcjet, and hydrogen is used as a propellant. The arc is modeled as a Gaussian distribution across the centerline of the constrictor. Temperature, pressure and velocity profiles for steady state converged solutions show both axial and radial changes in distributions resulting from their interaction with the arc energy source for specific input conditions. The temperature rise is largest at the centerline where there is a the greatest concentration arc energy. The solution does not converge for all initial inputs and the limitations in the range of obtainable solutions are discussed.

Micci, M. M.↗

A Blueprint for Demonstrating Quantum Supremacy with Superconducting Qubits

Long coherence times and high fidelity control recently achieved in scalable superconducting circuits paved the way for the growing number of experimental studies of many-qubit quantum coherent phenomena in these devices. Albeit full implementation of quantum error correction and fault tolerant quantum computation remains a challenge the near term pre-error correction devices could allow new fundamental experiments despite inevitable accumulation of errors. One such open question foundational for quantum computing is achieving the so called quantum supremacy, an experimental demonstration of a computational task that takes polynomial time on the quantum computer whereas the best classical algorithm would require exponential time and/or resources. It is possible to formulate such a task for a quantum computer consisting of less than a 100 qubits. The computational task we consider is to provide approximate samples from a non-trivial quantum distribution. This is a generalization for the case of superconducting circuits of ideas behind boson sampling protocol for quantum optics introduced by Arkhipov and Aaronson. In this presentation we discuss a proof-of-principle demonstration of such a sampling task on a 9-qubit chain of superconducting gmon qubits developed by Google. We discuss theoretical analysis of the driven evolution of the device resulting in output approximating samples from a uniform distribution in the Hilbert space, a quantum chaotic state. We analyze quantum chaotic characteristics of the output of the circuit and the time required to generate a sufficiently complex quantum distribution. We demonstrate that the classical simulation of the sampling output requires exponential resources by connecting the task of calculating the output amplitudes to the sign problem of the Quantum Monte Carlo method. We also discuss the detailed theoretical modeling required to achieve high fidelity control and calibration of the multi-qubit unitary evolution in the device. We use a novel cross-entropy statistical metric as a figure of merit to verify the output and calibrate the device controls. Finally, we demonstrate the statistics of the wave function amplitudes generated on the 9-gmon chain and verify the quantum chaotic nature of the generated quantum distribution. This verifies the implementation of the quantum supremacy protocol.

Kechedzhi, Kostyantyn↗

A search for chemical laser action in low pressure metal vapor flames

Optical emissions were studied from low pressure (approximately 1 torr) dilute diffusion flames of Ca and Mg vapor with O2, N2O and mixtures of CCl4 and O2. The Ca flames with O2 and N2O revealed high vibrational excitation of the product CaO molecule (up to v=30). The flames with CCl4 revealed extreme nonequilibrium metal atom electronic excitation, up to the metal atom ionization limit (6.1 eV for Ca, 7.6 eV for Mg). The metal atom excited electronic state populations did not follow a Boltzmann distribution, but the excitation rates ('pumping rate') were found to obey an Arrhenius-type expression, with the electronic excitation energy playing the role of activation energy and a temperature of about 5000 K for triplet excited states and 2500 K for singlets (vs. approximately 500 K translational temperature).

Zwillenberg, M. L.↗

Initial analysis of “stone” size Ryugu samples: current status

As a part of the initial analysis of the Ryugu samples, we perform a variety of analyses of millimeter-sized "stones". Our goals are to elucidate the entire formation process of C-type asteroid Ryugu from the viewpoint of petrology and mineralogy ande obtain necessary information by sample analysis, and then simulate the formation of Ryugu based on the evidencef obtained rom sample analysis. Eighteen stones (8 from Room A and 10 from Room C) were received from the ISAS curation facility on June 1, 2021, and brought into a fully nitrogen-displaced glove box at Tohoku University. At the same time, we also received the powder samples from Room A and Room C. To completely block the atmosphere from leaking into the container, all samples were put in the sample transport containers prepared by ISAS, transferred from the main chamber to the glovebox at ISAS, and then all containers were completely sealed in plastic bags with moisture and oxygen absorbers. No moisture or oxygen was detected when the bags were opened in the glove box at Tohoku University, so it was confirmed that there was no exposure to the atmosphere during transport. To date, a number of analyses have been carried out successfully and almost on schedule. The analysis started with the measurement of reflectance spectra, which are sensitive to atmospheric oxidation, hydroxylation, and adsorbed water. The visible, near-infrared, and mid-infrared reflectance spectra were measured while the samples were kept airtight. The spectra of powder samples and stone samples (as aggregates and as single stone) were successfully obtained. A major feature of the stone team's analysis is the use of synchrotron radiation facilities around the world. Since this analysis is non-destructive, stone samples whose reflectance spectra were measured were sent to KEK, SPring-8, ESRF (France), SOLEIL (France), DESY (Germany), and APS (USA). Using these synchrotron radiation facilities, high spatial resolution and sensitivity XRD, STXM, XANES, CT [1], IR-CT, FT-IR [2, 3], XRF, and Mössbauer analyses were performed. Most of the analyses were carried out under air-tight conditions on the stone samples and the particulates separated from the stone samples. These analyses allowed us to determine the three-dimensional distribution of minerals and elements, redox state, density and porosity of the stone samples. Furthermore, as a characteristic analysis of the stone team, light elemental analysis using negative Muon was performed at the MLF facility of J-PARC with an exceptionally long allocation of machine time. This is the only non-destructive method to measure the concentration of light elements in the whole (not the surface) of stone samples. Because the characteristic X-rays produced by muon irradiation are much higher in energy than the fluorescent X-rays produced by X-ray irradiation, there is little effect of self-absorption by the sample, and therefore, the concentration of light elements such as carbon, oxygen, and Na in the entire "stone" sample can be determined. Some stone samples are currently being measured for heat and strength physical properties in order to understand the physical properties of asteroid Ryugu. The data obtained from these measurements are useful for interpreting the remote-sensing data data taken from the surface layer of the asteroid Ryugu [4-7]. It is also important for understanding the behavior of the Ryugu material during impact events. The surfaces of many stone samples were observed by electron microscopy and other techniques, especially on natural “flat” surfaces formed on 5 stomes. As a result, characteristric mineral aggregates formed by the reactions with water and characteiristic impact features were observed on some of the samples. Based on the observations of surfaces and the synchrotron measurements of the whole stones, important objects such as characteristic structures and specific crystal aggregates for understanding the formation history of the asteroid were identified, and these parts were separated from the stone samples using an Xe beam (pFIB) and analyzed by various methods including transmission electron microscopy and synchrotron radiation analysis. Many stone samples, from which important objects have been separated, are embedded in epoxy resin and cut to produce many polished sections. Electron microscopy and spectroscopic measurements of the polished surfaces are being carried out to reveal the detailed mineralogical properties and elemental distribution inside the stone samples. In the fall, machine time for synchrotron radiation will begin, and we plan to analyze single crystals and characteristic objects separated from the stone samples.

Tomoki Nakamura↗