Results for “Statistical field theory” — Page 11

NASA NTRS1975

A survey of the interplanetary magnetic field

All interplanetary magnetic field measurements from near-earth satellites have been collected in the form of hourly averages for the years 1963-1974. Solar cycle variations of the field and the exceptional events of this cycle are discussed in relation to the associated cosmic ray events, and statistical properties of the field over this period are presented. The presence of magnetic blobs as a permanent feature of the field is demonstrated and explained. The implications of these results for the theory of cosmic ray modulation are presented and analyzed.

Barouch, E.↗

NASA NTRS2014

Statistical Mechanics of Turbulent Dynamos

Incompressible magnetohydrodynamic (MHD) turbulence and magnetic dynamos, which occur in magnetofluids with large fluid and magnetic Reynolds numbers, will be discussed. When Reynolds numbers are large and energy decays slowly, the distribution of energy with respect to length scale becomes quasi-stationary and MHD turbulence can be described statistically. In the limit of infinite Reynolds numbers, viscosity and resistivity become zero and if these values are used in the MHD equations ab initio, a model system called ideal MHD turbulence results. This model system is typically confined in simple geometries with some form of homogeneous boundary conditions, allowing for velocity and magnetic field to be represented by orthogonal function expansions. One advantage to this is that the coefficients of the expansions form a set of nonlinearly interacting variables whose behavior can be described by equilibrium statistical mechanics, i.e., by a canonical ensemble theory based on the global invariants (energy, cross helicity and magnetic helicity) of ideal MHD turbulence. Another advantage is that truncated expansions provide a finite dynamical system whose time evolution can be numerically simulated to test the predictions of the associated statistical mechanics. If ensemble predictions are the same as time averages, then the system is said to be ergodic; if not, the system is nonergodic. Although it had been implicitly assumed in the early days of ideal MHD statistical theory development that these finite dynamical systems were ergodic, numerical simulations provided sufficient evidence that they were, in fact, nonergodic. Specifically, while canonical ensemble theory predicted that expansion coefficients would be (i) zero-mean random variables with (ii) energy that decreased with length scale, it was found that although (ii) was correct, (i) was not and the expected ergodicity was broken. The exact cause of this broken ergodicity was explained, after much investigation, by greatly extending the statistical theory of ideal MHD turbulence. The mathematical details of broken ergodicity, in fact, give a quantitative explanation of how coherent structure, dynamic alignment and force-free states appear in turbulent magnetofluids. The relevance of these ideal results to real MHD turbulence occurs because broken ergodicity is most manifest in the ideal case at the largest length scales and it is in these largest scales that a real magnetofluid has the least dissipation, i.e., most closely approaches the behavior of an ideal magnetofluid. Furthermore, the effects grow stronger when cross and magnetic helicities grow large with respect to energy, and this is exactly what occurs with time in a real magnetofluid, where it is called selective decay. The relevance of these results found in ideal MHD turbulence theory to the real world is that they provide at least a qualitative explanation of why confined turbulent magnetofluids, such as the liquid iron that fills the Earth's outer core, produce stationary, large-scale magnetic fields, i.e., the geomagnetic field. These results should also apply to other planets as well as to plasma confinement devices on Earth and in space, and the effects should be manifest if Reynolds numbers are high enough and there is enough time for stationarity to occur, at least approximately. In the presentation, details will be given for both theoretical and numerical results, and references will be provided.

Shebalin, John V.↗

NASA NTRS2014

MHD Turbulence and Magnetic Dynamos

