Search NASA⌕ Search

SEARCH · Search NASA

Results for “Dynamical cluster approximation”

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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Temperature dependence of structural, dynamical, and electronic properties of amorphous Bi 2 Te 3 : an ab initio study

Bismuth telluride (Bi 2 Te 3 ) has garnered significant interest in thermoelectric applications and three-dimensional topological insulators due to its unique electronic, transport, and thermal properties. Bi 2 Te 3 and Sb 2 Te 3 chalcogenide compounds have the same crystal structure. While Sb 2 Te 3 has been shown to be a prototypical phase change memory (PCM) compound along the pseudobinary tie-line of Ge-Sb-Te alloys, whether Bi 2 Te 3 can also exhibit PCM functionality is still not well established. In this work, a systematic study on the structural, dynamical, and electronic properties of amorphous Bi 2 Te 3 during the quenching process has been performed by using ab initio molecular dynamics simulations. Pair correlation function, coordination number, bond-angle distribution functions, and a novel atomistic cluster alignment method are used to explore the structural characteristics of Bi 2 Te 3 as a function of temperature. Our study shows that there are many distorted octahedral clusters in amorphous Bi 2 Te 3 . In comparison with the local structures in Sb 2 Te 3 , we found that the degree of distortion of the octahedrons in the Bi 2 Te 3 system is smaller than that in Sb 2 Te 3 system. Moreover, the changes in the dynamical properties of Bi 2 Te 3 from liquid to glassy state are also explored. The approximate range of liquid-to-glass transition temperature is determined to be between 673 and 723 K. The electronic properties of Bi2Te 3 and Sb 2 Te 3 are also analysed by density-of-states and Bader charge calculations, both of them in glass state are semiconductors. Our studies provide useful insights into the local structure and dynamical properties of Bi 2 Te 3 at the atomistic level during the fast cooling process, and suggest that the compound can be a candidate for PCM materials.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Weyl rings and enhanced susceptibilities in pyrochlore iridates: k · p analysis of cluster dynamical mean-field theory results

We match analytic results to numerical calculations to provide a detailed picture of the metal-insulator and topological transitions found in density functional plus cluster dynamical mean-field calculations of pyrochlore iridates. We discuss the transition from Weyl metal to Weyl semimetal regimes, and then analyze in detail the properties of the Weyl semimetal phase and its evolution into the topologically trivial insulator. The energy scales in the Weyl semimetal phase are found to be very small, as are the anisotropy parameters. The electronic structure can, to a good approximation, be described as “Weyl rings” and one of the two branches that contributes to the Weyl bands is essentially flat, leading to enhanced susceptibilities. The optical longitudinal and Hall conductivities are determined; the frequency dependence includes pronounced features that reveal the basic energy scales of the Weyl semimetal phase.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Investigating the crust of neutron stars with neural-network quantum states

An accurate description of low-density nuclear matter is crucial for explaining the physics of neutron star crusts. In the density range between approximately 0.01 fm −3 and 0.1 fm −3 , matter transitions from neutron-rich nuclei to various higher-density pasta shapes, before ultimately reaching a uniform liquid. In this work, we introduce a variational Monte Carlo method based on a neural Pfaffian-Jastrow quantum state, which allows us to model the transition from the liquid phase to neutron-rich nuclei microscopically. At low densities, nuclear clusters dynamically emerge from the microscopic interactions among protons and neutrons, which we model based on pionless effective field theory. Our variational Monte Carlo approach represents a significant improvement over the state-of-the-art auxiliary-field diffusion Monte Carlo method, which is severely hindered by the fermion-sign problem in this low-density regime and cannot capture the onset of clusters. In addition to computing the energy per particle of symmetric nuclear matter and pure neutron matter, we analyze an intermediate isospin-asymmetry configuration to elucidate the formation of nuclear clusters. We also provide evidence that the presence of such nuclear clusters influences the amount of protons in the crust compared to protons in beta-equilibrated, neutrino-transparent matter.

