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

Investigation of the isoscalar response of 24 Mg to 6 Li scattering

24 Mg is a strongly deformed nucleus in the ground state. Deformation effects can be observed in the structure of the isoscalar giant monopole and quadrupole resonances. 24 Mg is also a nucleus that is well known to present different types of cluster-oscillation modes. Both giant resonances and cluster states are strongly populated by isoscalar transitions. To extract the E0, E1, and E2 transition strengths via 6 Li scattering. The 6 Li probe is a powerful tool for investigating the isoscalar nuclear response with a very favorable ratio of resonance-to-continuum background. Double-differential cross sections of 6 Li inelastic scattering, at the beam energy of 100 MeV/u, were measured in the excitation-energy range 10 – 40 MeV and scattering angles 0 - 3°. A multipole-decomposition analysis was performed for extracting the isoscalar E0 , E1, and E2 strength distributions. Results: The extracted multipole strengths were compared with predictions from consistent quasiparticle random phase approximation calculations. The theoretical predictions are in fair agreement with the experimental data. The E0 strength was also compared with results from antisymmetrized molecular dynamics calculations found in the literature. A few peaks in the experimental data might be associated with clustering in 24 Mg. Ground-state deformation effects were observed in the isoscalar giant monopole resonance (ISGMR) and isoscalar giant quadrupole resonance (ISGQR) distributions. The ISGMR strength is split in two peaks around 19 and 28 MeV. The ISGQR exhibits a pronounced peak at 20 MeV with a broadening at the low-energy region, similar to predictions from microscopic calculations. Signatures of excitation of cluster states were observed in the E0 response. Further studies including particle-decay measurements will be required to confirm the nature of the observed peaks.

20 ≤ A ≤ 38↗
DOE OSTI2021

Origins of Clustering of Metalate–Extractant Complexes in Liquid–Liquid Extraction

Effective and energy-efficient separation of precious and rare metals is very important for a variety of advanced technologies. Liquid-liquid extraction (LLE) is a relatively less energy intensive separation technique, widely used in separation of lanthanides, actinides, and platinum group metals (PGMs). In LLE, the distribution of an ion between an aqueous phase and an organic phase is determined by enthalpic (coordination interactions) and entropic (fluid reorganization) contributions. The molecular scale details of these contributions are not well understood. Preferential extraction of an ion from the aqueous phase is usually correlated with the resulting fluid organization in the organic phase, as the longer-range organization increases with metal loading. However, it is difficult to determine the extent to which organic phase fluid organization causes, or is caused by, metal loading. In this study, we demonstrate that two systems with the same metal loading may impart very different organic phase organizations and investigate the underlying molecular scale mechanism. Small-angle X-ray scattering shows that the structure of a quaternary ammonium extractant solution in toluene is affected differently by the extraction of two metalates (octahedral PtCl 6 2- and square-planar PtCl 6 2- ), although both are completely transferred into the organic phase. The aggregates formed by the metalate-extractant complexes (approximated as reverse micelles) exhibit a more long-range order (clustering) with PtCl 6 2- compared to that with PtCl 6 2- . Vibrational sum frequency generation spectroscopy and complementary atomistic molecular dynamics simulations on model Langmuir monolayers indicate that the two metalates affect the interfacial hydration structures differently. Furthermore, the interfacial hydration is correlated with water extraction into the organic phase. Overall, these results support a strong relationship between the organic phase organizational structure and the different local hydration present within the aggregates of metalate-extractant complexes, which is independent of metalate concentration.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1983

CCD photometry of Abell clusters. II - Surface photometry of 249 cluster galaxies

The luminosity of the brightest cluster galaxies remains the best standard candle for observational cosmology. Studies of over 250 relatively nearby clusters show that the luminosity of these objects within an aperture of constant metric diameter has a dispersion of only 0.35 mag. The observations considered by Schneider et al. (1983) were acquired to measure the aperture luminosities of the brightest cluster galaxies. However, the data base, because of its high signal-to-noise ratio and large spatial coverage, allows several other investigations to be performed. In the present investigation, surface photometry of the brightness three galaxies within approximately 250 kpc of the cluster center is presented, taking into account the cores of 83 Abell clusters. The relevance of the photometry to cosmological studies is also discussed. It is found that dynamical friction does play an extremely important role in determining properties of brightest cluster galaxies and cluster morphology.

