Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “CLUSTER MODEL”

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 415 records · Page 23

Shadow Molecular Dynamics and Atomic Cluster Expansions for Flexible Charge Models

Here, a shadow molecular dynamics scheme for flexible charge models is presented where the shadow Born–Oppenheimer potential is derived from a coarse-grained approximation of range-separated density functional theory. The interatomic potential, including the atomic electronegativities and the charge-independent short-range part of the potential and force terms, is modeled by the linear atomic cluster expansion (ACE), which provides a computationally efficient alternative to many machine learning methods. The shadow molecular dynamics scheme is based on extended Lagrangian (XL) Born–Oppenheimer molecular dynamics (BOMD). XL-BOMD provides stable dynamics while avoiding the costly computational overhead associated with solving an all-to-all system of equations, which normally is required to determine the relaxed electronic ground state prior to each force evaluation. To demonstrate the proposed shadow molecular dynamics scheme for flexible charge models using atomic cluster expansion, we emulate the dynamics generated from self-consistent charge density functional tight-binding (SCC-DFTB) theory using a second-order charge equilibration (QEq) model. The charge-independent potentials and electronegativities of the QEq model are trained for a supercell of uranium oxide (UO 2 ) and a molecular system of liquid water. The combined ACE+XL-QEq molecular dynamics simulations are stable over a wide range of temperatures both for the oxide and for the molecular systems and provide a precise sampling of the Born–Oppenheimer potential energy surfaces. Accurate ground Coulomb energies are produced by the ACE-based electronegativity model during an NVE simulation of UO 2 , predicted to be within 1 meV of those from SCC-DFTB on average during comparable simulations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Joint Analysis of X-Ray and Sunyaev-Zel'Dovich Observations of Galaxy Clusters Using an Analytic Model of the Intracluster Medium

We perform a joint analysis of X-ray and Sunyaev-Zel'dovich effect data using an analytic model that describes the gas properties of galaxy clusters. The joint analysis allows the measurement of the cluster gas mass fraction profile and Hubble constant independent of cosmological parameters. Weak cosmological priors are used to calculate the overdensity radius within which the gas mass fractions are reported. Such an analysis can provide direct constraints on the evolution of the cluster gas mass fraction with redshift. We validate the model and the joint analysis on high signal-to-noise data from the Chandra X-ray Observatory and the Sunyaev-Zel'dovich Array for two clusters, A2631 and A2204.

Galaxies↗

Production of recycled pulsars in globular clusters via two-body tidal capture

The numbers of binary systems containing a neutron star that are formed in the cores of globular clusters via two-body tidal capture are investigated. The subsequent evolution of the resulting binary systems is examined. Monte Carlo techniques are employed to trace the histories of many neutron stars in the core of a model globular cluster. Various scenarios for mass stripping during capture and orbital circularization, and the effects of magnetic braking on the circularized orbits are explored. Calculations are carried out for two model clusters: 47 Tuc and Omega Cen. Distributions of orbital periods for newly formed binaries containing a neutron star in orbit with either a main-sequence star or a giant are presented. Each binary is followed until contact between the captured star and the neutron star is established. It is argued that these systems should be regarded as potential progenitors of 'recycled' pulsars.

Di Stefano, R↗

A Search for CMB Decrements Towards Distant Cluster Candidates PC 1643+4631 and VLA 1312+4237 at 28.5 GHz

Recently, Jones et al. used the Ryle telescope, operating at a frequency of 15 GHz, to detect a flux decrement in the direction of the quasar pair PC 1643+461A,B. They interpreted this signal as the Sunyaev-Zel'dovich effect (SZE) produced by a distant cluster of galaxies. In the course of an effort to measure CMB anisotropies using the VLA at 8.4 GHz, Richards et al. (1997) detected a similar, but smaller, decrement which we refer to as VLA 1312+4237. They also proposed that this signal might be explained as the SZE signal of a distant galaxy cluster. We report observations in the direction of these claimed sources with the Berkeley Illinois Maryland Association (BIMA) interferometer operating at 28.5 GHz. We find no evidence for SZE emission in the direction of either of the claimed sources. In the case of PC 1643+4631, the BIMA data are inconsistent with the cluster emission model proposed by Jones et al. at greater than 99.99% confidence. Together with published x-ray and optical searches, these results make a compelling case against the existence of a massive cluster in the direction of PC 1643+4631. Because of the different scales to which the VLA and BIMA instruments are sensitive, the BIMA observations are not as constraining for the VLA 1312+4237 source. The BIMA data are inconsistent with the cluster model proposed by Richards et al. (1997) at approximately 80% confidence.

