Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “LIMITING VALUES”

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 613 records · Page 34
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

The rate of binary black hole mergers inferred from advanced LIGO observations surrounding GW150914

A transient gravitational-wave signal, GW150914, was identified in the twin Advanced LIGO detectors on 2015 September 2015 at 09:50:45 UTC. To assess the implications of this discovery, the detectors remained in operation with unchanged configurations over a period of 39 days around the time of the signal. At the detection statistic threshold corresponding to that observed for GW150914, our search of the 16 days of simultaneous two-detector observational data is estimated to have a false-alarm rate (FAR) of <4.9×10{sup −6} yr{sup −1}, yielding a p-value for GW150914 of <2×10{sup −7}. Parameter estimation follow-up on this trigger identifies its source as a binary black hole (BBH) merger with component masses (m{sub 1},m{sub 2})=(36{sub −4}{sup +5},29{sub −4}{sup +4}) M{sub ⊙} at redshift z=0.09{sub −0.04}{sup +0.03} (median and 90% credible range). Here, we report on the constraints these observations place on the rate of BBH coalescences. Considering only GW150914, assuming that all BBHs in the universe have the same masses and spins as this event, imposing a search FAR threshold of 1 per 100 years, and assuming that the BBH merger rate is constant in the comoving frame, we infer a 90% credible range of merger rates between 2--53 Gpc{sup −3} yr{sup −1} (comoving frame). Incorporating all search triggers that pass a much lower threshold while accounting for the uncertainty in the astrophysical origin of each trigger, we estimate a higher rate, ranging from 13--600 Gpc{sup −3} yr{sup −1} depending on assumptions about the BBH mass distribution. All together, our various rate estimates fall in the conservative range 2--600 Gpc{sup −3} yr{sup −1}.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

REVISED STELLAR PROPERTIES OF KEPLER TARGETS FOR THE QUARTER 1-16 TRANSIT DETECTION RUN

We present revised properties for 196,468 stars observed by the NASA Kepler mission and used in the analysis of Quarter 1-16 (Q1-Q16) data to detect and characterize transiting planets. The catalog is based on a compilation of literature values for atmospheric properties (temperature, surface gravity, and metallicity) derived from different observational techniques (photometry, spectroscopy, asteroseismology, and exoplanet transits), which were then homogeneously fitted to a grid of Dartmouth stellar isochrones. We use broadband photometry and asteroseismology to characterize 11,532 Kepler targets which were previously unclassified in the Kepler Input Catalog (KIC). We report the detection of oscillations in 2762 of these targets, classifying them as giant stars and increasing the number of known oscillating giant stars observed by Kepler by ∼20% to a total of ∼15,500 stars. Typical uncertainties in derived radii and masses are ∼40% and ∼20%, respectively, for stars with photometric constraints only, and 5%-15% and ∼10% for stars based on spectroscopy and/or asteroseismology, although these uncertainties vary strongly with spectral type and luminosity class. A comparison with the Q1-Q12 catalog shows a systematic decrease in radii of M dwarfs, while radii for K dwarfs decrease or increase depending on the Q1-Q12 provenance (KIC or Yonsei-Yale isochrones). Radii of F-G dwarfs are on average unchanged, with the exception of newly identified giants. The Q1-Q16 star properties catalog is a first step toward an improved characterization of all Kepler targets to support planet-occurrence studies.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

INITIAL PLANETESIMAL SIZES AND THE SIZE DISTRIBUTION OF SMALL KUIPER BELT OBJECTS

