Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “constrained dynamics”

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 703 records · Page 39

Large Exciton Diffusion Coefficients in Two Dimensional Covalent Organic Frameworks with Different Domain Sizes Revealed by Ultrafast Exciton Dynamics

Large singlet exciton diffusion lengths are a hallmark of high performance in organic-based devices such as photovoltaics, chemical sensors, and photodetectors. In this study, exciton dynamics of a two-dimensional covalent organic framework, 2D COF-S, is investigated using ultrafast spectroscopic techniques. After photoexcitation, the COF-S exciton decays via three pathways: (1) excimer formation (4 +/- 2 ps), (2) excimer relaxation (160 +/- 40 ps), and (3) excimer decay (>3 ns). Excitation fluence-dependent transient absorption studies suggest that COF-5 has a relatively large diffusion coefficient (0.08 cm(2)/s). Furthermore, exciton-exciton annihilation processes are characterized as a function of COF-S crystallite domain size in four different samples, which reveal domain-size-dependent exciton diffusion kinetics. These results reveal that exciton diffusion in COF-S is constrained by its crystalline domain size. These insights indicate the outstanding promise of delocalized excitonic processes available in 2D COFs, which motivate their continued design and implementation into optoelectronic devices.

Flanders, Nathan C.↗

Integrated flux and pool size analysis in plant central metabolism reveals unique roles of glycine and serine during photorespiration

Photorespiration is an essential process juxtaposed between plant carbon and nitrogen metabolism that responds to dynamic environments. Photorespiration recycles inhibitory intermediates arising from oxygenation reactions catalysed by Rubisco back into the C 3 cycle, but it is unclear what proportions of its nitrogen-containing intermediates (glycine and serine) are exported into other metabolisms in vivo and how these pool sizes affect net CO 2 gas exchange during photorespiratory transients. Here, to address this uncertainty, we measured rates of amino acid export from photorespiration using isotopically non-stationary metabolic flux analysis. This analysis revealed that ~23–41% of the photorespiratory carbon was exported from the pathway as serine under various photorespiratory conditions. Furthermore, we determined that the build-up and relaxation of glycine pools constrained a large portion of photosynthetic acclimation during photorespiratory transients. Furthermore, these results reveal the unique and important roles of glycine and serine in successfully maintaining various photorespiratory fluxes that occur under environmental fluctuations in nature and providing carbon and nitrogen for metabolism.

59 BASIC BIOLOGICAL SCIENCES↗

Comparative analysis of cosmological simulations via angular power spectra

Cosmology is the study of the universe's origin, large-scale structure, and dynamics. One such area of study is understanding dark energy and dark matter. Dark matter has not been observed to interact with luminous matter other than through gravitational attractions, forming large scale structures observed throughout the universe. In this area, the primary model is the Lambda Cold Dark Matter (ΛCDM) model. This standard model consists of cosmological dark energy (Λ), cold dark matter (CDM), and normal luminous matter. We use Angular Power Spectra (APS) to compare simulated ΛCDM models with observed large-scale galaxy distributions from the DECam Local Volume Exploration Survey (DELVE). After ensuring that the simulations can produce APS which match observational results, we vary cosmological parameters to constrain the allowed ranges of parameter values that are consistent with the observed data.

Feithen, William↗

Impact of Varying Dark Energy on Future Large-Scale Structure Studies