Holzapfel, W. L.↗

Nonpolytropic model for the Coma Cluster

In this article we demonstrate, for the first time, how a physically motivated static model for both the gas and galaxies in the Coma Cluster of galaxies can jointly fit all available X-ray and optical imaging and spectroscopic data. The principal assumption of this nonpolytropic model (Cavaliere & Fusco-Femiano 1981, hereafter CFF), is that the intracluster gas temperature is proportional to the square of the galaxy velocity dispersion everywhere throughout the cluster; no other assumption about the gas temperature distribution is required. After demonstrating that the CFF nonpolytropic model is an adequate representation of the gas and galaxy distributions, the radial velocity dispersion profile, and the gas temperature distribution, we derive the following information about the Coma Cluster: 1. The central temperature is about 9 keV and the central density is 2.8 x 10(exp -3)/cm(exp 3) for the X-ray emitting plasma; 2. The binding mass of the cluster is approximately 2 x 10(exp 15) solar mass within 5 Mpc for (H(sub 0) = 50 km/sec/Mpc), with a mass-to-light ratio of approximately 160 solar mass/solar luminosity; 3. The contribution of the gas to the total virial mass increases with distance from the cluster center, and we estimate that this ratio is no greater than approximately 50% within 5 Mpc. The ability of the CFF nonpolytropic model to describe the current X-ray and optical data for the Coma Cluster suggests that a significant fraction of the thermal energy contained in the hot gas in this as well as other rich galaxy clusters may have come from the interaction between the galaxies and the ambient cluster medium. interaction between the galaxies and the ambient cluster medium.

Fusco-Femiano, R.↗

Shadow molecular dynamics and atomic cluster expansions for flexible charge models

A shadow molecular dynamics scheme for flexible charge models is presented, where the shadow Born-Oppenheimer potential is derived from a coarse-grained approximation of range-separated density functional theory. The interatomic potential, including the atomic electronegativities and the charge-independent short-range part of the potential and force terms, are modeled by the linear atomic cluster expansion (ACE), which provides a computationally efficient alternative to many machine learning methods. The shadow molecular dynamics scheme is based on extended Lagrangian (XL) Born-Oppenheimer molecular dynamics (BOMD) [Eur. Phys. J. B 94, 164 (2021)]. XL-BOMD provides a stable dynamics, while avoiding the costly computational overhead associated with solving an all-to-all system of equations, which normally is required to determine the relaxed electronic ground state prior to each force evaluation. To demonstrate the proposed shadow molecular dynamics scheme for flexible charge models using the atomic cluster expansion, we emulate the dynamics generated from self-consistent charge density functional tight-binding (SCC-DFTB) theory using a second-order charge equilibration (QEq) model. The charge-independent potentials and electronegativities of the QEq model are trained for a supercell of uranium oxide (UO2) and a molecular system of liquid water. The combined ACE + XL-QEq dynamics are stable over a wide range of temperatures both for the oxide and the molecular systems, and provide a precise sampling of the Born-Oppenheimer potential energy surfaces. Accurate ground Coulomb energies are produced by the ACE-based electronegativity model during an NVE simulation of UO 2 , predicted to be within 1 meV of those from SCC-DFTB on average during comparable simulations.

74 ATOMIC AND MOLECULAR PHYSICS↗

Modeling Redshift-space Clustering with Abundance Matching

