Search NASA⌕ Search

SEARCH · Search NASA

Results for “set of small ordered structures”

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 OSTI2013

Neutron irradiation effects on domain wall mobility and reversibility in lead zirconate titanate thin films

The effects of neutron-induced damage on the ferroelectric properties of thin film lead zirconate titanate (PZT) were investigated. Two sets of PbZr{sub 0.52}Ti{sub 0.48}O{sub 3} films of varying initial quality were irradiated in a research nuclear reactor up to a maximum 1 MeV equivalent neutron fluence of (5.16 {+-} 0.03) Multiplication-Sign 10{sup 15} cm{sup -2}. Changes in domain wall mobility and reversibility were characterized by polarization-electric field measurements, Rayleigh analysis, and analysis of first order reversal curves (FORC). With increasing fluence, extrinsic contributions to the small-signal permittivity diminished. Additionally, redistribution of irreversible hysterons towards higher coercive fields was observed accompanied by the formation of a secondary hysteron peak following exposure to high fluence levels. The changes are attributed to the radiation-induced formation of defect dipoles and other charged defects, which serve as effective domain wall pinning sites. Differences in damage accumulation rates with initial film quality were observed between the film sets suggesting a dominance of pre-irradiation microstructure on changes in macroscopic switching behavior.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Development of neural network force fields for corrosion studies

To fully understand the chemistry and physics of corrosion, novel methods of simulation must be developed. One approach is designing machine learning (ML) algorithms integrated with density functional theory to develop adaptive force fields to gain insight into corrosion behavior namely at the surface of metal oxides. Current methods of modeling corrosion are slow due to the computational cost of resolving both reaction mechanics and mass transport processes. Machine learning methods can be implemented to obtain structure-activity relationships at both the molecular and bulk scale while still retaining the accuracy of density functional theory (DFT) and significantly decreasing the time needed for simulations of complex chemical processes in the various environments of corrosion. Multiscale models are needed for corrosion studies to fully understand its processes not only at the atomic length scale (chemical bonding, energies, and forces), but also at the nano and meso length scales (solid-state physics and material science processes). Current methods of study include DFT, molecular dynamics, and Monte Carlo. The limitation of DFT is that only a small number of atoms or molecules can be simulated at that level of theory. Density functional theory is used to study the electronic structure of atoms and molecules, and calculate the force component of each atom. However, these calculations are limited to about 1000 atoms. Custom periodic boundary conditions (PBC) can be used to describe the various environments and defects that affect the atomic forces to produce a large data set from which a training set can be derived. Machine learning can be utilized to overcome the barrier of modeling macroscopic and multi-scale processes from ab initio calculations through the development of adaptive force fields. Local environments determine the atomic forces of a given system, therefore adaptive force fields must be created to produce reliable quantum mechanical calculations. This can be achieved by developing a learning algorithm that uses the mapped atomic forces or fingerprint as an input to produce energies and magnetic moments as output. A systematic approach was used to begin to build a data set in order to accurately describe the atomic forces in various environments. In Figure 4 below, a simple PBC cell of Fe{sub 2}O{sub 3} was first optimized. A surface optimization was performed next, followed by a hydroxylated surface optimization. Once this calculation has converged, the adsorption of halide species to the hydroxylated surface will be investigated. TensorFlow is an open source platform for machine learning developed by Google. Using a high level application program interface (API) such as Keras allows for building and training ML models easily in a number of different environments and languages. For this project, a neural network was developed within Anaconda in Python. Future Work: Further development of reference data set; Refining neural network and learning algorithm; Fingerprinting atomic environment to enable mapping of atomic force components; Choosing appropriate training set from reference data; Learning from training set and enabling non-linear mapping of training set fingerprints and the atomic forces; Estimation of uncertainty to identify ranges of outside applicability; Testing and analysis of molecular dynamic simulations.

36 MATERIALS SCIENCE↗
NASA NTRS1999

A Discussion of Zero Spring Rate Mechanisms Used for the Active Isolation Mount Experiment