The Kuiper Belt is a remnant from the early solar system and its size distribution contains many important constraints that can be used to test models of planet formation and collisional evolution. We show, by comparing observations with theoretical models, that the observed Kuiper Belt size distribution is well matched by coagulation models, which start with an initial planetesimal population with radii of about 1 km, and subsequent collisional evolution. We find that the observed size distribution above R {approx} 30 km is primordial, i.e., it has not been modified by collisional evolution over the age of the solar system, and that the size distribution below R {approx} 30 km has been modified by collisions and that its slope is well matched by collisional evolution models that use published strength laws. We investigate in detail the resulting size distribution of bodies ranging from 0.01 km to 30 km and find that its slope changes several times as a function of radius before approaching the expected value for an equilibrium collisional cascade of material strength dominated bodies for R {approx}< 0.1 km. Compared to a single power-law size distribution that would span the whole range from 0.01 km to 30 km, we find in general a strong deficit of bodies around R {approx} 10 km and a strong excess of bodies around 2 km in radius. This deficit and excess of bodies are caused by the planetesimal size distribution left over from the runaway growth phase, which left most of the initial mass in small planetesimals while only a small fraction of the total mass is converted into large protoplanets. This excess mass in small planetesimals leaves a permanent signature in the size distribution of small bodies that is not erased after 4.5 Gyr of collisional evolution. Observations of the small Kuiper Belt Object (KBO) size distribution can therefore test if large KBOs grew as a result of runaway growth and constrained the initial planetesimal sizes. We find that results from recent KBO occultation surveys and the observed KBO size distribution can be best matched by an initial planetesimal population that contained about equal mass per logarithmic mass bin in bodies ranging from 0.4 km to 4 km in radius. We further find that we cannot match the observed KBO size distribution if most of the planetesimal mass was contained in bodies that were 10 km in radius or larger simply because their resulting size distribution cannot be sufficiently depleted over 4.5 Gyr to match observations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Surface magnetic field strengths: New tests of magnetoconvective models of M dwarfs

