SEARCH · Engineering Papers

Results for “PHASE SPACE”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 379 records · Page 21

All sources NASA NTRS DOE OSTI

DOE OSTI2017

Creating Binary Cu–Bi Compounds via High-Pressure Synthesis: A Combined Experimental and Theoretical Study

Exploration beyond the known phase space of thermodynamically stable compounds into the realm of metastable materials is a frontier of materials chemistry. The application of high pressure in experiment and theory provides a powerful vector by which to explore this uncharted phase space, allowing discovery of complex new structures and bonding in the solid state. We harnessed this approach for the Cu–Bi system, where the realization of new phases offers potential for exotic properties such as superconductivity. This potential is due to the presence of bismuth, which, by virtue of its status as one of the heaviest stable elements, forms a critical component in emergent materials such as superconductors and topological insulators. To fully investigate and understand the Cu–Bi system, we welded theoretical predictions with experiment to probe the Cu–Bi system under high pressures. By employing the powerful approach of in situ X-ray diffraction in a laser-heated diamond anvil cell (LHDAC), we thoroughly explored the high-pressure and high-temperature (high-PT) phase space to gain insight into the formation of intermetallic compounds at these conditions. We employed density functional theory (DFT) calculations to calculate a pressure versus temperature phase diagram, which correctly predicts that CuBi is stabilized at lower pressures than Cu11Bi7, and allows us to uncover the thermodynamic contributions responsible for the stability of each phase. Detailed comparisons between the NiAs structure type and the two high-pressure Cu–Bi phases, Cu11Bi7 and CuBi, reveal the preference for elemental segregation within the Cu–Bi phases, and highlight the unique channels and layers formed by ordered Cu vacancies. The electron localization function from DFT calculations account for the presence of these “voids” as a manifestation of the lone pair orientation on the Bi atoms. Our study demonstrates the power of joint experimental–computational work in exploring the chemistry occurring at high-PT conditions. The existence of multiple high-pressure-stabilized phases in the Cu–Bi binary system, which can be readily identified with in situ techniques, offers promise for other systems in which no ambient pressure phases are known to exist.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Galactic echoes

Gaia has revealed a variety of substructures in the phase space of stars in the Solar neighbourhood, including the vertical ‘Snail’ $(z, v_z)$ in space. Such substructures are often interpreted as the incompletely phase-mixed response of the disc stars to a single perturbation, such as an impulsive encounter with a satellite galaxy. In this paper, we consider the possibility that such structures contain manifestations of phase-space echoes. First established in plasma physics in the 1960s, echoes arise when a collisionless system is perturbed twice: the macroscopic responses to both perturbations mix to small scales in phase space, whereupon they couple non-linearly, producing a third macroscopic ‘echo’ response without the need for a third perturbation. We derive the galactic analogue of the plasma echo theory using angle-action variables and apply it to a one-dimensional model of vertical motion in the Milky Way. We verify the predicted echo behaviour using idealized test particle simulations, both with and without the inclusion of diffusion through orbital scattering off molecular clouds. While we conclude that the Gaia Snail itself is unlikely a (pure) echo effect, the basic physics we uncover is sufficiently generic that we expect phase-space echoes to be common in disc galaxies.

Galaxy: evolution↗

DOE OSTI2021

On the Effect of Beating during Nonlinear Frequency Chirping

