Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “density correlation”

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 235 records · Page 13

Mitigating imaging systematics for DESI 2024 emission Line Galaxies and beyond

Emission Line Galaxies (ELGs) are one of the main tracers that the Dark Energy Spectroscopic Instrument (DESI) uses to probe the universe. However, they are afflicted by strong spurious correlations between target density and observing conditions known as imaging systematics. In this paper, we present the imaging systematics mitigation applied to the DESI Data Release 1 (DR1) large-scale structure catalogs used in the DESI 2024 cosmological analyses. We also explore extensions of the fiducial treatment. This includes a combined approach, through forward image simulations (Obiwan) in conjunction with neural network-based regression, to obtain an angular selection function that mitigates the imaging systematics observed in the DESI DR1 ELGs target density. We further derive a line of sight selection function from the forward model that removes the strong redshift dependence between imaging systematics and low redshift ELGs. Combining both angular and redshift-dependent systematics, we construct a three-dimensional selection function and assess the impact of all selection functions on clustering statistics. We quantify differences between these extended treatments and the fiducial treatment in terms of the measured 2-point statistics. We find that the results are generally consistent with the fiducial treatment and conclude that the differences are far less than the imaging systematics uncertainty included in DESI 2024 full-shape measurements. We extend our investigation to the ELGs at 0.6 < z < 0.8, i.e., beyond the redshift range (0.8 < z < 1.6) adopted for the DESI clustering catalog, and demonstrate that determining the full three-dimensional selection function is necessary in this redshift range. Our tests showed that all changes are consistent with statistical noise for BAO analyses indicating they are robust to even severe imaging systematics. Specific tests for the full-shape analysis will be presented in a companion paper.

79 ASTRONOMY AND ASTROPHYSICS↗

Structure of Charge Density Waves in La1.875Ba0.125CuO4

Although charge-density wave (CDW) correlations exist in several families of cuprate superconductors, they exhibit substantial variation in CDW wavevector and correlation length, indicating a key role for CDW-lattice interactions. We investigated this interaction in La1.875Ba0.125CuO4 using single crystal x-ray diffraction to collect a large number of CDW peak intensities, and determined the Cu and La/Ba atomic distortions induced by the formation of CDW order. Within the CuO2 planes, the distortions involve a periodic modulation of the Cu-Cu spacing along the direction of the ordering wave vector. The charge ordering within the copper-oxygen layer induces an out-of-plane breathing modulation of the surrounding lanthanum layers, which leads to a related distortion on the adjacent copper-oxygen layer. Our result implies that the CDW-related structural distortions do not remain confined to a single layer but rather propagate an appreciable distance through the crystal. This leads to overlapping structural modulations, in which CuO2 planes exhibit distortions arising from the orthogonal CDWs in adjacent layers as well as distortions from the CDW within the layer itself. We attribute this striking effect to the weak c-axis charge screening in cuprates and suggest this effect could help couple the CDW between adjacent planes in the crystal.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Ground-state phase diagram and superconductivity of the doped Hubbard model on six-leg square cylinders

Here, we have studied the ground state properties of Hubbard model on long six-leg square cylinders with doped hole concentration per site 5.55%≤δ≤12.5% using density-matrix renormalization group. By keeping a large number of states for long system sizes, we find that the nature of the ground state is remarkably sensitive to the presence of next-nearest-neighbor electron hopping t'. In the positive t' side, we find a robust d-wave superconducting (SC) phase characterized by coexisting quasi-long-range SC and charge density wave (CDW) correlations. Without t' the ground state forms an insulating stripe phase with long-range CDW order but short-range spin-spin and SC correlations. In stark contrast to four-leg cylinders, our results show that the lightly doped Hubbard model on six-leg cylinders remains insulating in the negative t' side where the SC correlations decay exponentially with short correlation lengths. In the larger negative t' side, the doped holes form a novel holon Wigner crystal with one doped hole per emergent unit cell and short-range spin-spin correlations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Simultaneous velocity and density measurements of fully developed Rayleigh-Taylor mixing