The cosmos withholds multiple mysteries such as dark forms of matter and energy that are yet beyond human comprehension. In the standard cosmological model, $\Lambda$CDM, the cosmological constant $\Lambda$ is thought to be responsible for the late accelerated expansion of the Universe. However, recent results from the Dark Energy Spectroscopic Instrument (DESI) suggest the possibility of evolving dark energy, which warrants further exploration. Future surveys such as the Rubin Observatory Legacy Survey of Space and Time (LSST) will map the large-scale structure (LSS) with unprecedented precision, giving us valuable statistical information about the cosmos. The goal of this research is to investigate the potential impact of an evolving dark energy scenario on cosmological parameters constrained by LSS probes, as will be mapped by the LSST. For this, we examine the power spectra of lens galaxies, source galaxies, and cross-power spectra between lens galaxies and source galaxies. The combination of these statistics is commonly referred to as '3 $\times$ 2 points'. For this investigation, we created a set of simulations resembling LSST data and used them to perform cosmological parameter inference in two scenarios: one in which the simulated data is based on the fiducial model and another on evolving dark energy. We then examined the degeneracy between the cosmological parameters and checked for potential shifts in the parametric space when the data contains dynamical dark energy but the modeling assumes $\Lambda$CDM. Our findings indicate that mismodeling the dark energy equation of state can significantly impact parameter inference, particularly affecting the total matter density, $\Omega_m$, and the growth of structures, as represented by the $S_8$ parameter. These results highlight the importance of further exploring extensions of the $\Lambda$CDM model in future LSS studies.

Yaman Acharya, A.↗

COMBLE ARM Mobile Facility (AMF) Measurements of Ice Nucleating Particles Field Campaign Report

During the Cold-Air Outbreaks in the Marine Boundary Layer Experiment (COMBLE), a project with the overarching goal to shed light on the dynamics and microphysical properties of clouds and precipitation in the high-latitude marine boundary layer during cold-air outbreaks (CAOs), we were tasked with providing and assisting the collection of aerosol filter samples for measuring ice nucleating particle (INP) concentrations. Little is known about the properties of CAO cloud systems, how they vary with surface, environmental, and aerosol conditions, the role of cold-air outbreaks in the global atmospheric and ocean circulation, and the accuracy of the treatment of this atmospheric regime in climate models. Uncertainties in climate projections for the arctic and sub-arctic region are large, in part due to the lack of observational guidance to constrain the treatment of aerosol-cloud-precipitation linkages in climate models.

54 ENVIRONMENTAL SCIENCES↗

Comparative Analysis of Cosmological Simulations via Angular Power Spectra

Cosmology is the study of the universe's origin, large-scale structure, and dynamics. One such area of study is understanding dark energy and dark matter. Dark matter has not been observed to interact with luminous matter other than through gravitational attractions, forming large scale structures observed throughout the universe. In this area, the primary model is the Lambda Cold Dark Matter (ΛCDM) model. This standard model consists of cosmological dark energy (Λ), cold dark matter (CDM), and normal luminous matter. We use angular power spectra (APS) to compare simulated ΛCDM models with observed large-scale galaxy distributions from the DECam Local Volume Exploration Survey (DELVE) [1]. After ensuring that the simulations can produce APS that match observational results, we vary cosmological parameters to constrain the allowed ranges of parameter values that are consistent with the observed data.

79 ASTRONOMY AND ASTROPHYSICS↗

Dynamic, Adaptive, Systems and Materials: Complex, Simple and Emergent Behaviors