Spectroscopic analyses of energetic particle (EP) driven bursts of MHD fluctuations in magnetically confined plasmas often exhibit chirps that occur simultaneously in groups of two or more. While the superposition of oscillations at multiple frequencies necessarily causes beating in the signal acquired by a localized external probe, self-consistent hybrid simulations of chirping EP modes in a JT-60U tokamak plasma have demonstrated the possibility of global beating, where the mode's electromagnetic field vanishes globally between beats and reappears with opposite phase [Bierwage et al., Nucl. Fusion 57, 016036 (2017)]. This implies that there can be a single coherent field mode that oscillates at multiple frequencies simultaneously when it is resonantly driven by multiple density waves in EP phase space. Conversely, this means that the EP density waves are mutually coupled and interfere with each other via the jointly driven field, a mechanism ignored in some theories of chirping. In this thesis-style treatise, we study the role of field pulsations in general and beating in particular using the Hamiltonian guiding center orbit-following code ORBIT with a reduced wave-particle interaction model in realistic geometry. Beating is found to drive the evolution of EP phase space structures. Additionally, a key mechanism is the pulsation of effective phase space islands combined with the alternation of their effective O- and X-points due to phase jumps between each beat. Observations: (1) Beating causes density wave fronts to advance radially in a pulsed manner and the resulting chirps become staircase-like. (2) The pulsations facilitate convective transfer of material between neighboring layers of phase space density waves. On the one hand, this may inhibit the early detachment of solitary phase space vortices. On the other hand, it facilitates the accumulation of hole and clump fragments into larger structures. (3) Long-range chirping is observed when massive holes or clumps detach and drift away from the turbulent belt around the seed resonance. It is remarkable that the detached vortices remain robust and,on average, maintain their concentric nested layers while being visibly perturbed by the field's continued beating.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2015

Impact of the differential fluence distribution of brachytherapy sources on the spectroscopic dose-rate constant

Purpose: To investigate why dose-rate constants for {sup 125}I and {sup 103}Pd seeds computed using the spectroscopic technique, Λ{sub spec}, differ from those computed with standard Monte Carlo (MC) techniques. A potential cause of these discrepancies is the spectroscopic technique’s use of approximations of the true fluence distribution leaving the source, φ{sub full}. In particular, the fluence distribution used in the spectroscopic technique, φ{sub spec}, approximates the spatial, angular, and energy distributions of φ{sub full}. This work quantified the extent to which each of these approximations affects the accuracy of Λ{sub spec}. Additionally, this study investigated how the simplified water-only model used in the spectroscopic technique impacts the accuracy of Λ{sub spec}. Methods: Dose-rate constants as described in the AAPM TG-43U1 report, Λ{sub full}, were computed with MC simulations using the full source geometry for each of 14 different {sup 125}I and 6 different {sup 103}Pd source models. In addition, the spectrum emitted along the perpendicular bisector of each source was simulated in vacuum using the full source model and used to compute Λ{sub spec}. Λ{sub spec} was compared to Λ{sub full} to verify the discrepancy reported by Rodriguez and Rogers. Using MC simulations, a phase space of the fluence leaving the encapsulation of each full source model was created. The spatial and angular distributions of φ{sub full} were extracted from the phase spaces and were qualitatively compared to those used by φ{sub spec}. Additionally, each phase space was modified to reflect one of the approximated distributions (spatial, angular, or energy) used by φ{sub spec}. The dose-rate constant resulting from using approximated distribution i, Λ{sub approx,i}, was computed using the modified phase space and compared to Λ{sub full}. For each source, this process was repeated for each approximation in order to determine which approximations used in the spectroscopic technique affect the accuracy of Λ{sub spec}. Results: For all sources studied, the angular and spatial distributions of φ{sub full} were more complex than the distributions used in φ{sub spec}. Differences between Λ{sub spec} and Λ{sub full} ranged from −0.6% to +6.4%, confirming the discrepancies found by Rodriguez and Rogers. The largest contribution to the discrepancy was the assumption of isotropic emission in φ{sub spec}, which caused differences in Λ of up to +5.3% relative to Λ{sub full}. Use of the approximated spatial and energy distributions caused smaller average discrepancies in Λ of −0.4% and +0.1%, respectively. The water-only model introduced an average discrepancy in Λ of −0.4%. Conclusions: The approximations used in φ{sub spec} caused discrepancies between Λ{sub approx,i} and Λ{sub full} of up to 7.8%. With the exception of the energy distribution, the approximations used in φ{sub spec} contributed to this discrepancy for all source models studied. To improve the accuracy of Λ{sub spec}, the spatial and angular distributions of φ{sub full} could be measured, with the measurements replacing the approximated distributions. The methodology used in this work could be used to determine the resolution that such measurements would require by computing the dose-rate constants from phase spaces modified to reflect φ{sub full} binned at different spatial and angular resolutions.