Incompressible magnetohydrodynamic (MHD) turbulence and magnetic dynamos, which occur in magnetofluids with large fluid and magnetic Reynolds numbers, will be discussed. When Reynolds numbers are large and energy decays slowly, the distribution of energy with respect to length scale becomes quasi-stationary and MHD turbulence can be described statistically. In the limit of infinite Reynolds numbers, viscosity and resistivity become zero and if these values are used in the MHD equations ab initio, a model system called ideal MHD turbulence results. This model system is typically confined in simple geometries with some form of homogeneous boundary conditions, allowing for velocity and magnetic field to be represented by orthogonal function expansions. One advantage to this is that the coefficients of the expansions form a set of nonlinearly interacting variables whose behavior can be described by equilibrium statistical mechanics, i.e., by a canonical ensemble theory based on the global invariants (energy, cross helicity and magnetic helicity) of ideal MHD turbulence. Another advantage is that truncated expansions provide a finite dynamical system whose time evolution can be numerically simulated to test the predictions of the associated statistical mechanics. If ensemble predictions are the same as time averages, then the system is said to be ergodic; if not, the system is nonergodic. Although it had been implicitly assumed in the early days of ideal MHD statistical theory development that these finite dynamical systems were ergodic, numerical simulations provided sufficient evidence that they were, in fact, nonergodic. Specifically, while canonical ensemble theory predicted that expansion coefficients would be (i) zero-mean random variables with (ii) energy that decreased with length scale, it was found that although (ii) was correct, (i) was not and the expected ergodicity was broken. The exact cause of this broken ergodicity was explained, after much investigation, by greatly extending the statistical theory of ideal MHD turbulence. The mathematical details of broken ergodicity, in fact, give a quantitative explanation of how coherent structure, dynamic alignment and force-free states appear in turbulent magnetofluids. The relevance of these ideal results to real MHD turbulence occurs because broken ergodicity is most manifest in the ideal case at the largest length scales and it is in these largest scales that a real magnetofluid has the least dissipation, i.e., most closely approaches the behavior of an ideal magnetofluid. Furthermore, the effects grow stronger when cross and magnetic helicities grow large with respect to energy, and this is exactly what occurs with time in a real magnetofluid, where it is called selective decay. The relevance of these results found in ideal MHD turbulence theory to the real world is that they provide at least a qualitative explanation of why confined turbulent magnetofluids, such as the liquid iron that fills the Earth's outer core, produce stationary, large-scale magnetic fields, i.e., the geomagnetic field. These results should also apply to other planets as well as to plasma confinement devices on Earth and in space, and the effects should be manifest if Reynolds numbers are high enough and there is enough time for stationarity to occur, at least approximately. In the presentation, details will be given for both theoretical and numerical results, and references will be provided.

Shebalin, John V↗

DOE OSTI2024

Thermal Bekenstein-Hawking entropy from the worldsheet

We define and compute the leading sphere diagram contribution to the entropy of the BTZ black hole supported by Kalb-Ramond flux in bosonic string theory. In a winding condensate description, integrating exactly over the constant mode for the radial direction of AdS 3 reduces the problem to one of the correlation functions of winding operators in the free theory. The volume of the residual PSL(2,ℂ) gauge group of the sphere is canceled by the action of conformal transformations on the winding interaction insertions. We formulate a precise version of the replica trick in terms of (infinitesimally) non-integer winding condensates to produce the entropy of the BTZ black hole. The resulting entropy can be calculated from the one-point function of a non-local operator on the worldsheet.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2024

Coherent spin states and emergent de Sitter quasinormal modes

As a toy model for the microscopic description of matter in de Sitter space, we consider a Hamiltonian acting on the spin-j representation of SU(2). This is a model with a finite-dimensional Hilbert space, from which quasinormal modes emerge in the large-spin limit. The path integral over coherent spin states can be evaluated at the semiclassical level and from it we find the single-particle de Sitter density of states, including 1/j corrections. Along the way, we discuss the use of quasinormal modes in quantum mechanics, starting from the paradigmatic upside-down harmonic oscillator.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2023

The endpoint of partial deconfinement

We study the matrix quantum mechanics of two free hermitian N × N matrices subject to a singlet constraint in the microcanonical ensemble. This is the simplest example of a theory that at large N has a confinement/deconfinement transition. In the microcanonical ensemble, it also exhibits partial deconfinement with a Hagedorn density of states. We argue that the entropy of these configurations, based on a combinatorial counting of Young diagrams, are dominated by Young diagrams that have the VKLS shape. When the shape gets to the maximal depth allowed for a Young diagram of SU(N), namely N, we argue that the system stops exhibiting the Hagedorn behavior. The number of boxes (energy) at the transition is N 2 /4, independent of the charge of the state.

1/N Expansion↗

NASA NTRS1997

Foreshock Langmuir Waves for Unusually Constant Solar Wind Conditions: Data and Implications for Foreshock Structure