In the summer of 1995 the Structural Dynamics Branch at NASA Langley Research Center set out to conceive a small, lightweight, low frequency isolation mount that could be used for spaceflight experiments. The Engineering Design Branch undertook the task of developing the isolation mount. This report describes the engineering process that led to three phases of a study entitled "Active Isolation Mounts" (AIM). A zero spring rate mechanism was used to achieve low fundamental frequencies for a payloads in the 1 to 10 pound range. It worked by balancing both a positive and a negative stiffness so that the net result was a small positive stiffness. The study demonstrated devices that could reduce the initial corner frequency by a factor of six for brief periods and a factor of two for extended periods. The designs were relatively simple and minimized weight, volume, and power. They could be scaled down and they were made of spaceflight compatible materials. All designs offered the ability to continuously vary the fundamental frequency. Yet, the goal of reducing the frequency by an order of magnitude was not achieved because the systems were too unstable at low frequencies. There was a trade between performance and stability.

Teter, John E., Jr.↗
NASA NTRS2001

Diagnosing Tropical Divergent Circulations From Satellite-Derived Diabatic Forcing and Variational Integral Constraints

Large-scale divergent circulations are part of the atmospheric dynamic response to diabatic heating from condensation, radiative processes, and surface energy fluxes. Vertical motion and the associated divergent wind is thus intimately tied to the hydrologic cycle and the global beat balance. Despite its importance, the divergent circulation is too small in comparison to the rotational flow to measure directly with any accuracy. Vertical motions are recovered diagnostically from reanalyses and, as such, are subject to shortcomings in model physics, numerics, and data availability. While reanalysis estimates of tropical divergent circulations are much improved over those from the Global Weather Experiment, there are still substantial differences between products from operational centers. This is because these circulations are still forced largely by model physics and only secondarily by observations. In order to produce a refined estimate of tropical divergence and its interannual variability we have used a number of remotely-sensed data sets along with variational constraints to improve upon reanalysis estimates. Among these are: precipitation from SSM/I and GPCP, TOVS Path-A vertical cloud distributions; ISCCP radiative cooling rates; TOA radiative fluxes from ERBS, surface radiative fluxes from the SRB project, and surface latent and sensible flux estimates from SSM/I. The TOVS Path-A data constrain the divergent outflow in precipitating regions to have the same vertical structure as observed cloudiness. Using integral constraints for moisture, heat, and mass balance, we retrieve consistent divergent wind flows. We examine the ability of this type analysis to capture regional details of ENSO related perturbations to the divergent wind and associated tropical energy balance. Precipitation from satellite is found to be the major constraint in supplying this horizontal structure. We also consider the ability of this analysis to quantify integrated (land vs. ocean) fluctuations in the global monsoonal flow and energy balance. This depends, of course, on the error characteristics of the data constraints. Some estimate of these properties and their relative importance is provided.

Robertson, Franklin R.↗
DOE OSTI2022

Nondestructive Photo-Cross-Linking of Microphase-Separated Diblock Polymers through Coumarin Dimerization

The effect of long-wavelength ultraviolet photocrosslinking on microphase-separated coumarin-containing block polymers was studied by photo-rheometry and small-angle X-ray scattering. This model system consisted of three photocrosslinkable diblock polymers of poly(methoxyethyl acrylate)-b-poly(hexyl methacrylate-co-coumarin methacrylate) with different volume fractions of the crosslinkable coumarin-containing block, which microphase separated into lamellar and cylindrical morphologies. All polymers stiffened upon exposure to 365 nm light, with much greater relative increases in moduli recorded for lamellae-forming polymers (ca. 3200% increase) compared to the cylinder-forming polymer (ca. 550% increase). Disordering transitions that were evident in uncrosslinked samples were no longer observed after crosslinking in the ordered state, and domain sizes were found to remain stable to heating. The photocrosslinking reaction only proceeded under active irradiation, indicating a high degree of spatiotemporal control over curing in this system. Finally, at constant concentration of couamarin within the crosslinkable block, the cure rate was largely independent of polymer composition, suggesting a constant local concentration of coumarin moieties within the segregated crosslinkable domains. Lastly, these findings establish a set of specific structure-property relationships governing the phase-selective photocrosslinking of diblock polymers that can guide the design of robust nanostructured materials.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2016

