Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “quantum probability”

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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Identifying and tracking bubbles and drops in simulations: A toolbox for obtaining sizes, lineages, and breakup and coalescence statistics

Knowledge of bubble and drop size distributions in two-phase flows is important for characterizing a wide range of phenomena, including combustor ignition, sonar communication, and cloud formation. The physical mechanisms driving the background flow also drive the time evolution of these distributions. Accurate and robust identification and tracking algorithms for the dispersed phase are necessary to reliably measure this evolution and thereby quantify the underlying mechanisms in interface-resolving flow simulations. The identification of individual bubbles and drops traditionally relies on an algorithm used to identify connected regions. This traditional algorithm can be sensitive to the presence of spurious structures. A cost-effective refinement is proposed to maximize volume accuracy while minimizing the identification of spurious bubbles and drops. An accurate identification scheme is crucial for distinguishing bubble and drop pairs with large size ratios. The identified bubbles and drops need to be tracked in time to obtain breakup and coalescence statistics that characterize the evolution of the size distribution, including breakup and coalescence frequencies, and the probability distributions of parent and child bubble and drop sizes. An algorithm based on mass conservation is proposed to construct bubble and drop lineages using simulation snapshots that are not necessarily from consecutive time steps. These lineages are then used to detect breakup and coalescence events, and obtain the desired statistics. Accurate identification of large-size-ratio bubble and drop pairs enables accurate detection of breakup and coalescence events over a large size range. Accurate detection of successive breakup and coalescence events requires that the snapshot interval be an order of magnitude smaller than the characteristic breakup and coalescence times to capture these successive events while minimizing the identification of repeated confounding events. Together, these algorithms serve as a toolbox for detailed analysis of two-phase simulations, and enable insights into the mechanisms behind bubble and drop formation and evolution in flows of practical importance.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

A Deep Study of an Intermediate-age Open Cluster SAI 35 (Juchert 20) Using Ground-based Imaging and Gaia EDR3 Astrometry

We present a CCD UBVI photometric study of poorly studied intermediate-age open cluster SAI 35 (Juchert 20) for the first time. To accomplish this study, we also used LAMOST DR5, Two Micron All Sky Survey, and Gaia EDR3 databases. We identified 214 most probable cluster members with membership probability higher than 50%. The mean proper motion of the cluster is found as μ{sub α}cos δ=1.10±0.01 and μ {sub δ} = −1.66 ± 0.01 mas yr{sup −1}. We find the normal interstellar extinction law using the various two-color diagrams. The age, distance, reddening, and radial velocity of the cluster are estimated to be 360 ± 40 Myr, 2.9 ± 0.15 kpc, 0.72 ± 0.05 mag, and −91.62 ± 6.39 km s{sup −1}, respectively. The overall mass function slope for main-sequence stars is found to be 1.49 ± 0.16 within the mass range 1.1–3.1 M {sub ⊙}, which is in agreement with Salpeter’s value within uncertainty. The present study demonstrates that SAI 35 is dynamically relaxed. Galactic orbital parameters are determined using Galactic potential models. We found that this object follows a circular path around the Galactic center.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Entanglement phase transitions in random stabilizer tensor networks

Here we explore a class of random tensor network models with “stabilizer” local tensors which we name random stabilizer tensor networks (RSTNs). For RSTNs defined on a two-dimensional square lattice, we perform extensive numerical studies of entanglement phase transitions between volume-law and area-law entangled phases of the one-dimensional boundary states. These transitions occur when either (a) the bond dimension D of the constituent tensors is varied or (b) the tensor network is subject to random breaking of bulk bonds, implemented by forced measurements. In the absence of broken bonds, we find that the RSTN supports a volume-law entangled boundary state with bond dimension D ≥ 3 where D is a prime number, and an area-law entangled boundary state for D = 2. Upon breaking bonds at random in the bulk with probability p, there exists a critical measurement rate p c for each D ≥ 3 above which the boundary state becomes area-law entangled. To explore the conformal invariance at these entanglement transitions for different prime D, we consider tensor networks on a finite rectangular geometry with a variety of boundary conditions, and extract universal operator scaling dimensions via extensive numerical calculations of the entanglement entropy, mutual information, and mutual negativity at their respective critical points. Our results at large D approach known universal data of percolation conformal field theory, while showing clear discrepancies at smaller D, suggesting a distinct entanglement transition universality class for each prime D. We further study universal entanglement properties in the volume-law phase and demonstrate quantitative agreement with the recently proposed description in terms of a directed polymer in a random environment.