Abstract We explore the degrees of freedom required to jointly fit projected and redshift-space clustering of galaxies selected in three bins of stellar mass from the Sloan Digital Sky Survey Main Galaxy Sample (SDSS MGS) using a subhalo abundance matching (SHAM) model. We employ emulators for relevant clustering statistics in order to facilitate our analysis, leading to large speed gains with minimal loss of accuracy. We are able to simultaneously fit the projected and redshift-space clustering of the two most massive galaxy samples that we consider with just two free parameters: scatter in stellar mass at fixed SHAM proxy, and the dependence of the SHAM proxy on dark matter halo concentration. We find some evidence for models that include velocity bias, but including orphan galaxies improves our fits to the lower-mass samples significantly. We also model the clustering signals of specific star formation rate (sSFR) selected samples using conditional abundance matching (CAM). We obtain acceptable fits to projected and redshift-space clustering as a function of sSFR and stellar mass using two CAM variants, although the fits are worse than for stellar-mass-selected samples alone. By incorporating nonunity correlations between the CAM proxy and sSFR, we are able to resolve previously identified discrepancies between CAM predictions and SDSS observations of the environmental dependence of quenching for isolated central galaxies.

79 ASTRONOMY AND ASTROPHYSICS↗

Vibronic coupling in molecular crystals: A Franck-Condon Herzberg-Teller model of H-aggregate fluorescence based on quantum chemical cluster calculations

Here, we present a general approach to treating vibronic coupling in molecular crystals based on atomistic simulations of large clusters. Such clusters comprise model aggregates treated at the quantum chemical level embedded within a realistic environment treated at the molecular mechanics level. As we calculate ground and excited state equilibrium geometries and vibrational modes of model aggregates, our approach is able to capture effects arising from coupling to intermolecular degrees of freedom, absent from existing models relying on geometries and normal modes of single molecules. Using the geometries and vibrational modes of clusters, we are able to simulate the fluorescence spectra of aggregates for which the lowest excited state bears negligible oscillator strength (as is the case, e.g., ideal H-aggregates) by including both Franck-Condon (FC) and Herzberg-Teller (HT) vibronic transitions. The latter terms allow the adiabatic excited state of the cluster to couple with vibrations in a perturbative fashion via derivatives of the transition dipole moment along nuclear coordinates. While vibronic coupling simulations employing FC and HT terms are well established for single-molecules, to our knowledge this is the first time they are applied to molecular aggregates. Here, we apply this approach to the simulation of the low-temperature fluorescence spectrum of para-distyrylbenzene single-crystal H-aggregates and draw comparisons with coarse-grained Frenkel-Holstein approaches previously extensively applied to such systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

On the equilibrium distribution of the elements in the gas in the Coma cluster

The equilibrium distribution of the elements for the case of an isothermal plasma in a gravitational potential well is given governing equations and a rapidly converging iterative procedure for their solution. Results are applied to the intracluster gas in clusters of galaxies, and implications are identified for the controversy as to whether iron nuclei settle out in the cluster core, and with respect to the interpretation of distribution size and shape for continuum X-ray emission from clusters. A model for the Coma cluster is considered in which the intracluster gas is in equilibrium and partially processed. It is found that this model can account for the angular dependence of the continuum X-ray emission from Coma, as determined by Einstein Observatory observations, at least as well as the uniform composition model.

Abramopoulos, F.↗

Models Predicting Void Swelling Incubation Dose as a function of Irradiation Conditions

This report documents the development status of a cluster dynamics model to predict microstructural evolution in high-energy particle irradiated ferritic martensitic steels, in which we have recently implemented a model for the biased absorption of self-interstitial defects at small, nanometer-sized cavities that can provide an intrinsic nucleation barrier for the onset of void swelling. Following a brief review of the cluster dynamics modeling approach, this report describes the atomistic modeling results that demonstrate this phenomenon, along with an investigation on the influence of helium content on modifying the intrinsic cavity bias. Subsequently, the report describes the initial results obtained from our cluster dynamics modeling of the incorporation of this cavity bias on the predicted swelling incubation dose as a function of irradiation temperature. Finally, the report concludes with a discussion of the remaining modeling research required to fully validate the model for application to high dose neutron or ion beam irradiation exposure conditions.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Application of cluster-plus-glue-atom model to barrierless Cu–Ni–Ti and Cu–Ni–Ta films