A critical appraisal of the zero-multipole method: Structural, thermodynamic, dielectric, and dynamical properties of a water system

We performed extensive and strict tests for the reliability of the zero-multipole (summation) method (ZMM), which is a method for estimating the electrostatic interactions among charged particles in a classical physical system, by investigating a set of various physical quantities. This set covers a broad range of water properties, including the thermodynamic properties (pressure, excess chemical potential, constant volume/pressure heat capacity, isothermal compressibility, and thermal expansion coefficient), dielectric properties (dielectric constant and Kirkwood-G factor), dynamical properties (diffusion constant and viscosity), and the structural property (radial distribution function). We selected a bulk water system, the most important solvent, and applied the widely used TIP3P model to this test. In result, the ZMM works well for almost all cases, compared with the smooth particle mesh Ewald (SPME) method that was carefully optimized. In particular, at cut-off radius of 1.2 nm, the recommended choices of ZMM parameters for the TIP3P system are α ≤ 1 nm{sup −1} for the splitting parameter and l = 2 or l = 3 for the order of the multipole moment. We discussed the origin of the deviations of the ZMM and found that they are intimately related to the deviations of the equilibrated densities between the ZMM and SPME, while the magnitude of the density deviations is very small.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2014

Reconstruction of the primordial power spectrum of curvature perturbations using multiple data sets

Detailed knowledge of the primordial power spectrum of curvature perturbations is essential both in order to elucidate the physical mechanism ('inflation') which generated it, and for estimating the cosmological parameters from observations of the cosmic microwave background and large-scale structure. Hence it ought to be extracted from such data in a model-independent manner, however this is difficult because relevant cosmological observables are given by a convolution of the primordial perturbations with some smoothing kernel which depends on both the assumed world model and the matter content of the universe. Moreover the deconvolution problem is ill-conditioned so a regularisation scheme must be employed to control error propagation. We demonstrate that 'Tikhonov regularisation' can robustly reconstruct the primordial spectrum from multiple cosmological data sets, a significant advantage being that both its uncertainty and resolution are then quantified. Using Monte Carlo simulations we investigate several regularisation parameter selection methods and find that generalised cross-validation and Mallow's C{sub p} method give optimal results. We apply our inversion procedure to data from the Wilkinson Microwave Anisotropy Probe, other ground-based small angular scale CMB experiments, and the Sloan Digital Sky Survey. The reconstructed spectrum (assuming the standard ΛCDM cosmology) is not scale-free but has an infrared cutoff at k∼<5 × 10{sup −4} Mpc{sup −1} (due to the anomalously low CMB quadrupole) and several features with ∼ 2σ significance at k/Mpc{sup −1} ∼ 0.0013–0.0025, 0.0362–0.0402 and 0.051–0.056, reflecting the 'WMAP glitches'. To test whether these are indeed real will require more accurate data, such as from the Planck satellite and new ground-based experiments.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

TURBULENCE SETS THE INITIAL CONDITIONS FOR STAR FORMATION IN HIGH-PRESSURE ENVIRONMENTS