Nuclear astrophysics↗
DOE OSTI2021

Pressure Effects on the Relaxation of an Excited Ethane Molecule in High-Pressure Bath Gases

Here, we use molecular dynamics to calculate the rotational and vibrational energy relaxation of C 2 H 6 in Ar, Kr, and Xe bath gases over a pressure range of 10 to 400 atm and at temperatures of 300 K and 800 K. The C 2 H 6 is instantaneously excited by 80 kcal/mol randomly distributed into both vibrational and rotational modes. The computed relaxation rates show little sensitivity to the identity of the noble gas in the bath. Vibrational relaxation rates show a non-linear pressure dependence at 300 K. At 800 K the reduced range of bath gas densities covered by the range of pressures do not yet show any non-linearity in the pressure dependence. Rotational relaxation is characterized with two relaxation rates. The slower rate is comparable to the vibrational relaxation rate. The faster rate has a linear pressure dependence at 300 K but an irregular, nonlinear pressure dependence at 800 K. To understand this, a model was developed based on approximating the periodic box used in the molecular dynamics simulations by an equal-volume collection of cubes where each cube is sized to allow only single occupancy by the noble gas or the molecule. Combinatorial statistics then leads to a pressure and temperature dependent analytic distribution of the bath gas species the molecule encounters in a collision. This distribution, the dissociation energy of molecule/bath gas complexes and bath gas clusters, and the computed energy release per collision combine to show that only at 300 K is the energy release sufficient to dissociate likely complexes and clusters. This suggests that persistent and pressure-dependent clusters and complexes at 800 K may be responsible for the non-linear pressure dependence of rotational relaxation.

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

Dynamic Observation of Dendritic Quasicrystal Growth upon Laser-Induced Solid-State Transformation

In this paper, we report the laser-induced solid-state transformation between a periodic “approximant” and quasicrystal in the Al-Cr system during rapid quenching. Dynamic transmission electron microscopy allows us to capture in situ the dendritic growth of the metastable quasicrystals. The formation of dendrites during solid-state transformation is a rare phenomenon, which we attribute to the structural similarity between the two intermetallics. Through ab initio molecular dynamics simulations, we identify the dominant structural motif to be a 13-atom icosahedral cluster transcending the phases of matter.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Vibrational coherence transfer illuminates dark modes in models of the FeFe hydrogenase active site

Within the conceptual framework of Redfield theory, the optical response function arises from the dynamical evolution of the system’s density operator, where nonunitary relaxation is encoded in the Redfield relaxation superoperator. In the conventional approach, the so-called secular approximation neglects terms that induce transitions between distinct coherences and among coherences and populations. The rationale is that these nonsecular terms are small in comparison to the far more dominant population relaxation and coherence dephasing contributions. Since two-dimensional infrared (2D-IR) spectroscopy has significant contributions arising from population relaxation and transfer pathways, it can be challenging to isolate signatures of the nonsecular relaxation. We report om this paper that in three diiron dithiolate hexacarbonyl complexes that serve as small-molecule models of the [FeFe] hydrogenase H-cluster subsite, a fortuitous vibrational energy structure enables direct and clear signatures of vibrational coherence transfer in alkane solution. This finding holds promise towards developing a molecularly detailed understanding of the mechanism of vibrational coherence transfer processes, thanks to the ease of synthesizing derivatives based on the chemical modularity of these well studied diiron compounds. In addition to the fundamental need to characterize coherence transfer in molecular spectroscopy, we find in this set of molecules a practical utility for the nonsecular dynamics: the ability to determine the frequency of an IR-inactive mode. A coherence generated during the waiting time of the 2D-IR measurement transfers to a coherence involving the single dark CO stretching mode, which modulates some peak amplitudes in the 2D spectrum, revealing its transient excitation.

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

Near infrared overtone (v OH = 2 ← 0) spectroscopy of Ne–H 2 O clusters