This program has been funded by DoE/BES for twenty years. It has moved into and out of various subjects as it has developed, but it has retained its focus on complexity and complex systems. The project has evolved in the following way: Self-Assembly and Biomimetic Self-Assembly: All self-assembling systems depend upon a minimum of two types of interaction: a repulsion and an attraction. For the familiar molecular systems, attractive interactions are typically hydrogen bonds and electrostatic interactions. Repulsive interactions include steric effects, hydrophobic effects (in biological systems), and charge-charge repulsion. We have expanded this repertoire to include surface interactions, magnetic interactions, and others. I list these systems in the order in which we have explored them: i) A key emphasis in current work is in understanding how the movement of ions in a magnetic field (the Lorentz effect) interacts with catalytic systems. We have demonstrated that an acceleration in rate of reduction of CO 2 to CO can be accomplished by applying an external magnetic field. This acceleration is largely due to the application of the Lorentz effect on mass transport at the catalyst’s surface. ii) We have also extensively explored the influence of electrostatics, as exhibited in self-assembling systems, by tribocharging. iii) Another key system involves surface tension effects; examples include interactions between heavy particles floating at a liquid-air interface, and interacting by changes in surface area; interactions of bubbles and bubble rafts, behaviors of bubble trains in microfluidic networks, and behaviors of microorganisms in constraining environments. iv) This work has intentionally de-emphasized biological systems; but it does include some work on protein-ligand interactions and interactions among microorganisms. v) We have also explored applications of some of these effects, these explorations include bubble rafts as diffraction gratings, exploration of the structures that can be obtained by tribocharging and uses of these structures in exploring nucleation and melting of crystals. vi) Although not a major focus of this work, several other topics have emerged and offer opportunities for future work. These include the behavior of bubble trains and bubble rafts in microfluidic systems. A particularly interesting example is the formation of bubble trains that repeat in the alteration of large and small bubbles according to rules we do not presently understand, but are uniquely large-period oscillating systems. These systems offer a new route into understanding the instabilities of the type represented by oscillations. vii) We have also begun exploratory projects on magnetic levitation (especially to determine molecular density), and information storage (in molecules). Magnetic Levitation: Self-assembling and biomimetic systems require both attraction and repulsion. We have used electrostatics (tribocharging), interfacial free-energies (surface tension and related forces) and others. Potential uses include reconfigurable diffraction gratings and liquid lenses; exploration of mechanisms in tribocharging; tunneling in EGaIn junctions; and bubble trains (especially in micro-fluidic systems). Examples of systems representing these topics is included in the following papers: Complexity: Disks rotating at a water-air interface; Benard-Marangoni effects; Vortex-Crystals from spinning magnetic disks (Marangoni effects); EGaIn Electrode to study quantum tunneling; Self-Assembly of electrostatically-charged metallic spheres (electrets); Dynamically reconfigurable lens; Using computational designs of ligands for enzymes; Electrostatic self-assembly by tribocharging; Monodisperse bubble trains in microchannel systems; Inverted dripping faucet; Flames; Printing of micro-organisms to regenerate the “ink” of printing device; Using micro-organisms to move loads (“microoxen”); Motion of bacterial swarms near surfaces; Making monodisperse particles in microfluidic systems; Coding/decoding of information stored in droplet trains in microfluidic networks; Magnetic levitation; and Information storage. i) Tribocharging. The change in focus of this work on electrets from the fundamentals of charging to applications of these materials in studying self-assembly using electrostatic interactions. ii) Bubbles in Microchannels. The realization that systems of bubbles in microchannels represented a major opportunity to study complexity in a very tractable system, and the development of a semi-quantitative theory of this subject. iii) Flames. The growth of “flames” remains an exploratory subject for the research, although their currently relatively little active work involving it ongoing. iv) Systems with Microorganisms. The removal of work in biological systems from this project. Based on work supported in this program, we now have a significant project on the development of microfluidic tools for studying C. elegans (a nematode), but this work was not appropriate for a program focused on complexity, and we developed separate support for it. (It is, however, an example of successful seeding of a new area by BES.) The work on electrets has gone through a period in which a part of the program was the subject of a MURI; the focus of this work was to develop materials that did not charge electrostatically on friction or contact. The MURI is now over, and the work on dynamic self-assembly (supported by BES) is the major focus. “Flames” has also enjoyed synergistic support, in terms of a project supported by DARPA on flame suppression (in the absence of extinguishing agents, using acoustic and electrostatic interactions). This work was helpful in understanding some of the basics of flames, but is entirely distinct from the BES focus in complexity. A growing interest is in the Lorentz effect. The Lorentz effect is the force exerted on charged particles (electrons, ions, charged molecules) when they move through a perpendicular magnetic field. The Lorentz effect is almost ubiquitous in modern technology: examples of applications include electric motors, dynamos, cathode ray tubes, many batteries, and most systems that control electrical currents with magnetic forces. We have begun to explore the Lorentz effect in electrochemical systems and heterogeneous catalytic systems involving charged organic species and inorganic ions. This work is still at an early stage, but initial studies that Lorentz effects can be large when ions move through magnetic fields, or magnetic fields move in the presence of ions.