Despite the simplicity of theoretical models of supersonically turbulent, isothermal media, their predictions successfully match the observed gas structure and star formation activity within low-pressure (P/k < 10{sup 5} K cm{sup –3}) molecular clouds in the solar neighborhood. However, it is unknown whether or not these theories extend to clouds in high-pressure (P/k > 10{sup 7} K cm{sup –3}) environments, like those in the Galaxy's inner 200 pc central molecular zone (CMZ) and in the early universe. Here, we present Atacama Large Millimeter/submillimeter Array 3 mm dust continuum emission within a cloud, G0.253+0.016, which is immersed in the high-pressure environment of the CMZ. While the log-normal shape and dispersion of its column density probability distribution function (PDF) are strikingly similar to those of solar neighborhood clouds, there is one important quantitative difference: its mean column density is one to two orders of magnitude higher. Both the similarity and difference in the PDF compared to those derived from solar neighborhood clouds match predictions of turbulent cloud models given the high-pressure environment of the CMZ. The PDF shows a small deviation from log-normal at high column densities confirming the youth of G0.253+0.016. Its lack of star formation is consistent with the theoretically predicted, environmentally dependent volume density threshold for star formation which is orders of magnitude higher than that derived for solar neighborhood clouds. Our results provide the first empirical evidence that the current theoretical understanding of molecular cloud structure derived from the solar neighborhood also holds in high-pressure environments. We therefore suggest that these theories may be applicable to understand star formation in the early universe.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Stochastic embedding DFT: Theory and application to p -nitroaniline in water

Over this past decade, we combined the idea of stochastic resolution of identity with a variety of electronic structure methods. In our stochastic Kohn-Sham density functional theory (DFT) method, the density is an average over multiple stochastic samples, with stochastic errors that decrease as the inverse square root of the number of sampling orbitals. Here in this paper, we develop a stochastic embedding density functional theory method (se-DFT) that selectively reduces the stochastic error (specifically on the forces) for a selected subsystem(s). The motivation, similar to that of other quantum embedding methods, is that for many systems of practical interest, the properties are often determined by only a small subsystem. In stochastic embedding DFT, two sets of orbitals are used: a deterministic one associated with the embedded subspace and the rest, which is described by a stochastic set. The method agrees exactly with deterministic calculations in the limit of a large number of stochastic samples. We apply se-DFT to study a p-nitroaniline molecule in water, where the statistical errors in the forces on the system (the p-nitroaniline molecule) are reduced by an order of magnitude compared with nonembedding stochastic DFT.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

A Strong Gravitational Lens Is Worth a Thousand Dark Matter Halos: Inference on Small-scale Structure Using Sequential Methods

Strong gravitational lenses are a singular probe of the Universe’s small-scale structure—they are sensitive to the gravitational effects of low-mass (<10 10 M ⊙ ) halos even without a luminous counterpart. Recent strong-lensing analyses of dark matter structure rely on simulation-based inference (SBI). Modern SBI methods, which leverage neural networks as density estimators, have shown promise in extracting the halo-population signal. However, it is unclear whether the constraints from these models are limited by the methodology or the data. In this study, we introduce an accelerator-optimized simulation pipeline that can generate lens images with realistic subhalo populations in milliseconds. Leveraging this simulator, we identify the main limitation of our fiducial SBI analysis: training set size. We then adopt a sequential neural posterior estimation (SNPE) approach, allowing us to refine the training distribution to align with the observed data. Using only one-fifth as many mock Hubble Space Telescope images, SNPE matches the constraints on the low-mass halo population produced by our best nonsequential model. Our experiments suggest that an over 3 order-of-magnitude increase in training set size and GPU hours would be required to achieve an equivalent result without sequential methods. While the full potential of the existing lens sample remains to be explored, the notable improvement in constraining power enabled by our sequential approach highlights that current constraints are limited primarily by methodology and not the data itself. Moreover, our results emphasize the need to treat training set generation and model optimization as interconnected stages of any cosmological analysis using SBI.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

A numerical method for deriving shape functions of nanoparticles for pair distribution function refinements

We report that in the structural refinement of nanoparticles, discrete atomistic modeling can be used for small nanocrystals (< 15 nm), but becomes computationally unfeasible at larger sizes, where instead unit-cell-based small-box modeling is usually employed. However, the effect of the nanocrystal's shape is often ignored or accounted for with a spherical model regardless of the actual shape due to the complexities of solving and implementing accurate shape effects. Recent advancements have provided a way to determine the shape function directly from a pair distribution function calculated from a discrete atomistic model of any given shape, including both regular polyhedra (e.g. cubes, spheres, octahedra) and anisotropic shapes (e.g. rods, discs, ellipsoids) [Olds et al. (2015). J. Appl. Cryst. 48, 1651–1659], although this approach is still limited to small size regimes due to computational demands. In order to accurately account for the effects of nanoparticle size and shape in small-box refinements, a numerical or analytical description is needed. This article presents a methodology to derive numerical approximations of nanoparticle shape functions by fitting to a training set of known shape functions; the numerical approximations can then be employed on larger sizes yielding a more accurate and physically meaningful refined nanoparticle size. Lastly, the method is demonstrated on a series of simulated and real data sets, and a table of pre-calculated shape function expressions for a selection of common shapes is provided.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Electronic structure and lattice dynamics of 1⁢𝑇−VSe 2 : Origin of the three-dimensional charge density wave