Vibrationally state selective overtone spectroscopy and dynamics of weakly bound Ne–H 2 O complexes (D 0 (para) = 31.67 cm -1 , D 0 (ortho) = 34.66 cm -1 ) are reported for the first time, based on near infrared excitation of van der Waals cluster bands correlating with v OH = 2 ← 0 overtone transitions (|02 - $\rangle$←|00+$\rangle$ and |02 + $\rangle$←|00 + $\rangle$) out of the ortho (1 01 ) and para (0 00 ) internal rotor states of the H 2 O moiety. Quantum theoretical calculations for nuclear motion on a high level ab initio potential energy surface (CCSD(T)/VnZ-f12 (n = 3,4), corrected for basis set superposition error and extrapolated to the complete basis set limit) are employed for assignment of Σ←Σ,Π←Σ, and Σ←Π infrared bands in the overtone spectra, where Σ(K = 0) and Π (K = 1) represent approximate projections (K) of the body angular momentum along the Ne–H 2 O internuclear axis. End-over-end tumbling of the ortho Ne–H 2 O cluster is evident via rotational band contours observed, with band origins and rotational progressions in excellent agreement with ab initio frequency and intensity predictions. A clear Q branch in the corresponding |02 + $\rangle$ f Π(1 11 )← e Σ(0 00 ) para Ne–H 2 O spectrum provides evidence for a novel e/f parity-dependent metastability in these weakly bound clusters, in agreement with ab initio bound state calculations and attributable to the symmetry blocking of an energetically allowed channel for internal rotor predissociation. Lastly, Boltzmann analysis of the rotational spectra reveals anomalously low jet temperatures (T rot ≈ 4(1) K), which are attributed to “evaporative cooling” of weakly bound Ne–H 2 O clusters and provide support for similar cooling dynamics in rare gas-tagging studies.

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

OctoFAS: A Two-Level Fair Scheduler That Increases Fairness in Network-Based Key-Value Storage

We identified a fairness problem in a network-based key-value storage system using Intel Storage Performance Development Kit (SPDK) in a multitenant environment. In such an environment, each tenant’s I/O service rate is not fairly guaranteed compared to that of other tenants. To address the fairness problem, we propose OctoFAS, a two-level fair scheduler designed to improve overall throughput and fairness among tenants. The two-level scheduler of OctoFAS consists of (i) inter-core scheduling and (ii) intra-core scheduling. Through inter-core scheduling, OctoFAS addresses the load imbalance problem that is inherent in SPDK on the storage server by dynamically migrating I/O requests from overloaded cores to underloaded cores, thereby increasing overall throughput. Intra-core scheduling prioritizes handling requests from starving tenants over well-fed tenants within core-specific event queues to ensure fair I/O services among multiple tenants. OctoFAS is deployed on a Linux cluster with SPDK. Through extensive evaluations, we found that OctoFAS ensures that the total system throughput remains high and balanced, while enhancing fairness by approximately 10% compared to the baseline, when both scheduling levels operate in a hybrid fashion.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

A Generative Model for Realistic Galaxy Cluster X-Ray Morphologies

Abstract The X-ray morphologies of clusters of galaxies display significant variations, reflecting their dynamical histories and the nonlinear dependence of X-ray emissivity on the density of the intracluster gas. Qualitative and quantitative assessments of X-ray morphology have long been considered a proxy for determining whether clusters are dynamically active or “relaxed.” Conversely, the use of circularly or elliptically symmetric models for cluster emission can be complicated by the variety of complex features realized in nature, spanning scales from megaparsecs down to the resolution limit of current X-ray observatories. In this work, we use mock X-ray images from simulated clusters from The Three Hundred project to define a basis set of cluster image features. We take advantage of the clusters’ approximate self-similarity to minimize the differences between images before encoding the remaining diversity through a distribution of high-order polynomial coefficients. Principal component analysis then provides an orthogonal basis for this distribution, corresponding to natural perturbations from an average model. This representation allows novel, realistically complex X-ray cluster images to be easily generated, and we provide code to do so. The approach provides a simple way to generate training data for cluster image analysis algorithms and could be straightforwardly adapted to generate clusters displaying specific types of features or selected by physical characteristics available in the original simulations.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1994

