Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “quantum statistical mechanics”

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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Classical discrete time crystals

The spontaneous breaking of time-translation symmetry in periodically driven quantum systems leads to a new phase of matter: the discrete time crystal (DTC). This phase exhibits collective subharmonic oscillations that depend upon an interplay of non-equilibrium driving, many-body interactions and the breakdown of ergodicity. However, subharmonic responses are also a well-known feature of classical dynamical systems ranging from predator–prey models to Faraday waves and a.c.-driven charge density waves. This raises the question of whether these classical phenomena display the same rigidity characteristic of a quantum DTC. In this work, we explore this question in the context of periodically driven Hamiltonian dynamics coupled to a finite-temperature bath, which provides both friction and, crucially, noise. Focusing on one-dimensional chains, where in equilibrium any transition would be forbidden at finite temperature, we provide evidence that the combination of noise and interactions drives a sharp, first-order dynamical phase transition between a discrete time-translation invariant phase and an activated classical discrete time crystal (CDTC) in which time-translation symmetry is broken out to exponentially long timescales. Power-law correlations are present along a first-order line, which terminates at a critical point. We analyse the transition by mapping it to the locked-to-sliding transition of a d.c.-driven charge density wave. To conclude, building upon results from the field of probabilistic cellular automata, we conjecture the existence of classical time crystals with true long-range order, where time-translation symmetry is broken out to infinite times.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Sample-Efficient Adaptive Calibration of Quantum Networks Using Bayesian Optimization

All physical systems employed for quantum information tasks must act as unbiased carriers of encoded quantum states. Ensuring such indistinguishability of information carriers is a major challenge in many quantum information applications, including advanced quantum communication protocols. For photons, the workhorses of quantum communication networks, it is difficult to obtain and maintain their indistinguishability because of environment-induced transformations and loss imparted by communication channels, especially in noisy scenarios. Conventional strategies to mitigate these transformations often require hardware or software overhead that is restrictive (e.g., adding noise), infeasible (e.g., on a satellite), or time-consuming for deployed networks. In this work we propose and develop resource-efficient Bayesian optimization techniques to rapidly and adaptively calibrate the indistinguishability of individual photons for quantum networks using only information derived from their measurement. To experimentally validate our approach, we demonstrate the optimization of Hong-Ou-Mandel interference between two photons-a central task in quantum networking- finding rapid, efficient, and reliable convergence towards maximal photon indistinguishability in the presence of high loss and shot noise. We expect our resource-optimized and experimentally friendly methodology will allow fast and reliable calibration of indistinguishable quanta, a necessary task in distributed quantum computing, communications, and sensing, as well as for fundamental investigations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Does a Free Electron Laser Exhibit Non-Standard Statistics? [Slides and Report]

The originally proposed work under this award sought to develop a theoretical basis for previous experimental observations in free-electron laser (FEL) physics involving certain non-classical effects in the light produced by these devices. A sound and valid theoretical interpretation of these experiments would represent a transformational understanding of free-electron lasers, and of electron-photon interactions in general. The specific objective of the completed work, as described in the attached progress report by Jeongwan Park, was to perform a detailed analysis of an experiment by Chen and Madey (Physical Review Letters; volume 86 number 26, 2001) concerning the purported sub-Poissonian photon statistics in an FEL. Since that published work contradicted FEL theory as presently understood, the purpose of the work under the present subcontract was to obtain a better understanding of the quantum nature of the FEL, to investigate the possibility of theoretical and experimental evidence for the non-standard photon statistics of the FEL, and to critically examine the experimental evidence through a re-evaluation of the data analysis of Chen and Madey. The results of this study showed that there were numerous experimental conditions under which the observed photon statistics could be explained by combining both the photon clustering property and the dead-time effect, and consequently, that the observed photon statistics in the Chen-Madey experiment did not require sub-Poissonian photon statistics for their explanation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Determination of the critical manifold tangent space and curvature with Monte Carlo renormalization group