Schneider, D. P.↗
DOE OSTI2025

Simulations of Attosecond Metallization in Quartz and Diamond Probed with Inner-Shell Transient Absorption Spectroscopy

When dielectrics are hit with intense infrared (IR) laser pulses, transient metalization can occur. The initial attosecond dynamics behind this metallization are not entirely understood. Therefore, simulations are needed to understand this process and to help interpret experimental observations of it, such as with attosecond transient absorption (ATA). In this paper, we present first-principles simulations of ATA based on bulk-mimicking clusters and real-time time-dependent density functional theory (RT-TDDFT), with Koopmans-tuned range-separated hybrid functionals and Gaussian basis sets. Our method gives good agreement with the experiment for the breakdown threshold in silica and diamond. This breakdown voltage corresponds to a Keldysh parameter of approximately one and thus involves a transition to a regime where the dynamics are driven by tunneling. Pumping at an amplitude just below this value causes a mixture of multiphoton and tunneling excitations across the band gap to occur. The computed extreme ultraviolet and X-ray attosecond transient spectra also agree well with the experiment and show a decrease in optical density due to the transient population of the conduction band from the IR field. First-principles approaches such as this are valuable for interpreting the complicated modulations in a spectrum and for guiding future attosecond experiments on solids.

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

Cluster spin glass correlations and dynamics in Zn 0.5⁢ Mn 0.5⁢ Te

Here, we present a combined magnetometry, muon spin-relaxation (𝜇⁢SR), and neutron-scattering study of the insulating spin glass Zn 0.5 ⁢Mn 0.5 ⁢Te, for which magnetic Mn 2+ and nonmagnetic Zn 2+ ions are randomly distributed on a face-centered cubic lattice. The magnetometry and 𝜇⁢SR results confirm a spin freezing transition around 𝑇 𝑓 ≈ 23 K, with the spin-fluctuation rate decreasing gradually and somewhat inhomogeneously through the sample volume as the temperature decreases toward 𝑇 𝑓 . Characteristic spin-correlation times well above 𝑇 𝑓 are on the order of 10 −10 s, much slower than typically observed in canonical spin glasses but in line with expectations for a cluster spin glass. Using magnetic pair distribution function (mPDF) analysis and reverse Monte Carlo (RMC) modeling of the magnetic diffuse neutron-scattering data, we show that the spin-glass ground state consists of clusters of spins exhibiting short-range-ordered type-III antiferromagnetic correlations with a locally ordered moment of 3.1⁢(1)⁢𝜇 B between nearest-neighbor spins. The type-III correlations decay exponentially as a function of spin separation distance with a correlation length of approximately 5 Å. The diffuse magnetic scattering and corresponding mPDF show no significant changes across 𝑇 𝑓 , indicating that the dynamically fluctuating short-range spin correlations in the paramagnetic state retain the same basic type-III configuration that characterizes the spin-glass state; the only change apparent from the neutron-scattering data is a gradual reduction of the correlation length and locally ordered moment with increasing temperature. Taken together, these results paint a unique and detailed picture of the local magnetic structure and dynamics in Zn 0.5 ⁢Mn 0.5⁢ Te and provide strong evidence that this material is best described as a cluster spin glass. In addition, this work showcases a statistical method for extracting diffuse scattering signals from neutron powder diffraction data, which we developed to facilitate the mPDF and RMC analysis of the neutron data. This method has the potential to be broadly useful for neutron powder diffraction experiments on a variety of materials with short-range atomic or magnetic order.

magnetism↗
NASA NTRS2017

Time Dependence of Collision Probabilities During Satellite Conjunctions