36 MATERIALS SCIENCE↗
DOE OSTI2021

OSSOS. XXI. Collision Probabilities in the Edgeworth–Kuiper Belt

Here, we present results on the intrinsic collision probabilities, P {sub I}, and range of collision speeds, V {sub I}, as a function of the heliocentric distance, r, in the trans-Neptunian region. The collision speed is one of the parameters that serves as a proxy for a collisional outcome (e.g., disruption and scattering of fragments, or formation of a crater, as both processes are related to the impact energy). We utilize an improved and debiased model of the trans-Neptunian object (TNO) region from the “Outer Solar System Origins Survey” (OSSOS). It provides a well-defined model of TNO orbital distribution, based on multiple opposition observations of more than 1000 bodies. We compute collisional probabilities for the OSSOS models of the main classical, resonant, detached+outer, and scattering TNO populations. The intrinsic collision probabilities and collision speeds are computed using Öpik’s approach, as revised and modified by Wetherill for noncircular and inclined orbits. The calculations are carried out for each of the dynamical TNO groups, allowing for inter-population collisions as well as collisions within each TNO population, resulting in 28 combinations in total. Our results indicate that collisions in the trans-Neptunian region are possible over a wide range in (r, V {sub I}) phase space. Although collisions are calculated to happen within r ∼ 20–200 au and V {sub I} ∼ 0.1 km s{sup −1} to as high as V {sub I} ∼ 9 km s{sup −1}, most of the collisions are likely to happen at low relative velocities V {sub I} < 1 km s{sup −1} and are dominated by the main classical belt.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Distribution of centrality measures on undirected random networks via the cavity method

The Katz centrality of a node in a complex network is a measure of the node’s importance as far as the flow of information across the network is concerned. For ensembles of locally tree-like undirected random graphs, this observable is a random variable. Its full probability distribution is of interest but difficult to handle analytically because of its “global” character and its definition in terms of a matrix inverse. Leveraging a fast Gaussian Belief Propagation-Cavity algorithm to solve linear systems on tree-like structures, we show that i) the Katz centrality of a single instance can be computed recursively in a very fast way, and ii) the probability P ( K ) that a random node in the ensemble of undirected random graphs has centrality K satisfies a set of recursive distributional equations, which can be analytically characterized and efficiently solved using a population dynamics algorithm. We test our solution on ensembles of Erdős-Rényi and Scale Free networks in the locally tree-like regime, with excellent agreement. The analytical distribution of centrality for the configuration model conditioned on the degree of each node can be employed as a benchmark to identify nodes of empirical networks with over- and underexpressed centrality relative to a null baseline. We also provide an approximate formula based on a rank- 1 projection that works well if the network is not too sparse, and we argue that an extension of our method could be efficiently extended to tackle analytical distributions of other centrality measures such as PageRank for directed networks in a transparent and user-friendly way.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Social network structure and the spread of complex contagions from a population genetics perspective

Ideas, behaviors, and opinions spread through social networks. If the probability of spreading to a new individual is a non-linear function of the fraction of the individuals’ affected neighbors, such a spreading process becomes a “complex contagion”. This non-linearity does not typically appear with physically spreading infections, but instead can emerge when the concept that is spreading is subject to game theoretical considerations (e.g. for choices of strategy or behavior) or psychological effects such as social reinforcement and other forms of peer influence (e.g. for ideas, preferences, or opinions). Here we study how the stochastic dynamics of such complex contagions are affected by the underlying network structure. Motivated by simulations of complex contagions on real social networks, we present a framework for analyzing the statistics of contagions with arbitrary non-linear adoption probabilities based on the mathematical tools of population genetics. The central idea is to use an effective lower-dimensional diffusion process to approximate the statistics of the contagion. This leads to a tradeoff between the effects of ”selection” (microscopic tendencies for an idea to spread or die out), random drift, and network structure. Our framework illustrates intuitively several key properties of complex contagions: stronger community structure and network sparsity can significantly enhance the spread, while broad degree distributions dampen the effect of selection compared to random drift. Finally, we show that some structural features can exhibit critical values that demarcate regimes where global contagions become possible for networks of arbitrary size. Our results draw parallels between the competition of genes in a population and memes in a world of minds and ideas. Our tools provide insight into the spread of information, behaviors, and ideas via social influence, and highlight the role of macroscopic network structure in determining their fate.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Bootstrap-determined p values in lattice QCD