We show that the critical manifold of a statistical mechanical system in the vicinity of a critical point is locally accessible through correlation functions at that point. A practical numerical method is presented to determine the tangent space and the curvature to the critical manifold with variational Monte Carlo renormalization group. Because of the use of a variational bias potential of the coarse-grained variables, critical slowing down is greatly alleviated in the Monte Carlo simulation. In addition, this method is free of truncation error. Furthermore, we study the isotropic Ising model on square and cubic lattices, the anisotropic Ising model, and the tricritical Ising model on square lattices to illustrate the method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1977

Quantum statistical theory of semiconductor junctions in thermal equilibrium

Free carrier and electric field distributions of one-dimensional semiconductor junctions are evaluated using a quantum mechanical phase-space distribution and its corresponding Boltzmann equation. Attention is given to quantum and exchange corrections in cases of high doping concentrations when carrier densities become degenerate. Quantitative differences between degenerate and classical junction characteristics, e.g., maximum electric field and built-in voltage and carrier concentration within the transition region, are evaluated numerically.

Von Roos, O.↗
DOE OSTI2020

Equations of State in Computational Physics: What, Why, and How

Summary: 1. Equations of state furnish thermodynamic properties of a material. 2. Usually involve P, p (or V), and T, or P, r, and E 3. Necessary ingredient to solve compressible flow problems. 4. Varieties: i. Analytical expressions (ideal gas, Grüneisen) – fast, convenient for verification and test problems ii. Graphs and charts – good for hand calculations iii. Tables and databases – generally most accurate, require interpolation i. Steam tables ii. NIST database iii. SESAME EOS Library 5. Constructed by using semi-empirical models from Statistical Mechanics and fitting to available data

36 MATERIALS SCIENCE↗
DOE OSTI2022

Algebraic Bethe Circuits

The Algebraic Bethe Ansatz (ABA) is a highly successful analytical method used to exactly solve several physical models in both statistical mechanics and condensed-matter physics. Here we bring the ABA into unitary form, for its direct implementation on a quantum computer. This is achieved by distilling the non-unitary R matrices that make up the ABA into unitaries using the QR decomposition. Our algorithm is deterministic and works for both real and complex roots of the Bethe equations. We illustrate our method on the spin-$\frac{1}{2}$ XX and XXZ models. We show that using this approach one can efficiently prepare eigenstates of the XX model on a quantum computer with quantum resources that match previous state-of-the-art approaches. We run small-scale error mitigated implementations on the IBM quantum computers, including the preparation of the ground state for the XX and XXZ models on 4 sites. Finally, we derive a new form of the Yang-Baxter equation using unitary matrices, and also verify it on a quantum computer.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Obstacles to quantum annealing in a planar embedding of XORSAT

We introduce a planar embedding of the k-regular k-XORSAT problem, in which solutions are encoded in the ground state of a classical statistical mechanics model of reversible logic gates arranged on a square grid and acting on bits that represent the Boolean variables of the problem. The special feature of this embedding is that the resulting model lacks a finite-temperature phase transition, thus bypassing the first-order thermodynamic transition known to occur in the random graph representation of XORSAT. In spite of this attractive feature, the thermal relaxation into the ground state displays remarkably slow glassy behavior. The question addressed in this paper is whether this planar embedding can afford an efficient path to solution of k-regular k-XORSAT via quantum adiabatic annealing. We first show that our model bypasses an avoided level crossing and consequent exponentially small gap in the limit of small transverse fields. We then present quantum Monte Carlo results for our embedding of the k-regular k-XORSAT that strongly support a picture in which second-order and first-order transitions develop at a finite transverse field for k = 2 and k = 3, respectively. Furthermore, this translates into power-law and exponential dependences in the scaling of energy gaps with system size, corresponding to times-to-solution which are, respectively, polynomial and exponential in the number of variables. We conclude that neither classical nor quantum annealing can efficiently solve our reformulation of XORSAT, even though the original problem can be solved in polynomial time by Gaussian elimination.