The NASA Conjunction Assessment Risk Analysis (CARA) team has recently implemented updated software to calculate the probability of collision (P (sub c)) for Earth-orbiting satellites. The algorithm can employ complex dynamical models for orbital motion, and account for the effects of non-linear trajectories as well as both position and velocity uncertainties. This “3D P (sub c)” method entails computing a 3-dimensional numerical integral for each estimated probability. Our analysis indicates that the 3D method provides several new insights over the traditional “2D P (sub c)” method, even when approximating the orbital motion using the relatively simple Keplerian two-body dynamical model. First, the formulation provides the means to estimate variations in the time derivative of the collision probability, or the probability rate, R (sub c). For close-proximity satellites, such as those orbiting in formations or clusters, R (sub c) variations can show multiple peaks that repeat or blend with one another, providing insight into the ongoing temporal distribution of risk. For single, isolated conjunctions, R (sub c) analysis provides the means to identify and bound the times of peak collision risk. Additionally, analysis of multiple actual archived conjunctions demonstrates that the commonly used “2D P (sub c)” approximation can occasionally provide inaccurate estimates. These include cases in which the 2D method yields negligibly small probabilities (e.g., P (sub c)) is greater than 10 (sup -10)), but the 3D estimates are sufficiently large to prompt increased monitoring or collision mitigation (e.g., P (sub c) is greater than or equal to 10 (sup -5)). Finally, the archive analysis indicates that a relatively efficient calculation can be used to identify which conjunctions will have negligibly small probabilities. This small-P (sub c) screening test can significantly speed the overall risk analysis computation for large numbers of conjunctions.

Hall, Doyle T.↗
DOE OSTI2017

Model Order Reduction Algorithm for Estimating the Absorption Spectrum

The ab initio description of the spectral interior of the absorption spectrum poses both a theoretical and computational challenge for modern electronic structure theory. Due to the often spectrally dense character of this domain in the quantum propagator’s eigenspectrum for medium-to-large sized systems, traditional approaches based on the partial diagonalization of the propagator often encounter oscillatory and stagnating convergence. Electronic structure methods which solve the molecular response problem through the solution of spectrally shifted linear systems, such as the complex polarization propagator, offer an alternative approach which is agnostic to the underlying spectral density or domain location. This generality comes at a seemingly high computational cost associated with solving a large linear system for each spectral shift in some discretization of the spectral domain of interest. In this work, we present a novel, adaptive solution to this high computational overhead based on model order reduction techniques via interpolation. Model order reduction reduces the computational complexity of mathematical models and is ubiquitous in the simulation of dynamical systems and control theory. The efficiency and effectiveness of the proposed algorithm in the ab initio prediction of X-ray absorption spectra is demonstrated using a test set of challenging water clusters which are spectrally dense in the neighborhood of the oxygen K-edge. On the basis of a single, user defined tolerance we automatically determine the order of the reduced models and approximate the absorption spectrum up to the given tolerance. We also illustrate that, for the systems studied, the automatically determined model order increases logarithmically with the problem dimension, compared to a linear increase of the number of eigenvalues within the energy window. Furthermore, we observed that the computational cost of the proposed algorithm only scales quadratically with respect to the problem dimension.

IDREAM, K-Edge, theoretical, computational, electr↗
DOE OSTI2018

Dissection of the multichannel reaction of acetylene with atomic oxygen: from the global potential energy surface to rate coefficients and branching dynamics