07 ISOTOPE AND RADIATION SOURCES↗

DOE OSTI2014

Toward multi-differential cross sections: measuring two angularities on a single jet

The analytic study of differential cross sections in QCD has typically focused on individual observables, such as mass or thrust, to great success. Here, we present a first study of double differential jet cross sections considering two recoil-free angularities measured on a single jet. By analyzing the phase space defined by the two angularities and using methods from soft-collinear effective theory, we prove that the double differential cross section factorizes at the boundaries of the phase space. We also show that the cross section in the bulk of the phase space cannot be factorized using only soft and collinear modes, excluding the possibility of a global factorization theorem in soft-collinear effective theory. Nevertheless, we are able to define a simple interpolation procedure that smoothly connects the factorization theorem at one boundary to the other. We present an explicit example of this at next-to-leading logarithmic accuracy and show that the interpolation is unique up to α s 4 order in the exponent of the cross section, under reasonable assumptions. This is evidence that the interpolation is sufficiently robust to account for all logarithms in the bulk of phase space to the accuracy of the boundary factorization theorem. We compare our analytic calculation of the double differential cross section to Monte Carlo simulation and find qualitative agreement. Because our arguments rely on general structures of the phase space, we expect that much of our analysis would be relevant for the study of phenomenologically well-motivated observables, such as N -subjettiness, energy correlation functions, and planar flow.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2015

The nature and orbit of the Ophiuchus stream

The Ophiuchus stream is a recently discovered stellar tidal stream in the Milky Way. We present high-quality spectroscopic data for 14 stream member stars obtained using the Keck and MMT telescopes. We confirm the stream as a fast moving (v{sub los} ∼ 290 km s{sup −1}), kinematically cold group (σ{sub v{sub l{sub o{sub s}}}}≲1 km s{sup −1}) of α-enhanced and metal-poor stars ([α/Fe] ∼ 0.4 dex, [Fe/H] ∼ −2.0 dex). Using a probabilistic technique, we model the stream simultaneously in line-of-sight velocity, color–magnitude, coordinate, and proper motion space, and so determine its distribution in 6D phase-space. We find that the stream extends in distance from 7.5 to 9 kpc from the Sun; it is 50 times longer than wide, merely appearing highly foreshortened in projection. The analysis of the stellar population contained in the stream suggests that it is ∼12 Gyr old, and that its initial stellar mass was ∼2 × 10{sup 4} M{sub ⊙} (or at least ≳7 × 10{sup 3} M{sub ⊙}). Assuming a fiducial Milky Way potential, we fit an orbit to the stream that matches the observed phase-space distribution, except for some tension in the proper motions: the stream has an orbital period of ∼350 Myr, and is on a fairly eccentric orbit (e ∼ 0.66) with a pericenter of ∼3.5 kpc and an apocenter of ∼17 kpc. The phase-space structure and stellar population of the stream show that its progenitor must have been a globular cluster that was disrupted only ∼240 Myr ago. We do not detect any significant overdensity of stars along the stream that would indicate the presence of a progenitor, and conclude that the stream is all that is left of the progenitor.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2016

SU-F-T-145: Modeling of a Proton Pencil Beam Partially Blocked by a Collimator

Purpose: To develop a mathematical model of the fluence distribution of a proton pencil beam that is partially blocked by a collimator. Methods: In one dimension, an exact analytical solution of the proton distribution behind a collimator edge was derived, taking the beam’s phase space, the edge position, and the scattering downstream of the collimator as input. Extension to two dimensions was achieved using a Clarkson-like Method: determine the radial distances from each calculation point to the collimator outline, calculate the beam’s phase space in the direction of each radial line, apply the one-dimensional model along each radial line, and sum the contributions. A Matlab implementation of the full model produced high-resolution distributions with relatively short calculation times (of about 2 ms per calculation point). The model was validated by comparison to Monte Carlo simulations for select pencil beam phase space distributions and collimator shapes. Results: Agreement between model and simulation confirm the exact mathematical solution of the one-dimensional model. In two dimensions the agreement is reasonable, but not exact, due to the approximations introduced in the Clarkson-like method. Conclusion: A model has been developed and validated that gives the lateral fluence distribution of a pencil beam beyond a collimator. Notably, this model does not include free parameters and the final fluence is completely determined by the phase space parameters of the beam, already existing in the treatment planning system, as well as both axial and lateral position of the aperture edge. As such, it can be incorporated in a full pencil-beam dose calculation algorithm allowing for the modeling of field-specific collimators in proton pencil beam scanning and, potentially, the optimization of collimator shape and/or pencil phase space in addition to the customary spot weight optimization.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2018