36 MATERIALS SCIENCE↗

Predicting multiphase flow and tracer transport for an underground chemical explosive test

Detecting radionuclide gas seepage from clandestine underground nuclear tests is central to nonproliferation explosion monitoring research. Yet, early-time (<6 day) gas transport driven by the explosive pressure wave remains poorly constrained due to scarcity of field data. We simulate multi-phase gas transport in the vadose zone using pre-shot data from a recent chemical explosion in P-Tunnel at the Nevada National Security Site, USA. Despite using a simplified 2D-radial model, predictions of tracer arrival matched observations within one order-of-magnitude. Our results show how transient blast forcing rapidly mobilizes gases from the cavity into surrounding rock – critical for optimizing sensor placement and test planning. This unique integration of field data and modeling represents a significant improvement in our ability to predict gas migration from underground explosions. More broadly, it offers insights into the coupled dynamics of pressure waves and contaminant transport in the vadose zone, with implications for monitoring and hazard assessment.

54 ENVIRONMENTAL SCIENCES↗

Using quasar X-ray and UV flux measurements to constrain cosmological model parameters

ABSTRACT Risaliti and Lusso have compiled X-ray and UV flux measurements of 1598 quasars (QSOs) in the redshift range 0.036 ≤ z ≤ 5.1003, part of which, z ∼ 2.4 − 5.1, is largely cosmologically unprobed. In this paper we use these QSO measurements, alone and in conjunction with baryon acoustic oscillation (BAO) and Hubble parameter [H(z)] measurements, to constrain cosmological parameters in six different cosmological models, each with two different Hubble constant priors. In most of these models, given the larger uncertainties, the QSO cosmological parameter constraints are mostly consistent with those from the BAO + H(z) data. A somewhat significant exception is the non-relativistic matter density parameter Ωm0 where QSO data favour Ωm0 ∼ 0.5 − 0.6 in most models. As a result, in joint analyses of QSO data with H(z) + BAO data the 1D Ωm0 distributions shift slightly towards larger values. A joint analysis of the QSO + BAO + H(z) data is consistent with the current standard model, spatially-flat ΛCDM, but mildly favours closed spatial hypersurfaces and dynamical dark energy. Since the higher Ωm0 values favoured by QSO data appear to be associated with the z ∼ 2 − 5 part of these data, and conflict somewhat with strong indications for Ωm0 ∼ 0.3 from most z < 2.5 data as well as from the cosmic microwave background anisotropy data at z ∼ 1100, in most models, the larger QSO data Ωm0 is possibly more indicative of an issue with the z ∼ 2 − 5 QSO data than of an inadequacy of the standard flat ΛCDM model.

79 ASTRONOMY AND ASTROPHYSICS↗

Use of Shallow Ice Core Measurements to Evaluate and Constrain 1980–1990 Global Reanalyses of Ice Sheet Precipitation Rates

Abstract Sea‐level rise (SLR) projections by Earth System Models (ESMs) depend on ice sheet surface mass balances. Accurate, global atmosphere reanalyses would be ideal for providing equilibrated ice sheet model initial conditions in fully coupled ESM simulations. Here we present the first evaluation of 1980–1990 global reanalysis precipitation over Greenland and Antarctica that uses independent observations of net accumulation rates derived from shallow ice cores. Precipitation distributions from both the European Centre for Medium‐Range Weather Forecast's Reanalysis (ERA5) and the Modern‐Era Retrospective Analysis for Research and Applications (MERRA‐2) are highly correlated with contemporaneous co‐located net accumulation rates from Greenland ( r 2 > 0.95) and West Antarctica ( r 2 > 0.7). Three other commonly used reanalyses (WFDE5, CRUNCEP, and GSWP3) exhibit significantly weaker correlations on one or both ice sheets. Our findings imply that ESMs should use ERA5 or MERRA‐2 in data‐forced simulations to validate ice sheet model dynamics and precondition firn for SLR projections.