Gradients in giant branch morphology in the core of 47 Tucanae

I describe an algorithm which uses the high spatial resolution of the Hubble Space Telescope to complement the high spatial-to-noise, approximately symmetric point response function, relatively large spatial coverage, and standard filters available from ground based images of crowded fields. Applying this technique to the central regions of the globular cluster 47 Tucanae, I find that the morphology of the giant branch in the core is significantly different from that in more distant regions (r approximately equals 5 to 10 core radii) of the cluster. In particular, there appear to be fewer bright giants in the core, along with an enhanced `asymptotic giant branch' (AGB) sequence. Depletion of giants has been observed in the cores of other dense clusters, and may be due to `stripping' of large stars by stellar encounters and/or mass transfer in binary systems. Central concentrations of true asymptotic giant branch stars are not expected to result from dynamical processes; possibly some of these stars may be evolved blue stragglers.

Bailyn, Charles D.↗
DOE OSTI2023

Evolution of 1 / f Flux Noise in Superconducting Qubits with Weak Magnetic Fields

The microscopic description of 1/f magnetic flux noise in superconducting circuits has remained an open question for several decades despite extensive experimental and theoretical investigation. Recent progress in superconducting devices for quantum information has highlighted the need to mitigate sources of qubit decoherence, driving a renewed interest in understanding the underlying noise mechanism(s). Though a consensus has emerged attributing flux noise to surface spins, their identity and interaction mechanisms remain unclear, prompting further study. Here, we apply weak in-plane magnetic fields to a capacitively shunted flux qubit (where the Zeeman splitting of surface spins lies below the device temperature) and study the flux-noise-limited qubit dephasing, revealing previously unexplored trends that may shed light on the dynamics behind the emergent 1/f noise. Notably, we observe an enhancement (suppression) of the spin-echo (Ramsey) pure-dephasing time in fields up to B = 100 G. With direct noise spectroscopy, we further observe a transition from a 1/f to approximately Lorentzian frequency dependence below 10 Hz and a reduction of the noise above 1 MHz with increasing magnetic field. Here, we suggest that these trends are qualitatively consistent with an increase of spin cluster sizes with magnetic field. These results should help to inform a complete microscopic theory of 1/f flux noise in superconducting circuits.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Ultra-fast computation of electronic spectra for large systems by tight-binding based simplified Tamm-Dancoff approximation (sTDA-xTB)

The computational bottleneck of the extremely fast simplified Tamm-Dancoff approximated (sTDA) time-dependent density functional theory procedure [S. Grimme, J. Chem. Phys. 138, 244104 (2013)] for the computation of electronic spectra for large systems is the determination of the ground state Kohn-Sham orbitals and eigenvalues. This limits such treatments to single structures with a few hundred atoms and hence, e.g., sampling along molecular dynamics trajectories for flexible systems or the calculation of chromophore aggregates is often not possible. The aim of this work is to solve this problem by a specifically designed semi-empirical tight binding (TB) procedure similar to the well established self-consistent-charge density functional TB scheme. The new special purpose method provides orbitals and orbital energies of hybrid density functional character for a subsequent and basically unmodified sTDA procedure. Compared to many previous semi-empirical excited state methods, an advantage of the ansatz is that a general eigenvalue problem in a non-orthogonal, extended atomic orbital basis is solved and therefore correct occupied/virtual orbital energy splittings as well as Rydberg levels are obtained. A key idea for the success of the new model is that the determination of atomic charges (describing an effective electron-electron interaction) and the one-particle spectrum is decoupled and treated by two differently parametrized Hamiltonians/basis sets. The three-diagonalization-step composite procedure can routinely compute broad range electronic spectra (0-8 eV) within minutes of computation time for systems composed of 500-1000 atoms with an accuracy typical of standard time-dependent density functional theory (0.3-0.5 eV average error). An easily extendable parametrization based on coupled-cluster and density functional computed reference data for the elements H–Zn including transition metals is described. The accuracy of the method termed sTDA-xTB is first benchmarked for vertical excitation energies of open- and closed-shell systems in comparison to other semi-empirical methods and applied to exemplary problems in electronic spectroscopy. As side products of the development, a robust and efficient valence electron TB method for the accurate determination of atomic charges as well as a more accurate calculation scheme of dipole rotatory strengths within the Tamm-Dancoff approximation is proposed.

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