Quantum Representation of Classical Dynamics

A constructive proof is given to show that there is an exact embedding of a classical system of N ordinary differential equations (ODEs) within a quantum mechanical system corresponding to a classical phase space of dimension 2N. Doubling the classical phase space allows one to express any set of ODEs through a classical Hamiltonian that is a sum over a set of constraints. The Lagrange multipliers that enforce the constraints are canonically conjugate to the dependent variables of the original ODE and satisfy equations that ensure that the evolution of the system is symplectic. For the quantized system on extended phase space, Heisenberg’s uncertainty principle is satisfied by each variable and its canonically conjugate momentum, the Lagrange multiplier. However, there is no uncertainty in a simultaneous measurement for any of the variables of the original ODE. Hence,there is complete fidelity for the quantum representation of the classical system.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2015

Ergodic theory and visualization. II. Fourier mesochronic plots visualize (quasi)periodic sets

We present an application and analysis of a visualization method for measure-preserving dynamical systems introduced by I. Mezić and A. Banaszuk [Physica D 197, 101 (2004)], based on frequency analysis and Koopman operator theory. This extends our earlier work on visualization of ergodic partition [Z. Levnajić and I. Mezić, Chaos 20, 033114 (2010)]. Our method employs the concept of Fourier time average [I. Mezić and A. Banaszuk, Physica D 197, 101 (2004)], and is realized as a computational algorithms for visualization of periodic and quasi-periodic sets in the phase space. The complement of periodic phase space partition contains chaotic zone, and we show how to identify it. The range of method's applicability is illustrated using well-known Chirikov standard map, while its potential in illuminating higher-dimensional dynamics is presented by studying the Froeschlé map and the Extended Standard Map.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2017

Semiclassical approach to dynamics of interacting fermions

Understanding the behaviour of interacting fermions is of fundamental interest in many fields ranging from condensed matter to high energy physics. Developing numerically efficient and accurate simulation methods is an indispensable part of this. Already in equilibrium, fermions are notoriously hard to handle due to the sign problem. Out of equilibrium, an important outstanding problem is the efficient numerical simulation of the dynamics of these systems. In this work we develop a new semiclassical phase-space approach (a.k.a. the truncated Wigner approximation) for simulating the dynamics of interacting fermions in arbitrary dimensions. As fermions are essentially non-classical objects, a phase-space is constructed out of all fermionic bilinears. Classical phase-space is thus comprised of highly non-local (hidden) variables representing these bilinears, and the cost of the method is that it scales quadratic rather than linear with system size. We demonstrate the strength of the method by comparing the results to the exact quantum dynamics of fermion expansion in the Hubbard model and quantum thermalization in the Sachdev–Ye–Kitaev (SYK) model for small systems, where the semiclassics nearly perfectly reproduces correct results. We furthermore analyse fermion expansion in a larger, intractable by exact methods, 2D Hubbard model, which is directly relevant to recent cold atom experiments. - Highlights: • We propose a semiclassical method to simulate dynamics of interacting fermions. • A phase-space is constructed out of all fermionic bilinears. • We study quantum thermalization in the Sachdev–Ye–Kitaev model.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2020

Bounded motion design in the Earth zonal problem using differential algebra based normal form methods