To characterize in detail the charge density wave (CDW) transition of 1⁢𝑇−VSe 2 , its electronic structure and lattice dynamics are comprehensively studied by means of x-ray diffraction, muon spectroscopy, angle resolved photoemission (ARPES), diffuse and inelastic x-ray scattering, and state-of-the-art first-principles density functional theory calculations. Resonant elastic x-ray scattering does not show any resonant enhancement at either V or Se, indicating that the CDW peak at the 𝐾 edges describes a purely structural modulation of the electronic ordering. ARPES experiments identify (i) a pseudogap at 𝑇 > 𝑇⁢ CDW , which leads to a depletion of the density of states in the ML-M'L' plane at 𝑇 < 𝑇⁢ CDW , and (ii) anomalies in the electronic dispersion reflecting a sizable impact of phonons on it. A diffuse scattering precursor, characteristic of soft phonons, is observed at room temperature (RT) and leads to the full collapse of the low-energy phonon (𝜔 1 ) with propagation vector (0.25 0 −0.3) r.l.u. Here, we show that the frequency and linewidth of this mode are anisotropic in momentum space, reflecting the momentum dependence of the electron-phonon interaction (EPI), hence demonstrating that the origin of the CDW is, to a much larger extent, due to the momentum dependent EPI with a small contribution from nesting. The pressure dependence of the 𝜔 1 soft mode remains nearly constant up to 13 GPa at RT, with only a modest softening before the transition to the high-pressure monoclinic C2/m phase. The wide set of experimental data is well captured by our state-of-the art first-principles anharmonic calculations with the inclusion of van der Waals corrections in the exchange-correlation functional. The comprehensive description of the electronic and dynamical properties of VSe 2 reported here adds important pieces of information to the understanding of the electronic modulations in the family of transition-metal dichalcogenides.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Localizing Clinical Patterns of Blast Traumatic Brain Injury Through Computational Modeling and Simulation

Blast traumatic brain injury is ubiquitous in modern military conflict with significant morbidity and mortality. Yet the mechanism by which blast overpressure waves cause specific intracranial injury in humans remains unclear. Reviewing of both the clinical experience of neurointensivists and neurosurgeons who treated service members exposed to blast have revealed a pattern of injury to cerebral blood vessels, manifested as subarachnoid hemorrhage, pseudoaneurysm, and early diffuse cerebral edema. Additionally, a seminal neuropathologic case series of victims of blast traumatic brain injury (TBI) showed unique astroglial scarring patterns at the following tissue interfaces: subpial glial plate, perivascular, periventricular, and cerebral gray-white interface. The uniting feature of both the clinical and neuropathologic findings in blast TBI is the co-location of injury to material interfaces, be it solid-fluid or solid-solid interface. This motivates the hypothesis that blast TBI is an injury at the intracranial mechanical interfaces. In order to investigate the intracranial interface dynamics, we performed a novel set of computational simulations using a model human head simplified but containing models of gyri, sulci, cerebrospinal fluid (CSF), ventricles, and vasculature with high spatial resolution of the mechanical interfaces. Simulations were performed within a hybrid Eulerian—Lagrangian simulation suite (CTH coupled via Zapotec to Sierra Mechanics). Because of the large computational meshes, simulations required high performance computing resources. Twenty simulations were performed across multiple exposure scenarios—overpressures of 150, 250, and 500 kPa with 1 ms overpressure durations—for multiple blast exposures (front blast, side blast, and wall blast) across large variations in material model parameters (brain shear properties, skull elastic moduli). All simulations predict fluid cavitation within CSF (where intracerebral vasculature reside) with cavitation occurring deep and diffusely into cerebral sulci. These cavitation events are adjacent to high interface strain rates at the subpial glial plate. Larger overpressure simulations (250 and 500kPa) demonstrated intraventricular cavitation—also associated with adjacent high periventricular strain rates. Additionally, models of embedded intraparenchymal vascular structures—with diameters as small as 0.6 mm—predicted intravascular cavitation with adjacent high perivascular strain rates. The co-location of local maxima of strain rates near several of the regions that appear to be preferentially damaged in blast TBI (vascular structures, subpial glial plate, perivascular regions, and periventricular regions) suggest that intracranial interface dynamics may be important in understanding how blast overpressures leads to intracranial injury.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2013