Schneider, Adam↗

Nutrient limitation shapes functional traits of mycorrhizal fungi and phosphorus-cycling bacteria across an elevation gradient

In nutrient-limited high-elevation ecosystems, plants rely on arbuscular mycorrhizal (AM) fungi to provide mineral phosphorus (P) in the form of phosphate (PO43-). AM fungi gather these nutrients from phosphorus-cycling bacteria (PCBs) that can mineralize PO43- from organic matter and solubilize mineral-bound P. How climate, soil factors, and nutrient limitation influence AM fungi and PCB assembly remains unclear. We collected soil from montane meadows across a 1,000-m elevation gradient on three replicate mountainsides and analyzed AM fungal marker genes, P-cycling genes from shotgun metagenomes, and edaphic measurements. High-elevation soils had nearly 50-fold less soil PO₄³⁻ and 60% more AM fungal hyphae than low-elevation soils. AM fungal turnover was linked to changes in pH, organic carbon, and PO₄³-. The composition of 198 P-cycling genes was influenced by the AM fungal community structure. Drivers of individual PCB functional genes, including pH and organic carbon, varied with gene phylogeny. We found a trade-off in P-cycling strategies across elevation: P-rich, low-elevation soils supported root-colonizing AM fungi and organic P-mineralizing bacteria. P-poor, high-elevation soils were dominated by stress-tolerant AM fungi and mineral P-solubilizing bacteria. Our results suggest that AM fungi and PCB community turnover across elevation are both shaped by pH, organic carbon, and P availability. With continued climate warming, the structure and function of mountaintop ecosystems might shift to resemble lower elevations, disrupting long-established and specialized microbial assemblages, with consequences for P-cycling dynamics and the total P available to plant communities.IMPORTANCEPhosphorus (P) limits plant productivity in high-elevation ecosystems, yet the microbial networks that mobilize P, including arbuscular mycorrhizal (AM) fungi and phosphorus-cycling bacteria (PCBs), remain under-characterized in these nutrient-poor soils. We show that across a 10,00-m elevation gradient, AM fungi and P-cycling gene assemblages shift predictably with pH, organic carbon, and phosphate availability. Higher elevations, with less available P, select for stress-tolerant AM fungal taxa and PCB strategies geared toward mineral solubilization, while low-elevation sites favor root colonization by AM fungi and organic P mineralization. These results suggest that nutrient limitation can constrain microbial community assembly in consistent ways across landscapes. High mountain soils are low in P and rely on a network of underground AM fungi and PCB to deliver nutrients to plants. This study shows how those underground relationships reorganize with elevation and how climate change could collapse long-standing microbial strategies by pushing high-elevation ecosystems toward lowland conditions. As soils warm and dry, the microbial scaffolding that supports alpine plant life may become increasingly unstable.

arbuscular mycorrhizal fungi↗

A two-phase substrate model for enzymatic hydrolysis of lignocellulose: application to batch and continuous reactors