Plasma wave data are compared with ISEE 1's position in the electron foreshock for an interval with unusually constant (but otherwise typical) solar wind magnetic field and plasma characteristics. For this period, temporal variations in the wave characteristics can be confidently separated from sweeping of the spatially varying foreshock back and forth across the spacecraft. The spacecraft's location, particularly the coordinate D(sub f) downstream from the foreshock boundary (often termed DIFF), is calculated by using three shock models and the observed solar wind magnetometer and plasma data. Scatterplots of the wave field versus D(sub f) are used to constrain viable shock models, to investigate the observed scatter in the wave fields at constant D(sub f), and to test the theoretical predictions of linear instability theory. The scatterplots confirm the abrupt onset of the foreshock waves near the upstream boundary, the narrow width in D(sub f) of the region with high fields, and the relatively slow falloff of the fields at large D(sub f), as seen in earlier studies, but with much smaller statistical scatter. The plots also show an offset of the high-field region from the foreshock boundary. It is shown that an adaptive, time-varying shock model with no free parameters, determined by the observed solar wind data and published shock crossings, is viable but that two alternative models are not. Foreshock wave studies can therefore remotely constrain the bow shock's location. The observed scatter in wave field at constant D(sub f) is shown to be real and to correspond to real temporal variations, not to unresolved changes in D(sub f). By comparing the wave data with a linear instability theory based on a published model for the electron beam it is found that the theory can account qualitatively and semiquantitatively for the abrupt onset of the waves near D(sub f) = 0, for the narrow width and offset of the high-field region, and for the decrease in wave intensity with increasing D(sub f). Quantitative differences between observations and theory remain, including large overprediction of the wave fields and the slower than predicted falloff at large D(sub f) of the wave fields. These differences, as well as the unresolved issue of the electron beam speed in the high-field region of the foreshock, are discussed. The intrinsic temporal variability of the wave fields, as well as their overprediction based on homogeneous plasma theory, are indicative of stochastic growth physics, which causes wave growth to be random and varying in sign, rather than secular.

Cairns, Iver H.↗

DOE OSTI2020

Chirality transfer and chiral turbulence in gauge theories

Chirality transfer between fermions and gauge fields plays a crucial role for understanding the dynamics of anomalous transport phenomena such as the Chiral Magnetic Effect. In this proceeding we present a first principles study of these processes based on classical-statistical real-time lattice simulations of strongly coupled QED (e 2 N f = 64). Furthermore, our simulations demonstrate that a chirality imbalance in the fermion sector triggers chiral plasma instabilities in the gaugefield sector, which ultimately lead to the generation of long range helical magnetic fields via an self-similar turbulent cascade of the magnetic helicity.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2025

Production of neutron-rich heavy nuclei in deep-inelastic 208 Pb + 208 Pb collisions within the stochastic mean-field theory

In deep-inelastic collisions of heavy nuclei, reaction products with a wide range of mass and charge are produced. Such collisions been considered as a possible way to produce superheavy nuclei, as an alternative to fusion reactions. To provide reliable theoretical predictions, it is desired to develop microscopic approaches that correctly and accurately describe nucleon transfer processes in dissipative collisions of heavy nuclei. The purpose of the present work is (1) to investigate the mechanism of nucleon transfers in dissipative collisions of two heavy nuclei, and (2) to explore possible pathways to produce neutron-rich heavy nuclei, through detailed theoretical analyses of fluctuations and correlations in nucleon transfers in 208 Pb + 208 Pb reactions. Three-dimensional time-dependent Hartree-Fock (TDHF) calculations are performed for the collisions of 208 Pb + 208 Pb at 𝐸 c.m. = 832, 936, and 1040 MeV, using the Skyrme SLy4d energy density functional. To calculate fluctuations and correlations in nucleon transfers, we employ the stochastic mean-field (SMF) theory, and the results are compared with another theoretical framework currently available, the time-dependent random phase approximation (TDRPA). Primary and secondary production cross sections are calculated with the SMF theory combined with a statistical model, GEMINI ++ . Using information of nucleon flow across a neck of colliding nuclei in TDHF calculations, we solve quantal diffusion equations for fluctuations and correlations in nucleon transfers based on the SMF theory. From the SMF calculations, we obtain the time evolution of diffusion coefficients as well as fluctuations and correlations in nucleon transfers for a range of initial orbital angular momenta. We compare the results of the SMF calculations with those of TDRPA, showing that TDRPA tends to predict substantially larger fluctuations and correlations in strongly damped collisions of heavy nuclei, which exhibit complex initial angular momentum dependence, while the SMF results provide almost constant (stable) values. Using the obtained fluctuations and correlations, we calculate primary and secondary production cross sections for the 208 Pb + 208 Pb collisions. From the results, we find that both lighter and heavier reaction products as compared to 208 Pb are produced for a wide region in the 𝑁−𝑍 plane as primary products, thanks to the quantal diffusion mechanism in the dissipative collisions. However, we show that cross sections for production of heavy nuclei with 𝑍 ≳ 90 or 𝑁 ≳ 135 are washed out due to secondary particle evaporation and/or fission processes. On the other hand, we find that there remain sizable cross sections for production of neutron-rich nuclei along 𝑁 = 126 with 𝑍< 82, even after secondary disintegration processes. We demonstrate that the secondary production cross sections depend weakly on incident energies, but lower (higher) energy is slightly preferred for production of nuclei with smaller (larger) atomic numbers as compared to 𝑍 = 82. Here, based on the microscopic SMF calculations, it has been shown that deep-inelastic collisions of heavy nuclei, such as 208 Pb + 208 Pb examined in this study, can be a promising means to produce neutron-rich heavy nuclei along 𝑁=126. Discrepancies between the SMF and TDRPA approaches are left unsolved for future investigations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2022