We present a general method to determine the probability that stochastic Monte Carlo data, in particular those generated in a lattice QCD calculation, would have been obtained were that data drawn from the distribution predicted by a given theoretical hypothesis. Such a probability, or p -value, is often used as an important heuristic measure of the validity of that hypothesis. The proposed method offers the benefit that it remains usable in cases where the standard Hotelling T 2 methods based on the conventional χ 2 statistic do not apply, such as for uncorrelated fits. Specifically, we analyze q 2 , defined as the correlated χ 2 statistic obtained using an arbitrary covariance matrix estimator, and show how to use the bootstrap as a data-driven method to determine the expected distribution of q 2 for a given hypothesis with minimal assumptions. This distribution can then be used to determine the p -value for a fit to the data. We also describe a bootstrap approach for quantifying the impact upon this p -value of estimating population parameters from a single ensemble of N samples. The overall method is accurate up to a 1 / N bias which we do not attempt to quantify. Published by the American Physical Society 2025

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

A Stellar Activity F-statistic for Exoplanet Surveys (SAFE)

In the search for planets orbiting distant stars, the presence of stellar activity in the atmospheres of observed stars can obscure the radial velocity signal used to detect such planets. Furthermore, this stellar activity contamination is set by the star itself and cannot simply be avoided with better instrumentation. Various stellar activity indicators have been developed that may correlate with this contamination. We introduce a new stellar activity indicator called the Stellar Activity F-statistic for Exoplanet surveys (SAFE) that has higher statistical power (i.e., probability of detecting a true stellar activity signal) than many traditional stellar activity indicators in a simulation study of an active region on a Sun-like star with a moderate-to-high signal-to-noise ratio. Also through simulation, the SAFE is demonstrated to be associated with the projected area on the visible side of the star covered by active regions. We also demonstrate that the SAFE detects statistically significant stellar activity in most of the spectra for HD 22049, a star known to have high stellar variability. Additionally, the SAFE is calculated for recent observations of the three low-variability stars HD 34411, HD 10700, and HD 3651, the latter of which is known to have a planetary companion. As expected, the SAFE for these three only occasionally detects activity. Furthermore, initial exploration appears to indicate that the SAFE may be useful for disentangling stellar activity signals from planet-induced Doppler shifts.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Continuum radiated power density $(P^{\textrm{C}}_{\textrm{rad}})$ and effective charge $(Z_{\textrm{eff}})$ estimates from multi-energy photon-counting measurements

Multi-energy soft x-ray pinhole cameras have been designed, built, calibrated, and operated at Madison symmetric Torus, Alcator C-Mod, and more recently at Tokamak a configuration variable and Tungsten Environment in Steady-State Tokamak (WEST), to measure plasma emission across multiple energy ranges. Here we describe a new methodology to estimate the local continuum radiated power density and the plasma effective charge (Z eff ) directly from photon-counting measurements of the line-free continuum emission (Bremsstrahlung and Radiative Recombination) in several energy bands between 11 and 18 keV. This capability is particularly valuable for confinement systems using metal plasma-facing components, where x-ray losses from interactions with the sputtered wall can represent a significant fraction of the total radiated power (P rad ). The approach leverages a well-characterized detector responsivity, modeled by a complementary error function, and interprets the differential multi-energy measurements between adjacent energy levels through the probability density function of a Gaussian distribution. The implementation of this diagnostic technique is currently under development on the WEST tokamak, aiming at the goal of providing real-time P rad and Z eff measurements during long-pulse operation (up to 1000 s) in the 2026 campaign.

WEST↗
DOE OSTI2023

Transport across interfaces in symmetric orbifolds