Establishing long-term relative bounded motion between orbits in perturbed dynamics is a key challenge in astrodynamics to enable cluster flight with minimum propellant expenditure. Herein, we present an approach that allows for the design of long-term relative bounded motion considering a zonal gravitational model. Entire sets of orbits are obtained via high-order Taylor expansions of Poincarè return maps about reference fixed points. The high-order normal form algorithm is used to determine a change in expansion variables of the map into normal form space, in which the phase space behavior is circular and can be easily parameterized by action–angle coordinates. The action–angle representation of the normal form coordinates is then used to parameterize the original Poincarè return map and average it over a full phase space revolution by a path integral along the angle parameterization. As a result, the averaged nodal period and drift in the ascending node are obtained, for which the bounded motion conditions are straightforwardly imposed. Sets of highly accurate bounded orbits are obtained, extending over several thousand kilometers, and valid for decades.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2014

Influence of polarization and a source model for dose calculation in MRT

Purpose: Microbeam Radiation Therapy (MRT), an alternative preclinical treatment strategy using spatially modulated synchrotron radiation on a micrometer scale, has the great potential to cure malignant tumors (e.g., brain tumors) while having low side effects on normal tissue. Dose measurement and calculation in MRT is challenging because of the spatial accuracy required and the arising high dose differences. Dose calculation with Monte Carlo simulations is time consuming and their accuracy is still a matter of debate. In particular, the influence of photon polarization has been discussed in the literature. Moreover, it is controversial whether a complete knowledge of phase space trajectories, i.e., the simulation of the machine from the wiggler to the collimator, is necessary in order to accurately calculate the dose. Methods: With Monte Carlo simulations in the Geant4 toolkit, the authors investigate the influence of polarization on the dose distribution and the therapeutically important peak to valley dose ratios (PVDRs). Furthermore, the authors analyze in detail phase space information provided byMartínez-Rovira et al. [“Development and commissioning of a Monte Carlo photon model for the forthcoming clinical trials in microbeam radiation therapy,” Med. Phys. 39(1), 119–131 (2012)] and examine its influence on peak and valley doses. A simple source model is developed using parallel beams and its applicability is shown in a semiadjoint Monte Carlo simulation. Results are compared to measurements and previously published data. Results: Polarization has a significant influence on the scattered dose outside the microbeam field. In the radiation field, however, dose and PVDRs deduced from calculations without polarization and with polarization differ by less than 3%. The authors show that the key consequences from the phase space information for dose calculations are inhomogeneous primary photon flux, partial absorption due to inclined beam incidence outside the field center, increased beam width and center to center distance due to the beam propagation from the collimator to the phantom surface and imperfect absorption in the absorber material of the Multislit Collimator. These corrections have an effect of approximately 10% on the valley dose and suffice to describe doses in MRT within the measurement uncertainties of currently available dosimetry techniques. Conclusions: The source for the first clinical pet trials in MRT is characterized with respect to its phase space and the photon polarization. The results suggest the use of a presented simplified phase space model in dose calculations and hence pave the way for alternative and fast dose calculation algorithms. They also show that the polarization is of minor importance for the clinical important peak and valley doses inside the microbeam field.

61 RADIATION PROTECTION AND DOSIMETRY↗

DOE OSTI2022

Modelling the Alfvén eigenmode induced fast-ion flow measured by an imaging neutral particle analyzer

An imaging neutral particle analyzer (INPA) provides energy and radially resolved measurements of the confined fast-ion population ranging from the high-field side to the edge on the midplane of the DIII-D tokamak. In recent experiments, it was used to diagnose fast-ion flow in the INPA-interrogated phase-space driven by multiple, marginally unstable Alfvén eigenmodes (AEs). Here, the key features of this measured fast-ion flow are: (I) a fast-ion flow from q min and the injection energy (81 keV) towards lower energies and plasma periphery.(II) A flow from the same location towards higher energies and the plasma core, (III) a phase-space 'hole' at the injected energy and plasma core and (IV) a pile-up at the plasma core at lower energies (~60 keV). Ad hoc energetic particle diffusivity modelling of TRANSP significantly deviates from the observation. Comparably, a reduced modelling, i.e. a combination of NOVA-K and ASCOT5 code with the measured mode structure and amplitude, generally reproduce some key features of the observed phase-space flow, but largely failed to interpret fast ion depletion near the plasma axis. At last, self-consistent, first-principle multi-phase hybrid simulations that include realistic neutral beam injection and collisions are able to reproduce most features of the time-resolved phase-space flow. During consecutive hybrid phases, an RSAE consistent with the experiment grows and saturates, redistributing the injected fast ions. The resulting synthetic INPA images are in good agreement with the measurement near the injection energy. The simulations track the fast-ion redistribution within the INPA range, confirming that the measured fast-ion flow follows streamlines defined by the intersection of phase-space surfaces of constant magnetic moment μ and constant E' = nE + ωP φ , where n and ω are the instability toroidal mode number and frequency, and E and P φ the ion energy and toroidal canonical momentum. Nonperturbative effects are required to reproduce the depletion of fast ions near the magnetic axis at the injection energy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