2-dimensional systems↗
DOE OSTI2020

Normal modes, rotational inertia, and thermal fluctuations of trapped ion crystals

The normal modes of a trapped ion crystal are derived using an approach based on the Hermitian properties of the system’s dynamical matrix. This method is equivalent to the standard Bogoliubov method, but for classical systems it is arguably simpler and more general in that canonical coordinates are not necessary. The theory is developed for stable, unstable, and neutrally-stable systems. The method is then applied to ion crystals in a Penning trap. Reduced eigenvalue problems for the case of large applied magnetic field are developed, for which the spectrum breaks into E x B drift modes, axial modes, and cyclotron modes. Thermal fluctuation levels in these modes are analyzed and shown to be consistent with the Bohr-van-Leeuwen theorem, provided that neutrally-stable modes associated with crystal rotations are included in the analysis. An expression for the rotational inertia of the crystal is derived, and a magnetic contribution to this inertia, which dominates in large magnetic fields, is described. An unusual limit is discovered for the special case of spherically-symmetric confinement, in which the rotational inertia does not exist and changes in angular momentum leave the rotation frequency unaffected.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Calculating adsorption isotherms using the two-phase thermodynamic method and molecular dynamics simulations

We describe the calculation of adsorption isotherms from molecular dynamics simulations based on the two-phase thermodynamic (2PT) model. The 2PT model developed for bulk fluid phases treats the gas-like components as hard spheres (HSs), which correctly recovers the limiting behaviors of unconfined fluids. We showed that this treatment, however, does not always lead to the correct zero-loading behavior in strongly confining systems. For methane adsorption into zeolite MFI, the HS reference state underestimates entropy by up to 20% at low loadings and leads to an order-of-magnitude increase in the adsorption onset pressure. To fix these issues, we propose the use of ideal adsorbed gas (IAG) as the gas reference model, the properties of which can be computed using the Widom insertion method on an empty adsorbent. We further describe three routes to compute adsorption isotherms from the Helmholtz free energy at different loadings. Comparing against established Monte Carlo (MC) methods, we found that the adsorption isotherms obtained using the IAG reference state agrees to within 40%, which corresponds to deviations of <5% in adsorption free energy. The isotherms calculated using the HS reference state underestimate the adsorption uptake at low to medium loadings in strongly confining systems, but its accuracy improves at higher loadings and as the pore size increases relative to the sorbate diameter. The methods described here provide an alternative approach for computing adsorption isotherms when MC simulations in an open ensemble are undesirable and enable a direct comparison of computed adsorption thermodynamics with experiments.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

A dissipation-induced superradiant transition in a strontium cavity-QED system

Driven-dissipative many-body systems are ubiquitous in nature and a fundamental resource for quantum technologies. However, they are also complex and hard to model because they cannot be described by the standard tools in equilibrium statistical mechanics. Probing nonequilibrium critical phenomena in pristine setups can illuminate fresh perspectives on these systems. Here, we use an ensemble of cold 88 Sr atoms coupled to a driven high-finesse cavity to study the cooperative resonance fluorescence (CRF) model, a classic driven-dissipative model describing coherently driven dipoles superradiantly emitting light. We observe its nonequilibrium phase diagram characterized by a second-order phase transition. Below a critical drive strength, the atoms quickly reach the so-called superradiant steady state featuring a macroscopic dipole moment; above the critical point, the atoms undergo persistent Rabi-like oscillations. At longer times, spontaneous emission transforms the second-order transition into a discontinuous first-order transition. Our observations pave the way for harnessing robust entangled states and exploring boundary time crystals in driven-dissipative systems.

Science & Technology - Other Topics↗
DOE OSTI2020

Statistical properties of homogeneous and isotropic turbulence in He II measured via particle tracking velocimetry