We examine how conformal boundaries encode energy transport coefficients — namely transmission and reflection probabilities — of corresponding conformal interfaces in symmetric orbifold theories. These constitute a large class of irrational theories and are closely related to holographic setups. Our central goal is to compare such coefficients at the orbifold point (a field theory calculation) against their values when the orbifold is highly deformed (a gravity calculation) — an approach akin to past AdS/CFT-guided comparisons of physical quantities at strong versus weak coupling. At the orbifold point, we find that the (weighted-average) transport coefficients are simply averages of coefficients in the underlying seed theory. We then focus on the symmetric orbifold of the T 4 sigma model interface CFT dual to type IIB supergravity on the 3d Janus solution. We compare the holographic transmission coefficient, which was found by [1], to that of the orbifold point. We find that the profile of the transmission coefficient substantially increases with the coupling, in contrast to boundary entropy. We also present some related ideas about twisted-sector data encoded by boundary states.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Dynamics of disordered mechanical systems with large connectivity, free probability theory, and quasi-Hermitian random matrices

Disordered mechanical systems with high connectivity represent a limit opposite to the more familiar case of disordered crystals. Individual ions in a crystal are subjected essentially to nearest-neighbor interactions. In contrast, the systems studied in this paper have all their degrees of freedom coupled to each other. Thus, the problem of linearized small oscillations of such systems involves two full positive-definite and non-commuting matrices, as opposed to the sparse matrices associated with disordered crystals. Consequently, the familiar methods for determining the averaged vibrational spectra of disordered crystals, introduced many years ago by Dyson and Schmidt, are inapplicable for highly connected disordered systems. In this paper we apply random matrix theory (RMT) to calculate the averaged vibrational spectra of such systems, in the limit of infinitely large system size. At the heart of our analysis lies a calculation of the average spectrum of the product of two positive definite random matrices by means of free probability theory techniques. We also show that this problem is intimately related with quasi-hermitian random matrix theory (QHRMT), which means that the ‘hamiltonian’ matrix is hermitian with respect to a non-trivial metric. This extends ordinary hermitian matrices, for which the metric is simply the unit matrix. The analytical results we obtain for the spectrum agree well with our numerical results. The latter also exhibit oscillations at the high-frequency band edge, which fit well the Airy kernel pattern. We also compute inverse participation ratios of the corresponding amplitude eigenvectors and demonstrate that they are all extended, in contrast with conventional disordered crystals. Finally, we compute the thermodynamic properties of the system from its spectrum of vibrations. In addition to matrix model analysis, we also study the vibrational spectra of various multi-segmented disordered pendula, as concrete realizations of highly connected mechanical systems. A universal feature of the density of vibration modes, common to both pendula and the matrix model, is that it tends to a non-zero constant at vanishing frequency.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Early-time measure in eternal inflation

In a situation like eternal inflation, where our data is replicated at infinitely-many other space-time events, it is necessary to make a prior assumption about our location to extract predictions. The principle of mediocrity entails that we live at asymptotic late times, when the occupational probabilities of vacua has settled to a near-equilibrium distribution. In this paper we further develop the idea that we instead exist during the approach to equilibrium, much earlier than the exponentially-long mixing time. In this case we are most likely to reside in vacua that are easily accessed dynamically. Using first-passage statistics, we prove that vacua that maximize their space-time volume at early times have: 1. maximal ever-hitting probability; 2. minimal mean first-passage time; and 3. minimal decay rate. These requirements are succinctly captured by an early-time measure. The idea that we live at early times is a predictive guiding principle, with many phenomenological implications. First, our vacuum should lie deep in a funneled region, akin to folding energy landscapes of proteins. Second, optimal landscape regions are characterized by relatively short-lived vacua, with lifetime of order the de Sitter Page time. For our vacuum, this lifetime is ~ 10 130 years, which is consistent with the Standard Model estimate due to Higgs metastability. Third, the measure favors vacua with small, positive vacuum energy. This can address the cosmological constant problem, provided there are sufficiently many vacua in the entire ensemble of funnels. As a concrete example, we study the Bousso-Polchinski lattice of flux vacua, and find that the early-time measure favors lattices with the fewest number of flux dimensions. This favors compactifications with a large hierarchy between the lightest modulus and all other Kähler and complex structure moduli.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Three-State Majority-vote Model on Scale-Free Networks and the Unitary Relation for Critical Exponents