High-throughput ab initio design of atomic interfaces using InterMatch

Forming a hetero-interface is a materials-design strategy that can access an astronomically large phase space. However, the immense phase space necessitates a high-throughput approach for an optimal interface design. Here we introduce a high-throughput computational framework, InterMatch, for efficiently predicting charge transfer, strain, and superlattice structure of an interface by leveraging the databases of individual bulk materials. Specifically, the algorithm reads in the lattice vectors, density of states, and the stiffness tensors for each material in their isolated form from the Materials Project. From these bulk properties, InterMatch estimates the interfacial properties. We benchmark InterMatch predictions for the charge transfer against experimental measurements and supercell density-functional theory calculations. We then use InterMatch to predict promising interface candidates for doping transition metal dichalcogenide MoSe 2 . Finally, we explain experimental observation of factor of 10 variation in the supercell periodicity within a few microns in graphene/α-RuCl 3 by exploring low energy superlattice structures as a function of twist angle using InterMatch. We anticipate our open-source InterMatch algorithm accelerating and guiding ever-growing interfacial design efforts. Moreover, the interface database resulting from the InterMatch searches presented in this paper can be readily accessed online.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2023

Intrinsic thermal conductivity of ZrC from low to ultrahigh temperatures: A critical revisit

Current phonon transport theory based on ground-state calculations has been successful in predicting thermal conductivity at room and medium temperatures but may misrepresent behavior at high temperatures. Here, in this work, we predict the thermal conductivity (κ) of ZrC including electronic and phonon contributions from 300 to 3500 K, by including high-order phonon scattering; lattice expansion; temperature-dependent (TD) second-, third-, and fourth-order force constants (2FC, 3FC, and 4FC); and interband phonon conduction by using first principles. For the phonon transport, we find that four-phonon scattering (4ph) significantly reduces the phonon thermal conductivity (κ ph ), by as much as ~75% at 3500 K. After including 4ph scattering and all other factors, κ ph shows a ~T –1.5 rather than ~T –1 dependence. TD 2FC decreases three-phonon scattering rates but increases 4ph rates by decreasing and increasing the scattering phase spaces, respectively. For 4ph phase space, the TD 2FC flattens phonon bands, and allows more redistribution-4ph processes (1 + 2 → 3 + 4) to happen. The combination effect of TD 2FC and TD 4FC reduces 4ph rates of acoustic modes but increases those of optical modes. The TD 3FC and 4FC decrease the phonon scattering cross section and increase the κ ph significantly (by 52% at 3500 K). The contribution from interband (Wigner) phonon conduction is small, even at ultrahigh temperatures. For electronic thermal transport, we find that it is sensitive to and can be changed by 20% by the TD lattice constants. The Lorenz number varies from 1.6 to 3.3 × 10 –8 WΩK –2 at different temperatures. The theoretical prediction in the literature overpredicts κph (e.g., ~28%) and underpredicts the κ el (e.g., ~38%), resulting in an overall underprediction of κ (~26% at 1500 K). The impacts of grain size and defects are found to be strong, leading to the lower observed thermal conductivity in experiments.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2021

Distant Relatives: The Chemical Homogeneity of Comoving Pairs Identified in Gaia