Uncertainty-aware molecular dynamics from Bayesian active learning for phase transformations and thermal transport in SiC

Abstract Machine learning interatomic force fields are promising for combining high computational efficiency and accuracy in modeling quantum interactions and simulating atomistic dynamics. Active learning methods have been recently developed to train force fields efficiently and automatically. Among them, Bayesian active learning utilizes principled uncertainty quantification to make data acquisition decisions. In this work, we present a general Bayesian active learning workflow, where the force field is constructed from a sparse Gaussian process regression model based on atomic cluster expansion descriptors. To circumvent the high computational cost of the sparse Gaussian process uncertainty calculation, we formulate a high-performance approximate mapping of the uncertainty and demonstrate a speedup of several orders of magnitude. We demonstrate the autonomous active learning workflow by training a Bayesian force field model for silicon carbide (SiC) polymorphs in only a few days of computer time and show that pressure-induced phase transformations are accurately captured. The resulting model exhibits close agreement with both ab initio calculations and experimental measurements, and outperforms existing empirical models on vibrational and thermal properties. The active learning workflow readily generalizes to a wide range of material systems and accelerates their computational understanding.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Uncertainty quantification for incident helium flux in plasma-exposed tungsten

In this work, the surface response of a tungsten plasma-facing component was simulated by a cluster-dynamics code, Xolotl, with a focus on quantifying the impact of uncertainty in one of the input parameters to Xolotl, namely, the incident helium flux. The simulated conditions involve a tungsten surface exposed to 100 eV helium ion implantations with a flux of either 4 × 10 22 or 4 × 10 25 He m -2 s -1 . Furthermore, two sources were used to describe the implanted helium depth distribution in tungsten, either molecular dynamics (MD) or a binary collision approximation code, the stopping and range of ions in matter (SRIM). The aim of this work is to evaluate and examine uncertain predictions on the helium retention based on these two different modeling methodologies that either neglect electronic energy loss or the crystalline structure of the solid, respectively. An embedded model-form error approach was pursued here in order to arrive at predictions that account for variability due to the two different data sources, and the impact of this model-form uncertainty in incident helium flux on Xolotl output was presented for the two implantation fluxes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Hetero-site Double Core Ionization Energies with Sub-electronvolt Accuracy from Delta-Coupled-Cluster Calculations

Benchmark scalar-relativistic delta-coupled-cluster calculations of hetero-site double core ionization energies of small molecules containing second-row elements are reported. The present study has focused on the high-spin triplet components of two-site double core-ionized states, which are single reference in character and consistent with the use of standard coupledcluster methods. Contributions to computed double core ionization energies from electron-correlation and basis-set effects as well as corrections to the core-valence separation approximation have been analyzed. Furthermore, on the basis of systematic convergence of computational results with respect to these effects, delta-coupled-cluster calculations have been shown to be capable of providing accurate double core ionization energies with remaining errors estimated to be below 0.3 eV, and thus they are recommended for use to facilitate experimental studies of two-site double core-ionized states that are involved in X-ray pump/X-ray probe studies of electronic and molecular dynamics following inner shell ionization or excitation.

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

Communication: Isotopic effects on tunneling motions in the water trimer