SASI ACTIVITY IN THREE-DIMENSIONAL NEUTRINO-HYDRODYNAMICS SIMULATIONS OF SUPERNOVA CORES

The relevance of the standing accretion shock instability (SASI) compared to neutrino-driven convection in three-dimensional (3D) supernova-core environments is still highly controversial. Studying a 27 M{sub Sun} progenitor, we demonstrate, for the first time, that violent SASI activity can develop in 3D simulations with detailed neutrino transport despite the presence of convection. This result was obtained with the PROMETHEUS-VERTEX code with the same sophisticated neutrino treatment so far used only in one-dimensional and two-dimensional (2D) models. While buoyant plumes initially determine the nonradial mass motions in the postshock layer, bipolar shock sloshing with growing amplitude sets in during a phase of shock retraction and turns into a violent spiral mode whose growth is only quenched when the infall of the Si/SiO interface leads to strong shock expansion in response to a dramatic decrease of the mass accretion rate. In the phase of large-amplitude SASI sloshing and spiral motions, the postshock layer exhibits nonradial deformation dominated by the lowest-order spherical harmonics (l = 1, m = 0, {+-}1) in distinct contrast to the higher multipole structures associated with neutrino-driven convection. We find that the SASI amplitudes, shock asymmetry, and nonradial kinetic energy in three dimensions can exceed those of the corresponding 2D case during extended periods of the evolution. We also perform parameterized 3D simulations of a 25 M{sub Sun} progenitor, using a simplified, gray neutrino transport scheme, an axis-free Yin-Yang grid, and different amplitudes of random seed perturbations. They confirm the importance of the SASI for another progenitor, its independence of the choice of spherical grid, and its preferred growth for fast accretion flows connected to small shock radii and compact proto-neutron stars as previously found in 2D setups.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1994

Steps toward determination of the size and structure of the broad-line region in active galactic nuclei. 6: Variability of NGC 3783 from ground-based data

The Seyfert 1 galaxy NGC 3783 was intensely monitored in several bands between 1991 December and 1992 August. This paper presents the results from the ground-based observations in the optical and near-IR bands, which complement the data set formed by the International Ultraviolet Explorer (IUE) spectra, discussed elsewhere. Spectroscopic and photometric data from several observatories were combined in order to obtain well-sampled light curves of the continuum and of H(beta). During the campaign the source underwent significant variability. The light curves of the optical continuum and of H(beta) display strong similarities to those obtained with the IUE. The near-IR flux did not vary significantly except for a slight increase at the end of the campaign. The cross-correlation analysis shows that the variations of the optical continuum have a lag of 1 day or less with respect to those of the UV continuum, with an uncertainty of is less than or equal to 4 days. The integrated flux of H(beta) varies with a delay of about 8 days. These results confirm that (1) the continuum variations occur simultaneously or with a very small lag across the entire UV-optical range, as in the Seyfert galaxy NGC 5548; and (2) the emission lines of NGC 3783 respond to ionizing continuum variations with less delay than those of NGC 5548. As observed in NGC 5548, the lag of H(beta) with respect to the continuum is greater than those of the high-ionization lines.