Comoving pairs, even at the separations of $\mathcal{O}$(10 6 ) are a predicted reservoir of conatal stars. Here we present detailed chemical abundances of 62 stars in 31 comoving pairs with separations of 10 2 –10 7 au and 3D velocity differences <2 km s -1 . This sample includes both bound comoving pairs/wide binaries and unbound comoving pairs. Observations were taken using the Magellan Inamori Kyocera Echelle (MIKE) spectrograph on board the Magellan/Clay Telescope at high resolution (R~ 45,000) with a typical signal-to-noise ratio of 150 pixel -1 . With these spectra, we measure surface abundances for 24 elements, including Li, C, Na, Mg, Al, Si, Ca, Sc, Ti, V, Cr, Mn, Fe, Co, Ni, Cu, Zn, Sr, Y, Zr, Ba, La, Nd, and Eu. Taking iron as the representative element, our sample of wide binaries is chemically homogeneous at the level of 0.05 dex, which agrees with prior studies on wide binaries. Importantly, even systems at separations 2 × 10 5 –10 7 au are homogeneous to 0.09 dex, as opposed to the random pairs, which have a dispersion of 0.23 dex. Assuming a mixture model of the wide binaries and random pairs, we find that 73 ± 22% of the comoving pairs at separations 2 × 10 5 –10 7 au are conatal. Our results imply that a much larger parameter space of phase space may be used to find conatal stars, to study M-dwarfs, star cluster evolution, exoplanets, chemical tagging, and beyond.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2013

LUMINOUS SATELLITES VERSUS DARK SUBHALOS: CLUSTERING IN THE MILKY WAY

The observed population of the Milky Way satellite galaxies offers a unique testing ground for galaxy formation theory on small scales. Our novel approach was to investigate the clustering of the known Milky Way satellite galaxies and to quantify the amount of substructure within their distribution using a two-point correlation function statistic in each of three spaces: configuration space, line-of-sight velocity space, and four-dimensional (4D) phase space. These results were compared to those for three sets of subhalos in the Via Lactea II cold dark matter (CDM) simulation defined to represent the luminous dwarfs. We found no evidence at a significance level above 2{sigma} of substructure within the distribution of the Milky Way satellite galaxies in any of the three spaces. The 'luminous' subhalo sets are more strongly clustered than are the Milky Way satellites in all three spaces and over a broader range of scales in 4D phase space. Each of the 'luminous' subhalo sets are clustered as a result of substructure within their line-of-sight velocity space distributions at greater than 3{sigma} significance, whereas the Milky Way satellite galaxies are randomly distributed in line-of-sight velocity space. While our comparison is with only one CDM simulation, the inconsistencies between the Milky Way satellite galaxies and the Via Lactea II subhalo sets for all clustering methods suggest a potential new 'small-scale' tension between CDM theory and the observed Milky Way satellites. Future work will obtain a more robust comparison between the observed Milky Way satellites and CDM theory by studying additional simulations.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2019

Predicting Reaction Products and Automating Reactive Trajectory Characterization in Molecular Simulations with Support Vector Machines

A machine learning-based methodology for the prediction of chemical reaction products, along with automated elucidation of mechanistic details via phase space analysis of reactive trajectories, is introduced using low-dimensional heuristic models and then applied to ab initio computer simulations of the photodissociation of acetaldehyde, an important chemical system in atmospheric chemistry. Our method is centered around training Support Vector Machines (SVMs) to identify optimal separatrices that delineate the regions of phase space that lead to different chemical reaction products. In contrast to the more common “black box” type machine learning methodologies for analyzing chemical simulation data, this SVM-based methodology allows for mechanistic insight to be gleaned from further analysis of the positioning of the phase space points used to train the SVM with respect to the separatrices. For example, a pair of phase space points that are in close proximity to each other but on opposite sides of a separatrix may be situated on opposite sides of a transition state, while phase space points occurring early in a simulation that are distant from a separatrix are likely to belong to trajectories that are highly biased toward the product state associated with the basin of attraction to which they belong. In addition to inferring mechanistic details about multiple-pathway chemical reactions, our method can also be used to increase reactive trajectory sampling efficiency in molecular simulations via rejection sampling, with trajectories leading to undesired product states being identified and terminated early in the simulation rather than being carried to completion.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