The O( 3 P) + C 2 H 2 reaction is the first step in acetylene oxidation. The accurate kinetic data and the understanding of the reaction dynamics is of great importance. To this end, in this study a full-dimensional global potential energy surface (PES) for the ground triplet state of the O( 3 P) + C 2 H 2 reaction is constructed based on approximately 85,000 ab initio points calculated at the level of explicitly correlated unrestricted coupled cluster single, double, and perturbative triple excitations with the explicitly correlated polarized valence triple zeta basis set (UCCSD(T)-F12b/VTZ-F12). The PES is fit using the permutation invariant polynomial-neural network (PIP-NN) approach with a total root mean square error of 0.21 kcal mol –1 . The key topographic features of the PES, including multiple potential wells and saddle points along different reaction pathways, are well represented by this fit PES. The kinetics and dynamics of the O( 3 P) + C 2 H 2 reaction are investigated using the quasi-classical trajectory (QCT) method. The calculated rate coefficients are in good agreement with experimental data over a wide temperature range, especially when the temperature is lower than 1500 K. The product branch ratio has also been determined, which indicates the H + HCCO channel as the dominant reaction pathway at 298–3000 K, accounting for 80–90% of the overall rate coefficient, in agreement with experimental observations. The dynamics of the reaction is analyzed in detail.

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

Thermofield Theory for Finite-Temperature Electronic Structure

Wave function methods have offered a robust, systematically improvable means to study ground-state properties in quantum many-body systems. Theories like coupled cluster and their derivatives provide highly accurate approximations to the energy landscape at a reasonable computational cost. Analogues of such methods to study thermal properties, though highly desirable, have been lacking because evaluating thermal properties involve a trace over the entire Hilbert space, which is a formidable task. Besides, excited-state theories are generally not as well studied as ground-state ones. In this mini-review, we present an overview of a finite-temperature wave function formalism based on thermofield dynamics to overcome these difficulties. Thermofield dynamics allows us to map the equilibrium thermal density matrix to a pure state, i.e., a single wave function, albeit in an expanded Hilbert space. Ensemble averages become expectation values over this so-called thermal state. Around this thermal state, we have developed a procedure to generalize ground-state wave function theories to finite temperatures. As explicit examples, we highlight formulations of mean-field, configuration interaction, and coupled cluster theories for thermal properties of Fermions in the grand-canonical ensemble. To assess the quality of these approximations, we also show benchmark studies for the one-dimensional Hubbard model, while comparing against exact results. We will see that the thermal methods perform similarly to their ground-state counterparts, while merely adding a prefactor to the asymptotic computational cost. Furthermore, they also inherit all the properties, good or bad, from the ground-state methods, signifying the robustness of our formalism and the scope for future development.

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

CLASH-VLT: DISSECTING THE FRONTIER FIELDS GALAXY CLUSTER MACS J0416.1-2403 WITH ∼800 SPECTRA OF MEMBER GALAXIES

We present VIMOS-Very Large Telescope (VLT) spectroscopy of the Frontier Fields cluster MACS J0416.1-2403 ( z = 0.397). Taken as part of the CLASH-VLT survey, the large spectroscopic campaign provided more than 4000 reliable redshifts over ∼600 arcmin{sup 2}, including ∼800 cluster member galaxies. The unprecedented sample of cluster members at this redshift allows us to perform a highly detailed dynamical and structural analysis of the cluster out to ∼2.2 r {sub 200} (∼4 Mpc). Our analysis of substructures reveals a complex system composed of a main massive cluster ( M {sub 200} ∼ 0.9 × 10{sup 15} M {sub ⊙} and σ{sub V,r200} ∼ 1000 km s{sup −1}) presenting two major features: (i) a bimodal velocity distribution, showing two central peaks separated by Δ V {sub rf} ∼ 1100 km s{sup −1} with comparable galaxy content and velocity dispersion, and (ii) a projected elongation of the main substructures along the NE–SW direction, with a prominent sub-clump ∼600 kpc SW of the center and an isolated BCG approximately halfway between the center and the SW clump. We also detect a low-mass structure at z ∼ 0.390, ∼10′ south of the cluster center, projected at ∼3 Mpc, with a relative line-of-sight velocity of Δ V{sub rf} ∼ −1700 km s{sup −1}. The cluster mass profile that we obtain through our dynamical analysis deviates significantly from the “universal” NFW, being best fit by a Softened Isothermal Sphere model instead. The mass profile measured from the galaxy dynamics is found to be in relatively good agreement with those obtained from strong and weak lensing, as well as with that from the X-rays, despite the clearly unrelaxed nature of the cluster. Our results reveal an overall complex dynamical state of this massive cluster and support the hypothesis that the two main subclusters are being observed in a pre-collisional phase, in agreement with recent findings from radio and deep X-ray data. In this article, we also release the entire redshift catalog of 4386 sources in the field of this cluster, which includes 60 identified Chandra X-ray sources and 105 JVLA radio sources.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Molecular Dynamics Driven by the Many-Body Expansion (MBE-MD)