Abstract We investigate the three-state majority-vote model for opinion dynamics on scale-free and regular networks. In this model, an individual selects an opinion equal to the opinion of the majority of its neighbors with probability 1 − q , and different to it with probability q . The parameter q is called the noise parameter of the model. We build a network of interactions where z neighbors are selected by each added site in the system, a preferential attachment network with degree distribution k − λ , where λ = 3 for a large number of nodes N . In this work, z is called the growth parameter. Using finite-size scaling analysis, we obtain that the critical exponents $$\beta /\bar{\nu }$$ β / ν ¯ and $$\gamma /\bar{\nu }$$ γ / ν ¯ associated with the magnetization and the susceptibility, respectively. Using Monte Carlo simulations, we calculate the critical noise parameter q c as a function of z for the scale-free networks and obtain the phase diagram of the model. We find that the critical exponents add up to unity when using a special volumetric scaling, regardless of the dimension of the network of interactions. We verify this result by obtaining the critical noise and the critical exponents for the two and three-state majority-vote model on cubic lattice networks.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Universal energy-speed-accuracy trade-offs in driven nonequilibrium systems

The connection between measure theoretic optimal transport and dissipative nonequilibrium dynamics provides a language for quantifying nonequilibrium control costs, leading to a collection of thermodynamic speed limits, which rely on the assumption that the target probability distribution is perfectly realized. This is almost never the case in experiments or numerical simulations, so here we address the situation in which the external controller is imperfect. We obtain a lower bound for the dissipated work in generic nonequilibrium control problems that (1) is asymptotically tight and (2) matches the thermodynamic speed limit in the case of optimal driving. Along with analytically solvable examples, we refine this imperfect driving notion to systems in which the controlled degrees of freedom are slow relative to the nonequilibrium relaxation rate, and identify independent energy contributions from fast and slow degrees of freedom. Furthermore, we develop a strategy for optimizing minimally dissipative protocols based on optimal transport flow matching, a generative machine learning technique. Furthermore, this latter approach ensures the scalability of both the theoretical and computational framework we put forth. Crucially, we demonstrate that we can compute the terms in our bound numerically using efficient algorithms from the computational optimal transport literature and that the protocols we learn saturate the bound.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Top Loading Helium Cryostat Integrated with High-Pressure Cell with Fast Remote Pressure Control (Final Report for SBIR Phase 1)

A large number of physical phenomena, such as superconductivity and quantum critical phenomena, often appear only at very low temperatures below 5 K. There is an immense interest to investigate these phenomena at high pressure as a means of tuning interatomic distances, and thus the interaction parameters controlling these phenomena, in a continuous and controlled fashion. The current P-T condition for neutron scattering experiments are limited to either relatively low pressures of about 2 GPa at temperatures below 5 K, or to relatively high temperatures at pressures or tens of GPa. The purpose of this research is to develop neutron sample environment instrumentation for reaching 10-20 GPa at 1-2K with rapid and reliable online pressure and temperature control (i.e. without having to interrupt the experiment). The ultimate goal of the project is to design an integrated fast-cooling low-temperature sample environment cryogenic system compatible with state of the art neutron diamond anvil cells suitable for single-crystal neutron scattering experiments for temperatures down to 2K and pressures of several tens of GPa. The integrated system will consist of top-loading Helium flow cryostat with in-situ sample alignment mechanisms, large-volume diamond anvil cells (DAC) made from novel superalloy Pascalloy and optimized for fast cooling and heating, and a compact remote pressure control mechanism for the DAC based on a novel concept of inflatable bellows integrated with a lever-arm based force amplifier. In Phase I, we have designed, manufactured, and tested prototypes of the novel compact force-amplified pneumatic pressure control mechanism for Neutron Diamond Anvil Cell (nDAC), which allows to use pneumatic bellows system for smooth remote pressure control in the nDAC inside a top-loading cryostat with bore size of 70 mm or larger. This allows significant minimum temperature decrease in remotely controlled nDAC from 5-10 K down to 2 K. We also studied mechanical properties of a novel non-magnetic superalloy Pascalloy with different heat-treatment conditions and preliminary results indicate that up to date this is probably the strongest and most suitable material for making more compact and lightweight cryogenic nDACs for neuron scattering experiments at extreme conditions, allowing much faster cooling and heating. Potential applications: The new developments will allow to create a wide range of compact Diamond Anvil Cells for neutron diffraction which can reach several tens of GPA pressure at 2-4 K while preserving accurate remote pressure control capabilities. The new concepts can be used for developing sample environment instrumentation outside the neutron scattering field. The new development have very strong potential for expanding experimental capabilities in materials sciences and have high commercialization potential. The proposed new cryogenic high-pressure system or any of its individual components will be in high demand not only in neutron scattering facilities, but also at synchrotron beamlines and other high-pressure research facilities around the world.