Stirpe, G. M.↗
DOE OSTI2020

TuckerMPI: A Parallel C++/MPI Software Package for Large-scale Data Compression via the Tucker Tensor Decomposition

With this study, our goal is compression of massive-scale grid-structured data, such as the multi-terabyte output of a high-fidelity computational simulation. For such data sets, we have developed a new software package called TuckerMPI, a parallel C++/MPI software package for compressing distributed data. The approach is based on treating the data as a tensor, i.e., a multidimensional array, and computing its truncated Tucker decomposition, a higher-order analogue to the truncated singular value decomposition of a matrix. The result is a low-rank approximation of the original tensor-structured data. Compression efficiency is achieved by detecting latent global structure within the data, which we contrast to most compression methods that are focused on local structure. In this work, we describe TuckerMPI, our implementation of the truncated Tucker decomposition, including details of the data distribution and in-memory layouts, the parallel and serial implementations of the key kernels, and analysis of the storage, communication, and computational costs. We test the software on 4.5 and 6.7 terabyte data sets distributed across 100 s of nodes (1,000 s of MPI processes), achieving compression ratios between 100 and 200,000×, which equates to 99--99.999% compression (depending on the desired accuracy) in substantially less time than it would take to even read the same dataset from a parallel file system. Moreover, we show that our method also allows for reconstruction of partial or down-sampled data on a single node, without a parallel computer so long as the reconstructed portion is small enough to fit on a single machine, e.g., in the instance of reconstructing/visualizing a single down-sampled time step or computing summary statistics. The code is available at https://gitlab.com/tensors/TuckerMPI.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1996

On the Dynamics of Small-Scale Solar Magnetic Elements

We report on the dynamics of the small-scale solar magnetic field, based on analysis of very high resolution images of the solar photosphere obtained at the Swedish Vacuum Solar Telescope. The data sets are movies from 1 to 4 hr in length, taken in several wavelength bands with a typical time between frames of 20 s. The primary method of tracking small-scale magnetic elements is with very high contrast images of photospheric bright points, taken through a 12 A bandpass filter centered at 4305 A in the Fraunhofer 'G band.' Previous studies have established that such bright points are unambiguously associated with sites of small-scale magnetic flux in the photosphere, although the details of the mechanism responsible for the brightening of the flux elements remain uncertain. The G band bright points move in the intergranular lanes at speeds from 0.5 to 5 km/s. The motions appear to be constrained to the intergranular lanes and are primarily driven by the evolution of the local granular convection flow field. Continual fragmentation and merging of flux is the fundamental evolutionary mode of small-scale magnetic structures in the solar photosphere. Rotation and folding of chains or groups of bright points are also observed. The timescale for magnetic flux evolution in active region plage is on the order of the correlation time of granulation (typically 6-8 minutes), but significant morphological changes can occur on timescales as short as 100 S. Smaller fragments are occasionally seen to fade beyond observable contrast. The concept of a stable, isolated subarcsecond magnetic 'flux tube' in the solar photosphere is inconsistent with the observations presented here.

Berger, T. E.↗
DOE OSTI2021

Magnetic topological quantum chemistry

For over 100 years, the group-theoretic characterization of crystalline solids has provided the foundational language for diverse problems in physics and chemistry. However, the group theory of crystals with commensurate magnetic order has remained incomplete for the past 70 years, due to the complicated symmetries of magnetic crystals. In this work, we complete the 100-year-old problem of crystalline group theory by deriving the small corepresentations, momentum stars, compatibility relations, and magnetic elementary band corepresentations of the 1,421 magnetic space groups (MSGs), which we have made freely accessible through tools on the Bilbao Crystallographic Server. We extend Topological Quantum Chemistry to the MSGs to form a complete, real-space theory of band topology in magnetic and nonmagnetic crystalline solids – Magnetic Topological Quantum Chemistry (MTQC). Using MTQC, we derive the complete set of symmetry-based indicators of electronic band topology, for which we identify symmetry-respecting bulk and anomalous surface and hinge states.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