Here, we present a protocol for classical and nuclear quantum dynamics, in which the energies and forces are generated by the many-body expansion (MBE), and apply it to water clusters using the TTM2.1-F and MB-Pol interaction potentials at various temperatures. We carry out MBE-MD classical and nuclear quantum dynamical simulations in which the energies and forces of the full system are approximated by the 2-, 3-, and 4-body terms of the MBE expansion and compare the average potential and the vibrational density of states (VDOS) with the full simulation, i.e., the one for which no MBE is used. Our results indicate that the thermally averaged potential energy from the MBE up to the 4-body term converges with near identical behavior to the one from the full simulation. The 3-body makes a substantial contribution (~20%) to the energy, whereas the 4- body is necessary for obtaining quantitatively accurate energetics and forces, albeit making a small contribution to each (~2%). We further show that the harmonic frequencies are reproduced to within a few wavenumbers (cm -1 ) at the 4-body level, and that the slowest modes to converge with the MBE rank are those involving the strongest hydrogen bonds. Anharmonicity exacerbates this effect, so that a 4-body description of the energies and forces is needed to achieve accurate anharmonic vibrational frequencies in the hydrogen-bonded OH-stretching region. We also discuss the asymptotic scaling of the MBE-MD protocol with respect to the cost of the underlying potential energy evaluation, suggesting that electronic structure methods that scale at least as N 4 ,N being the size of the system, are needed to result in savings over the traditional full MD simulation. We anticipate that MBE-MD protocol can evolve into a powerful and practical method, which will allow for highly accurate ab initio MD simulations on a much broader range of molecular systems than can be currently handled.

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

Optical properties of Platonic clusters of plasmonic nanoparticles

In the framework of the dipole approximation, we have developed a model of optical properties of a meta-atom consisting of spherical nanoparticles located at the vertices of Platonic solids. Based on the model, we have found and analysed the dynamics of changes in the optical spectra with a change in the length of the edge of a polyhedron. We have observed strong hybridisation and splitting of initially degenerate modes of individual nanoparticles. The obtained results can be used as the basis for the development of an optical nanosensor, which can determine the change in the chemical and biological composition of the environment. (paper)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Subhalos in Galaxy Clusters: Coherent Accretion and Internal Orbits

Subhalo dynamics in galaxy cluster host halos govern the observed distribution and properties of cluster member galaxies. We use the IllustrisTNG simulation to investigate the accretion and orbits of subhalos found in cluster-size halos. We find that the median change in the major axis direction of cluster-size host halos is approximately 80° between a ∼ 0.1 and the present day. We identify coherent regions in the angular distribution of subhalo accretion, and ∼68% of accreted subhalos enter their host halo through ∼38% of the surface area at the virial radius. The majority of galaxy clusters in the sample have ∼2 such coherent regions. We further measure angular orbits of subhalos with respect to the host major axis and use a clustering algorithm to identify distinct orbit modes with varying oscillation timescales. The orbit modes correlate with subhalo accretion conditions. Subhalos in orbit modes with shorter oscillations tend to have lower peak masses and accretion directions somewhat more aligned with the major axis. One orbit mode, exhibiting the least oscillatory behavior, largely consists of subhalos that accrete near the plane perpendicular to the host halo major axis. Our findings are consistent with expectations from inflow from major filament structures and internal dynamical friction: most subhalos accrete through coherent regions, and more massive subhalos experience fewer orbits after accretion. Our work offers a unique quantification of subhalo dynamics that can be connected to how the intracluster medium strips and quenches cluster galaxies.