Abstract Background Enzymatic hydrolysis continues to have a significant projected production cost for the biological conversion of biomass to fuels and chemicals, motivating research into improved enzyme and reactor technologies in order to reduce enzyme usage and equipment costs. However, technology development is stymied by a lack of accurate and computationally accessible enzymatic-hydrolysis reaction models. Enzymatic deconstruction of cellulosic materials is an exceedingly complex physico-chemical process. Models which elucidate specific mechanisms of deconstruction are often too computationally intensive to be accessible in process or multi-physics simulations, and empirical models are often too inflexible to be effectively applied outside of their batch contexts. In this paper, we employ a phenomenological modeling approach to represent rate slowdown due to substrate structure (implemented as two substrate phases) and feedback inhibition, and apply the model to a continuous reactor system. Results A phenomenological model was developed in order to predict glucose and solids concentrations in batch and continuous enzymatic-hydrolysis reactors from which liquor is continuously removed by ultrafiltration. A series of batch experiments were performed, varying initial conditions (solids, enzyme, and sugar concentrations), and best-fit model parameters were determined using constrained nonlinear least-squares methods. The model achieved a good fit for overall sugar yield and insoluble solids concentration, as well as for the reduced rate of sugar production over time. Additionally, without refitting model coefficients, good quantitative agreement was observed between results from continuous enzymatic-hydrolysis experiments and model predictions. Finally, the sensitivity of the model to its parameters is explored and discussed. Conclusions Although the phenomena represented by the model correspond to behaviors that emerge from clusters of mechanisms, and hence a set of model coefficients are unique to the substrate and the enzyme system, the model is efficient to solve and may be applied to novel reactor schema and implemented in computational fluid dynamics (CFD) simulations. Hence, this modeling approach finds the right balance between model complexity and computational efficiency. These capabilities have broad application to reactor design, scale-up, and process optimization.

09 BIOMASS FUELS↗

Characterizing the initial state and dynamical evolution in XeXe and PbPb collisions using multiparticle cumulants

For the first time, correlations among mixed-order moments of two or three flow harmonics —( v n k , v m l ) and ( v n k , v m l , v p q ), with k, l, and q denoting the respective orders—are measured in xenon–xenon (XeXe) collisions and compared with lead–lead (PbPb) results, providing a novel probe of collective behavior in heavy ion collisions. These measurements compare a nearly spherical, doubly-magic 208 Pb nucleus to a triaxially deformed 129 Xe nucleus, emphasizing the sensitivity to initial-state geometry fluctuations arising from nuclear deformation. The dependence of these results (v n , n = 2 , 3 , 4 ) on the shape and size of the nuclear overlap region is studied. Comparisons between v 2 , v 3 , and v 4 demonstrate the importance of v 3 and v 4 in exploring the nonlinear hydrodynamic response of the quark-gluon plasma (QGP) to the initial spatial anisotropy. The results constrain initial-state model parameters that influence the evolution of the QGP. The CMS detector was used to collect XeXe and PbPb data at nucleon-nucleon center-of-mass energies of S N N = 5.44 and 5.36 TeV, respectively. Correlations are extracted using multiparticle mixed-harmonic cumulants (up to eight-particle cumulants) with charged particles in the pseudorapidity range |η| < 2.4 and transverse momentum range 0.5 < p T < 3 GeV/c.

CMS↗

PET-FBA: A lightweight enzyme allocation and thermodynamics-constrained flux analysis approach to explore Escherichia coli metabolic adaptation to intracellular acidification

Escherichia coli employs diverse strategies to adapt to acidic environments that disrupt enzyme activity and the thermodynamic feasibility of essential reactions. To understand the impact of pH stress on cell metabolism, we present the PET-FBA (pH-, Enzyme protein allocation-, and Thermodynamics-constrained Flux Balance Analysis) framework. PET-FBA extends genome-scale modeling by integrating enzyme protein costs and reaction Gibbs free energy changes. Additionally, by incorporating pH-dependent enzyme kinetics in response to intracellular acidification, this framework enables the simulation of E. coli's metabolic adjustments across varying external pH levels. The model's accuracy is validated by comparing in silico growth simulations with experimental measurements under both anaerobic and aerobic conditions, as well as in silico gene knockouts of essential genes. By explicitly incorporating pH effects, our model accurately replicates the metabolic shift towards lactate production as the primary fermentation product at low pH in anaerobic conditions. This shift is only predicted when enzyme kinetics are dynamically adjusted as a function of pH. Further analysis revealed that this shift can be attributed to the reduced protein efficiency of the acetyl-CoA branch compared to lactate dehydrogenase under acidic stress, which then becomes crucial for maintaining NAD regeneration and cell growth at low pH. Furthermore, we identified strategies for enhancing cell growth under acidic anaerobic conditions by improving the enzyme activity of lactate dehydrogenase and pyruvate formate lyase, which increases NAD production efficiency and reduces enzyme protein allocation costs. Designed as a lightweight yet versatile framework, PET-FBA enables efficient genome-scale metabolic analysis. Using E. coli as a model system, our framework provides a systematic approach to understanding metabolic responses to environmental stress, pinpointing key metabolic bottlenecks, and identifying potential targets for strain optimization.