Precision modeling of M dwarfs has become worthwhile in recent years due to the increasingly precise values of masses and radii which can be obtained from eclipsing binary studies. In a recent paper, Torres has identified four prime M dwarf pairs with the most precise empirical determinations of masses and radii. The measured radii are consistently larger than standard stellar models predict by several percent. These four systems potentially provide the most challenging tests of precision evolutionary models of cool dwarfs at the present time. We have previously modeled M dwarfs in the context of a criterion due to Gough and Tayler in which magnetic fields inhibit the onset of convection according to a physics-based prescription. In the present paper, we apply our magnetoconvective approach to the four prime systems in the Torres list. Going a step beyond what we have already modeled in CM Dra (one of the four Torres systems), we note that new constraints on magnetoconvective models of M dwarfs are now available from empirical estimates of magnetic field strengths on the surfaces of these stars. In the present paper, we consider how well our magnetoconvective models succeed when confronted with this new test of surface magnetic field strengths. Among the systems listed by Torres, we find that plausible magnetic models work well for CM Dra, YY Gem, and CU Cnc. (The fourth system in Torres's list does not yet have enough information to warrant magnetic modeling.) Our magnetoconvection models of CM Dra, YY Gem, and CU Cnc yield predictions of the magnetic fluxes on the stellar surface which are consistent with the observed correlation between magnetic flux and X-ray luminosity.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Constraints on cosmology from the cosmic microwave background power spectrum of the 2500 deg{sup 2} SPT-SZ survey

We explore extensions to the ΛCDM cosmology using measurements of the cosmic microwave background (CMB) from the recent SPT-SZ survey, along with data from WMAP7 and measurements of H {sub 0} and baryon acoustic oscillation (BAO). We check for consistency within ΛCDM between these data sets, and find some tension. The CMB alone gives weak support to physics beyond ΛCDM, due to a slight trend relative to ΛCDM of decreasing power toward smaller angular scales. While it may be due to statistical fluctuation, this trend could also be explained by several extensions. We consider running of the primordial spectral index (dn{sub s} /dln k), as well as two extensions that modify the damping tail power (the primordial helium abundance Y{sub p} and the effective number of neutrino species N {sub eff}) and one that modifies the large-scale power due to the integrated Sachs-Wolfe effect (the sum of neutrino masses ∑m {sub ν}). These extensions have similar observational consequences and are partially degenerate when considered simultaneously. Of the six one-parameter extensions considered, we find CMB to have the largest preference for dn{sub s} /dln k with –0.046 < dn{sub s} /dln k < –0.003 at 95% confidence, which strengthens to a 2.7σ indication of dn{sub s} /dln k < 0 from CMB+BAO+H {sub 0}. Detectable dn{sub s} /dln k ≠ 0 is difficult to explain in the context of single-field, slow-roll inflation models. We find N {sub eff} = 3.62 ± 0.48 for the CMB, which tightens to N {sub eff} = 3.71 ± 0.35 from CMB+BAO+H {sub 0}. Larger values of N {sub eff} relieve the mild tension between CMB, BAO, and H {sub 0}. When the Sunyaev-Zel'dovich selected galaxy cluster abundances (SPT{sub CL}) data are also included, we obtain N {sub eff} = 3.29 ± 0.31. Allowing for ∑m {sub ν} gives a 3.0σ detection of ∑m {sub ν} > 0 from CMB+BAO+H {sub 0} +SPT{sub CL}. The median value is (0.32 ± 0.11) eV, a factor of six above the lower bound set by neutrino oscillation observations. All data sets except H {sub 0} show some preference for massive neutrinos; data combinations including H {sub 0} favor nonzero masses only if BAO data are also included. We also constrain the two-parameter extensions N {sub eff} + ∑m {sub ν} and N {sub eff} + Y{sub p} to explore constraints on additional light species and big bang nucleosynthesis, respectively.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

THE AGES OF 55 GLOBULAR CLUSTERS AS DETERMINED USING AN IMPROVED ΔV{sup HB}{sub TO} METHOD ALONG WITH COLOR-MAGNITUDE DIAGRAM CONSTRAINTS, AND THEIR IMPLICATIONS FOR BROADER ISSUES

Ages have been derived for 55 globular clusters (GCs) for which Hubble Space Telescope Advanced Camera for Surveys photometry is publicly available. For most of them, the assumed distances are based on fits of theoretical zero-age horizontal-branch (ZAHB) loci to the lower bound of the observed distributions of HB stars, assuming reddenings from empirical dust maps and metallicities from the latest spectroscopic analyses. The age of the isochrone that provides the best fit to the stars in the vicinity of the turnoff (TO) is taken to be the best estimate of the cluster age. The morphology of isochrones between the TO and the beginning part of the subgiant branch (SGB) is shown to be nearly independent of age and chemical abundances. For well-defined color-magnitude diagrams (CMDs), the error bar arising just from the 'fitting' of ZAHBs and isochrones is ≈ ± 0.25 Gyr, while that associated with distance and chemical abundance uncertainties is ∼ ± 1.5-2 Gyr. The oldest GCs in our sample are predicted to have ages of ≈13.0 Gyr (subject to the aforementioned uncertainties). However, the main focus of this investigation is on relative GC ages. In conflict with recent findings based on the relative main-sequence fitting method, which have been studied in some detail and reconciled with our results, ages are found to vary from mean values of ≈12.5 Gyr at [Fe/H] ∼< – 1.7 to ≈11 Gyr at [Fe/H] ∼> –1. At intermediate metallicities, the age-metallicity relation (AMR) appears to be bifurcated: one branch apparently contains clusters with disk-like kinematics, whereas the other branch, which is displaced to lower [Fe/H] values by ≈0.6 dex at a fixed age, is populated by clusters with halo-type orbits. The dispersion in age about each component of the AMR is ∼ ± 0.5 Gyr. There is no apparent dependence of age on Galactocentric distance (R{sub G}) nor is there a clear correlation of HB type with age. As previously discovered in the case of M3 and M13, subtle variations have been found in the slope of the SGB in the CMDs of other metal-poor ([Fe/H] ∼< – 1.5) GCs. They have been tentatively attributed to cluster-to-cluster differences in the abundance of helium. Curiously, GCs that have relatively steep 'M13-like' SGBs tend to be massive systems, located at small R{sub G}, that show the strongest evidence of in situ formation of multiple stellar populations. The clusters in the other group are typically low-mass systems (with 2-3 exceptions, including M3) that, at the present time, should not be able to retain the matter lost by mass-losing stars due either to the development of GC winds or to ram-pressure stripping by the halo interstellar medium. The apparent separation of the two groups in terms of their present-day gas retention properties is difficult to understand if all GCs were initially ∼20 times their current masses. The lowest-mass systems, in particular, may have never been massive enough to retain enough gas to produce a significant population of second-generation stars. In this case, the observed light element abundance variations, which are characteristic of all GCs, were presumably present in the gas out of which the observed cluster stars formed.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

The radio relics and halo of El Gordo, a massive z = 0.870 cluster merger

We present 610 MHz and 2.1 GHz imaging of the massive Sunyaev-Zel'dovich Effect selected z = 0.870 cluster merger ACT-CL J0102–4915 ({sup E}l Gordo{sup )}, obtained with the Giant Metrewave Radio Telescope and the Australia Telescope Compact Array (ATCA), respectively. We detect two complexes of radio relics separated by 3.'4 (1.6 Mpc) along the system's northwest-to-southeast collision axis that have high integrated polarization fractions (33%) and steep spectral indices (α between 1 and 2; S {sub ν}∝ν{sup –α}), consistent with creation via Fermi acceleration by shocks in the intracluster medium triggered by the cluster collision. From the spectral index of the relics, we compute a Mach number M=2.5{sub −0.3}{sup +0.7} and shock speed of 2500{sub −300}{sup +400} km s{sup −1}. With our wide-bandwidth, full-polarization ATCA data, we compute the Faraday depth φ across the northwest relic and find a range of values spanning Δφ = 30 rad m{sup –2}, with a mean value of (φ) = 11 rad m{sup –2} and standard deviation σ{sub φ} = 6 rad m{sup –2}. With the integrated line-of-sight gas density derived from new Chandra X-ray observations, our Faraday depth measurement implies B {sub ∥} ∼ 0.01 μG in the cluster outskirts. The extremely narrow shock widths in the relics (d {sub shock} ≤ 23 kpc), caused by the short synchrotron cooling timescale of relativistic electrons at z = 0.870, prevent us from placing a meaningful constraint on the magnetic field strength B using cooling time arguments. In addition to the relics, we detect a large (r {sub H} ≅ 1.1 Mpc radius), powerful (log (L {sub 1.4}/W Hz{sup –1}) = 25.66 ± 0.12) radio halo with a shape similar to El Gordo's 'bullet'-like X-ray morphology. The spatially resolved spectral-index map of the halo shows the synchrotron spectrum is flattest near the relics, along the system's collision axis, and in regions of high T {sub gas}, all locations associated with recent energy injection. The spatial and spectral correlation between the halo emission and cluster X-ray properties supports primary-electron processes like turbulent reacceleration as the halo production mechanism. The halo's integrated 610 MHz to 2.1 GHz spectral index is a relatively flat α = 1.2 ± 0.1, consistent with the cluster's high T {sub gas} in view of previously established global scaling relations. El Gordo is the highest-redshift cluster known to host a radio halo and/or radio relics, and provides new constraints on the non-thermal physics in clusters at z > 0.6.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

SU-F-T-128: Dose-Volume Constraints for Particle Therapy Treatment Planning

Purpose: Determine equivalent Organ at Risk (OAR) tolerance dose (TD) constraints for MV x-rays and particle therapy. Methods: Equivalent TD estimates for MV x-rays are determined from an isoeffect, regression-analysis of published and in-house constraints for various fractionation schedules (n fractions). The analysis yields an estimate of (α/β) for an OAR. To determine equivalent particle therapy constraints, the MV x-ray TD(n) values are divided by the RBE for DSB induction (RBE{sub DSB}) or cell survival (RBE{sub S}). Estimates of (RBE{sub DSB}) are computed using the Monte Carlo Damage Simulation, and estimates of RBES are computed using the Repair-Misrepair-Fixation (RMF) model. A research build of the RayStation™ treatment planning system implementing the above model is used to estimate (RBE{sub DSB}) for OARs of interest in 16 proton therapy patient plans (head and neck, thorax, prostate and brain). Results: The analysis gives an (α/β) estimate of about 20 Gy for the trachea and heart and 2–4 Gy for the esophagus, spine, and brachial plexus. Extrapolation of MV x-ray constraints (n = 1) to fast neutrons using RBE{sub DSB} = 2.7 are in excellent agreement with clinical experience (n = 10 to 20). When conventional (n > 30) x-ray treatments are used as the reference radiation, fast neutron RBE increased to a maximum of 6. For comparison to a constant RBE of 1.1, the RayStation™ analysis gave estimates of proton RBE{sub DSB} from 1.03 to 1.33 for OARs of interest. Conclusion: The presented system of models is a convenient formalism to synthesize from multiple sources of information a set of self-consistent plan constraints for MV x-ray and hadron therapy treatments. Estimates of RBE{sub DSB} from the RayStation™ analysis differ substantially from 1.1 and vary among patients and treatment sites. A treatment planning system that incorporates patient and anatomy-specific corrections in proton RBE would create opportunities to increase the therapeutic ratio. The research build of the RayStation used in the study was made available to the University of Washington free of charge. RaySearch Laboratories did not provide any monetary support for the reported studies.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

CLASH-VLT: testing the nature of gravity with galaxy cluster mass profiles

We use high-precision kinematic and lensing measurements of the total mass profile of the dynamically relaxed galaxy cluster MACS J1206.2-0847 at z=0.44 to estimate the value of the ratio η=Ψ/Φ between the two scalar potentials in the linear perturbed Friedmann-Lemaitre-Robertson-Walker metric. An accurate measurement of this ratio, called anisotropic stress, could show possible, interesting deviations from the predictions of the theory of General Relativity, according to which Ψ should be equal to Φ. Complementary kinematic and lensing mass profiles were derived from exhaustive analyses using the data from the Cluster Lensing And Supernova survey with Hubble (CLASH) and the spectroscopic follow-up with the Very Large Telescope (CLASH-VLT). Whereas the kinematic mass profile tracks only the time-time part of the perturbed metric (i.e. only Φ), the lensing mass profile reflects the contribution of both time-time and space-space components (i.e. the sum Φ+Ψ). We thus express η as a function of the mass profiles and perform our analysis over the radial range 0.5 Mpc≤ r≤ r{sub 200}=1.96 Mpc. Using a spherical Navarro-Frenk-White mass profile, which well fits the data, we obtain η(r{sub 200})=1.01 {sub −0.28}{sup +0.31} at the 68% C.L. We discuss the effect of assuming different functional forms for mass profiles and of the orbit anisotropy in the kinematic reconstruction. Interpreting this result within the well-studied f(R) modified gravity model, the constraint on η translates into an upper bound to the interaction length (inverse of the scalaron mass) smaller than 2 Mpc. This tight constraint on the f(R) interaction range is however substantially relaxed when systematic uncertainties in the analysis are considered. Our analysis highlights the potential of this method to detect deviations from general relativity, while calling for the need of further high-quality data on the total mass distribution of clusters and improved control on systematic effects.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

A Polymer Canvas with the Stiffness of the Bone Matrix to Study and Control Mesenchymal Stem Cell Response

Reproducing in vitro the complex multiscale physical features of human tissues creates novel biomedical opportunities and fundamental understanding of cell–environment interfaces and interactions. While stiffness has been recognized as a key driver of cell behavior, systematic studies on the role of stiffness have been limited to values in the KPa–MPa range, significantly below the stiffness of bone. Here, in this work, a platform enabling the tuning of the stiffness of a biocompatible polymeric interface up to values characteristic of human bone is reported, which are in the GPa range, by using extremely thin polymer films on glass and cross-linking the films using ultraviolet (UV) light irradiation. It is shown that a higher stiffness is related to better adhesion, proliferation, and osteogenic differentiation, and that it is possible to switch on/off cell attachment and growth by solely tuning the stiffness of the interface, without any surface chemistry or topography modification. Since the stiffness is tuned directly by UV irradiation, this platform is ideal for rapid and simple fabrication of stiffness patterns and gradients, thus representing an innovative tool for combinatorial studies of the synergistic effect of tissue environmental cues on cell behavior, and creates new opportunities for next-generation biosensors, single-cell patterning, and lab-on-a-chip devices.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

A Deep Quench Approach to the Optimal Control of an Allen–Cahn Equation with Dynamic Boundary Conditions and Double Obstacles

In this paper, we investigate optimal control problems for Allen-Cahn variational inequalities with a dynamic boundary condition involving double obstacle potentials and the Laplace-Beltrami operator. The approach covers both the cases of distributed controls and of boundary controls. The cost functional is of standard tracking type, and box constraints for the controls are prescribed. We prove existence of optimal controls and derive first-order necessary conditions of optimality. The general strategy is the following: we use the results that were recently established by two of the authors for the case of (differentiable) logarithmic potentials and perform a so-called “deep quench limit”. Using compactness and monotonicity arguments, it is shown that this strategy leads to the desired first-order necessary optimality conditions for the case of (non-differentiable) double obstacle potentials.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Weak Necessary and Sufficient Stochastic Maximum Principle for Markovian Regime-Switching Diffusion Models

In this paper we prove a weak necessary and sufficient maximum principle for Markovian regime switching stochastic optimal control problems. Instead of insisting on the maximum condition of the Hamiltonian, we show that 0 belongs to the sum of Clarke’s generalized gradient of the Hamiltonian and Clarke’s normal cone of the control constraint set at the optimal control. Under a joint concavity condition on the Hamiltonian and a convexity condition on the terminal objective function, the necessary condition becomes sufficient. We give four examples to demonstrate the weak stochastic maximum principle.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Robust Utility Maximization Under Convex Portfolio Constraints

We study a robust maximization problem from terminal wealth and consumption under a convex constraints on the portfolio. We state the existence and the uniqueness of the consumption–investment strategy by studying the associated quadratic backward stochastic differential equation. We characterize the optimal control by using the duality method and deriving a dynamic maximum principle.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Hydrochemical constraints between the karst Tabular Middle Atlas Causses and the Saïs basin (Morocco): implications of groundwater circulation

The karst Tabular Middle Atlas Causses reservoir is the main drinking-water supply of Fez-Meknes region (Saïs Basin) in Morocco. Recent analyses showed a decline in associated groundwater chemical quality and increased turbidity. To understand this hydrosystem, four surveys were undertaken during fall and spring, 2009–2011. Hydrogeochemical studies coupled with isotopic analyses (δ{sup 18}O, δD and {sup 222}Rn) showed that the aquifers between the causses (mountains) and the Saïs Basin are of Liassic origin and at the southern extremities are of Triassic origin. Five recharge zones of different altitudes have been defined, including two main mixing zones in the south. Deuterium excess results suggest local recharge, while a plot of δ{sup 18}O versus δD characterizes a confined aquifer in the eastern sector. {sup 222}Rn results reveal areas of rapid exchanges with an upwelling time of less than 2 weeks. A schematic conceptual model is presented to explain the groundwater circulation system and the behavior of this karst system.

58 GEOSCIENCES↗
DOE OSTI2018

Optimization of groundwater artificial recharge systems using a genetic algorithm: a case study in Beijing, China

An optimization approach is used for the operation of groundwater artificial recharge systems in an alluvial fan in Beijing, China. The optimization model incorporates a transient groundwater flow model, which allows for simulation of the groundwater response to artificial recharge. The facilities’ operation with regard to recharge rates is formulated as a nonlinear programming problem to maximize the volume of surface water recharged into the aquifers under specific constraints. This optimization problem is solved by the parallel genetic algorithm (PGA) based on OpenMP, which could substantially reduce the computation time. To solve the PGA with constraints, the multiplicative penalty method is applied. In addition, the facilities’ locations are implicitly determined on the basis of the results of the recharge-rate optimizations. Two scenarios are optimized and the optimal results indicate that the amount of water recharged into the aquifers will increase without exceeding the upper limits of the groundwater levels. Optimal operation of this artificial recharge system can also contribute to the more effective recovery of the groundwater storage capacity.

58 GEOSCIENCES↗
DOE OSTI2018

On the Origin of a Small Hole Pocket in the Fermi Surface of Underdoped YBa{sub 2}Cu{sub 3}O{sub y}

The comprehensive experimental studies of the quantum oscillations (QOs) in YBa{sub 2}Cu{sub 3}O{sub y} (YBCO) lead N. Doiron-Leyraud et al. (Nature Communications 6:6034, 2015) to detection of small hole pockets below the transition into the charge-ordered (CO) phase. Their interpretation of the observed energy spectrum imposes serious constraints for the possible origin of electron and hole pockets—both of them should be due to the Fermi surface (FS) reconstruction at charge ordering with biaxial order parameter. However, the results of recent X-ray experiments which reveal the uniaxial CO in the underdoped YBCO make such interpretation questionable. We show that uniaxial CO parameter in combination with preexisting electron pocket at the Γ-point of Brillouin zone suggests a plausible interpretation of QOs associated with small hole pockets in YBCO.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Statistics of Point Vortex Turbulence in Non-neutral Flows and in Flows with Translational and Rotational Symmetries

A theory (Esler and Ashbee in J Fluid Mech 779:275–308, 2015) describing the statistics of N freely-evolving point vortices in a bounded two-dimensional domain is extended. First, the case of a non-neutral vortex gas is addressed, and it is shown that the density of states function can be identified with the probability density function of an infinite sum of independent non-central chi-squared random variables, the details of which depend only on the shape of the domain. Equations for the equilibrium energy spectrum and other statistical quantities follow, the validity of which are verified against direct numerical simulations of the equations of motion. Second, domains with additional conserved quantities associated with a symmetry (e.g., circle, periodic channel) are investigated, and it is shown that the treatment of the non-neutral case can be modified to account for the additional constraint.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Exact Maximum-Entropy Estimation with Feynman Diagrams

A longstanding open problem in statistics is finding an explicit expression for the probability measure which maximizes entropy with respect to given constraints. In this paper a solution to this problem is found, using perturbative Feynman calculus. The explicit expression is given as a sum over weighted trees.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