galaxy clusters↗
DOE OSTI2016

Multiscale modelling of evolving foams

We present a set of multi-scale interlinked algorithms to model the dynamics of evolving foams. These algorithms couple the key effects of macroscopic bubble rearrangement, thin film drainage, and membrane rupture. For each of the mechanisms, we construct consistent and accurate algorithms, and couple them together to work across the wide range of space and time scales that occur in foam dynamics. These algorithms include second order finite difference projection methods for computing incompressible fluid flow on the macroscale, second order finite element methods to solve thin film drainage equations in the lamellae and Plateau borders, multiphase Voronoi Implicit Interface Methods to track interconnected membrane boundaries and capture topological changes, and Lagrangian particle methods for conservative liquid redistribution during rearrangement and rupture. We derive a full set of numerical approximations that are coupled via interface jump conditions and flux boundary conditions, and show convergence for the individual mechanisms. We demonstrate our approach by computing a variety of foam dynamics, including coupled evolution of three-dimensional bubble clusters attached to an anchored membrane and collapse of a foam cluster.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1995

Fast Multipole Methods for Three-Dimensional N-body Problems

We are developing computational tools for the simulations of three-dimensional flows past bodies undergoing arbitrary motions. High resolution viscous vortex methods have been developed that allow for extended simulations of two-dimensional configurations such as vortex generators. Our objective is to extend this methodology to three dimensions and develop a robust computational scheme for the simulation of such flows. A fundamental issue in the use of vortex methods is the ability of employing efficiently large numbers of computational elements to resolve the large range of scales that exist in complex flows. The traditional cost of the method scales as Omicron (N(sup 2)) as the N computational elements/particles induce velocities at each other, making the method unacceptable for simulations involving more than a few tens of thousands of particles. In the last decade fast methods have been developed that have operation counts of Omicron (N log N) or Omicron (N) (referred to as BH and GR respectively) depending on the details of the algorithm. These methods are based on the observation that the effect of a cluster of particles at a certain distance may be approximated by a finite series expansion. In order to exploit this observation we need to decompose the element population spatially into clusters of particles and build a hierarchy of clusters (a tree data structure) - smaller neighboring clusters combine to form a cluster of the next size up in the hierarchy and so on. This hierarchy of clusters allows one to determine efficiently when the approximation is valid. This algorithm is an N-body solver that appears in many fields of engineering and science. Some examples of its diverse use are in astrophysics, molecular dynamics, micro-magnetics, boundary element simulations of electromagnetic problems, and computer animation. More recently these N-body solvers have been implemented and applied in simulations involving vortex methods. Koumoutsakos and Leonard (1995) implemented the GR scheme in two dimensions for vector computer architectures allowing for simulations of bluff body flows using millions of particles. Winckelmans presented three-dimensional, viscous simulations of interacting vortex rings, using vortons and an implementation of a BH scheme for parallel computer architectures. Bhatt presented a vortex filament method to perform inviscid vortex ring interactions, with an alternative implementation of a BH scheme for a Connection Machine parallel computer architecture.

Koumoutsakos, P.↗
NASA NTRS1991

The small-scale clustering properties of dwarf galaxies

Two results on the small-scale clustering properties of dwarf galaxies are reported, which were identified in the vicinity of early-type Shapley-Ames galaxies on high-resolution photographic plates. The first result indicates that dwarf galaxies display the same trend of stronger clustering toward earlier morphological type on small scales as their giant counterparts on larger scales. It is suggested that early-type dwarfs can be used as dynamical probes of dark halos around early-type giant galaxies and as tracers of the dynamical evolution of such halos in dense environments. The second result pertains to the trend of increasing early-type dwarf frequency per early-type giant with environment richness previously established for rich groups. It is found that a minimum value of isolated early-type galaxies is approximately 0.25, as compared to a maximum of approximately 8 in rich environments like the Virgo Cluster.

Vader, J. P.↗