42 ENGINEERING↗

A data-driven peridynamic continuum model for upscaling molecular dynamics

Nonlocal models, including peridynamics, often use integral operators that embed lengthscales in their definition. However, the integrands in these operators are difficult to define from the data that are typically available for a given physical system, such as laboratory mechanical property tests. In contrast, molecular dynamics (MD) does not require these integrands, but it suffers from computational limitations in the length and time scales it can address. To combine the strengths of both methods and to obtain a coarse-grained, homogenized continuum model that efficiently and accurately captures materials’ behavior, we propose a learning framework to extract, from MD data, an optimal Linear Peridynamic Solid (LPS) model as a surrogate for MD displacements. To maximize the accuracy of the learnt model we allow the peridynamic influence function to be partially negative, while preserving the well-posedness of the resulting model. To achieve this, we provide sufficient well-posedness conditions for discretized LPS models with sign-changing influence functions and develop a constrained optimization algorithm that minimizes the equation residual while enforcing such solvability conditions. This framework guarantees that the resulting model is mathematically well-posed, physically consistent, and that it generalizes well to settings that are different from the ones used during training. We illustrate the efficacy of the proposed approach with several numerical tests for single layer graphene. Our two-dimensional tests show the robustness of the proposed algorithm on validation data sets that include thermal noise, different domain shapes and external loadings, and discretizations substantially different from the ones used for training.

homogenization↗

Domain Aware Deep-learning Algorithms Integrated with Scientific-computing Technologies (DADAIST)

This technical report summarized the contribution of the DADAIST project funded by the Data Model Convergence Initiative via the Laboratory Directed Research and Development (LDRD) investments at Pacific Northwest National Laboratory (PNNL). Specifically, we report the development of the NeuroMANCER (Neural Modules with Adaptive Nonlinear Constraints and Efficient Regularizations), a new open-source Scientific Machine Learning library for formulating and solving parametric constrained optimization problems, physics-informed system identification, and parametric optimal control problems. NeuroMANCER is using differentiable programming to combine modern data-driven models and optimization modeling language into a coherent algorithmic and software framework. NeuroMANCER is a Pytorch-based framework and adopts much of its philosophy focused on research and development, rapid prototyping, and streamlined deployment. Strong emphasis is given to extensibility, interoperability with the PyTorch ecosystem, and quick adaptability to custom domain problems. Neuromancer repository contains a comprehensive library of differentiable modules, including custom activation functions, matrix factorizations, deep learning architectures, neural differential equations, differential equation solvers, implicit layers such as iterative solvers, high-level API for symbolic expressions, API for modeling and control of dynamical systems, and extensive set of tutorial code examples in the form of python scripts and jupyter notebooks.

97 MATHEMATICS AND COMPUTING↗

Planck 2018 results