To improve the thermal stability of copper and avoid its diffusion into surrounding dielectrics or interfacial reactions with them, the authors applied the cluster-plus-glue-atom model to investigate barrierless Cu–Ni–M (M = Ti or Ta) seed layers. The dissolution of the third element (Ti or Ta) in the Cu lattice with the aid of Ni significantly improved the thermal stability of the Cu seed layer. The appropriate M/Ni (M = Ti or Ta) ratio was selected to obtain a low resistivity: the resistivity was as low as 2.5 μΩ cm for the (Ti{sub 1.5/13.5}Ni{sub 12/13.5}){sub 0.3}Cu{sub 99.7} film and 2.8 μΩ cm for the (Ta{sub 1.1/13.1}Ni{sub 12/13.1}){sub 0.4}Cu{sub 99.6} film after annealing at 500 °C for 1 h. After annealing at 500 °C for 40 h, the two films remained stable without forming a Cu{sub 3}Si compound. The authors confirmed that the range of applications of the cluster-plus-glue-atom model could be extended. Therefore, a third element M with negative enthalpies of mixing with both Cu and Ni could be selected, under the premise that the mixing enthalpy of M–Ni is more negative than that of M–Cu.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Cluster-continuum quantum mechanical models to guide the choice of anions for Li{sup +}-conducting ionomers

A quantum-mechanical investigation on Li poly(ethylene oxide)-based ionomers was performed in the cluster-continuum solvation model (CCM) that includes specific solvation in the first shell surrounding the cation, all surrounded by a polarizable continuum. A four-state model, including a free Li cation, Li{sup +}-anion pair, triple ion, and quadrupole was used to represent the states of Li{sup +} within the ionomer in the CCM. The relative energy of each state was calculated for Li{sup +} with various anions, with dimethyl ether representing the ether oxygen solvation. The population distribution of Li{sup +} ions among states was estimated by applying Boltzmann statistics to the CCM energies. Entropy difference estimates are needed for populations to better match the true ionomer system. The total entropy change is considered to consist of four contributions: translational, rotational, electrostatic, and solvent immobilization entropies. The population of ion states is reported as a function of Bjerrum length divided by ion-pair separation with/without entropy considered to investigate the transition between states. Predicted concentrations of Li{sup +}-conducting states (free Li{sup +} and positive triple ions) are compared among a series of anions to indicate favorable features for design of an optimal Li{sup +}-conducting ionomer; the perfluorotetraphenylborate anion maximizes the conducting positive triple ion population among the series of anions considered.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Constraining quenching time-scales in galaxy clusters by forward-modelling stellar ages and quiescent fractions in projected phase space

We forward-model mass-weighted stellar ages (MWAs) and quiescent fractions (f Q ) in projected phase space (PPS), using data from the Sloan Digital Sky Survey, to jointly constrain an infall quenching model for galaxies in log (M vir /M ⊙ ) > 14 galaxy clusters at z ∼ 0. We find the average deviation in MWA from the MWA–M⋆ relation depends on position in PPS, with a maximum difference between the inner cluster and infalling interloper galaxies of ∼1 Gyr. Our model employs infall information from N-body simulations and stochastic star-formation histories from the universemachine model. We find total quenching times of t Q = 3.7 ± 0.4 Gyr and t Q = 4.0 ± 0.2 Gyr after first pericentre, for 9 < log (M ⋆ /M ⊙ ) < 10 and 10 < log (M ⋆ /M ⊙ ) < 10.5 galaxies, respectively. By using MWAs, we break the degeneracy in time of quenching onset and time-scale of star formation rate (SFR) decline. We find that time of quenching onset relative to pericentre is $t_{\mathrm{delay}}=3.5^{+0.6}_{-0.9}$ Gyr and $t_{\mathrm{delay}}=-0.3^{+0.8}_{-1.0}$ Gyr for 9 < log (M ⋆ /M ⊙ ) < 10 and 10 < log (M ⋆ /M ⊙ ) < 10.5 galaxies, respectively, and exponential SFR suppression time-scales are τenv ≤ 1.0 Gyr for 9 < log (M ⋆ /M ⊙ ) < 10 galaxies and τenv ∼ 2.3 Gyr for 10 < log (M ⋆ /M ⊙ ) < 10.5 galaxies. Stochastic star formation histories remove the need for rapid infall quenching to maintain the bimodality in the SFR of cluster galaxies; the depth of the green valley prefers quenching onsets close to first pericentre and a longer quenching envelope, in slight tension with the MWA-driven results. Taken together these results suggest that quenching begins close to, or just after pericentre, but the time-scale for quenching to be fully complete is much longer and therefore ram-pressure stripping is not complete on first pericentric passage.