36 MATERIALS SCIENCE↗
DOE OSTI2021

A Comprehensive Analysis of NGC 2158 in the Gaia Era: Photometric Parameters, Apex, and Orbit

We present an investigation of NGC 2158 using Gaia DR2 data. We identified 800 most likely cluster members with membership probability higher than 90%. The mean proper motions of this object are determined as (μ {sub x} = −0.203 ± 0.003, μ {sub y} = −1.99 ± 0.004) mas yr{sup −1}. The limiting radius, log(age), and distance of the cluster are obtained as 23.5 arcmin, 9.38 ± 0.04 Gyr, and 4.69 ± 0.22 kpc, respectively. The overall mass function slope (0.93 ± 0.14) is flatter than the Salpeter value (1.35) within the mass range of 1.17–1.44 M {sub ⊙}. This cluster also shows the mass-segregation effect and our study demonstrates that NGC 2158 is a dynamically relaxed open cluster. Using the apex diagram (AD), the apex coordinates of the cluster are obtained in different ways and are examined using the (μ {sub U},μ {sub T}) diagram. The best value of the apex coordinates is determined as A = 8724 ± 160, D = −3661 ± 530. We also determined the orbit of the cluster and found that NGC 2158 moves almost in the solar antapex direction. The resulting spatial velocity of NGC 2158 is 51 km s{sup −1}. A significant oscillation along the Z-coordinate up to 529 pc is detected. Various scenarios regarding the origin of this cluster are also discussed.

42 ENGINEERING↗
DOE OSTI2021

Adversarial sampling of unknown and high-dimensional conditional distributions

Many engineering problems require the prediction of realization-to-realization variability or a refined description of modeled quantities. In that case, it is necessary to sample elements from unknown high-dimensional spaces with possibly millions of degrees of freedom. While there exist methods able to sample elements from probability density functions (PDF) with known shapes, several approximations need to be made when the distribution is unknown. In this paper the sampling method, as well as the inference of the underlying distribution, are both handled with a data-driven method known as generative adversarial networks (GAN), which trains two competing neural networks to produce a network that can effectively generate samples from the training set distribution. In practice, it is often necessary to draw samples from conditional distributions. When the conditional variables are continuous, only one (if any) data point corresponding to a particular value of a conditioning variable may be available, which is not sufficient to estimate the conditional distribution. This work handles this problem using an a priori estimation of the conditional moments of a PDF. Herein, two approaches, stochastic estimation, and an external neural network are compared for computing these moments; however, any preferred method can be used. The algorithm is demonstrated in the case of the deconvolution of a filtered turbulent flow field. It is shown that all the versions of the proposed algorithm effectively sample the target conditional distribution with minimal impact on the quality of the samples compared to state-of-the-art methods. Additionally, the procedure can be used as a metric for the diversity of samples generated by a conditional GAN (cGAN) conditioned with continuous variables.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Spin Change of Asteroid 2012 TC4 Probably by Radiation Torques

Asteroid 2012 TC4 is a small (∼10 m) near-Earth object that was observed during its Earth close approaches in 2012 and 2017. Earlier analyses of light curves revealed its excited rotation state. We collected all available photometric data from the two apparitions to reconstruct its rotation state and convex shape model. We show that light curves from 2012 and 2017 cannot be fitted with a single set of model parameters; the rotation and precession periods are significantly different for these two data sets, and they must have changed between or during the two apparitions. Nevertheless, we could fit all light curves with a dynamically self-consistent model assuming that the spin states of 2012 TC4 in 2012 and 2017 were different. To interpret our results, we developed a numerical model of its spin evolution in which we included two potentially relevant perturbations: (i) gravitational torque due to the Sun and Earth and (ii) radiation torque, known as the Yarkovsky–O’Keefe–Radzievskii–Paddack (YORP) effect. Despite our model simplicity, we found that the role of gravitational torques is negligible. Instead, we argue that the observed change of its spin state may be plausibly explained as a result of the YORP torque. To strengthen this interpretation, we verify that (i) the internal energy dissipation due to material inelasticity and (ii) an impact with a sufficiently large interplanetary particle are both highly unlikely causes of its observed spin state change. If true, this is the first case where the YORP effect has been detected for a tumbling body.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