Fractional vortices, Z 2 gauge theory, and the confinement-deconfinement transition

Here, in this Letter we discuss the classical three-dimensional XY model whose nearest-neighbor interaction is a mixture of cos (θ i -θ j ) (ferromagnetic) and cos 2 (θ i -θ j ) (nematic). This model is dual to a theory with integer and half-integer vortices. While both types of vortices interact with a noncompact U(1) gauge field, the half-integer vortices interact with an extra interaction mediated by a Z 2 gauge field. We shall discuss the confinement-deconfinement transition of the half-integer vortices, the Wilson and the 't Hooft loops, and their mutual statistics in path integral language. In addition, we shall present a quantum version of the classical model which exhibits these physics.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

NASA NTRS1977

Computer simulation of the velocity diffusion of cosmic rays

Monte Carlo simulation experiments were performed in order to study the velocity diffusion of charged particles in a static turbulent magnetic field. By following orbits of particles moving in a large ensemble of random magnetic field realizations with suitable chosen statistical properties, a pitch-angle diffusion coefficient is derived. Results are presented for a variety of particle rigidities and rms random field strengths and compared with the predictions of standard quasi-linear theory and the nonlinear partially averaged field theory.

Kaiser, T. B.↗

NASA NTRS1978

Computer simulation of the velocity diffusion of cosmic rays

Monte Carlo simulation experiments have been performed in order to study the velocity diffusion of charged particles in a static turbulent magnetic field. By following orbits of particles moving in a large ensemble of random magnetic field realizations with suitably chosen statistical properties, a pitch-angle diffusion coefficient is derived. Results are presented for a variety of particle rigidities and rms random field strengths and compared with the predictions of standard quasi-linear theory and the nonlinear partially averaged field theory.

Kaiser, T. B.↗

NASA NTRS1991

Neutron stars and the distance to gamma-ray bursters

Assuming that gamma-ray bursts originate from galactic neutron stars, an analytic method for studying their statistical properties is outlined. If a significant fraction of all neutron stars are born with space velocities of less than approximately 100 km/s, as suggested by studies of pulsar statistics, then the sampling distance to gamma-ray burst sources should be less than about several hundred pc. These results have important implications on theories of radio-pulsar evolution and magnetic-field decay.

Dermer, Charles D.↗

NASA NTRS2012

Statistical Theory of the Ideal MHD Geodynamo

A statistical theory of geodynamo action is developed, using a mathematical model of the geodynamo as a rotating outer core containing an ideal (i.e., no dissipation), incompressible, turbulent, convecting magnetofluid. On the concentric inner and outer spherical bounding surfaces the normal components of the velocity, magnetic field, vorticity and electric current are zero, as is the temperature fluctuation. This allows the use of a set of Galerkin expansion functions that are common to both velocity and magnetic field, as well as vorticity, current and the temperature fluctuation. The resulting dynamical system, based on the Boussinesq form of the magnetohydrodynamic (MHD) equations, represents MHD turbulence in a spherical domain. These basic equations (minus the temperature equation) and boundary conditions have been used previously in numerical simulations of forced, decaying MHD turbulence inside a sphere [1,2]. Here, the ideal case is studied through statistical analysis and leads to a prediction that an ideal coherent structure will be found in the form of a large-scale quasistationary magnetic field that results from broken ergodicity, an effect that has been previously studied both analytically and numerically for homogeneous MHD turbulence [3,4]. The axial dipole component becomes prominent when there is a relatively large magnetic helicity (proportional to the global correlation of magnetic vector potential and magnetic field) and a stationary, nonzero cross helicity (proportional to the global correlation of velocity and magnetic field). The expected angle of the dipole moment vector with respect to the rotation axis is found to decrease to a minimum as the average cross helicity increases for a fixed value of magnetic helicity and then to increase again when average cross helicity approaches its maximum possible value. Only a relatively small value of cross helicity is needed to produce a dipole moment vector that is aligned at approx.10deg with the rotation axis.