79 ASTRONOMY AND ASTROPHYSICS↗

Role of Defects in Swelling and Creep of Irradiated SiC

Silicon carbide is a promising cladding material because of its high strength and relatively good corrosion resistance. However, SiC is brittle and therefore SiC-based components need to be carefully designed to avoid cracking and failure by fracture. In design of SiC-based composites for nuclear reactor applications it is essential to take into account how mechanical properties are affected by radiation and temperature, or in other words, what strains and stresses develop in this material due to environmental conditions. While thermal strains in SiC can be predicted using classical theories, radiation-induced strains are much less understood. In particular, it is critical to correctly account for radiation swelling and radiation creep, which contribute significantly to dimensional instability of SiC under radiation. Swelling typically increases logarithmically with radiation dose and saturates at relatively low doses (damage levels of a few dpa). Consequently, swelling-induced stresses are likely to develop within a few months of operation of a reactor. Radiation-induced volume swelling in SiC can be as high as 2%, which is significantly higher than the cracking strain of 0.1% in SiC. Swelling-induced strains will lead to enormous stresses and fracture, unless these stresses can be relaxed via some other mechanism. An effective way to achieve stress relaxation is via radiation creep. Although it has been hypothesized that both radiation swelling and radiation creep are driven by formation of defect clusters, existing models for swelling and creep in SiC are limited by the lack of understanding of specific defects that form due to radiation in the range of temperatures relevant to fuel cladding in light water reactors (LWRs) (<1000°C). For example, defects that can be detected with traditional transmission electron microscopy (TEM) techniques account only for 10-45% of the swelling measured in irradiated SiC. Here, we have undertaken an integrated experimental and modeling effort to discover the previously invisible defects in irradiated SiC and to determine the contributions of these defects to radiation swelling. Knowledge of the most stable defect structures and the rate controlling processes during defect evolution is essential for development of predictive models for swelling and creep as a function of temperature and radiation dose. This research has been enabled by state-of-the-art imaging techniques, such as the aberration corrected scanning transmission electron microscopy (STEM) (FEI TITAN) closely coupled with multi-scale models of stable defect clusters and their evolution.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

An Assessment of Antarctic Circumpolar Current and Southern Ocean Meridional Overturning Circulation During 1958-2007 in a Suite of Interannual CORE-2 Simulations

In the framework of the second phase of the Coordinated Ocean-ice Reference Experiments (CORE-II), we present an analysis of the representation of the Antarctic Circumpolar Current (ACC) and Southern Ocean Meridional Overturning Circulation (MOC) in a suite of seventeen global ocean-sea ice models. We focus on the mean, variability and trends of both the ACC and MOC over the 19582007 period, and discuss their relationship with the surface forcing. We aim to quantify the degree of eddy saturation and eddy compensation in the models participating in CORE-II, and compare our results with available observations, previous fine-resolution numerical studies and theoretical constraints. Most models show weak ACC transport sensitivity to changes in forcing during the past five decades, and they can be considered to be in an eddy saturated regime. Larger contrasts arise when considering MOC trends, with a majority of models exhibiting significant strengthening of the MOC during the late 20th and early 21st century. Only a few models show a relatively small sensitivity to forcing changes, responding with an intensified eddy-induced circulation that provides some degree of eddy compensation, while still showing considerable decadal trends. Both ACC and MOC interannual variability are largely controlled by the Southern Annular Mode (SAM). Based on these results, models are clustered into two groups. Models with constant or two-dimensional (horizontal) specification of the eddy-induced advection coefficient show larger ocean interior decadal trends, larger ACC transport decadal trends and no eddy compensation in the MOC. Eddy-permitting models or models with a three-dimensional time varying show smaller changes in isopycnal slopes and associated ACC trends, and partial eddy compensation. As previously argued, a constant in time or space is responsible for a poor representation of mesoscale eddy effects and cannot properly simulate the sensitivity of the ACC and MOC to changing surface forcing. Evidence is given for a larger sensitivity of the MOC as compared to the ACC transport, even when approaching eddy saturation. Future process studies designed for disentangling the role of momentum and buoyancy forcing in driving the ACC and MOC are proposed.