We report on the implications for cosmic inflation of the 2018 release of the Planck cosmic microwave background (CMB) anisotropy measurements. The results are fully consistent with those reported using the data from the two previous Planck cosmological releases, but have smaller uncertainties thanks to improvements in the characterization of polarization at low and high multipoles. Planck temperature, polarization, and lensing data determine the spectral index of scalar perturbations to be n s = 0.9649 ± 0.0042 at 68% CL. We find no evidence for a scale dependence of n s , either as a running or as a running of the running. The Universe is found to be consistent with spatial flatness with a precision of 0.4% at 95% CL by combining Planck with a compilation of baryon acoustic oscillation data. The Planck 95% CL upper limit on the tensor-to-scalar ratio, r0.002 < 0.10, is further tightened by combining with the BICEP2/Keck Array BK15 data to obtain r 0.002 < 0.056. In the framework of standard single-field inflationary models with Einstein gravity, these results imply that: (a) the predictions of slow-roll models with a concave potential, V"(Φ) < 0, are increasingly favoured by the data; and (b) based on two different methods for reconstructing the inflaton potential, we find no evidence for dynamics beyond slow roll. Three different methods for the non-parametric reconstruction of the primordial power spectrum consistently confirm a pure power law in the range of comoving scales 0.005 Mpc -1 ≲ k ≲ 0.2 Mpc -1 . A complementary analysis also finds no evidence for theoretically motivated parameterized features in the Planck power spectra. For the case of oscillatory features that are logarithmic or linear in k, this result is further strengthened by a new combined analysis including the Planck bispectrum data. The new Planck polarization data provide a stringent test of the adiabaticity of the initial conditions for the cosmological fluctuations. In correlated, mixed adiabatic and isocurvature models, the non-adiabatic contribution to the observed CMB temperature variance is constrained to 1.3%, 1.7%, and 1.7% at 95% CL for cold dark matter, neutrino density, and neutrino velocity, respectively. Planck power spectra plus lensing set constraints on the amplitude of compensated cold dark matter-baryon isocurvature perturbations that are consistent with current complementary measurements. The polarization data also provide improved constraints on inflationary models that predict a small statistically anisotropic quadupolar modulation of the primordial fluctuations. However, the polarization data do not support physical models for a scale-dependent dipolar modulation. All these findings support the key predictions of the standard single-field inflationary models, which will be further tested by future cosmological observations.

79 ASTRONOMY AND ASTROPHYSICS↗

The Gemini Planet Imager Exoplanet Survey: Dynamical Mass of the Exoplanet β Pictoris b from Combined Direct Imaging and Astrometry

Herein, we present new observations of the planet β Pictoris b from 2018 with the Gemini Planet Imager (GPI), the first GPI observations following conjunction. Based on these new measurements, we perform a joint orbit fit to the available relative astrometry from ground-based imaging, the Hipparcos Intermediate Astrometric Data (IAD), and the Gaia DR2 position, and demonstrate how to incorporate the IAD into direct imaging orbit fits. We find a mass consistent with predictions of hot-start evolutionary models and previous works following similar methods, though with larger uncertainties: $12.8^{+5.3}_{–3.2} M_{Jup}$. Our eccentricity determination of ${0.12}_{-0.03}^{+0.04}$ disfavors circular orbits. We consider orbit fits to several different imaging data sets, and find generally similar posteriors on the mass for each combination of imaging data. Our analysis underscores the importance of performing joint fits to the absolute and relative astrometry simultaneously, given the strong covariance between orbital elements. Time of conjunction is well-constrained within 2.8 days of 2017 September 13, with the star behind the planet's Hill sphere between 2017 April 11 and 2018 February 16 (±18 days). Following the recent radial velocity detection of a second planet in the system, β Pic c, we perform additional two-planet fits combining relative astrometry, absolute astrometry, and stellar radial velocities. These joint fits find a significantly smaller mass (8.0 ± 2.6 M Jup) for the imaged planet β Pic b, in a somewhat more circular orbit. We expect future ground-based observations to further constrain the visual orbit and mass of the planet in advance of the release of Gaia DR4.

79 ASTRONOMY AND ASTROPHYSICS↗