Shebalin, J. V.↗

DOE OSTI2024

Machine Learned Force Field Modeling of Metal Organic Frameworks for CO2 Direct Air Capture

Metal organic frameworks (MOFs) are a large class of porous materials and have garnered significant interest due to their large surface areas and their tunable physical and chemical properties. Numerous prior studies have been performed to screen large databases of this material class for promising DAC sorbent materials. These studies have often relied on classical model potentials. While density functional theory (DFT) calculations have been shown to be very accurate for modeling the interaction of CO2 with MOFs, such calculations are too computationally demanding for statistically significant adsorption predictions. To overcome this barrier, we developed methods for training models to achieve DFT-level accuracy for the forces and energies associated with MOF flexibility and CO2 adsorption using machine learned force fields (MLFFs). These methods were parametrized based on DFT calculations of CO2 in a flexible MOF and used to predict MOF structural properties as well as CO2 adsorption in several MOFs.

Findley, John↗

DOE OSTI2021

Quantum-classical hybrid algorithm for the simulation of all-electron correlation

While chemical systems containing hundreds to thousands of electrons remain beyond the reach of quantum devices, hybrid quantum-classical algorithms present a promising pathway toward a quantum advantage. Hybrid algorithms treat the exponentially scaling part of the calculation-the static correlation-on the quantum computer and the non-exponentially scaling part-the dynamic correlation-on the classical computer. While a variety of algorithms have been proposed, the dependence of many methods on the total wave function limits the development of easy-to-use classical post-processing implementations. Here, we present a novel combination of quantum and classical algorithms, which computes the all-electron energy of a strongly correlated molecular system on the classical computer from the 2-electron reduced density matrix (2-RDM) evaluated on the quantum device. Significantly, we circumvent the wave function in the all-electron calculations by using density matrix methods that only require input of the statically correlated 2-RDM. Although the algorithm is completely general, we test it with two classical density matrix methods, the anti-Hermitian contracted Schrödinger equation (ACSE) and multiconfiguration pair-density functional theories, using the recently developed quantum ACSE method for simulating the statically correlated 2-RDM. Furthermore, we obtain experimental accuracy for the relative energies of all three benzyne isomers and thereby demonstrate the ability of the developed algorithm to achieve chemically relevant and accurate results on noisy intermediate-scale quantum devices.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS1979

Optical communication with two-photon coherent states. II - Photoemissive detection and structured receiver performance

In a previous paper (1978), the authors developed a method of analyzing the performance of two-photon coherent state (TCS) systems for free-space optical communications. General theorems permitting application of classical point process results to detection and estimation of signals in arbitrary quantum states were derived. The present paper examines the general problem of photoemissive detection statistics. On the basis of the photocounting theory of Kelley and Kleiner (1964) it is shown that for arbitrary pure state illumination, the resulting photocurrent is in general a self-exciting point process. The photocount statistics for first-order coherent fields reduce to those of a special class of Markov birth processes, which the authors term single-mode birth processes. These general results are applied to the structure of TCS radiation, and it is shown that the use of TCS radiation with direct or heterodyne detection results in minimal performance increments over comparable coherent-state systems. However, significant performance advantages are offered by use of TCS radiation with homodyne detection. The abstract quantum descriptions of homodyne and heterodyne detection are derived and a synthesis procedure for obtaining quantum measurements described by arbitrary TCS is given.

Shapiro, J. H.↗

NASA NTRS1980

Guiding center plasma with gravitational or gradient drifts

It is noted that since its introduction in 1971, the subject of the statistical dynamics of the electrostatic guiding center plasma model has received considerable theoretical attention. The paper presents the theory and a simulation for the two-dimensional electrostatic guiding center plasma in a uniform gravitational field. Finally, it is shown that the gravitational field leads to large electrostatic energies at long wavelengths and associated vortex motion in the plasma.

Joyce, G.↗