Southern Ocean dynamics↗

Gravitational lens models of arcs in clusters

It is now well established that the luminous arcs discovered in clusters of galaxies, in particular those in Abell 370 and Cluster 2244-02, are produced by gravitational lensing of background sources. The arcs are modeled and constraints are placed on the distribution of the mass in the clusters and the shape and size of the sources. The models require, as expected, a large amount of dark matter in the clusters and a mass-to blue-light ratio for the cluster which exceeds 100 solar mass/solar luminosity and could be as high as 1000 solar mass/solar luminosity depending on cosmological parameters and the distribution of the dark matter. Furthermore, it is found that in the case of the arc in A370 the dark matter must have a different distribution than the luminous galaxies, while for the arc in Cl 2244 the dark matter can have a distribution similar to that of the light matter (galaxies) or a separate distribution.

Bergmann, Anton G.↗

IR emission and UV extinction in two open clusters

Recent models of interstellar extinction have shown the importance of understanding both the UV and IR properties of interstellar dust grains. IRAS data have shown variations in 60 and 100 micron emissions presumably due to the presence of IR cirrus, while recent observations in the UV by Fitzpatrick and Massa have identified components in the UV extinction curve which vary in different star regions. A Draine and Anderson model connects these results by proposing that different size variations in interstellar grains would cause distinct changes in both the IR emission and the UV extinction. In order to test this model it is necessary to make observations in well defined locations away from peculiar extinction regions. In the infrared this means looking away from the galactic plane so as to limit non-local sources of IR radiation. Two open clusters that are out of the galactic plane and which contain a number of late B and early A stars suitable for UV extinction studies, and whose IRAS data show variations in the 60/100 micron ratio were studied. Based on the Drain and Anderson model, variations were expected in their UV extinction curves that correlate with the IR cirrus emission.

Hackwell, James A.↗

Coupled cluster and dislocation dynamics modeling of microstructure evolution in irradiated materials

We develop here a coupled cluster and dislocation dynamics framework to study the microstructure evolution of irradiated materials. The framework not only accounts for the three dimensional diffusion of radiation-generated clusters, but also their interaction with dislocation networks and the resultant climb motion of discrete dislocations within finite crystals. The framework is solved with a superposition solution scheme, and is applied to investigate the evolution of the irradiation-induced dislocation loops in zirconium (Zr), considering the effects of various bias factors including the diffusion anisotropy difference (DAD) of interstitials and interstitial clusters, the dislocation bias of defects to discrete dislocation segments, and the production bias of defects from the radiation cascade. We find that the DAD is the most critical factor influencing the kinetics of the loop evolution in Zr, while the recombination/interaction of mobile defects can induce a strong spatial dependence of the loop evolution together with the DAD. Here, the method is also adopted to study the evolution of interstitial $\langle$a$\rangle$ and vacancy $\langle$c$\rangle$ dislocation loop ensembles consistent with the microstructure observed during irradiation-induced growth of Zr. Our findings not only reveal the spatial dependence of the size and ellipticity of the dislocation loops, but also suggest a limit on the anisotropy factor of interstitials to reproduce the co-growth of $\langle$a$\rangle$ and $\langle$c$\rangle$ loops in zirconium, in good agreement with experimental observations and other simulation results.

Bias factors↗