Here, the dynamics of molecular mixing and the energy transfer process in the Rayleigh-Taylor instability (RTI) are studied through the collection of simultaneous density-velocity measurements. These experiments provide simultaneous density-velocity field measurements, in contrast to previous point measurements. Statistically stationary experiments are performed in a “convective-type” gas tunnel facility, with density contrast achieved through the injection of helium into the bottom stream. Three experiments at Atwood number ≈ 0.1 are captured at three outer scale Reynolds numbers Re = 520, 2260, and 4050. Particle image velocimetry and laser induced fluorescence are employed simultaneously. Statistics of the density and velocity show self-similar collapse of RTI profiles at large Reynolds number Re >2000. Flat velocity profiles indicate homogeneous turbulence characteristics in the core of the mixing region. Significant anisotropy develops in the flow, with horizontal velocity fluctuations being only 60% of the vertical velocity fluctuations. The turbulent mass flux is found to be asymmetric about the centerline, with increased peak towards the spike. Measurements of the molecular mixing show that mixing is maximized at the core of the flow and increases with increased Reynolds number. The transport equation of density-specific-volume correlation b shows that it is mostly produced in the core of the mixing region, and that the spatial evolution of its profile is the result of transport by bulk motion of the bubble and spike. Energy transfer from gravitational potential energy to turbulent kinetic energy and viscous dissipation is observed to occur in the experiment with a ratio of dissipated energy to potential energy released of 38%. The analysis of the turbulent kinetic energy transport equation budget reveals that production is the primary mechanism towards the growth of turbulent kinetic energy in the core of the flow, and is asymmetrically slightly skewed towards the spike. However, it is through the transport that the strong advection at the edges of the mixing region is maintained.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Anomalous Resistivity and Electron Heating by Lower Hybrid Drift Waves during Magnetic Reconnection with a Guide Field

The lower hybrid drift wave (LHDW) has been a candidate for anomalous resistivity and electron heating inside the electron diffusion region of magnetic reconnection. In a laboratory reconnection layer with a finite guide field, quasielectrostatic LHDW (ES-LHDW) propagating along the direction nearly perpendicular to the local magnetic field is excited in the electron diffusion region. ES-LHDW generates large density fluctuations ($δn_{\text e}$, about 25% of the mean density) that are correlated with fluctuations in the out-of-plane electric field ($δE_Y$, about twice larger than the mean reconnection electric field). With a small phase difference (~ 30°) between two fluctuating quantities, the anomalous resistivity associated with the observed ES-LHDW is twice larger than the classical resistivity and accounts for 20% of the mean reconnection electric field. After we verify the linear relationship between $δn_{\text e}$ and $δE_Y$, anomalous electron heating by LHDW is estimated by a quasilinear analysis. The estimated electron heating is about 2.6 ± 0.3 MW / m 3 , which exceeds the classical Ohmic heating of about 2.0 ± 0.2 MW / m 3 . This LHDW-driven heating is consistent with the observed trend of higher electron temperatures when the wave amplitude is larger. Finally, presented results provide the first direct estimate of anomalous resistivity and electron heating power by LHDW, which demonstrates the importance of wave-particle interactions in magnetic reconnection.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Isotropy of Hubble Expansion in the Early and Late Universe

We test the isotropy of Hubble expansion by combining several probes for the first time, constructing full-sky maps of expansion rate variation using Type Ia supernovae, fundamental plane galaxies, and cosmic microwave background (CMB) temperature fluctuations. We find no hint of anisotropy or correlation between early- and late-Universe expansion across all systematic models. The 99% confidence upper limits on expansion rate anisotropy are 0.39% for low-redshift supernovae, 0.95% for high-redshift CMB, and 0.37% when combined at a 60-degree smoothing scale. A significant anomaly in the fundamental plane residual map may reflect systematics in the current dark energy spectroscopic instrument dataset, as evidenced by the absence of cross-correlation with other tracers and its correlation with spatial density variations.