We present results of ring polymer molecular dynamics simulations that shed light on the effects of nuclear quantum fluctuations on tunneling motions in cyclic [H{sub 2}O]{sub 3} and [D{sub 2}O]{sub 3}, at the representative temperature of T = 75 K. In particular, we focus attention on free energies associated with two key isomerization processes: The first one corresponds to flipping transitions of dangling OH bonds, between up and down positions with respect to the O–O–O plane of the cluster; the second involves the interchange between connecting and dangling hydrogen bond character of the H-atoms in a tagged water molecule. Zero point energy and tunneling effects lead to sensible reductions of the free energy barriers. Due to the lighter nature of the H nuclei, these modifications are more marked in [H{sub 2}O]{sub 3} than in [D{sub 2}O]{sub 3}. Estimates of the characteristic time scales describing the flipping transitions are consistent with those predicted based on standard transition-state-approximation arguments.

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

Investigating the Relativistic Motion of the Stars Near the Supermassive Black Hole in the Galactic Center

The S-star cluster in the Galactic center allows us to study the physics close to a supermassive black hole, including distinctive dynamical tests of general relativity. Our best estimates for the mass of and the distance to Sgr A* using the three stars with the shortest period (S2, S38, and S55/S0-102) and Newtonian models are M {sub BH} = (4.15 ± 0.13 ± 0.57) × 10{sup 6} M {sub ⊙} and R {sub 0} = 8.19 ± 0.11 ± 0.34 kpc. Additionally, we aim at a new and practical method to investigate the relativistic orbits of stars in the gravitational field near Sgr A*. We use a first-order post-Newtonian approximation to calculate the stellar orbits with a broad range of periapse distance r {sub p} . We present a method that employs the changes in orbital elements derived from elliptical fits to different sections of the orbit. These changes are correlated with the relativistic parameter defined as ϒ ≡ r {sub s} / r {sub p} (with r {sub s} being the Schwarzschild radius) and can be used to derive ϒ from observational data. For S2 we find a value of ϒ = 0.00088 ± 0.00080, which is consistent, within the uncertainty, with the expected value of ϒ = 0.00065 derived from M {sub BH} and the orbit of S2. We argue that the derived quantity is unlikely to be dominated by perturbing influences such as noise on the derived stellar positions, field rotation, and drifts in black hole mass.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Effect of stoichiometry and hydration level on water domain size and transport in poly(aryl piperidinium) alkaline anion-exchange membranes

Alkaline water electrolysis holds promise in decarbonizing the global economy by enabling renewable hydrogen production with non-precious group metal catalysts. Anion exchange membranes are an important component of alkaline water electrolyzers and would ideally be durable while allowing for high hydroxide conductivity. The poly(aryl piperidinium) (PAP) class of polymers has attracted recent interest due to their good mechanical robustness and high ionic conductivity. Here, in this work, we perform atomistic molecular dynamics (MD) simulations of several PAP polymers at experimentally relevant hydration levels and polymer ion exchange capacities (IECs) to gain nanoscale insight into their properties and to help elucidate the trade-offs that result from tuning the IECs through the polymer stoichiometry. Our MD-predicted macroscopic polymer properties were found to be in good agreement with experimentally available polymer swelling ratios, water-occupied volumes, X-ray scattering, and ionic conductivities. The models show that for hydration levels greater than 8H 2 O per cation a single water cluster will form that percolates through the system. The growth in water cluster size results in large polymer swelling, the creation of larger channels with widths of 7 Å or larger, and nanophase separation between the hydrophilic domains and the polymer with characteristic length scales of approximately 20–30 Å. The experimentally observed lack of a strong X-ray scattering peak at low wavevectors can be explained by a cancellation between the polymer-polymer/water-water and polymer-water correlations and not a loss in nanophase separation. The overlap in coordination environments of the hydroxide oxygen and polymer nitrogen atoms implies that vehicular diffusion between cationic groups could play a role in hydroxide transport. The polymers' hydroxide and water diffusion constants increase by approximately an order of magnitude between hydration levels of 8 and 20H 2 O per cation. However, there are diminishing returns in hydroxide diffusion constant once the IEC exceeds 2.4 meq/g.

08 HYDROGEN↗