Despite being a quantum two-fluid system, superfluid helium-4 (He II) is observed to behave similarly to classical fluids when a flow is generated by mechanical forcing. This similarity has brought up the feasibility of utilizing He II for high Reynolds number classical turbulence research, considering the small kinematic viscosity of He II. However, it has been suggested that the nonclassical dissipation mechanism in He II at small scales may alter its turbulent statistics and intermittency. In this work, we report our study of a nearly homogeneous and isotropic turbulence (HIT) generated by a towed grid in He II. We measure the velocity field using particle tracking velocimetry with solidified deuterium particles as the tracers. By correlating the velocities measured simultaneously on different particle trajectories or at different times along the same particle trajectory, we are able to conduct both Eulerian and Lagrangian flow analyses. Spatial velocity structure functions obtained through the Eulerian analysis show scaling behaviors in the inertial subrange similar to that for classical HIT but with enhanced intermittency. The Lagrangian analysis allows us to examine the flow statistics down to below the dissipation length scale. Interestingly, abnormal deviations from the classical scaling behaviors are observed in this regime. Lastly, we discuss how these deviations may relate to the motion of quantized vortices in the superfluid component in He II.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Measurement-induced criticality in random quantum circuits

We investigate the critical behavior of the entanglement transition induced by projective measurements in (Haar) random unitary quantum circuits. Using a replica approach, we map the calculation of the entanglement entropies in such circuits onto a two-dimensional statistical-mechanics model. In this language, the area- to volume-law entanglement transition can be interpreted as an ordering transition in the statistical-mechanics model. We derive the general scaling properties of the entanglement entropies and mutual information near the transition using conformal invariance. We analyze in detail the limit of infinite on-site Hilbert space dimension in which the statistical-mechanics model maps onto percolation. In particular, we compute the exact value of the universal coefficient of the logarithm of subsystem size in the n th Rényi entropies for n ≥ 1 in this limit using relatively recent results for the conformal field theory describing the critical theory of two-dimensional (2D) percolation, and we discuss how to access the generic transition at finite on-site Hilbert space dimension from this limit, which is in a universality class different from 2D percolation. Here, we also comment on the relation to the entanglement transition in random tensor networks, studied previously in Vasseur et al. [Phys. Rev. B 100, 134203 (2019)].

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Disentangling the impact of quasiparticles and two-level systems on the statistics of superconducting-qubit lifetime

Temporal fluctuations in the superconducting qubit lifetime, T 1 , present additional challenges in the pursuit of fault-tolerant quantum computing. Although the exact mechanisms remain unclear, T 1 fluctuations are generally attributed to strong coupling between the qubit and a few near-resonant two-level systems (TLSs), which can exchange energy with an ensemble of thermally fluctuating two-level fluctuators (TLFs) at low frequencies. Here, we report T 1 measurements of qubits with varying geometrical footprints and surface dielectrics as a function of temperature. By analyzing the noise spectrum of the qubit depolarization rate, Γ 1 = 1 / T 1 , we disentangle the contributions of TLSs, nonequilibrium quasiparticles (QPs), and equilibrium (thermally excited) QPs to the variance in Γ 1 . We find that the Γ 1 variance in qubits with smaller footprints is more susceptible to QP and TLS fluctuations than that in larger-footprint qubits. Furthermore, the QP-induced variances in all qubits align with the theoretical framework of QP diffusion and fluctuation. These findings offer valuable insights for future qubit design and engineering optimization.

Zhu, Shaojiang [Fermilab] (ORCID:0000000293180092)↗
DOE OSTI2021

Filtering, averaging, and scale dependency in homogeneous variable density turbulence