Zhou, Alan Junzhe [Chicago U.; Chicago U., KICP] (↗

Reliability and Lifetime Prediction Model of Sintered Silver Under High-Temperature Cycling

Although excellent reliability has been reported for sintered silver as a die-attach material under both thermal and power cycling loads in power electronics applications, the promise of this material as a large-area attachment at temperatures beyond 200 degrees C needs to be investigated. This paper presents insights into the thermomechanical behavior and reliability of sintered silver under extreme thermal cycling conditions. In this study, we bonded sintered silver samples and subjected it to a thermal cycling profile of -40 °C to 200 °C with high ramp rates. We periodically monitored samples under thermal cycling to detect the presence of any failure mechanisms using a scanning acoustic microscope. We also included 95Pb5Sn solder in the study to obtain reference data. Results show the occurrence of cracks in sintered silver followed by a rapid rate of crack growth that exceeded the failure criterion in just 50 cycles. The predominant failure mechanism we observed was adhesive failure. As a large-area attachment, solder exhibited a higher reliability than sintered silver but failed within 100 cycles. Finally, we performed thermomechanical modeling to compute strain energy density values and correlated these with the experimentally observed crack growth rates to formulate a lifetime prediction model for sintered silver.

30 DIRECT ENERGY CONVERSION↗

The catalog-to-cosmology framework for weak lensing and galaxy clustering for LSST

We present TXPipe, a modular, automated and reproducible pipeline for ingesting catalog data and performing all the calculations required to obtain quality-assured two-point measurements of lensing and clustering, and their covariances, with the metadata necessary for parameter estimation. The pipeline is developed within the Rubin Observatory Legacy Survey of Space and Time (LSST) Dark Energy Science Collaboration (DESC), and designed for cosmology analyses using LSST data. In this paper, we present the pipeline for the so-called 3x2pt analysis -- a combination of three two-point functions that measure the auto- and cross-correlation between galaxy density and shapes. We perform the analysis both in real and harmonic space using TXPipe and other LSST-DESC tools. We validate the pipeline using Gaussian simulations and show that it accurately measures data vectors and recovers the input cosmology to the accuracy level required for the first year of LSST data under this simplified scenario. We also apply the pipeline to a realistic mock galaxy sample extracted from the CosmoDC2 simulation suite (Korytov et al. 2019). TXPipe establishes a baseline framework that can be built upon as the LSST survey proceeds. Furthermore, the pipeline is designed to be easily extended to science probes beyond the 3x2pt analysis.

79 ASTRONOMY AND ASTROPHYSICS↗

UT-GOM2-2 Preliminary Report: Terrebonne Basin Northern Gulf of Mexico, 30 July-28 September 2023

In the summer and fall of 2023, the Gulf of Mexico Deepwater Hydrate Coring Expedition (UT-GOM2-2) drilled, cored, made downhole measurements, and analyzed samples from the seafloor to the base of the gas hydrate stability zone in one location (Site H, WR313) in the Terrebonne basin, deepwater Gulf of Mexico. Analyses of data and samples from the expedition will inform biological, geochemical, and geomechanical models to constrain the role of gas hydrates in the carbon cycle and the potential for gas hydrates as an energy resource. Pressure and conventional cores were collected continuously to a depth of 155.1 meters below the seafloor (mbsf). At deeper depths, cores were taken periodically from hydrate-bearing sands and their bounding muds to a total depth of 861.3 mbsf. 162.6 m of conventional core and 54.8 m of pressure core were obtained. Twelve temperature measurements were made between 27.1 and 144.5 mbsf to determine the geothermal gradient. At the seafloor, more than 4 m of sandy silt of unknown origin was encountered. Beneath this sand, to a depth of ~200 mbsf, the section was composed of interbedded mud and biogenic carbonate ooze. The biogenic ooze correlated to low density and high porosity intervals observed in the previously acquired logging while drilling (LWD) data and as measured. Calcareous nannofossil biostratigraphy constrains the entire record to the Pleistocene (< 0.91 million years) with a pronounced increase in sedimentation rate with depth. Beneath 200 mbsf, the section was predominantly composed of mud with two thicker, hydrate-bearing coarse-grained intervals, which are commonly known as the Blue and Orange sands. The dissolved gas concentration was quantified from pressure cores. In the shallow section, dissolved methane concentration increased below the sulfate-methane transition zone (SMTZ) and reaches saturation (the limit of solubility for methane) at 147 mbsf. Gas expansion was very common in conventional and depressurized pressure (conventionalized) cores below the SMTZ. At deeper depths, the methane concentration within muds bounding the Blue and Orange reservoirs was generally found to be less than saturation. The dissolved and hydrate gas composition is consistent with a microbial source, containing greater than 99.99% methane and only trace concentrations of ethane, propane, and butane. The methane to ethane ratio (C 1 /C 2 ) and the methane to ethane plus propane (C 1 /(C 2 +C 3 )) decrease with depth down to at least 678 mbsf, mainly driven by the increase in ethane with depth. It is unclear if this trend continues through the Orange sand interval. The δ 13 C isotopic signature of methane ranges between -69.9 and -78.5 ‰ Vienna Pee Dee Belemnite (VPDB). Pressure core recovery of all sandy intervals was poor. However, pressure core logs of the Orange sand show intervals of low density and high velocity, which are indicative of high hydrate saturation. One core from within the Orange sand was composed of interbedded graded sandy silt and mud. The sandy silts from this core are composed of mainly quartz and feldspar with some lithics. Most of the recovered pressure core samples are maintained at near in-situ pressure and temperature (within the hydrate stability field) at the University of Texas Pressure Core Center awaiting analysis. In the shallow section, samples will be used to determine the flux of organic carbon through the basin system, find the rate at which that carbon was consumed, and understand the microbial population responsible for these processes. In the deeper section, samples from in and around the hydrate reservoirs will be used to determine the petrophysical properties of the reservoir and bounding seals in these systems.

03 NATURAL GAS↗

Noble Gas Transport in the MSRE

This study explores the relationship between the physicochemical parameter known as Henry’s gas constant and gas transport across a two-layer film interface. The investigation utilized the Gibbs free energy, incorporating surface and volume terms to elucidate trends in enthalpy and entropy. Notably, our findings align with experimental data and offer predictive insights into the Henry’s gas constants for helium and krypton, which hold significance for future experiments and theoretical developments. Furthermore, this study enhances the Gibbs free energy theory pertaining to the liquid–gas interface. It underscores the substantial contribution of noble gases in this region to volumetric energy as temperature increases. Additionally, we employed Monte Carlo simulations to analyze the effective thermal neutron multiplication factor, denoted as k eff . Our analysis reveals a linear correlation between graphite density and uniform density as a function of temperature. For the 1D Molten Salt Reactor Experiment (MSRE) system, we employed the Mole code to conduct heat and mass transfer calculations. These computations enable us to ascertain the distribution of fuel temperature based on coefficients and thermal properties. We also studied delayed neutron precursors during fuel cycling, taking into account the drift of cycling fuel through Mole–Griffin coupling. Our model represents k eff and β eff across various volume flow rates and salinity compositions. Finally, this study leveraged xenon-135 for continuous on-line monitoring of fuel salts to investigate the impact of steady-state xenon-135 on the MSRE and to better understand its distribution. These efforts build upon previous research related to removal processes.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Expedition UT-GOM2-2 Site H

Pressure and conventional cores were collected at Site H of the Walker Ridge Protracted Area Block 313 in the Terrebonne Basin, deepwater Gulf of America (Gulf of Mexico) during the University of Texas (UT) Deepwater Hydrate Coring Expedition (UT-GOM2-2). Pressure and conventional cores were collected continuously to a depth of 155.1 meters below the seafloor (mbsf). At deeper depths, cores were taken periodically from hydrate-bearing sands and their bounding muds to a total depth of 861.3 mbsf. 162.6 m of conventional core and 54.8 m of pressure core were recovered. Twelve temperature measurements were made between 27.1 and 144.5 mbsf to determine the geothermal gradient. At the seafloor, more than 4 m of sandy silt of unknown origin was encountered. Beneath this sand, to a depth of about 200 mbsf, the section was composed of interbedded mud and biogenic carbonate ooze. The ooze correlated to low density and high porosity intervals observed in the previously acquired logging while drilling (LWD) data and as measured. These ooze intervals also correspond to lighter sediment color, increased Ca content based on X-ray florescence (XRF) core scanning, and increased calcareous nannofossil abundance. Calcareous nannofossil biostratigraphy constrains the entire record to the Pleistocene (< 0.91 million years), with a pronounced increase in sedimentation rate with depth. Below 200 mbsf, the section was predominantly composed of mud with two thicker, hydrate-bearing coarse-grained intervals, which are commonly known as the Blue and Orange sands. The dissolved gas concentration was quantified from pressure cores. In the shallow section, dissolved methane concentration increased below the sulfate-methane transition zone (SMTZ) and reaches saturation (the limit of solubility for methane) at 147 mbsf. Gas expansion was very common in conventional and depressurized pressure (conventionalized) cores below the SMTZ. At deeper depths, the methane concentration within muds bounding the Blue and Orange reservoirs was generally found to be less than saturation. The dissolved and hydrate gas composition is consistent with a microbial source, containing greater than 99.99% methane and only trace concentrations of ethane, propane, and butane. The methane to ethane ratio (C1 /C2 ) and the methane to ethane plus propane (C1 /(C2 +C3 )) decrease with depth down to at least 678 mbsf, mainly driven by the increase in ethane with depth. It is unclear if this trend continues through the Orange sand interval. The δ13C isotopic signature of methane ranges between -69.9 and -78.5 ‰ relative to the Vienna Pee Dee Belemnite (VPDB) standard. Pressure core recovery in sandy intervals was poor. However, pressure core logs of the Orange sand show intervals of low density and high velocity, which are indicative of high hydrate saturation. One pressure core was degassed and the average hydrate saturation in the core was determined to be 24%. One core from within the Orange sand was composed of interbedded graded sandy silt and mud. The sandy silts from this core are composed of mainly quartz and feldspar with some lithics. Most of the recovered pressure core samples are maintained at near in-situ pressure and temperature (within the hydrate stability field) at the University of Texas Pressure Core Center awaiting analysis. In the shallow section, samples will be used to determine the flux of organic carbon through the basin system, find the rate at which that carbon was consumed, and understand the microbial population responsible for these processes. In the deeper section, samples from in and around the hydrate reservoirs will be used to determine the petrophysical properties of the reservoir and bounding seals in these systems.

03 NATURAL GAS↗

Kaon Femtoscopy with Lévy-Stable Sources from $\sqrt{s_{NN}}$ = 200 GeV Au+Au Collisions at RHIC

Femtoscopy has the capacity to probe the space-time geometry of the particle-emitting source in heavy-ion collisions. In particular, femtoscopy of like-sign kaon pairs may shed light on the origin of non-Gaussianity of the spatial emission probability density. The momentum correlations between like-sign kaon pairs are measured in data recorded by the STAR experiment, from $\sqrt{s_{NN}}$ = 200 GeV Au+Au collisions at RHIC, BNL. Preliminary results hint at the possible existence of non-Gaussian, Lévy-stable sources, with the likely presence of an anomalous diffusion process in the signal for the identically charged kaon pairs so produced. More statistically significant studies at lower centre-of-mass energies may contribute to the search for the critical end point of QCD.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Standardized Luminosity of the Tip of the Red Giant Branch Utilizing Multiple Fields in NGC 4258 and the CATs Algorithm

Abstract The tip of the red giant branch provides a luminous standard candle for calibrating distance ladders that reach Type Ia supernova (SN Ia) hosts. However, recent work reveals that tip measurements vary at the ∼0.1 mag level for different stellar populations and locations within a host, which may lead to inconsistencies along the distance ladder. We pursue a calibration of the tip using 11 Hubble Space Telescope fields around the maser host, NGC 4258, that is consistent with SN Ia hosts by standardizing tip measurements via their contrast ratios. We find F814W-band tips that exhibit a full 0.3 mag range and 0.1 mag dispersion. We do not find any correlation between Hicolumn density and the apparent tip to 0.04 ± 0.03 mag/cm −2 . We search for a tip–contrast relation (TCR) and measure the TCR within the fields of NGC 4258 of −0.015 ± 0.008 mag/R, whereRis the contrast ratio. This value is consistent with the TCR originally discovered in the GHOSTS sample of −0.023 ± 0.005 mag/R. Combining these measurements, we find a global TCR of −0.021 ± 0.004 mag/Rand a calibration of M I TRGB = − 4.025 ± 0.035 − ( R − 4 ) × 0.021 mag. We also use stellar models to simulate single age and metallicity stellar populations with [Fe/H] from −2.0 to −0.7 and ages from 3 to 12 Gyr and reconstruct the global TCR found here to a factor of ∼2. This work is combined in a companion analysis with tip measurements of nearby SN Ia hosts to measureH 0 .

Astronomy & Astrophysics↗

Modeling the Far-infrared Polarization Spectrum of a High-mass Star-forming Cloud

The polarization spectrum, or wavelength dependence of the polarization fraction, of interstellar dust emission provides important insights into the grain alignment mechanism of interstellar dust grains. We investigate the far-infrared polarization spectrum of a realistic simulated high-mass star-forming cloud under various models of grain alignment and emission. We find that neither a homogeneous grain alignment model nor a grain alignment model that includes collisional dealignment is able to produce the falling spectrum seen in observations. On the other hand, we find that a grain alignment model with grain alignment efficiency dependent on local temperature is capable of producing a falling spectrum that is in qualitative agreement with observations of OMC-1. For the model most in agreement with OMC-1, we find no correlation between the temperature and the slope of the polarization spectrum. However, we do find a positive correlation between the column density and the slope of the polarization spectrum. We suggest this latter correlation to be the result of wavelength-dependent polarization by absorption.

79 ASTRONOMY AND ASTROPHYSICS↗

Subsurface Characterization and Machine Learning Predictions at Brady Hot Springs: Preprint

Subsurface data analysis, reservoir modeling, and machine learning (ML) techniques have been applied to the Brady Hot Springs (BHS) geothermal field in Nevada, USA to further characterize the subsurface and assist with optimizing reservoir management. Hundreds of reservoir simulations have been conducted in TETRAD-G and CMG STARS to explore different injection and production fluid flow rates and allocations and to develop a training data set for ML. This process included simulating the historical injection and production since 1979 and prediction of future performance through 2040. ML networks were created and trained using TensorFlow based on multilayer perceptron (MLP), long short-term memory (LSTM), and convolutional neural network (CNN) architectures. These networks took as input selected flow rates, injection temperatures, and historical field operation data and produced estimates of future production temperatures. This approach was first successfully tested on a simplified single fracture doublet system, followed by the application to the BHS reservoir. Using an initial BHS dataset with 37 simulated scenarios, the trained and validated network predicted the production temperature for 6 production wells with the mean absolute percentage error of less than 8%. In a complementary analysis effort, the principal component analysis applied to 13 BHS geological parameters revealed that vertical fracture permeability shows the strongest correlation with fault density and fault intersection density. A new BHS reservoir model was developed considering the fault intersection density as proxy for permeability. This new reservoir model helps to explore under-exploited zones in the reservoir. A data gathering plan to obtain additional subsurface data was developed; it includes temperature surveying for three idle injection wells, at which the reservoir simulations indicate high bottom-hole temperatures. The collected data assist with calibrating the reservoir model and may lead to converting these wells to producers to access under-exploited zones in the reservoir. Data gathering activities are planned for the first quarter of 2021.

40 EE - Geothermal Technologies Office (EE-4G)↗

Identification of Suitable Vacuum Gas Oils as Plasticizers Using HT-GC × GC-HRMS

Vacuum gas oils (VGOs) have long been heralded as effective plasticizers due to their high boiling-points and lubricating properties. One of these VGOs, HyVac Oil 93050, is a paraffinic plasticizer that is currently utilized in polymer–plasticizer explosive formulations, while this oil effectively increases the elasticity and decreases the sensitivity for current formulations. We seek to identify alternative VGOs with similar density, viscosity, molecular composition, and impurities for future formulations. In this study, 18 VGOs, including HyVac Oil 93050, were initially evaluated for their density and viscosity. The oils were then ranked based on physical characteristics and analyzed for molecular composition using high-temperature comprehensive two-dimensional gas chromatography with high-resolution time-of-flight mass spectrometry (HT-GC × GC-HRMS). Further, the HT-GC × GC-HRMS chromatograms of the top 10 most similar VGOs to HyVac Oil 93050, according to density, were compared utilizing the Pearson correlation coefficient. Three of the oils with densities similar to those of HyVac Oil 93050 had a Pearson correlation within the lot-to-lot variation of HyVac Oil 93050. For the top 3 candidates, Pearson correlation was then utilized as a feature selection technique to discover significant chemical differences. Positive chemical ionization (PCI) and negative chemical ionization (NCI) HT-GC × GC-HRMS chromatograms were also evaluated and aided in the discovery and identification of several key compounds including additives that acted as stabilizers for the VGOs. After this further chemical analysis, two of the original 17 potential VGOs were determined to be physically and chemically similar to HyVac Oil 93050.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

ConKer: An algorithm for evaluating correlations of arbitrary order

Context. High order correlations in the cosmic matter density have become increasingly valuable in cosmological analyses. However, computing these correlation functions is computationally expensive. Aims. We aim to circumvent these challenges by developing a new algorithm called ConKer for estimating correlation functions. Methods. This algorithm performs convolutions of matter distributions with spherical kernels using FFT. Since matter distributions and kernels are defined on a grid, it results in some loss of accuracy in the distance and angle definitions. We study the algorithm setting at which these limitations become critical and suggest ways to minimize them. Results. ConKer is applied to the CMASS sample of the SDSS DR12 galaxy survey and corresponding mock catalogs, and is used to compute the correlation functions up to correlation order n = 5. We compare the n = 2 and n = 3 cases to traditional algorithms to verify the accuracy of the new algorithm. We perform a timing study of the algorithm and find that three of the four distinct processes within the algorithm are nearly independent of the catalog size N , while one subdominant component scales as O ( N ). The dominant portion of the calculation has complexity of O ( N c 4/3 log N c ), where N c is the of cells in a three-dimensional grid corresponding to the matter density. Conclusions. We find ConKer to be a fast and accurate method of probing high order correlations in the cosmic matter density, then discuss its application to upcoming surveys of large-scale structure.

79 ASTRONOMY AND ASTROPHYSICS↗

Symmetry breaking and self-interaction correction in the chromium atom and dimer

Density functional approximations to the exchange–correlation energy can often identify strongly correlated systems and estimate their energetics through energy-minimizing symmetry-breaking. In particular, the binding energy curve of the strongly correlated chromium dimer is described qualitatively by the local spin density approximation (LSDA) and almost quantitatively by the Perdew–Burke–Ernzerhof generalized gradient approximation (PBE-GGA), where the symmetry breaking is antiferromagnetic for both. Here, we show that a full Perdew–Zunger self-interaction-correction (SIC) to LSDA seems to go too far by creating an unphysical symmetry-broken state, with effectively zero magnetic moment but non-zero spin density on each atom, which lies ~4 eV below the antiferromagnetic solution. A similar symmetry-breaking, observed in the atom, better corresponds to the 3d↑↑4s↑3d↓↓4s↓ configuration than to the standard 3d↑↑↑↑↑4s↑. For this new solution, the total energy of the dimer at its observed bond length is higher than that of the separated atoms. In conclusion, these results can be regarded as qualitative evidence that the SIC needs to be scaled down in many-electron regions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