We investigate relationships between statistics obtained from filtering and from ensemble or Reynolds-averaging turbulence flow fields as a function of length scale. Generalized central moments in the filtering approach are expressed as inner products of generalized fluctuating quantities, q ' ( ξ , x ) = q ( ξ ) - q ¯ ( x ) , representing fluctuations of a field q ( ξ ) , at any point ξ, with respect to its filtered value at x. For positive-definite filter kernels, these expressions provide a scale-resolving framework, with statistics and realizability conditions at any length scale. In the small-scale limit, scale-resolving statistics become zero. In the large-scale limit, scale-resolving statistics and realizability conditions are the same as in the Reynolds-averaged description. Using direct numerical simulations (DNS) of homogeneous variable density turbulence, we diagnose Reynolds stresses, T i j , resolved kinetic energy, kr, turbulent mass-flux velocity, a i , and density-specific volume covariance, b, defined in the scale-resolving framework. These variables, and terms in their governing equations, vary smoothly between zero and their Reynolds-averaged definitions at the small and large scale limits, respectively. At intermediate scales, the governing equations exhibit interactions between terms that are not active in the Reynolds-averaged limit. For example, in the Reynolds-averaged limit, b follows a decaying process driven by a destruction term; at intermediate length scales, it is a balance between production, redistribution, destruction, and transport, where b grows as the density spectrum develops and then decays when mixing becomes strong enough. This work supports the notion of a generalized, length-scale adaptive model that converges to DNS at high resolutions and to Reynolds-averaged statistics at coarse resolutions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

The metamorphosis of semi-classical mechanisms of confinement: from monopoles on ℝ 3 × S 1 to center-vortices on ℝ 2 × T 2

There are two distinct regimes of Yang-Mills theory where we can demonstrate confinement, the existence of a mass gap, and the multi-branch structure of the effective potential as a function of the theta angle using a reliable semi-classical calculation. The two regimes are deformed Yang-Mills theory on ℝ 3 × S 1 , and Yang-Mills theory on ℝ 2 × T 2 where the torus is threaded by a ’t Hooft flux. The weak coupling regime is ensured by the small size of the circle or torus. In the first case the confinement mechanism is related to self-dual monopoles, whereas in the second case self-dual center-vortices play a crucial role. These two topological objects are distinct. In particular, they have different mutual statistics with Wilson loops. On the other hand, they carry the same topological charge and action. We consider the theory on ℝ × T 2 × S 1 and extrapolate both the monopole and vortex regimes to a quantum mechanical domain, where a cross-over takes place. Both sides of the cross-over are described by a deformed ℤ N TQFT. On ℝ 2 × S 1 × S 1 , we derive an effective field theory (EFT) of vortices from the EFT of monopoles in the presence of a ’t Hooft flux. This construction is based on a two-stage Higgs mechanism, reducing SU(N) to U(1) N−1 in 3d first, followed by reduction to a ℤ N EFT in 2d in the second step. This result shows how monopoles transmute into center-vortices, and suggests adiabatic continuity between the two confinement mechanisms. The basic mechanism is flux fractionalization: the magnetic flux of the monopoles splits up and is collimated in such a way that 2d Wilson loops detect it as a center vortex.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Machine learning to accelerate screening for Marcus reorganization energies

Understanding and predicting the charge transport properties of π-conjugated materials is an important challenge for designing new organic electronic devices, such as solar cells, plastic transistors, light-emitting devices, and chemical sensors. A key component of the hopping mechanism of charge transfer in these materials is the Marcus reorganization energy which serves as an activation barrier to hole or electron transfer. While modern density functional methods have proven to accurately predict trends in intramolecular reorganization energy, such calculations are computationally expensive. In this work, we outline active machine learning methods to predict computed intramolecular reorganization energies of a wide range of polythiophenes and their use toward screening new compounds with low internal reorganization energies. Our models have an overall root mean square error (RMSE) of ±0.113 eV, but a much smaller RMSE of only ±0.036 eV on the new screening set. Since the larger error derives from high-reorganization energy compounds, the new method is highly effective to screen for compounds with potentially efficient charge transport parameters.

14 SOLAR ENERGY↗